geometry

Module

Description

geometry

References:

None

Owner:

Daniel Price

Runtime parameters:

None

Dependencies:

None

Quick access

Variables:

cosa, cosb, hp_on_rp, ierr_invalid_dims, ierr_invalid_dimsin, ierr_invalid_dimsout, ierr_not_implemented, ierr_warning_assuming_cartesian, igeom_cartesian, igeom_cylindrical, igeom_flaredcyl, igeom_logflared, igeom_planetwake, igeom_rotated, igeom_spherical, igeom_toroidal, integrate_trap, labelcoord, labelcoordsys, large_number, maxcoordsys, p_index, phi_p, planet_t_integrand, planet_wake_phiw, q_index, rp, rtorus, sina, sinb, small_number, xref

Routines:

coord_is_length(), coord_is_periodic(), coord_transform(), coord_transform_limits(), get_coord_limits(), planet_wake_eta(), planet_wake_phi(), planet_wake_r(), planet_wake_t(), print_error(), set_flaring_index(), set_planet_wake(), set_rotation_angles(), vector_transform()

Variables

  • geometry/beta [real,private/optional/default=1.5]

    dzdash/dz

  • geometry/cosa [real,private/optional/default=sqrt(5.)/3.]
  • geometry/cosb [real,private/optional/default=1./sqrt(5.)]
  • geometry/hp_on_rp [real,private/optional/default=0.05]
  • geometry/ierr_invalid_dims [integer,public/parameter/optional/default=3]
  • geometry/ierr_invalid_dimsin [integer,public/parameter/optional/default=1]
  • geometry/ierr_invalid_dimsout [integer,public/parameter/optional/default=2]
  • geometry/ierr_not_implemented [integer,public/parameter/optional/default=-2]
  • geometry/ierr_warning_assuming_cartesian [integer,public/parameter/optional/default=-1]
  • geometry/igeom_cartesian [integer,public/parameter/optional/default=1]
  • geometry/igeom_cylindrical [integer,public/parameter/optional/default=2]
  • geometry/igeom_flaredcyl [integer,public/parameter/optional/default=6]
  • geometry/igeom_logflared [integer,public/parameter/optional/default=7]
  • geometry/igeom_planetwake [integer,public/parameter/optional/default=8]
  • geometry/igeom_rotated [integer,public/parameter/optional/default=5]
  • geometry/igeom_spherical [integer,public/parameter/optional/default=3]
  • geometry/igeom_toroidal [integer,public/parameter/optional/default=4]
  • geometry/integrate_trap [private]
  • geometry/labelcoord (3,8) [character,public/parameter/optional/default=reshape((/'x     ','y     ','z     ','r     ','\phi  ','z     ','r     ','\phi  ','\theta','r_t   ','\theta','\phi  ','x_1   ','x_2   ','x_3   ','r     ','\phi  ','zdash ','log r ','\phi  ','zdash ','t     ','\eta  ','z     '/),shape=(/3,maxcoordsys/))]
  • geometry/labelcoordsys (8) [character,public/parameter/optional/default=(/'cartesian   x,y,z      ','cylindrical r,phi,z    ','spherical   r,phi,theta','toroidal    r,theta,phi','rotated     x_1,x_2,x_3','flared cyl  r,phi,zdash','log flared  logr,phi,zd','planet wake t,eta,z    '/)]
  • geometry/large_number [real,public/parameter/optional/default=1e+30]
  • geometry/maxcoordsys [integer,public/parameter/optional/default=8]
  • geometry/p_index [real,private/optional/default=1.0]
  • geometry/phi_p [real,private/optional/default=0.0]
  • geometry/pi [real,private/parameter/optional/default=4.*atan(1.)]
  • geometry/planet_t_integrand [private]
  • geometry/planet_wake_phiw [private]
  • geometry/q_index [real,private/optional/default=0.25]
  • geometry/rp [real,private/optional/default=1.0]
  • geometry/rtorus [real,private/parameter/optional/default=1.0]
  • geometry/sina [real,private/optional/default=0.6666666666666666]
  • geometry/sinb [real,private/optional/default=2./sqrt(5.)]
  • geometry/small_number [real,public/parameter/optional/default=1e-15]
  • geometry/unknown_abstract_interface [private]
  • geometry/xref [real,private/optional/default=1.0]

Subroutines and functions

function  geometry/coord_is_length(ix, igeom)
Parameters:

  • ix [integer,in]

  • igeom [integer,in]

Return:

coord_is_length [logical,pure]

function  geometry/coord_is_periodic(ix, igeom)
Parameters:

  • ix [integer,in]

  • igeom [integer,in]

Return:

coord_is_periodic [logical,pure]

subroutine  geometry/print_error(ierr)
Parameters:

ierr [integer,in]

subroutine  geometry/set_rotation_angles([a[, b[, sin_a[, sin_b[, cos_a[, cos_b]]]]]])
Parameters:

  • a [real,in,]

  • b [real,in,]

  • sin_a [real,in,]

  • sin_b [real,in,]

  • cos_a [real,in,]

  • cos_b [real,in,]

subroutine  geometry/set_flaring_index(r_ref, findex)
Parameters:

  • r_ref [real,in]

  • findex [real,in]

subroutine  geometry/set_planet_wake(r, phi, p, q, honr)
Parameters:

  • r [real,in]

  • phi [real,in]

  • p [real,in]

  • q [real,in]

  • honr [real,in]

function  geometry/planet_wake_t(r)
Parameters:

r [real,in]

Return:

t [real]

Called from:

get_coord_limits()

function  geometry/planet_wake_r(t)
Parameters:

t [real,in]

Return:

r [real]

function  geometry/planet_wake_eta(r, phi)
Parameters:

  • r [real,in]

  • phi [real,in]

Return:

eta [real]

Called from:

get_coord_limits()

function  geometry/planet_wake_phi(r, eta)
Parameters:

  • r [real,in]

  • eta [real,in]

Return:

phi [real]

subroutine  geometry/coord_transform(xin, ndimin, itypein, xout, ndimout, itypeout[, err])
Parameters:

  • xin (ndimin) [real,in]

  • ndimin [integer,in,]

  • itypein [integer,in]

  • xout (ndimout) [real,out]

  • ndimout [integer,in]

  • itypeout [integer,in]

  • err [integer,out,]

Called from:

set_density_profile(), test_geometry()

subroutine  geometry/vector_transform(xin, vecin, ndimin, itypein, vecout, ndimout, itypeout[, err])
Parameters:

  • xin (ndimin) [real,in]

  • vecin (ndimin) [real,in]

  • ndimin [integer,in,]

  • itypein [integer,in]

  • vecout (ndimout) [real,out]

  • ndimout [integer,in]

  • itypeout [integer,in]

  • err [integer,out,]

Called from:

test_geometry()

subroutine  geometry/coord_transform_limits(xmin, xmax, itypein, itypeout, ndim)
Parameters:

  • xmin (ndim) [real,inout]

  • xmax (ndim) [real,inout]

  • itypein [integer,in]

  • itypeout [integer,in]

  • ndim [integer,in,]

subroutine  geometry/get_coord_limits(rad, xin, xout, xmin, xmax, itypein)
Parameters:

  • rad [real,in]

  • xin (3) [real,in]

  • xout (3) [real,out]

  • xmin (3) [real,out]

  • xmax (3) [real,out]

  • itypein [integer,in]

Called from:

test_geometry()

Call to:

planet_wake_t(), planet_wake_eta()