definit.h

Go to the documentation of this file.

#define ORDERSTORAGE 5

#define NEWMETRICSTORAGE 1

#define ASYMDIAGCHECK 0

#define ACCURATEDIAG 1
#define FLUXDIAGPARTS 0 // whether to compute parts of fluxes

// 0 : not accurate
// 1 : cumulative source accurate
// 2 : each component accurate
#define ACCURATESOURCEDIAG 1
#define DIAGSOURCECOMPSTEP 4 // how many fullsteps to wait to compute diag_source_comp() that is somewhat expensive but just diagnostics and not incredibly needed compared to diag_source_all()


#define ACCURATESINCOS 1

#define HOT2ENTROPY 1

#define HOT2COLD 1

#define HOT2FFDE 1

#define ENTROPY2COLD 1

#define ENTROPY2HOT 1

#define ENTROPY2FFDE 1




#define BADENERGYMAC(ugenergy,ugentropy) ((ugenergy) < 0.5*(ugentropy))
#define BADENERGY2MAC(ugenergy,ugentropy) ((ugenergy) < (ugentropy))

#define DIVMETHODPREFLUX 0
#define DIVMETHODPOSTFLUX 1

#define DIVERGENCEMETHOD 1 // choose, but normally 1 for best results.


#define COMPDIM 3

#define GETGEOMUSEPOINTER 0

#define SPLITNPR 0

#define ANALYTICMEMORY 1

#define WENOMEMORY 0

#define HIGHERORDERMEM 1

#define SPLITMAEMMEM (0) // CHANGINGMARK

#define DOMONOINTERP SMONOINTERP

#define SPLITPRESSURETERMINFLUXMA 0 // CHANGINGMARK

#define SPLITPRESSURETERMINFLUXEM 0 // CHANGINGMARK


#define MERGEDC2EA2CMETHODMA 0
#define MERGEDC2EA2CMETHODEM 0
#define MERGEDC2EA2CMETHOD (MERGEDC2EA2CMETHODMA || MERGEDC2EA2CMETHODEM)


#define HIGHERORDERROUGH 0
#define HIGHERORDERSMOOTH 1
#define HIGHERORDERVERYSMOOTH 2

#define HIGHERORDERTYPE HIGHERORDERROUGH // original method




#if(HIGHERORDERTYPE==HIGHERORDERROUGH)// non-smooth version:





#define DO_SPLITC2E NOSPLITA2C
#define DO_SPLITA2C ENERGY_CONTROLS_ALL_WEIGHTS // c2a and a2c for conserved quants
#define DO_SPLITA2CSMOOTH DO_SPLITA2C
#define DO_SPLITA2C4FLUX DO_SPLITA2C // more generally
#define DO_SPLITA2C4MAFLUX MASSENERGYMOMENTUM_IS_COUPLED_WEIGHTS // couples weights together
#define DO_SPLITA2C4SMOOTHFLUX DO_SPLITA2C
#define DO_SPLITSOURCE NOSPLITA2C // c2a for ENOSOURCETERM

#define EMFLUXFIELDTERMTYPE ENOFLUX
#define FIELDTERMTYPE ENOFLUX
#define FIELDFVTYPE ENOCONSERVED



#elif(HIGHERORDERTYPE==HIGHERORDERSMOOTH)  // smooth version



#define DO_SPLITC2E NOSPLITA2C
#define DO_SPLITA2C ENERGY_CONTROLS_ALL_WEIGHTS // c2a and a2c for conserved quants
#define DO_SPLITA2CSMOOTH CONSTANT_ALL_WEIGHTS






#define DO_SPLITA2C4FLUX MASSENERGYMOMENTUM_IS_COUPLED_WEIGHTS
#define DO_SPLITA2C4MAFLUX DO_SPLITA2C4FLUX

#define DO_SPLITA2C4SMOOTHFLUX CONSTANT_ALL_WEIGHTS
#define DO_SPLITSOURCE NOSPLITA2C // c2a for ENOSOURCETERM

#define EMFLUXFIELDTERMTYPE ENOSMOOTHFLUX
#define FIELDTERMTYPE ENOSMOOTHFLUX
#define FIELDFVTYPE ENOSMOOTHCONSERVED



#elif(HIGHERORDERTYPE==HIGHERORDERVERYSMOOTH) // supersmooth version




#define DO_SPLITC2E NOSPLITA2C
#define DO_SPLITA2C CONSTANT_ALL_WEIGHTS
#define DO_SPLITA2CSMOOTH CONSTANT_ALL_WEIGHTS
#define DO_SPLITA2C4FLUX CONSTANT_ALL_WEIGHTS
#define DO_SPLITA2C4MAFLUX CONSTANT_ALL_WEIGHTS
#define DO_SPLITA2C4SMOOTHFLUX CONSTANT_ALL_WEIGHTS
#define DO_SPLITSOURCE NOSPLITA2C

#define EMFLUXFIELDTERMTYPE ENOSMOOTHFLUX
#define FIELDTERMTYPE ENOSMOOTHFLUX
#define FIELDFVTYPE ENOSMOOTHCONSERVED

#endif





#define VELTERMSMINIMIZE(pl) (pl==RHO || pl==UU || pl==UU+dir)
#define ORTHOVEL1TERMSMINIMIZE(pl) (pl==UU+odir1)
#define ORTHOVEL2TERMSMINIMIZE(pl) (pl==UU+odir2)
#define PRESSUREMINIMIZE(pl) (pl==FLUXSPLITPMA(dir)) // pressure term when splitmaem==1 and doing MA term
#define EMFTERMS(pl) (pl>=B1 && pl<=B3) // only used if splitmaem==0

#define ALLOTHERSMINIMIZE(pl) (! (VELTERMSMINIMIZE(pl) || ORTHOVEL1TERMSMINIMIZE(pl) || ORTHOVEL2TERMSMINIMIZE(pl) || PRESSUREMINIMIZE(pl) || ((emffixedstencil==1)&&EMFTERMS(pl))  ) )




#define FIELDSTAGMEM 1

#define STOREFLUXSTATE 1

#define STORESHOCKINDICATOR 1

#define FIELDTOTHMEM 1 // CHANGINGMARK


#define TRACKVPOT 0

#define EVOLVEWITHVPOT TRACKVPOT // choice

#define MODIFYEMFORVPOT MODIFYVPOT

#define CORNGDETVERSION 1




#define N1 64  // number of zones in 1-direction
#define N2 64        // number of zones in 2-direction
#define N3      1               // number of zones in 3-direction

#define MAXBND 5





#define MCOORD KSCOORDS


#define MAXWELL PRIMMAXWELL


#define UD0PLUS1FIX 1


#define DO_VORTICITY_IMAGE 0


#define PRODUCTION 0







#define FULLOUTPUT 0


#define MAILWHENDONE 1
#define MAILFROMREMOTE 0
#define REMOTEHOST "ki-rh42.slac.stanford.edu"
#define REMOTEUSER "jon"
#define EMAILADDRESS "jmckinne@stanford.edu"
#define EMAILMESSAGE "Done with GRMHD run DEFAULT"



#define SUPERLONGDOUBLE 0


#define PERFTEST 0



#define DOAVG 0
#define DOAVGDUMP (DOAVG)



#define DOJETDIAG 1



#define DOAVG2 0 // only make 1 if DOAVG 1 above




#define DOGDUMP 1



#define DORDUMP 1



#define DOMAINDUMP 1
#define DOMAINDUMPDIAG (DOMAINDUMP)



#define DODEBUGDUMP 1



#define DODEBUG 1



#define DOFLOORDUMP 1
#define DOFLOORDIAG (DOFLOORDUMP)



#define DODISSMEASUREDUMP 0
#define DODISSMEASUREDIAG (DODISSMEASUREDUMP)



#define DOEOSDUMP 1
#define DOEOSDIAG (DOEOSDUMP)



#define DORADDUMP 1
#define DORADDIAG (DORADDUMP)




#define DOVPOTDUMP (1 && TRACKVPOT)

#define DOIMAGEDUMP (1)

#define DOAREAMAPDUMP (1)


#define DOENODEBUG 0





#define DOGRIDSECTIONING 0





#define WENO_USE_PRIM_REDUCTION 1





#define WENO_EXTRA_A2C_MINIMIZATION 0 // CHANGINGMARK

#define WENO_REDUCE_A2C_LOOK_OTHER_DIRECTIONS 1




#define LIMIT_FLUXC2A_PRIM_CHANGE 0


 
#define DO_WENO_DEBUG 0


#define DOSUPERDEBUG 0



#define ALLOWMETRICROT 0





#if(ALLOWMETRICROT)
#define DOMIXTHETAPHI 1 // for g_{\theta\phi} // no choice
#else
#define DOMIXTHETAPHI 0 // choice
#endif





#define DOEVOLVEMETRIC 0

#define EVOLVEMETRICSUBSTEP 0

#define LIMITDTWITHSOURCETERM 0

#define LIMITSOURCES 0


#define USEGRAVITYDTINDTLIMIT 0


#define RESTRICTDTSETTINGINSIDEHORIZON 1


#define DODISS 0

#define DODISSMEASURE 0



#define DOLUMVSR 0

#define DODISSVSR 0

#define WHICHENTROPYEVOLVE EVOLVESIMPLEENTROPY



#define DOSELFGRAVVSR 0

#define DOFIELDLINE 1



#define ROEAVERAGEDWAVESPEED 0

#define ATHENAROE 0

#define STOREWAVESPEEDS 0

#define PERCELLDT 1

#define USESTOREDSPEEDSFORFLUX (STOREWAVESPEEDS) // choice really independent of STOREWAVESPEEDS, but generall normally want to couple them

#define VCHARTYPE VERYLOCALVCHAR

#define TRUEFAST 0

#define FORCESOLVEL 0

#define FORCESOLVELFLUX 0


#define CHECKONINVERSION 1
#define CHECKONINVERSIONRAD 0 // check on radiation inversion.  Default 0 because generally limiters on |v| mean not the same, so only use for debugging.

#define PRECISEINVERSION 1


#define WHICHVEL VELREL4



#define REMOVERESTMASSFROMUU 1

#define RELTYPE RELEOM

#define EOMTYPE EOMGRMHD


#define EOMRADTYPE EOMRADNONE // EOMRADM1CLOSURE
#define EVOLVENRAD 0 // whether to evolve number density of radiation

#define WHICHRADSOURCEMETHOD SOURCEMETHODNONE




#define UTOPRIMTRYAGAIN 0


#define WHICHEOS IDEALGAS






#define DOENTROPY DONOENTROPY // normal total energy equation

#define FIXUPAFTERINIT 1

#define FIXUPAFTERRESTART 1

#define CHECKSOLUTION 0



#define GAMMAPERCDIFFMAX 2.0
#define UPERCDIFFMAX 10.0


#define DOEXTRAINTERP 0


#define DOYFL 0

#define DOYL 0

#define DOYNU DOYL


#define LIMADJUST LIMITERFIXED


#define HYDROLIMADJUSTONLY 1

#define FLUXADJUST FLUXFIXED


#define HYDROFLUXADJUSTONLY 1

#define STEPOVERNEGU NEGDENSITY_FIXONFULLSTEP


#define STEPOVERNEGRHO NEGDENSITY_FIXONFULLSTEP

#define STEPOVERNEGURHO NEGDENSITY_FIXONFULLSTEP




#define UTOPRIMADJUST UTOPRIMAVG

#define UTOPRIMFIXMPICONSISTENT 1

#define UTOPRIMFAILRETURNTYPE UTOPRIMRETURNADJUSTED

#define SMOOTHSING 1

#define COORDSINGFIXCYL 0 //whether perform the same fix for CYLMINKMETRIC's z axis

#define COORDSINGFIX 1




#define SINGSMALL (1000*NUMEPSILON) // must be larger than machine precision to work for outer M_PI boundary!.  1E-14 works, but need some insurance so use 1E-13


#define VOLUMEDIFF 0


#define DOSTOREPOSITIONDATA 1


#define STORETLAB2ORTHO 0

#define CONNDERTYPE DIFFNUMREC // improved now and much more accurate then DIFFGAMMIE in general. However this is too slow to be used when time-dependent metric is cycling near substeps.

#define CONNMACHINEBODY 1

#define CONNAXISYMM 1


#define WHICHEOM WITHGDET



#define GDETVOLDIFF 0

#define FIXGDETSPC_WHEN_1DRADIAL 1

#define FORCEGDETPOSITIVE 0

#define MINDT 1.e-20 // minimum dt

#define JONCHECKS 1 // for vel=3 extra checks (standard things)
#define JONCHECKS2 1 // check_pr() crazy thing for vel3 (crazy u^t check and fix -- 2D only)


#define FLOORDIAGS 1 // whether to compute floor diagnostics

#define ANALYTICGCON 0

#define ANALYTICCONNECTION 0 // whether to use analytic connection


#define ANALYTICSOURCE 0 // whether to use analytic source function



#define OUTFLOWAVOIDBC 0

#define FLUXDIMENSPLIT PERFECTUNSPLIT
#define A2CDIMENSPLIT PERFECTUNSPLIT

#define BOUNDFLUXRECON 0

#define BOUNDUNEW 0

#define DODQMEMORY 1

#define BOUNDARYINTERPADJUST 0

#define COMPUTEFRDOT 0

#define CALCFARADAYANDCURRENTS 0 // FUCK: Some memory leak maybe because Kraken crashes if this is 1.  Didn't crash before, so unsure what's wrong.  core dump gives some info, but not much.

#define WHICHCURRENTCALC CURRENTCALC1


#define FARADAYT0 1

#define CURRENTST0 1


#define BOUNDPLPR 0
#define NOFLUXCTONX1DN 0


#define OUTFLOW 0
#define SYMM 1
#define ASYMM 2
#define FIXED 3
#define POLARAXIS 4
#define FIXEDOUTFLOW 5 // means fixed inflow but allows outflow -- basically outflow if no inflow, but if inflow then set values to fixed quantities
#define NSSURFACE 6 // whatever in bounds.c for NS surface
#define PERIODIC 7 // periodic boundary conditions
#define OUTFLOWNOINFLOW 8  //copy if velocity directed outward; force velocity to zero and copy if directed inward (into the grid from the boundary)
#define RAMESHOUTFLOW 9 //OUTFLOW quantities according the Ramesh's power-law solution
#define HORIZONOUTFLOW 10

#define BCEXTRAP 100  //atch: extrapolation into the boundary with high order (can be specified)
#define CYLAXIS 101   //atch: cylindrical axis BC
#define BCEXTRAP_VEL3 102  //atch: the same as BCEXTRAP but extrapolates 3-velocity rather than 4-velocity; important for Hubble flow
#define JETINJECTION 103 //Fixed boundary condition for jet injection surrounded by an outflow condition
#define BCU1EXTRAPOTHERFIXED 104
#define BCEXTRAPCONSTRAINED 105
#define RESCALEOUTFLOW 106
#define RESCALEFIXEDALLOUTFLOWU1 107
#define FIXED_RESCALEOUTFLOWU1 108
#define BONDIMDOTOUTFLOW 109
#define BONDIINTOUTFLOW 110
#define DISKSURFACE 111
#define FREEOUTFLOW 112
#define R0SING 113
#define FIXEDUSEPANALYTIC 200


#define EVOLVECHECKS 0


#define FIXUPZONES FIXUP1ZONE

#define HLLBOUNDARY 0



#define FIXUPFLUX 0





#define ZEROOUTFLOWFLUX 0









#define ZEROPOLEFLUX 0







#define FLIPGDETAXIS 1

#if(0)
#define FLIPU3AXIS 1
#define FLIPB3AXIS 1

#define FLIPU2AXIS 1
#define FLIPB2AXIS 1

#define FLIPU1AXIS 0
#define FLIPB1AXIS 0

#else

#define FLIPU3AXIS 0
#define FLIPB3AXIS 0
#define FLIPU2AXIS 0
#define FLIPB2AXIS 0
#define FLIPU1AXIS 0
#define FLIPB1AXIS 0
#endif

#define DOPOLEDEATH 0
#define DOPOLESMOOTH 0 // can choose 1, but probably not necessary and generally won't treat flow correctly near pole even if possibly more robust.
#define DOPOLEGAMMADEATH 0



#define IF3DSPCTHENMPITRANSFERATPOLE 1 // working fine now that wavespeed bug in fluxctstag.c was fixed, extrapfunc B1,B2 bug fixed, and extrap gdet B3 instead of Bd3 that exaggerates extrapolation near poles and inconsistent with interpolation.  Also using VARTOINTERPFIELD GDETVERSION.




#define RESCALEINTERP 0



#define RESCALEINTERPFLUXCTSTAG 0



#define BDIRCONT 1



#define HYPERHLL 0


#define HORIZONSUPERFAST 0





#define VARTOINTERP PRIMTOINTERP

#define VARTOINTERPFIELD GDETVERSION // most consistent with fluxctstag.c and standard extrapfunc in bounds.tools.c




#define VLINEWITHGDETRHO 1

#define PLINEWITHFIELD 1

#define PARALINEUSESMONO 0


#define USEAVGPRIMITIVEFORWENOFLAT 1


#define USEPRIMITIVEFROMAVGCONSERVED 0


#define CONTACTINDICATOR 0

#define COMPUTEDRHODP 1


#define SUPERFASTDIVREDUCE 0 //atch

#define MINPREFORDER 3

#define SHOCKINDICATOR 1

#define WHICHPARA PARA4


#define BONDI_BOUNDARY_SET_PL_PR 0  //do not analytically set p_l & p_r at the outer boundary for the Bondi problem




#define NUMPANALYTICOTHER 0
#define DODUMPOTHER 0 // whether to dump other stuff

#define FLUXDUMP 0

#define NUMPOTHER 0







#if((REALTYPE==DOUBLETYPE)||(REALTYPE==FLOATTYPE))

#define CONNDELTA (MY1EM5) // default -- seems to work pretty good generally to reduce max error


#elif(REALTYPE==LONGDOUBLETYPE)

#define CONNDELTA (MY1EM5) // based on min of error for specific case

#endif



#define USERERSETREGION 0


#define BCSIGMACONSTATPOLE 0

#define DOONESTEPDUACCOUNTING 1

#define FIELDLINEGDETB 0

#define DO_ASSERTS 0










#define ADJUSTFLUXCT 0 // whether to adjust fluxCT

#define ADJUSTFLUX 0 // whether to adjust normal flux


#define DEBUGINOUTLOOPS 0

#define OUTEREXTRAP 3





#define HORIZONEXTRAP 3



#define NUMITERVPHI 5


#define POLEDEATH (MIN(DOPOLEDEATH,N2BND)) // with expansion by 1 point if detects jumps in densities or Lorentz factor (see poldeath())

#define MAXPOLEDEATH (MIN(DOPOLEDEATH+1,N2BND)) // can't be larger than N2BND
#define DEATHEXPANDAMOUNT 0

#define POLEINTERPTYPE 3 // 0=set uu2=bu2=0, 1=linearly interpolate uu2,bu2  2=interpolate B_\phi into pole  3 =linearly for uu2 unless sucking on pole





#define POLEGAMMADEATH (MIN(DOPOLEGAMMADEATH,N2BND))



#define GAMMAPOLEOUTGOING 1.1 // keep low
#define GAMMAPOLEINGOINGOUT 1.3 // keep low
#define GAMMAPOLEOUTGOINGPOWER 0.1 // keep low unless know will grow alot
#define GAMMAPOLEOUTGOINGRADIUS 6.0 // very model dependent
#define GAMMAPOLEOUTGOINGRADIUSIN 4.0 // very model dependent
#define GAMMAPOLEOUTGOINGMAX (5.0)
#define GAMMAPOLEINGOING GAMMAMAX

#define GAMMARADPOLEOUTGOING 1.3 // keep low
#define GAMMARADPOLEINGOINGOUT 1.7 // keep low
#define GAMMARADPOLEOUTGOINGPOWER 0.5
#define GAMMARADPOLEOUTGOINGRADIUS 6.0 // very model dependent
#define GAMMARADPOLEOUTGOINGRADIUSIN 4.0 // very model dependent
#define GAMMARADPOLEOUTGOINGMAX (3.0)
#define GAMMARADPOLEINGOING GAMMAMAXRAD

#define POLEDENSITYDROPFACTOR 5.0
#define POLEGAMMAJUMPFACTOR 2.0

#define RADIUSMOREDEATH (300.0)

#define AVERAGEINRADIUS 0 // not correct  across MPI boundaries since have to shift near boundary yet need that last cell to be consistent with as if no MPI boundary // OPENNPMARK: Also not correct for OpenMP
#define RADIUSTOSTARTAVERAGING 7 // should be beyond horizon so doesn't diffuse across horizon
#define RADIUSTOAVOIDRADIALSUCK (2.0*Rhor)


#define IFLIMITPOLEDEATH 0

// whether to limit poledeath to inside horizon and inside death region
#define LIMITEDPOLEDEATHINRADIUS 1


#define RADIUSLIMITPOLEDEATHIN (3.0) // choose r=3M since always close to BH but always slightly outside horizon to help control stability.

#define IFLIMITPOLEDEATHIOUT (-100)



#define POLESMOOTH (MIN(DOPOLESMOOTH,N2BND))

// defaults for bounds.tools.c
#define OUTERDEATH 0
#define OUTERDEATHRADIUS (1E3)
#define OUTERDEATHGAMMAMAX (6.0)
#define OUTERDEATHGAMMAMAXRAD (GAMMAMAXRAD)

// whether to use old non-continuous method for Bi along i interpolation in staggered field method
#define OLDNONCONT 0

#define OUTERRADIALSUPERFAST 0

#define DOCOMPTON 0
// whether to (during utoprimgen()) to do mhd inversion, recompute umhd, use energy conservation to get urad, and then invert that urad->prad.  Doesn't work if mhd>>rad when non-QTYPMHD methods used.
#define ENFORCEMHDCONS2RADCONS 0
#define DOPERF 0 // whether to enable performance opts, including enforcing only QTYPMHD method.
#define ENSURECONS 0 // when done with full step, whether to ensure energy conservation by modifying rad result if possible
#define BORROWENTROPY 0 // whether to try borrowing energy/entropy from radiation to get total energy conservation with entropy method