HARM
harm and utilities
 All Data Structures Files Functions Variables Typedefs Macros Pages
Macros | Functions | Variables
init.c File Reference

USER initial conditions (usually links to code in initbasecode directory) More...

#include "decs.h"
Include dependency graph for init.c:

Go to the source code of this file.

Macros

#define SLOWFAC   1.0 /* reduce u_phi by this amount */
 
#define MAXPASSPARMS   10
 
#define USER_THETAROTMETRIC   (0.0)
 
#define USER_THETAROTPRIMITIVES   (0.0)
 
#define NOFIELD   -1
 
#define DISK1FIELD   0
 
#define DISK2FIELD   1
 
#define VERTFIELD   2
 
#define DISK1VERT   3
 
#define DISK2VERT   4
 
#define BLANDFORDQUAD   5
 
#define TOROIDALFIELD   6
 
#define OHSUGAFIELD   7
 
#define MONOPOLAR   8
 
#define OLEKFIELD   9
 
#define FIELDJONMAD   10
 
#define FIELDWALD   11
 
#define MONOPOLE   12
 
#define SPLITMONOPOLE   13
 
#define DIMVARLIST   GGG,CCCTRUE,MSUNCM,MPERSUN,LBAR,TBAR,VBAR,RHOBAR,MBAR,ENBAR,UBAR,TEMPBAR,ARAD_CODE_DEF,XFACT,YFACT,ZFACT,MUMEAN,MUMEAN,OPACITYBAR,MASSCM,KORAL2HARMRHO(1.0),TEMPMIN
 
#define DIMTYPELIST   "%21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g\n"
 
#define NTUBE   31
 
#define WHICHRADATM   0
 
#define fsx(x)   ((x)>0 ? exp(-1.0/(x)) : 0.0)
 
#define gsx(x)   (fsx(x)/(fsx(x) + fsx(1.0-(x))))
 
#define stepfunction(x)   (1.0 - gsx((x)-0.5))
 
#define stepfunctionab(x, a, b)   (((a)-(b))*stepfunction(x) + (b))
 
#define stepfunctionab2(x, a, b)   (((a)-(b))*stepfunction(x)*stepfunction(x) + (b))
 
#define NUMTgasITERS   4
 
#define JONMADHPOW   (4.0)
 
#define JONMADR0   (0.0)
 
#define JONMADROUT   (300.0)
 
#define QPOWER   (1.0)
 
#define POWERNU   (2.0)
 
#define FCOVDERTYPE   DIFFGAMMIE
 
#define FCOVDXDELTA   1E-5
 
#define TAUADJUSTATM   (10.0)
 
#define KAPPA_ES_FERMICORR(rhocode, Tcode)   (1.0/(1.0+2.7E11*((rhocode)*RHOBAR)/prpow((Tcode)*TEMPBAR,2.0)))
 PURE ELASTIC SCATTERING. More...
 
#define KAPPA_ES_KNCORR(rhocode, Tcode)   (1.0/(1.0+prpow((Tcode)*TEMPBAR/4.5E8,0.86)))
 
#define KAPPA_ES_CODE(rhocode, Tcode)   (0.2*(1.0+XFACT)*KAPPA_ES_FERMICORR(rhocode,Tcode)*KAPPA_ES_KNCORR(rhocode,Tcode)/OPACITYBAR)
 kappaes = sigma_T n_e = sigma_T n_b (n_e/n_b) = sigma_T rho/mb (ne/nb) More...
 
#define KAPPA_ES_BASIC_CODE(rhocode, Tcode)   (0.2*(1.0+XFACT)/OPACITYBAR)
 
#define KAPPA_FORCOMPT_RELCORREP(ep)   ((1.0 + 3.683*(ep)+4.0*(ep)*(ep))/(1.0 + (ep)))
 
#define KAPPA_FORCOMPT_RELCORR(rhocode, Tcode)   (KAPPA_FORCOMPT_RELCORREP(K_BOLTZ*(Tcode)*TEMPBAR/(MELE*CCCTRUE*CCCTRUE)))
 
#define KAPPA_FORCOMPT_CODE(rhocode, Tcode)   (0.2*(1.0+XFACT)*KAPPA_ES_FERMICORR(rhocode,Tcode)*KAPPA_FORCOMPT_RELCORR(rhocode,Tcode)/OPACITYBAR)
 
#define KAPPA_ZETA(Tgcode, Trcode)   ((TEMPMIN+Trcode)/(TEMPMIN+Tgcode))
 EMISSION (Tr=Tg) or ABSORBPTION (Tr different from Tg) More...
 
#define KAPPA_FF_CODE(rhocode, Tgcode, Trcode)   (4.0E22*(1.0+XFACT)*(1.0-ZFACT)*((rhocode)*RHOBAR)*prpow((Tgcode)*TEMPBAR,-0.5)*prpow((Trcode)*TEMPBAR,-3.0)*prlog(1.0+1.6*KAPPA_ZETA(Tgcode,Trcode))*(1.0+4.4E-10*(Tgcode*TEMPBAR))/OPACITYBAR)
 
#define KAPPAN_FF_CODE(rhocode, Tgcode, Trcode)   KAPPA_FF_CODE(rhocode,Tgcode,Trcode)
 
#define KAPPA_BF_CODE(rhocode, Tgcode, Trcode)   (3.0E25*(ZFACT)*(1.0+XFACT+0.75*YFACT)*((rhocode)*RHOBAR)*prpow((Tgcode)*TEMPBAR,-0.5)*prpow((Trcode)*TEMPBAR,-3.0)*prlog(1.0+1.6*KAPPA_ZETA(Tgcode,Trcode))/OPACITYBAR)
 
#define KAPPA_CHIANTIBF_CODE(rhocode, Tgcode, Trcode)   (4.0E34*((rhocode*RHOBAR))*(ZFACT/ZSOLAR)*YELE*prpow((Tgcode)*TEMPBAR,-1.7)*prpow((Trcode)*TEMPBAR,-3.0)/OPACITYBAR)
 
#define KAPPA_HN_CODE(rhocode, Tgcode, Trcode)   (1.1E-25*prpow(ZFACT,0.5)*prpow((rhocode)*RHOBAR,0.5)*prpow((Tgcode)*TEMPBAR,7.7)/OPACITYBAR)
 
#define KAPPA_MOL_CODE(rhocode, Tgcode, Trcode)   (0.1*ZFACT/OPACITYBAR)
 
#define KAPPA_GENFF_CODE(rhocode, Tgcode, Trcode)   (1.0/(1.0/(KAPPA_MOL_CODE(rhocode,Tgcode,Trcode)+KAPPA_HN_CODE(rhocode,Tgcode,Trcode)) + 1.0/(KAPPA_CHIANTIBF_CODE(rhocode,Tgcode,Trcode)+KAPPA_BF_CODE(rhocode,Tgcode,Trcode)+KAPPA_FF_CODE(rhocode,Tgcode,Trcode))))
 
#define KAPPA   0.
 SYNCH STUFF. More...
 
#define KAPPAES   0.
 
#define KAPPAUSER(rho, B, Tg, Tr)   (rho*KAPPA*KAPPA_FF_CODE(rho,Tg,Tr))
 
#define KAPPAESUSER(rho, T)   (rho*KAPPAES*KAPPA_ES_BASIC_CODE(rho,T))
 
#define KAPPA   0.
 SYNCH STUFF. More...
 
#define KAPPAES   (1E3)
 
#define KAPPAUSER(rho, B, Tg, Tr)   (rho*KAPPA)
 
#define KAPPAESUSER(rho, T)   (rho*KAPPAES)
 
#define KAPPA   0.
 SYNCH STUFF. More...
 
#define KAPPAES   0.
 
#define KAPPAUSER(rho, B, Tg, Tr)   (rho*KAPPA*KAPPA_FF_CODE(rho,Tg,Tr))
 
#define KAPPAESUSER(rho, T)   (rho*KAPPAES*KAPPA_ES_BASIC_CODE(rho,T))
 
#define KAPPAESUSER(rho, T)   (0.0)
 
#define KAPPAUSER(rho, B, Tg, Tr)   (25*rho)
 
#define KAPPAUSER(rho, B, Tg, Tr)   (rho*1.0)
 
#define KAPPAESUSER(rho, T)   (rho*0.0)
 
#define KAPPA   0.
 SYNCH STUFF. More...
 
#define KAPPAES   0.
 
#define KAPPAUSER(rho, B, Tg, Tr)   (rho*KAPPA*KAPPA_FF_CODE(rho,Tg,Tr))
 
#define KAPPAESUSER(rho, T)   (rho*KAPPAES*KAPPA_ES_BASIC_CODE(rho,T))
 
#define KAPPA   0.
 SYNCH STUFF. More...
 
#define KAPPAES   0.
 
#define KAPPAUSER(rho, B, Tg, Tr)   (rho*KAPPA*KAPPA_FF_CODE(rho,Tg,Tr))
 
#define KAPPAESUSER(rho, T)   (rho*KAPPAES*KAPPA_ES_BASIC_CODE(rho,T))
 
#define KAPPA   0.
 SYNCH STUFF. More...
 
#define KAPPAES   1.
 
#define KAPPAUSER(rho, B, Tg, Tr)   (rho*KAPPA*KAPPA_FF_CODE(rho,Tg,Tr))
 
#define KAPPAESUSER(rho, T)   (rho*KAPPAES*KAPPA_ES_BASIC_CODE(rho,T))
 
#define KAPPA   0.
 SYNCH STUFF. More...
 
#define KAPPAES   0.
 
#define KAPPAUSER(rho, B, Tg, Tr)   (rho*KAPPA*KAPPA_FF_CODE(rho,Tg,Tr))
 
#define KAPPAESUSER(rho, T)   (rho*KAPPAES*KAPPA_ES_BASIC_CODE(rho,T))
 
#define KAPPAUSER(rho, B, Tg, Tr)   (rho*RADWAVE_KAPPA)
 
#define KAPPAESUSER(rho, T)   (rho*RADWAVE_KAPPAES)
 
#define KAPPAUSER(rho, B, Tg, Tr)   (0.)
 
#define KAPPAESUSER(rho, T)   (0.)
 
#define KAPPA   1.0
 SYNCH STUFF. More...
 
#define KAPPAES   1.0
 
#define KAPPAUSER(rho, B, Tg, Tr)   (rho*KAPPA*KAPPA_FF_CODE(rho,Tg,Tr))
 
#define KAPPAESUSER(rho, T)   (rho*KAPPAES*KAPPA_ES_BASIC_CODE(rho,T))
 
#define KAPPA   0.
 SYNCH STUFF. More...
 
#define KAPPAES   0.
 
#define KAPPAUSER(rho, B, Tg, Tr)   (rho*KAPPA*KAPPA_FF_CODE(rho,Tg,Tr))
 
#define KAPPAESUSER(rho, T)   (rho*KAPPAES*KAPPA_ES_BASIC_CODE(rho,T))
 
#define KAPPAUSER(rho, B, Tg, Tr)   (rho*KAPPA_ES_CODE(rho,Tg)/1E14*0.1)
 
#define KAPPAESUSER(rho, T)   (0.0)
 
#define KAPPA   1.0
 SYNCH STUFF. More...
 
#define KAPPAES   1.0
 
#define KAPPAUSER(rho, B, Tg, Tr)   (rho*KAPPA*(KAPPA_GENFF_CODE(SMALL+rho,Tg+TEMPMIN,Tr+TEMPMIN)))
 
#define KAPPANUSER(rho, B, Tg, Tr)   (rho*KAPPA*(KAPPA_GENFF_CODE(SMALL+rho,Tg+TEMPMIN,Tr+TEMPMIN)))
 
#define KAPPAESUSER(rho, T)   (rho*KAPPAES*KAPPA_ES_CODE(rho,T))
 
#define KAPPAUSER(rho, B, Tg, Tr)   (rho*KAPPA_ES_BASIC_CODE(rho,Tg)/1E14*0.0)
 
#define KAPPAESUSER(rho, T)   (0.0)
 
#define KAPPAUSER(rho, B, Tg, Tr)   (rho*(KAPPA_FF_CODE(SMALL+rho,Tg+TEMPMIN,Tr+TEMPMIN)))
 
#define KAPPAESUSER(rho, Tg)   (rho*KAPPA_ES_BASIC_CODE(rho,Tg)/100.0)
 
#define ISKAPPAEABS   0
 
#define ISKAPPANABS   1
 
#define ISKAPPAEEMIT   2
 
#define ISKAPPANEMIT   3
 
#define ISKAPPAES   4
 
#define E   2.718281828459045
 
#define JET6LIKEUSERCOORD   0
 
#define UNIHALFUSERCOORD   1
 
#define ORIGWALD   2
 
#define WHICHUSERCOORD   JET6LIKEUSERCOORD
 
#define cr(x)   (exp(-1.0/(x)))
 
#define tr(x)   (cr(x)/(cr(x) + cr(1.0-(x))))
 
#define trans(x, L, R)   ((x)<=(L) ? 0.0 : ((x)>=(R) ? 1.0 : tr(((x)-(L))/((R)-(L)))) )
 
#define transR(x, center, width)   ( 0.5*(1.0-tanh(+((x)-(center))/(width))))
 
#define transL(x, center, width)   ( 0.5*(1.0-tanh(-((x)-(center))/(width))))
 
#define transM(x, center, width)   ( exp(-pow(((x)-(center))/((width)*0.5),2.0) ) )
 
#define plateau(x, L, R, W)   (trans(x,(L)-0.5*(W),(L)+0.5*(W))*(1.0-trans(x,(R)-0.5*(W),(R)+0.5*(W))))
 
#define line1r(x, w)   (Rout)
 
#define line2r(x, w)   (Routeq)
 
#define line3r(x, w)   (Rout))
 
#define thetajon(x, w, xp1, xp2)   (line1r(x,w)*(1.0-trans(x,xp1,xp2)) + line2r(x,w)*trans(x,xp1,xp2))
 
#define lineeqr(x, w)   (R0 + exp(npow*(X[1]-mysx1)*0.70 + npow*mysx1 ) )
 
#define linepoler(x, w)   (R0 + exp(npow*X[1]))
 
#define thetaLr(x, wp, weq, xp1, xp2)   ( linepoler(x,wp)*(1.0-trans(x,xp1,xp2)) + lineeqr(x,weq)*trans(x,xp1,xp2) )
 
#define thetajon2(x, wp, weq, xp1, xp2)   ( x<0.5 ? thetaLr(x,wp,weq,xp1,xp2) : thetaLr(1.0-x,wp,weq,xp1,xp2) )
 
#define line1(x, w)   ((x)*(w))
 
#define line2(x, w)   ((x)*(w)+M_PI-(w))
 
#define line3(x, w)   ((x)*(w))
 
#define wparsam(x, r)   (h0 + pow(0.15 + ((r)-0.0)/10.0 , -njet))
 
#define thetasam(x, r, w, xp1, xp2)   (line1(x,w)*(1.0-trans(x,xp1,xp2)) + line2(x,w)*trans(x,xp1,xp2))
 
#define lineeq(x, w)   ((x)*(w)+(0.5*M_PI)-(0.5*w))
 
#define linepole(x, w)   (line1(x,w))
 
#define thetaL(x, wp, weq, xp1, xp2)   ( linepole(x,wp)*(1.0-trans(x,xp1,xp2)) + lineeq(x,weq)*trans(x,xp1,xp2) )
 
#define thetasam2(x, wp, weq, xp1, xp2)   ( x<0.5 ? thetaL(x,wp,weq,xp1,xp2) : -thetaL(1.0-x,wp,weq,xp1,xp2)+M_PI )
 
#define MAXIHOR   10
 
#define FRACN1   (0.1)
 
#define ADJUSTFRACT   (0.25)
 

Functions

static FTYPE nz_func (FTYPE R)
 
static FTYPE taper_func2 (FTYPE R, FTYPE rin, FTYPE rpow)
 
static int fieldprim (int whichmethod, int whichinversion, int *whichvel, int *whichcoord, int ii, int jj, int kk, FTYPE *pr)
 
static int get_full_rtsolution (int *whichvel, int *whichcoord, int opticallythick, FTYPE *pp, FTYPE *X, FTYPE *V, struct of_geom **ptrptrgeom)
 
static int make_nonrt2rt_solution (int *whichvel, int *whichcoord, int opticallythick, FTYPE *pp, FTYPE *X, FTYPE *V, struct of_geom **ptrptrgeom)
 
static int donut_analytical_solution (int *whichvel, int *whichcoord, int opticallythick, FTYPE *pp, FTYPE *X, FTYPE *V, struct of_geom **ptrptrgeom, FTYPE *ptptr)
 
static int process_solution (int *whichvel, int *whichcoord, int opticallythick, FTYPE *pp, FTYPE *X, FTYPE *V, struct of_geom **ptrptrgeom, FTYPE *ptptr)
 
int prepre_init_specific_init (void)
 
int pre_init_specific_init (void)
 
int init_conservatives (FTYPE(*prim)[NSTORE2][NSTORE3][NPR], FTYPE(*pstag)[NSTORE2][NSTORE3][NPR], FTYPE(*Utemp)[NSTORE2][NSTORE3][NPR], FTYPE(*U)[NSTORE2][NSTORE3][NPR])
 
int post_init_specific_init (void)
 
int init_consts (void)
 
int init_global (void)
 
int init_defcoord (void)
 
int init_grid (void)
 
int init_grid_post_set_grid (FTYPE(*prim)[NSTORE2][NSTORE3][NPR], FTYPE(*pstag)[NSTORE2][NSTORE3][NPR], FTYPE(*ucons)[NSTORE2][NSTORE3][NPR], FTYPE(*vpot)[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3], FTYPE(*Bhat)[NSTORE2][NSTORE3][NPR], FTYPE(*panalytic)[NSTORE2][NSTORE3][NPR], FTYPE(*pstaganalytic)[NSTORE2][NSTORE3][NPR], FTYPE(*vpotanalytic)[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3], FTYPE(*Bhatanalytic)[NSTORE2][NSTORE3][NPR], FTYPE(*F1)[NSTORE2][NSTORE3][NPR+NSPECIAL], FTYPE(*F2)[NSTORE2][NSTORE3][NPR+NSPECIAL], FTYPE(*F3)[NSTORE2][NSTORE3][NPR+NSPECIAL], FTYPE(*Atemp)[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3])
 
int init_primitives (FTYPE(*prim)[NSTORE2][NSTORE3][NPR], FTYPE(*pstag)[NSTORE2][NSTORE3][NPR], FTYPE(*ucons)[NSTORE2][NSTORE3][NPR], FTYPE(*vpot)[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3], FTYPE(*Bhat)[NSTORE2][NSTORE3][NPR], FTYPE(*panalytic)[NSTORE2][NSTORE3][NPR], FTYPE(*pstaganalytic)[NSTORE2][NSTORE3][NPR], FTYPE(*vpotanalytic)[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3], FTYPE(*Bhatanalytic)[NSTORE2][NSTORE3][NPR], FTYPE(*F1)[NSTORE2][NSTORE3][NPR+NSPECIAL], FTYPE(*F2)[NSTORE2][NSTORE3][NPR+NSPECIAL], FTYPE(*F3)[NSTORE2][NSTORE3][NPR+NSPECIAL], FTYPE(*Atemp)[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3])
 
int init_dsandvels (int inittype, int pos, int *whichvel, int *whichcoord, SFTYPE time, int i, int j, int k, FTYPE *pr, FTYPE *pstag)
 
int init_dsandvels_koral (int *whichvel, int *whichcoord, int i, int j, int k, FTYPE *pr, FTYPE *pstag)
 
int process_restart_toget_radiation (void)
 
static FTYPE taper_func_exp (FTYPE R, FTYPE rin, FTYPE POTENTIALorPRESSURE)
 
FTYPE setgpara (FTYPE myr, FTYPE th, FTYPE thpower)
 
FTYPE setblandfordfield (FTYPE r, FTYPE th)
 
int init_vpot_user (int *whichcoord, int l, SFTYPE time, int i, int j, int k, int loc, FTYPE(*prim)[NSTORE2][NSTORE3][NPR], FTYPE *V, FTYPE *A)
 
int init_postvpot (int i, int j, int k, FTYPE *pr, FTYPE *pstag, FTYPE *ucons)
 
void shortout_Bd3 (FTYPE(*prim)[NSTORE2][NSTORE3][NPR])
 
void Fcov_numerical (int whichcoord, FTYPE *X, FTYPE(*Fcov)[NDIM])
 
FTYPE mcov_func_mcoord (struct of_geom *ptrgeom, FTYPE *X, int ii, int jj)
 
FTYPE kcov_func_mcoord (struct of_geom *ptrgeom, FTYPE *X, int ii, int jj)
 
void Jcon_numerical (int whichcoord, FTYPE *X, FTYPE *Jcon)
 
FTYPE Fcon_func_mcoord (struct of_geom *ptrgeom, FTYPE *X, int ii, int jj)
 
int init_vpot2field_user (SFTYPE time, FTYPE(*A)[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3], FTYPE(*prim)[NSTORE2][NSTORE3][NPR], FTYPE(*pstag)[NSTORE2][NSTORE3][NPR], FTYPE(*ucons)[NSTORE2][NSTORE3][NPR], FTYPE(*Bhat)[NSTORE2][NSTORE3][NPR])
 
int normalize_densities (FTYPE(*prim)[NSTORE2][NSTORE3][NPR])
 
int getmax_densities_full (FTYPE(*prim)[NSTORE2][NSTORE3][NPR], SFTYPE *rhomax, SFTYPE *umax, SFTYPE *uradmax, SFTYPE *utotmax, SFTYPE *pmax, SFTYPE *pradmax, SFTYPE *ptotmax)
 
int get_maxes (FTYPE(*prim)[NSTORE2][NSTORE3][NPR], FTYPE *bsq_max, FTYPE *ptot_max, FTYPE *beta_min)
 
int normalize_field (FTYPE(*prim)[NSTORE2][NSTORE3][NPR], FTYPE(*pstag)[NSTORE2][NSTORE3][NPR], FTYPE(*ucons)[NSTORE2][NSTORE3][NPR], FTYPE(*vpot)[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3], FTYPE(*Bhat)[NSTORE2][NSTORE3][NPR])
 
int jetbound (int i, int j, int k, int loc, FTYPE *prin, FTYPE *prflux, FTYPE(*prim)[NSTORE2][NSTORE3][NPR])
 
static FTYPE kappa_func_fits (int which, FTYPE rho, FTYPE B, FTYPE Tg, FTYPE Tr, FTYPE varexpf)
 
int kappa_func_fits_all (FTYPE rho, FTYPE B, FTYPE Tg, FTYPE Tr, FTYPE varexpf, FTYPE *kappa, FTYPE *kappaemit, FTYPE *kappan, FTYPE *kappanemit, FTYPE *kappaes)
 
FTYPE Gcompt (FTYPE rho0, FTYPE Tgas, FTYPE Tradff, FTYPE Ruu)
 
int coolfunc_user (FTYPE h_over_r, FTYPE *pr, struct of_geom *geom, struct of_state *q, FTYPE(*dUcomp)[NPR])
 
void set_coord_parms_nodeps_user (int defcoordlocal)
 
void set_coord_parms_deps_user (int defcoordlocal)
 
void write_coord_parms_user (int defcoordlocal, FILE *out)
 
void read_coord_parms_user (int defcoordlocal, FILE *in)
 
void blcoord_user (FTYPE *X, FTYPE *V)
 
void dxdxp_analytic_user (FTYPE *X, FTYPE *V, FTYPE(*dxdxp)[NDIM])
 
FTYPE setRin_user (int ihor, FTYPE ihoradjust)
 
int setihor_user (void)
 

Variables

static SFTYPE rhomax =0
 
static SFTYPE umax =0
 
static SFTYPE uradmax =0
 
static SFTYPE utotmax =0
 
static SFTYPE pmax =0
 
static SFTYPE pradmax =0
 
static SFTYPE ptotmax =0
 
static SFTYPE bsq_max =0
 
static SFTYPE beta
 
static SFTYPE randfact
 
static SFTYPE rin
 
static SFTYPE rinfield
 
static SFTYPE routfield
 
static FTYPE rhodisk
 
FTYPE B0WALD
 
int DOWALDDEN =0
 
int WALDWHICHACOV =3
 
FTYPE BSQORHOWALD =100.0
 
FTYPE aforwald
 
FTYPE TILTWALD
 
FTYPE thindiskrhopow =-0.6
 
FTYPE normglobal
 
int inittypeglobal
 
int RADBEAM2D_BEAMNO
 
int RADBEAM2D_FLATBACKGROUND
 
FTYPE RADBEAM2D_RHOAMB
 
FTYPE RADBEAM2D_TAMB
 
int RADBEAM2D_BLOB
 
FTYPE RADBEAM2D_BLOBW
 
FTYPE RADBEAM2D_BLOBP
 
FTYPE RADBEAM2D_BLOBX
 
FTYPE RADBEAM2D_BLOBZ
 
FTYPE RADBEAM2D_PAR_D
 
FTYPE RADBEAM2D_PAR_E
 
int RADBEAM2D_IFBEAM
 
FTYPE RADBEAM2D_TLEFT
 
FTYPE RADBEAM2D_NLEFT
 
FTYPE RADBEAM2D_BEAML
 
FTYPE RADBEAM2D_BEAMR
 
int RADBEAM2DKSVERT_BEAMNO
 
FTYPE RADBEAMFLAT_FRATIO
 
FTYPE RADBEAMFLAT_ERAD
 
FTYPE RADBEAMFLAT_RHO
 
FTYPE RADBEAMFLAT_UU
 
FTYPE RADATM_MDOTEDD
 
FTYPE RADATM_LUMEDD
 
int RADATM_THINRADATM
 
FTYPE RADATM_FERATIO
 
FTYPE RADATM_FRATIO
 
FTYPE RADATM_RHOAMB
 
FTYPE RADATM_TAMB
 
FTYPE RADSHADOW_NLEFT
 
FTYPE RADSHADOW_ANGLE
 
FTYPE RADSHADOW_TLEFTOTAMB
 
FTYPE RADSHADOW_BEAMY
 
FTYPE RADDBLSHADOW_NLEFT
 
FTYPE RADDBLSHADOW_ANGLE
 
FTYPE RADDBLSHADOW_TLEFTOTAMB
 
FTYPE RADDBLSHADOW_BEAMY
 
int RADWAVE_NWAVE
 
int RADWAVE_NUMERO
 
int RADWAVE_WAVETYPE
 
FTYPE RADWAVE_PP
 
FTYPE RADWAVE_CC
 
FTYPE RADWAVE_KAPPA
 
FTYPE RADWAVE_RHOFAC
 
FTYPE RADWAVE_B0
 
FTYPE RADWAVE_DRRE
 
FTYPE RADWAVE_DRIM
 
FTYPE RADWAVE_DVRE
 
FTYPE RADWAVE_DVIM
 
FTYPE RADWAVE_DV2RE
 
FTYPE RADWAVE_DV2IM
 
FTYPE RADWAVE_DB2RE
 
FTYPE RADWAVE_DB2IM
 
FTYPE RADWAVE_DURE
 
FTYPE RADWAVE_DUIM
 
FTYPE RADWAVE_DERE
 
FTYPE RADWAVE_DEIM
 
FTYPE RADWAVE_DFRE
 
FTYPE RADWAVE_DFIM
 
FTYPE RADWAVE_DF2RE
 
FTYPE RADWAVE_DF2IM
 
FTYPE RADWAVE_OMRE
 
FTYPE RADWAVE_OMIM
 
FTYPE RADWAVE_DTOUT1
 
FTYPE RADWAVE_RHOZERO
 
FTYPE RADWAVE_AAA
 
FTYPE RADWAVE_KK
 
FTYPE RADWAVE_UINT
 
FTYPE RADWAVE_TEMP
 
FTYPE RADWAVE_ERAD
 
FTYPE RADWAVE_VX
 
FTYPE RADWAVE_KAPPAES
 
FTYPE RADWAVE_ERADFACTOR
 
FTYPE RADWAVE_GASFACTOR
 
FTYPE RADBONDI_TESTNO
 
FTYPE RADBONDI_PRADGAS
 
FTYPE RADBONDI_TGAS0
 
FTYPE RADBONDI_MDOTPEREDD
 
FTYPE RADBONDI_MDOTEDD
 
FTYPE RADBONDI_MINX
 
FTYPE RADBONDI_MAXX
 
FTYPE RADDOT_XDOT
 
FTYPE RADDOT_YDOT
 
FTYPE RADDOT_ZDOT
 
int RADDOT_IDOT
 
int RADDOT_JDOT
 
int RADDOT_KDOT
 
FTYPE RADDOT_FYDOT
 
FTYPE RADDOT_LTEFACTOR
 
FTYPE RADDOT_URFX
 
FTYPE RADDOT_F1
 
FTYPE RADDOT_F2
 
FTYPE RADNT_MINX
 
FTYPE RADNT_MAXX
 
FTYPE RADNT_KKK
 
FTYPE RADNT_ELL
 
FTYPE RADNT_UTPOT
 
FTYPE RADNT_RHOATMMIN
 
FTYPE RADNT_RHODONUT
 
FTYPE RADNT_UINTATMMIN
 
FTYPE RADNT_ERADATMMIN
 
FTYPE RADNT_DONUTTYPE
 
FTYPE RADNT_INFLOWING
 
FTYPE RADNT_TGASATMMIN
 
FTYPE RADNT_TRADATMMIN
 
FTYPE RADNT_ROUT
 
FTYPE RADNT_OMSCALE
 
FTYPE RADNT_FULLPHI
 
FTYPE RADNT_DONUTRADPMAX
 
FTYPE RADNT_HOVERR
 
FTYPE RADNT_LPOW
 
int RADCYLJET_TYPE
 
FTYPE RADCYLJET_VRSTAR
 
FTYPE RADCYLJET_EHATJET
 
FTYPE RADCYLJET_RHOJET
 
FTYPE RADCYLJET_TEMPJET
 
FTYPE RADCYLJET_UJET
 
int RADDONUT_OPTICALLYTHICKTORUS
 
static FTYPE npow
 
static FTYPE r1jet
 
static FTYPE njet1
 
static FTYPE njet
 
static FTYPE r0jet
 
static FTYPE rsjet
 
static FTYPE Qjet
 
static FTYPE ntheta
 
static FTYPE htheta
 
static FTYPE rsjet2
 
static FTYPE r0jet2
 
static FTYPE rsjet3
 
static FTYPE r0jet3
 
static FTYPE rs
 
static FTYPE r0
 
static FTYPE npow2
 
static FTYPE cpow2
 
static FTYPE rbr
 
static FTYPE x1br
 
static FTYPE h0
 
static FTYPE cpow3
 

Detailed Description

USER initial conditions (usually links to code in initbasecode directory)

Definition in file init.c.

Macro Definition Documentation

#define ADJUSTFRACT   (0.25)

Definition at line 10819 of file init.c.

#define BLANDFORDQUAD   5

Definition at line 76 of file init.c.

#define cr (   x)    (exp(-1.0/(x)))
#define DIMTYPELIST   "%21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g %21.15g\n"
#define DISK1FIELD   0

Definition at line 71 of file init.c.

#define DISK1VERT   3

Definition at line 74 of file init.c.

#define DISK2FIELD   1

Definition at line 72 of file init.c.

#define DISK2VERT   4

Definition at line 75 of file init.c.

#define E   2.718281828459045

Definition at line 9611 of file init.c.

#define FCOVDERTYPE   DIFFGAMMIE

Definition at line 8112 of file init.c.

#define FCOVDXDELTA   1E-5

Definition at line 8116 of file init.c.

#define FIELDJONMAD   10

Definition at line 81 of file init.c.

#define FIELDWALD   11

Definition at line 82 of file init.c.

#define FRACN1   (0.1)

Definition at line 10818 of file init.c.

#define fsx (   x)    ((x)>0 ? exp(-1.0/(x)) : 0.0)
#define gsx (   x)    (fsx(x)/(fsx(x) + fsx(1.0-(x))))
#define ISKAPPAEABS   0

Definition at line 9605 of file init.c.

#define ISKAPPAEEMIT   2

Definition at line 9607 of file init.c.

#define ISKAPPAES   4

Definition at line 9609 of file init.c.

#define ISKAPPANABS   1

Definition at line 9606 of file init.c.

#define ISKAPPANEMIT   3

Definition at line 9608 of file init.c.

#define JET6LIKEUSERCOORD   0

Definition at line 10199 of file init.c.

#define JONMADHPOW   (4.0)
#define JONMADR0   (0.0)
#define JONMADROUT   (300.0)
#define KAPPA   0.

SYNCH STUFF.

Synchrotron energy-opacity Synchrotron number-opacity

Definition at line 9559 of file init.c.

#define KAPPA   0.

SYNCH STUFF.

Synchrotron energy-opacity Synchrotron number-opacity

Definition at line 9559 of file init.c.

#define KAPPA   0.

SYNCH STUFF.

Synchrotron energy-opacity Synchrotron number-opacity

Definition at line 9559 of file init.c.

#define KAPPA   0.

SYNCH STUFF.

Synchrotron energy-opacity Synchrotron number-opacity

Definition at line 9559 of file init.c.

#define KAPPA   0.

SYNCH STUFF.

Synchrotron energy-opacity Synchrotron number-opacity

Definition at line 9559 of file init.c.

#define KAPPA   0.

SYNCH STUFF.

Synchrotron energy-opacity Synchrotron number-opacity

Definition at line 9559 of file init.c.

#define KAPPA   0.

SYNCH STUFF.

Synchrotron energy-opacity Synchrotron number-opacity

Definition at line 9559 of file init.c.

#define KAPPA   1.0

SYNCH STUFF.

Synchrotron energy-opacity Synchrotron number-opacity

Definition at line 9559 of file init.c.

#define KAPPA   0.

SYNCH STUFF.

Synchrotron energy-opacity Synchrotron number-opacity

Definition at line 9559 of file init.c.

#define KAPPA   1.0

SYNCH STUFF.

Synchrotron energy-opacity Synchrotron number-opacity

Definition at line 9559 of file init.c.

#define KAPPA_BF_CODE (   rhocode,
  Tgcode,
  Trcode 
)    (3.0E25*(ZFACT)*(1.0+XFACT+0.75*YFACT)*((rhocode)*RHOBAR)*prpow((Tgcode)*TEMPBAR,-0.5)*prpow((Trcode)*TEMPBAR,-3.0)*prlog(1.0+1.6*KAPPA_ZETA(Tgcode,Trcode))/OPACITYBAR)

Definition at line 9258 of file init.c.

#define KAPPA_CHIANTIBF_CODE (   rhocode,
  Tgcode,
  Trcode 
)    (4.0E34*((rhocode*RHOBAR))*(ZFACT/ZSOLAR)*YELE*prpow((Tgcode)*TEMPBAR,-1.7)*prpow((Trcode)*TEMPBAR,-3.0)/OPACITYBAR)

Definition at line 9259 of file init.c.

#define KAPPA_ES_BASIC_CODE (   rhocode,
  Tcode 
)    (0.2*(1.0+XFACT)/OPACITYBAR)

Definition at line 9210 of file init.c.

#define KAPPA_ES_CODE (   rhocode,
  Tcode 
)    (0.2*(1.0+XFACT)*KAPPA_ES_FERMICORR(rhocode,Tcode)*KAPPA_ES_KNCORR(rhocode,Tcode)/OPACITYBAR)

kappaes = sigma_T n_e = sigma_T n_b (n_e/n_b) = sigma_T rho/mb (ne/nb)

Definition at line 9209 of file init.c.

#define KAPPA_ES_FERMICORR (   rhocode,
  Tcode 
)    (1.0/(1.0+2.7E11*((rhocode)*RHOBAR)/prpow((Tcode)*TEMPBAR,2.0)))

PURE ELASTIC SCATTERING.

Definition at line 9206 of file init.c.

#define KAPPA_ES_KNCORR (   rhocode,
  Tcode 
)    (1.0/(1.0+prpow((Tcode)*TEMPBAR/4.5E8,0.86)))

Definition at line 9207 of file init.c.

#define KAPPA_FF_CODE (   rhocode,
  Tgcode,
  Trcode 
)    (4.0E22*(1.0+XFACT)*(1.0-ZFACT)*((rhocode)*RHOBAR)*prpow((Tgcode)*TEMPBAR,-0.5)*prpow((Trcode)*TEMPBAR,-3.0)*prlog(1.0+1.6*KAPPA_ZETA(Tgcode,Trcode))*(1.0+4.4E-10*(Tgcode*TEMPBAR))/OPACITYBAR)

Definition at line 9236 of file init.c.

#define KAPPA_FORCOMPT_CODE (   rhocode,
  Tcode 
)    (0.2*(1.0+XFACT)*KAPPA_ES_FERMICORR(rhocode,Tcode)*KAPPA_FORCOMPT_RELCORR(rhocode,Tcode)/OPACITYBAR)

Definition at line 9222 of file init.c.

#define KAPPA_FORCOMPT_RELCORR (   rhocode,
  Tcode 
)    (KAPPA_FORCOMPT_RELCORREP(K_BOLTZ*(Tcode)*TEMPBAR/(MELE*CCCTRUE*CCCTRUE)))

Definition at line 9221 of file init.c.

#define KAPPA_FORCOMPT_RELCORREP (   ep)    ((1.0 + 3.683*(ep)+4.0*(ep)*(ep))/(1.0 + (ep)))

Definition at line 9220 of file init.c.

#define KAPPA_GENFF_CODE (   rhocode,
  Tgcode,
  Trcode 
)    (1.0/(1.0/(KAPPA_MOL_CODE(rhocode,Tgcode,Trcode)+KAPPA_HN_CODE(rhocode,Tgcode,Trcode)) + 1.0/(KAPPA_CHIANTIBF_CODE(rhocode,Tgcode,Trcode)+KAPPA_BF_CODE(rhocode,Tgcode,Trcode)+KAPPA_FF_CODE(rhocode,Tgcode,Trcode))))

Definition at line 9263 of file init.c.

#define KAPPA_HN_CODE (   rhocode,
  Tgcode,
  Trcode 
)    (1.1E-25*prpow(ZFACT,0.5)*prpow((rhocode)*RHOBAR,0.5)*prpow((Tgcode)*TEMPBAR,7.7)/OPACITYBAR)

Definition at line 9260 of file init.c.

#define KAPPA_MOL_CODE (   rhocode,
  Tgcode,
  Trcode 
)    (0.1*ZFACT/OPACITYBAR)

Definition at line 9261 of file init.c.

#define KAPPA_ZETA (   Tgcode,
  Trcode 
)    ((TEMPMIN+Trcode)/(TEMPMIN+Tgcode))

EMISSION (Tr=Tg) or ABSORBPTION (Tr different from Tg)

Definition at line 9234 of file init.c.

#define KAPPAES   0.

Definition at line 9560 of file init.c.

#define KAPPAES   (1E3)

Definition at line 9560 of file init.c.

#define KAPPAES   0.

Definition at line 9560 of file init.c.

#define KAPPAES   0.

Definition at line 9560 of file init.c.

#define KAPPAES   0.

Definition at line 9560 of file init.c.

#define KAPPAES   1.

Definition at line 9560 of file init.c.

#define KAPPAES   0.

Definition at line 9560 of file init.c.

#define KAPPAES   1.0

Definition at line 9560 of file init.c.

#define KAPPAES   0.

Definition at line 9560 of file init.c.

#define KAPPAES   1.0

Definition at line 9560 of file init.c.

#define KAPPAESUSER (   rho,
 
)    (rho*KAPPAES*KAPPA_ES_BASIC_CODE(rho,T))

Definition at line 9588 of file init.c.

#define KAPPAESUSER (   rho,
 
)    (rho*KAPPAES)

Definition at line 9588 of file init.c.

#define KAPPAESUSER (   rho,
 
)    (rho*KAPPAES*KAPPA_ES_BASIC_CODE(rho,T))

Definition at line 9588 of file init.c.

#define KAPPAESUSER (   rho,
 
)    (0.0)

Definition at line 9588 of file init.c.

#define KAPPAESUSER (   rho,
 
)    (rho*0.0)

Definition at line 9588 of file init.c.

#define KAPPAESUSER (   rho,
 
)    (rho*KAPPAES*KAPPA_ES_BASIC_CODE(rho,T))

Definition at line 9588 of file init.c.

#define KAPPAESUSER (   rho,
 
)    (rho*KAPPAES*KAPPA_ES_BASIC_CODE(rho,T))

Definition at line 9588 of file init.c.

#define KAPPAESUSER (   rho,
 
)    (rho*KAPPAES*KAPPA_ES_BASIC_CODE(rho,T))

Definition at line 9588 of file init.c.

#define KAPPAESUSER (   rho,
 
)    (rho*KAPPAES*KAPPA_ES_BASIC_CODE(rho,T))

Definition at line 9588 of file init.c.

#define KAPPAESUSER (   rho,
 
)    (rho*RADWAVE_KAPPAES)

Definition at line 9588 of file init.c.

#define KAPPAESUSER (   rho,
 
)    (0.)

Definition at line 9588 of file init.c.

#define KAPPAESUSER (   rho,
 
)    (rho*KAPPAES*KAPPA_ES_BASIC_CODE(rho,T))

Definition at line 9588 of file init.c.

#define KAPPAESUSER (   rho,
 
)    (rho*KAPPAES*KAPPA_ES_BASIC_CODE(rho,T))

Definition at line 9588 of file init.c.

#define KAPPAESUSER (   rho,
 
)    (0.0)

Definition at line 9588 of file init.c.

#define KAPPAESUSER (   rho,
 
)    (rho*KAPPAES*KAPPA_ES_CODE(rho,T))

Definition at line 9588 of file init.c.

#define KAPPAESUSER (   rho,
 
)    (0.0)

Definition at line 9588 of file init.c.

#define KAPPAESUSER (   rho,
  Tg 
)    (rho*KAPPA_ES_BASIC_CODE(rho,Tg)/100.0)

Definition at line 9588 of file init.c.

#define KAPPAN_FF_CODE (   rhocode,
  Tgcode,
  Trcode 
)    KAPPA_FF_CODE(rhocode,Tgcode,Trcode)

Definition at line 9254 of file init.c.

#define KAPPANUSER (   rho,
  B,
  Tg,
  Tr 
)    (rho*KAPPA*(KAPPA_GENFF_CODE(SMALL+rho,Tg+TEMPMIN,Tr+TEMPMIN)))

Definition at line 9570 of file init.c.

#define KAPPAUSER (   rho,
  B,
  Tg,
  Tr 
)    (rho*KAPPA*KAPPA_FF_CODE(rho,Tg,Tr))

Definition at line 9587 of file init.c.

#define KAPPAUSER (   rho,
  B,
  Tg,
  Tr 
)    (rho*KAPPA)

Definition at line 9587 of file init.c.

#define KAPPAUSER (   rho,
  B,
  Tg,
  Tr 
)    (rho*KAPPA*KAPPA_FF_CODE(rho,Tg,Tr))

Definition at line 9587 of file init.c.

#define KAPPAUSER (   rho,
  B,
  Tg,
  Tr 
)    (25*rho)

Definition at line 9587 of file init.c.

#define KAPPAUSER (   rho,
  B,
  Tg,
  Tr 
)    (rho*1.0)

Definition at line 9587 of file init.c.

#define KAPPAUSER (   rho,
  B,
  Tg,
  Tr 
)    (rho*KAPPA*KAPPA_FF_CODE(rho,Tg,Tr))

Definition at line 9587 of file init.c.

#define KAPPAUSER (   rho,
  B,
  Tg,
  Tr 
)    (rho*KAPPA*KAPPA_FF_CODE(rho,Tg,Tr))

Definition at line 9587 of file init.c.

#define KAPPAUSER (   rho,
  B,
  Tg,
  Tr 
)    (rho*KAPPA*KAPPA_FF_CODE(rho,Tg,Tr))

Definition at line 9587 of file init.c.

#define KAPPAUSER (   rho,
  B,
  Tg,
  Tr 
)    (rho*KAPPA*KAPPA_FF_CODE(rho,Tg,Tr))

Definition at line 9587 of file init.c.

#define KAPPAUSER (   rho,
  B,
  Tg,
  Tr 
)    (rho*RADWAVE_KAPPA)

Definition at line 9587 of file init.c.

#define KAPPAUSER (   rho,
  B,
  Tg,
  Tr 
)    (0.)

Definition at line 9587 of file init.c.

#define KAPPAUSER (   rho,
  B,
  Tg,
  Tr 
)    (rho*KAPPA*KAPPA_FF_CODE(rho,Tg,Tr))

Definition at line 9587 of file init.c.

#define KAPPAUSER (   rho,
  B,
  Tg,
  Tr 
)    (rho*KAPPA*KAPPA_FF_CODE(rho,Tg,Tr))

Definition at line 9587 of file init.c.

#define KAPPAUSER (   rho,
  B,
  Tg,
  Tr 
)    (rho*KAPPA_ES_CODE(rho,Tg)/1E14*0.1)

Definition at line 9587 of file init.c.

#define KAPPAUSER (   rho,
  B,
  Tg,
  Tr 
)    (rho*KAPPA*(KAPPA_GENFF_CODE(SMALL+rho,Tg+TEMPMIN,Tr+TEMPMIN)))

Definition at line 9587 of file init.c.

#define KAPPAUSER (   rho,
  B,
  Tg,
  Tr 
)    (rho*KAPPA_ES_BASIC_CODE(rho,Tg)/1E14*0.0)

Definition at line 9587 of file init.c.

#define KAPPAUSER (   rho,
  B,
  Tg,
  Tr 
)    (rho*(KAPPA_FF_CODE(SMALL+rho,Tg+TEMPMIN,Tr+TEMPMIN)))

Definition at line 9587 of file init.c.

#define line1 (   x,
 
)    ((x)*(w))
#define line1r (   x,
 
)    (Rout)
#define line2 (   x,
 
)    ((x)*(w)+M_PI-(w))
#define line2r (   x,
 
)    (Routeq)
#define line3 (   x,
 
)    ((x)*(w))
#define line3r (   x,
 
)    (Rout))
#define lineeq (   x,
 
)    ((x)*(w)+(0.5*M_PI)-(0.5*w))
#define lineeqr (   x,
 
)    (R0 + exp(npow*(X[1]-mysx1)*0.70 + npow*mysx1 ) )
#define linepole (   x,
 
)    (line1(x,w))
#define linepoler (   x,
 
)    (R0 + exp(npow*X[1]))
#define MAXIHOR   10

Definition at line 10817 of file init.c.

#define MAXPASSPARMS   10

Definition at line 59 of file init.c.

#define MONOPOLAR   8

Definition at line 79 of file init.c.

#define MONOPOLE   12

Definition at line 83 of file init.c.

#define NOFIELD   -1

Definition at line 70 of file init.c.

#define NTUBE   31
#define NUMTgasITERS   4
#define OHSUGAFIELD   7

Definition at line 78 of file init.c.

#define OLEKFIELD   9

Definition at line 80 of file init.c.

#define ORIGWALD   2

Definition at line 10201 of file init.c.

#define plateau (   x,
  L,
  R,
 
)    (trans(x,(L)-0.5*(W),(L)+0.5*(W))*(1.0-trans(x,(R)-0.5*(W),(R)+0.5*(W))))
#define POWERNU   (2.0)
#define QPOWER   (1.0)
#define SLOWFAC   1.0 /* reduce u_phi by this amount */

Definition at line 58 of file init.c.

#define SPLITMONOPOLE   13

Definition at line 84 of file init.c.

#define stepfunction (   x)    (1.0 - gsx((x)-0.5))
#define stepfunctionab (   x,
  a,
 
)    (((a)-(b))*stepfunction(x) + (b))
#define stepfunctionab2 (   x,
  a,
 
)    (((a)-(b))*stepfunction(x)*stepfunction(x) + (b))
#define TAUADJUSTATM   (10.0)

Definition at line 8565 of file init.c.

#define thetajon (   x,
  w,
  xp1,
  xp2 
)    (line1r(x,w)*(1.0-trans(x,xp1,xp2)) + line2r(x,w)*trans(x,xp1,xp2))
#define thetajon2 (   x,
  wp,
  weq,
  xp1,
  xp2 
)    ( x<0.5 ? thetaLr(x,wp,weq,xp1,xp2) : thetaLr(1.0-x,wp,weq,xp1,xp2) )
#define thetaL (   x,
  wp,
  weq,
  xp1,
  xp2 
)    ( linepole(x,wp)*(1.0-trans(x,xp1,xp2)) + lineeq(x,weq)*trans(x,xp1,xp2) )
#define thetaLr (   x,
  wp,
  weq,
  xp1,
  xp2 
)    ( linepoler(x,wp)*(1.0-trans(x,xp1,xp2)) + lineeqr(x,weq)*trans(x,xp1,xp2) )
#define thetasam (   x,
  r,
  w,
  xp1,
  xp2 
)    (line1(x,w)*(1.0-trans(x,xp1,xp2)) + line2(x,w)*trans(x,xp1,xp2))
#define thetasam2 (   x,
  wp,
  weq,
  xp1,
  xp2 
)    ( x<0.5 ? thetaL(x,wp,weq,xp1,xp2) : -thetaL(1.0-x,wp,weq,xp1,xp2)+M_PI )
#define TOROIDALFIELD   6

Definition at line 77 of file init.c.

#define tr (   x)    (cr(x)/(cr(x) + cr(1.0-(x))))
#define trans (   x,
  L,
  R 
)    ((x)<=(L) ? 0.0 : ((x)>=(R) ? 1.0 : tr(((x)-(L))/((R)-(L)))) )
#define transL (   x,
  center,
  width 
)    ( 0.5*(1.0-tanh(-((x)-(center))/(width))))
#define transM (   x,
  center,
  width 
)    ( exp(-pow(((x)-(center))/((width)*0.5),2.0) ) )
#define transR (   x,
  center,
  width 
)    ( 0.5*(1.0-tanh(+((x)-(center))/(width))))
#define UNIHALFUSERCOORD   1

Definition at line 10200 of file init.c.

#define USER_THETAROTMETRIC   (0.0)

Definition at line 65 of file init.c.

#define USER_THETAROTPRIMITIVES   (0.0)

Definition at line 66 of file init.c.

#define VERTFIELD   2

Definition at line 73 of file init.c.

#define WHICHRADATM   0
#define WHICHUSERCOORD   JET6LIKEUSERCOORD

Definition at line 10203 of file init.c.

#define wparsam (   x,
 
)    (h0 + pow(0.15 + ((r)-0.0)/10.0 , -njet))

Function Documentation

void blcoord_user ( FTYPE X,
FTYPE V 
)

Definition at line 10356 of file init.c.

Here is the caller graph for this function:

int coolfunc_user ( FTYPE  h_over_r,
FTYPE pr,
struct of_geom geom,
struct of_state q,
FTYPE(*)  dUcomp[NPR] 
)

Definition at line 10192 of file init.c.

static int donut_analytical_solution ( int *  whichvel,
int *  whichcoord,
int  opticallythick,
FTYPE pp,
FTYPE X,
FTYPE V,
struct of_geom **  ptrptrgeom,
FTYPE ptptr 
)
static

Definition at line 6356 of file init.c.

Here is the call graph for this function:

Here is the caller graph for this function:

void dxdxp_analytic_user ( FTYPE X,
FTYPE V,
FTYPE(*)  dxdxp[NDIM] 
)

Definition at line 10687 of file init.c.

FTYPE Fcon_func_mcoord ( struct of_geom ptrgeom,
FTYPE X,
int  ii,
int  jj 
)

Definition at line 8420 of file init.c.

Here is the call graph for this function:

Here is the caller graph for this function:

void Fcov_numerical ( int  whichcoord,
FTYPE X,
FTYPE(*)  Fcov[NDIM] 
)

Definition at line 8121 of file init.c.

Here is the call graph for this function:

Here is the caller graph for this function:

static int fieldprim ( int  whichmethod,
int  whichinversion,
int *  whichvel,
int *  whichcoord,
int  ii,
int  jj,
int  kk,
FTYPE pr 
)
static

Definition at line 7542 of file init.c.

Here is the call graph for this function:

Here is the caller graph for this function:

FTYPE Gcompt ( FTYPE  rho0,
FTYPE  Tgas,
FTYPE  Tradff,
FTYPE  Ruu 
)

Definition at line 10072 of file init.c.

Here is the caller graph for this function:

static int get_full_rtsolution ( int *  whichvel,
int *  whichcoord,
int  opticallythick,
FTYPE pp,
FTYPE X,
FTYPE V,
struct of_geom **  ptrptrgeom 
)
static

Definition at line 5975 of file init.c.

Here is the call graph for this function:

Here is the caller graph for this function:

int get_maxes ( FTYPE(*)  prim[NSTORE2][NSTORE3][NPR],
FTYPE bsq_max,
FTYPE ptot_max,
FTYPE beta_min 
)

Definition at line 8500 of file init.c.

Here is the call graph for this function:

Here is the caller graph for this function:

int getmax_densities_full ( FTYPE(*)  prim[NSTORE2][NSTORE3][NPR],
SFTYPE rhomax,
SFTYPE umax,
SFTYPE uradmax,
SFTYPE utotmax,
SFTYPE pmax,
SFTYPE pradmax,
SFTYPE ptotmax 
)

Definition at line 8488 of file init.c.

Here is the call graph for this function:

Here is the caller graph for this function:

int init_conservatives ( FTYPE(*)  prim[NSTORE2][NSTORE3][NPR],
FTYPE(*)  pstag[NSTORE2][NSTORE3][NPR],
FTYPE(*)  Utemp[NSTORE2][NSTORE3][NPR],
FTYPE(*)  U[NSTORE2][NSTORE3][NPR] 
)

Definition at line 390 of file init.c.

Here is the caller graph for this function:

int init_consts ( void  )

Definition at line 465 of file init.c.

Here is the caller graph for this function:

int init_defcoord ( void  )

Definition at line 2612 of file init.c.

Here is the caller graph for this function:

int init_dsandvels ( int  inittype,
int  pos,
int *  whichvel,
int *  whichcoord,
SFTYPE  time,
int  i,
int  j,
int  k,
FTYPE pr,
FTYPE pstag 
)

Definition at line 3669 of file init.c.

Here is the caller graph for this function:

int init_dsandvels_koral ( int *  whichvel,
int *  whichcoord,
int  i,
int  j,
int  k,
FTYPE pr,
FTYPE pstag 
)

Definition at line 3706 of file init.c.

Here is the call graph for this function:

Here is the caller graph for this function:

int init_global ( void  )

Definition at line 482 of file init.c.

Here is the caller graph for this function:

int init_grid ( void  )

Definition at line 3453 of file init.c.

Here is the caller graph for this function:

int init_grid_post_set_grid ( FTYPE(*)  prim[NSTORE2][NSTORE3][NPR],
FTYPE(*)  pstag[NSTORE2][NSTORE3][NPR],
FTYPE(*)  ucons[NSTORE2][NSTORE3][NPR],
FTYPE(*)  vpot[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3],
FTYPE(*)  Bhat[NSTORE2][NSTORE3][NPR],
FTYPE(*)  panalytic[NSTORE2][NSTORE3][NPR],
FTYPE(*)  pstaganalytic[NSTORE2][NSTORE3][NPR],
FTYPE(*)  vpotanalytic[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3],
FTYPE(*)  Bhatanalytic[NSTORE2][NSTORE3][NPR],
FTYPE(*)  F1[NSTORE2][NSTORE3][NPR+NSPECIAL],
FTYPE(*)  F2[NSTORE2][NSTORE3][NPR+NSPECIAL],
FTYPE(*)  F3[NSTORE2][NSTORE3][NPR+NSPECIAL],
FTYPE(*)  Atemp[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3] 
)

Definition at line 3495 of file init.c.

Here is the call graph for this function:

int init_postvpot ( int  i,
int  j,
int  k,
FTYPE pr,
FTYPE pstag,
FTYPE ucons 
)

Definition at line 8011 of file init.c.

Here is the caller graph for this function:

int init_primitives ( FTYPE(*)  prim[NSTORE2][NSTORE3][NPR],
FTYPE(*)  pstag[NSTORE2][NSTORE3][NPR],
FTYPE(*)  ucons[NSTORE2][NSTORE3][NPR],
FTYPE(*)  vpot[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3],
FTYPE(*)  Bhat[NSTORE2][NSTORE3][NPR],
FTYPE(*)  panalytic[NSTORE2][NSTORE3][NPR],
FTYPE(*)  pstaganalytic[NSTORE2][NSTORE3][NPR],
FTYPE(*)  vpotanalytic[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3],
FTYPE(*)  Bhatanalytic[NSTORE2][NSTORE3][NPR],
FTYPE(*)  F1[NSTORE2][NSTORE3][NPR+NSPECIAL],
FTYPE(*)  F2[NSTORE2][NSTORE3][NPR+NSPECIAL],
FTYPE(*)  F3[NSTORE2][NSTORE3][NPR+NSPECIAL],
FTYPE(*)  Atemp[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3] 
)

Definition at line 3633 of file init.c.

Here is the caller graph for this function:

int init_vpot2field_user ( SFTYPE  time,
FTYPE(*)  A[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3],
FTYPE(*)  prim[NSTORE2][NSTORE3][NPR],
FTYPE(*)  pstag[NSTORE2][NSTORE3][NPR],
FTYPE(*)  ucons[NSTORE2][NSTORE3][NPR],
FTYPE(*)  Bhat[NSTORE2][NSTORE3][NPR] 
)

Definition at line 8448 of file init.c.

Here is the caller graph for this function:

int init_vpot_user ( int *  whichcoord,
int  l,
SFTYPE  time,
int  i,
int  j,
int  k,
int  loc,
FTYPE(*)  prim[NSTORE2][NSTORE3][NPR],
FTYPE V,
FTYPE A 
)

Definition at line 7152 of file init.c.

Here is the caller graph for this function:

void Jcon_numerical ( int  whichcoord,
FTYPE X,
FTYPE Jcon 
)

Definition at line 8315 of file init.c.

Here is the call graph for this function:

Here is the caller graph for this function:

int jetbound ( int  i,
int  j,
int  k,
int  loc,
FTYPE prin,
FTYPE prflux,
FTYPE(*)  prim[NSTORE2][NSTORE3][NPR] 
)

Definition at line 8972 of file init.c.

Here is the call graph for this function:

Here is the caller graph for this function:

static FTYPE kappa_func_fits ( int  which,
FTYPE  rho,
FTYPE  B,
FTYPE  Tg,
FTYPE  Tr,
FTYPE  varexpf 
)
static

Definition at line 9614 of file init.c.

Here is the call graph for this function:

int kappa_func_fits_all ( FTYPE  rho,
FTYPE  B,
FTYPE  Tg,
FTYPE  Tr,
FTYPE  varexpf,
FTYPE kappa,
FTYPE kappaemit,
FTYPE kappan,
FTYPE kappanemit,
FTYPE kappaes 
)

kappaes = sigma_T n_e = sigma_T n_b (n_e/n_b) = sigma_T rho/mb (ne/nb)

Definition at line 9646 of file init.c.

FTYPE kcov_func_mcoord ( struct of_geom ptrgeom,
FTYPE X,
int  ii,
int  jj 
)

Definition at line 8295 of file init.c.

Here is the caller graph for this function:

static int make_nonrt2rt_solution ( int *  whichvel,
int *  whichcoord,
int  opticallythick,
FTYPE pp,
FTYPE X,
FTYPE V,
struct of_geom **  ptrptrgeom 
)
static

Definition at line 6202 of file init.c.

Here is the call graph for this function:

Here is the caller graph for this function:

FTYPE mcov_func_mcoord ( struct of_geom ptrgeom,
FTYPE X,
int  ii,
int  jj 
)

Definition at line 8271 of file init.c.

Here is the caller graph for this function:

int normalize_densities ( FTYPE(*)  prim[NSTORE2][NSTORE3][NPR])

Definition at line 8464 of file init.c.

Here is the call graph for this function:

Here is the caller graph for this function:

int normalize_field ( FTYPE(*)  prim[NSTORE2][NSTORE3][NPR],
FTYPE(*)  pstag[NSTORE2][NSTORE3][NPR],
FTYPE(*)  ucons[NSTORE2][NSTORE3][NPR],
FTYPE(*)  vpot[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3],
FTYPE(*)  Bhat[NSTORE2][NSTORE3][NPR] 
)

Definition at line 8535 of file init.c.

Here is the call graph for this function:

Here is the caller graph for this function:

static FTYPE nz_func ( FTYPE  R)
static

Definition at line 7010 of file init.c.

Here is the caller graph for this function:

int post_init_specific_init ( void  )

Definition at line 406 of file init.c.

Here is the caller graph for this function:

int pre_init_specific_init ( void  )

Definition at line 316 of file init.c.

Here is the caller graph for this function:

int prepre_init_specific_init ( void  )

Definition at line 245 of file init.c.

Here is the caller graph for this function:

int process_restart_toget_radiation ( void  )

Definition at line 6976 of file init.c.

Here is the call graph for this function:

static int process_solution ( int *  whichvel,
int *  whichcoord,
int  opticallythick,
FTYPE pp,
FTYPE X,
FTYPE V,
struct of_geom **  ptrptrgeom,
FTYPE ptptr 
)
static

Definition at line 6959 of file init.c.

void read_coord_parms_user ( int  defcoordlocal,
FILE *  in 
)

Definition at line 10318 of file init.c.

Here is the caller graph for this function:

void set_coord_parms_deps_user ( int  defcoordlocal)

Definition at line 10303 of file init.c.

Here is the caller graph for this function:

void set_coord_parms_nodeps_user ( int  defcoordlocal)

Definition at line 10210 of file init.c.

Here is the caller graph for this function:

FTYPE setblandfordfield ( FTYPE  r,
FTYPE  th 
)

Definition at line 7106 of file init.c.

Here is the call graph for this function:

Here is the caller graph for this function:

FTYPE setgpara ( FTYPE  myr,
FTYPE  th,
FTYPE  thpower 
)

Definition at line 7090 of file init.c.

Here is the caller graph for this function:

int setihor_user ( void  )

Definition at line 10821 of file init.c.

FTYPE setRin_user ( int  ihor,
FTYPE  ihoradjust 
)

Definition at line 10790 of file init.c.

Here is the caller graph for this function:

void shortout_Bd3 ( FTYPE(*)  prim[NSTORE2][NSTORE3][NPR])

Definition at line 8050 of file init.c.

Here is the call graph for this function:

static FTYPE taper_func2 ( FTYPE  R,
FTYPE  rin,
FTYPE  rpow 
)
static

Definition at line 7033 of file init.c.

static FTYPE taper_func_exp ( FTYPE  R,
FTYPE  rin,
FTYPE  POTENTIALorPRESSURE 
)
static

Definition at line 7022 of file init.c.

void write_coord_parms_user ( int  defcoordlocal,
FILE *  out 
)

Definition at line 10311 of file init.c.

Here is the caller graph for this function:

Variable Documentation

FTYPE aforwald

Definition at line 46 of file init.c.

FTYPE B0WALD

Definition at line 38 of file init.c.

SFTYPE beta
static

Definition at line 30 of file init.c.

SFTYPE bsq_max =0
static

Definition at line 29 of file init.c.

FTYPE BSQORHOWALD =100.0

Definition at line 45 of file init.c.

FTYPE cpow2
static

Definition at line 10207 of file init.c.

FTYPE cpow3
static

Definition at line 10207 of file init.c.

int DOWALDDEN =0

Definition at line 42 of file init.c.

FTYPE h0
static

Definition at line 10207 of file init.c.

FTYPE htheta
static

Definition at line 10207 of file init.c.

int inittypeglobal

Definition at line 56 of file init.c.

FTYPE njet
static

Definition at line 10207 of file init.c.

FTYPE njet1
static

Definition at line 10207 of file init.c.

FTYPE normglobal

Definition at line 55 of file init.c.

FTYPE npow
static

Definition at line 10207 of file init.c.

FTYPE npow2
static

Definition at line 10207 of file init.c.

FTYPE ntheta
static

Definition at line 10207 of file init.c.

SFTYPE pmax =0
static

Definition at line 29 of file init.c.

SFTYPE pradmax =0
static

Definition at line 29 of file init.c.

SFTYPE ptotmax =0
static

Definition at line 29 of file init.c.

FTYPE Qjet
static

Definition at line 10207 of file init.c.

FTYPE r0
static

Definition at line 10207 of file init.c.

FTYPE r0jet
static

Definition at line 10207 of file init.c.

FTYPE r0jet2
static

Definition at line 10207 of file init.c.

FTYPE r0jet3
static

Definition at line 10207 of file init.c.

FTYPE r1jet
static

Definition at line 10207 of file init.c.

FTYPE RADATM_FERATIO

Definition at line 125 of file init.c.

FTYPE RADATM_FRATIO

Definition at line 126 of file init.c.

FTYPE RADATM_LUMEDD

Definition at line 123 of file init.c.

FTYPE RADATM_MDOTEDD

Definition at line 122 of file init.c.

FTYPE RADATM_RHOAMB

Definition at line 127 of file init.c.

FTYPE RADATM_TAMB

Definition at line 128 of file init.c.

int RADATM_THINRADATM

Definition at line 124 of file init.c.

FTYPE RADBEAM2D_BEAML

Definition at line 110 of file init.c.

int RADBEAM2D_BEAMNO

Definition at line 96 of file init.c.

FTYPE RADBEAM2D_BEAMR

Definition at line 111 of file init.c.

int RADBEAM2D_BLOB

Definition at line 100 of file init.c.

FTYPE RADBEAM2D_BLOBP

Definition at line 102 of file init.c.

FTYPE RADBEAM2D_BLOBW

Definition at line 101 of file init.c.

FTYPE RADBEAM2D_BLOBX

Definition at line 103 of file init.c.

FTYPE RADBEAM2D_BLOBZ

Definition at line 104 of file init.c.

int RADBEAM2D_FLATBACKGROUND

Definition at line 97 of file init.c.

int RADBEAM2D_IFBEAM

Definition at line 107 of file init.c.

FTYPE RADBEAM2D_NLEFT

Definition at line 109 of file init.c.

FTYPE RADBEAM2D_PAR_D

Definition at line 105 of file init.c.

FTYPE RADBEAM2D_PAR_E

Definition at line 106 of file init.c.

FTYPE RADBEAM2D_RHOAMB

Definition at line 98 of file init.c.

FTYPE RADBEAM2D_TAMB

Definition at line 99 of file init.c.

FTYPE RADBEAM2D_TLEFT

Definition at line 108 of file init.c.

int RADBEAM2DKSVERT_BEAMNO

Definition at line 115 of file init.c.

FTYPE RADBEAMFLAT_ERAD

Definition at line 118 of file init.c.

FTYPE RADBEAMFLAT_FRATIO

Definition at line 117 of file init.c.

FTYPE RADBEAMFLAT_RHO

Definition at line 119 of file init.c.

FTYPE RADBEAMFLAT_UU

Definition at line 120 of file init.c.

FTYPE RADBONDI_MAXX

Definition at line 189 of file init.c.

FTYPE RADBONDI_MDOTEDD

Definition at line 187 of file init.c.

FTYPE RADBONDI_MDOTPEREDD

Definition at line 186 of file init.c.

FTYPE RADBONDI_MINX

Definition at line 188 of file init.c.

FTYPE RADBONDI_PRADGAS

Definition at line 184 of file init.c.

FTYPE RADBONDI_TESTNO

Definition at line 183 of file init.c.

FTYPE RADBONDI_TGAS0

Definition at line 185 of file init.c.

FTYPE RADCYLJET_EHATJET

Definition at line 228 of file init.c.

FTYPE RADCYLJET_RHOJET

Definition at line 229 of file init.c.

FTYPE RADCYLJET_TEMPJET

Definition at line 230 of file init.c.

int RADCYLJET_TYPE

Definition at line 226 of file init.c.

FTYPE RADCYLJET_UJET

Definition at line 231 of file init.c.

FTYPE RADCYLJET_VRSTAR

Definition at line 227 of file init.c.

FTYPE RADDBLSHADOW_ANGLE

Definition at line 136 of file init.c.

FTYPE RADDBLSHADOW_BEAMY

Definition at line 138 of file init.c.

FTYPE RADDBLSHADOW_NLEFT

Definition at line 135 of file init.c.

FTYPE RADDBLSHADOW_TLEFTOTAMB

Definition at line 137 of file init.c.

int RADDONUT_OPTICALLYTHICKTORUS

Definition at line 233 of file init.c.

FTYPE RADDOT_F1

Definition at line 201 of file init.c.

FTYPE RADDOT_F2

Definition at line 202 of file init.c.

FTYPE RADDOT_FYDOT

Definition at line 198 of file init.c.

int RADDOT_IDOT

Definition at line 195 of file init.c.

int RADDOT_JDOT

Definition at line 196 of file init.c.

int RADDOT_KDOT

Definition at line 197 of file init.c.

FTYPE RADDOT_LTEFACTOR

Definition at line 199 of file init.c.

FTYPE RADDOT_URFX

Definition at line 200 of file init.c.

FTYPE RADDOT_XDOT

Definition at line 192 of file init.c.

FTYPE RADDOT_YDOT

Definition at line 193 of file init.c.

FTYPE RADDOT_ZDOT

Definition at line 194 of file init.c.

FTYPE RADNT_DONUTRADPMAX

Definition at line 222 of file init.c.

FTYPE RADNT_DONUTTYPE

Definition at line 215 of file init.c.

FTYPE RADNT_ELL

Definition at line 209 of file init.c.

FTYPE RADNT_ERADATMMIN

Definition at line 214 of file init.c.

FTYPE RADNT_FULLPHI

Definition at line 221 of file init.c.

FTYPE RADNT_HOVERR

Definition at line 223 of file init.c.

FTYPE RADNT_INFLOWING

Definition at line 216 of file init.c.

FTYPE RADNT_KKK

Definition at line 208 of file init.c.

FTYPE RADNT_LPOW

Definition at line 224 of file init.c.

FTYPE RADNT_MAXX

Definition at line 207 of file init.c.

FTYPE RADNT_MINX

Definition at line 206 of file init.c.

FTYPE RADNT_OMSCALE

Definition at line 220 of file init.c.

FTYPE RADNT_RHOATMMIN

Definition at line 211 of file init.c.

FTYPE RADNT_RHODONUT

Definition at line 212 of file init.c.

FTYPE RADNT_ROUT

Definition at line 219 of file init.c.

FTYPE RADNT_TGASATMMIN

Definition at line 217 of file init.c.

FTYPE RADNT_TRADATMMIN

Definition at line 218 of file init.c.

FTYPE RADNT_UINTATMMIN

Definition at line 213 of file init.c.

FTYPE RADNT_UTPOT

Definition at line 210 of file init.c.

FTYPE RADSHADOW_ANGLE

Definition at line 131 of file init.c.

FTYPE RADSHADOW_BEAMY

Definition at line 133 of file init.c.

FTYPE RADSHADOW_NLEFT

Definition at line 130 of file init.c.

FTYPE RADSHADOW_TLEFTOTAMB

Definition at line 132 of file init.c.

FTYPE RADWAVE_AAA

Definition at line 171 of file init.c.

FTYPE RADWAVE_B0

Definition at line 150 of file init.c.

FTYPE RADWAVE_CC

Definition at line 147 of file init.c.

FTYPE RADWAVE_DB2IM

Definition at line 158 of file init.c.

FTYPE RADWAVE_DB2RE

Definition at line 157 of file init.c.

FTYPE RADWAVE_DEIM

Definition at line 162 of file init.c.

FTYPE RADWAVE_DERE

Definition at line 161 of file init.c.

FTYPE RADWAVE_DF2IM

Definition at line 166 of file init.c.

FTYPE RADWAVE_DF2RE

Definition at line 165 of file init.c.

FTYPE RADWAVE_DFIM

Definition at line 164 of file init.c.

FTYPE RADWAVE_DFRE

Definition at line 163 of file init.c.

FTYPE RADWAVE_DRIM

Definition at line 152 of file init.c.

FTYPE RADWAVE_DRRE

Definition at line 151 of file init.c.

FTYPE RADWAVE_DTOUT1

Definition at line 169 of file init.c.

FTYPE RADWAVE_DUIM

Definition at line 160 of file init.c.

FTYPE RADWAVE_DURE

Definition at line 159 of file init.c.

FTYPE RADWAVE_DV2IM

Definition at line 156 of file init.c.

FTYPE RADWAVE_DV2RE

Definition at line 155 of file init.c.

FTYPE RADWAVE_DVIM

Definition at line 154 of file init.c.

FTYPE RADWAVE_DVRE

Definition at line 153 of file init.c.

FTYPE RADWAVE_ERAD

Definition at line 175 of file init.c.

FTYPE RADWAVE_ERADFACTOR

Definition at line 178 of file init.c.

FTYPE RADWAVE_GASFACTOR

Definition at line 179 of file init.c.

FTYPE RADWAVE_KAPPA

Definition at line 148 of file init.c.

FTYPE RADWAVE_KAPPAES

Definition at line 177 of file init.c.

FTYPE RADWAVE_KK

Definition at line 172 of file init.c.

int RADWAVE_NUMERO

Definition at line 144 of file init.c.

int RADWAVE_NWAVE

Definition at line 143 of file init.c.

FTYPE RADWAVE_OMIM

Definition at line 168 of file init.c.

FTYPE RADWAVE_OMRE

Definition at line 167 of file init.c.

FTYPE RADWAVE_PP

Definition at line 146 of file init.c.

FTYPE RADWAVE_RHOFAC

Definition at line 149 of file init.c.

FTYPE RADWAVE_RHOZERO

Definition at line 170 of file init.c.

FTYPE RADWAVE_TEMP

Definition at line 174 of file init.c.

FTYPE RADWAVE_UINT

Definition at line 173 of file init.c.

FTYPE RADWAVE_VX

Definition at line 176 of file init.c.

int RADWAVE_WAVETYPE

Definition at line 145 of file init.c.

SFTYPE randfact
static

Definition at line 30 of file init.c.

FTYPE rbr
static

Definition at line 10207 of file init.c.

FTYPE rhodisk
static

Definition at line 31 of file init.c.

SFTYPE rhomax =0
static

Definition at line 29 of file init.c.

SFTYPE rin
static

Definition at line 30 of file init.c.

SFTYPE rinfield
static

Definition at line 30 of file init.c.

SFTYPE routfield
static

Definition at line 30 of file init.c.

FTYPE rs
static

Definition at line 10207 of file init.c.

FTYPE rsjet
static

Definition at line 10207 of file init.c.

FTYPE rsjet2
static

Definition at line 10207 of file init.c.

FTYPE rsjet3
static

Definition at line 10207 of file init.c.

FTYPE thindiskrhopow =-0.6

Definition at line 53 of file init.c.

FTYPE TILTWALD

Definition at line 47 of file init.c.

SFTYPE umax =0
static

Definition at line 29 of file init.c.

SFTYPE uradmax =0
static

Definition at line 29 of file init.c.

SFTYPE utotmax =0
static

Definition at line 29 of file init.c.

int WALDWHICHACOV =3

Definition at line 43 of file init.c.

FTYPE x1br
static

Definition at line 10207 of file init.c.