HARM
harm and utilities
 All Data Structures Files Functions Variables Typedefs Macros Pages
defs.general.h
Go to the documentation of this file.
1 
14 
15 int BEGINOPENMPSHAREDLIST;
16 
17 
18 
19 #include "mpidefs.h"
20 
21 
22 #include "kazfulleos.defsglobalprivate.h" // put here so OpenMP private globals are defined before global.nondepmnemonics.h sets up thread private pragma's
23 
24 
25 
26 
27 #include "rancdefs.h"
28 
29 FTYPE Xmetricnew[NDIM],Xmetricold[NDIM]; // used to store time of latest and oldest metric
30 
31 SFTYPE *lumvsr,*lumvsr_tot;
32 
33 SFTYPE *dissvsr[NUMDISSVERSIONS],*dissvsr_tot[NUMDISSVERSIONS];
34 
35 FTYPE *rcent,*rcent_tot;
36 
37 SFTYPE *dVvsr,*dVvsr_tot,*vrsqvsr,*vrsqvsr_tot,*dMvsr,*dMvsr_tot,*dTrrvsr, *dTrrvsr_tot,*Mvsr_tot,*Mvsrface1_tot,*MOrvsr_tot,*phivsr_tot,*dJvsr,*dJvsr_tot,*Jvsr_tot,*Jvsrface1_tot;
38 
39 
42 /* physics parameters */
43 FTYPE gam,gamideal,ARAD_CODE;
44 
45 /* numerical parameters */
46 int defcoord;
47 FTYPE Rin, R0, Rout, hslope, Zin, Zout;
48 FTYPE Rin_array[NDIM], Rout_array[NDIM]; //atch -- arrays for a more general way of handling the grid dimensions
49 FTYPE Risco,Rhor;
50 FTYPE cour;
51 FTYPE dV, dVF, dx[NDIM], startx[NDIM], endx[NDIM], Diffx[NDIM];
52 SFTYPE dt,t,tf,tstepparti,tsteppartf;
53 SFTYPE fluxdumpdt;
54 long fluxdumprealnstep;
55 FTYPE TDYNORYEglobal,Hglobal;
56 FTYPE rcurr, hcurr;
57 FTYPE drsing;
58 
60 FTYPE FLUXDISSIPATION,EMFDISSIPATION;
61 
62 
64 #if(SIMULBCCALC!=-1)
65 int isc,iec,jsc,jec;
66 int isf1,ief1,jsf1,jef1,ksf1,kef1;
67 int isf2,ief2,jsf2,jef2,ksf2,kef2;
68 int isf3,ief3,jsf3,jef3,ksf3,kef3;
69 int ise,iee,jse,jee;
70 int isf1ct,ief1ct,jsf1ct,jef1ct;// GODMARK: other stage type requires more
71 int isf2ct,ief2ct,jsf2ct,jef2ct;
72 int isf3ct,ief3ct,jsf3ct,jef3ct;
73 int isdq,iedq,jsdq,jedq;
74 int ispdq,iepdq,jspdq,jepdq;
75 #endif
76 
77 FTYPE mydminarg1, mydminarg2;
78 long nstep;
79 int specialstep;
80 
81 int steppart,numstepparts;
82 
83 int gocont; // used to continue running(runtype, directinput, timereenter)
84 int runtype;
85 
86 /* output parameters */
87 FILE *log_file;
88 FILE *fail_file;
89 FILE *logfull_file;
90 FILE *logdt_file;
91 FILE *logdtfull_file;
92 FILE *logstep_file;
93 FILE *logperf_file;
94 SFTYPE DTstep,DTstepdot,DTperf,DTgocheck,DTtimecheck,DTperfdump;
95 int reallaststep,onemorestep;
96 
98 int PVER,GRIDVER,DVER,FLVER,NPVER,AVG1DVER,AVG2DVER,ENERVER,MODEVER,LOSSVER,SPVER,TSVER,LOGDTVER,STEPVER,PERFVER,ADVER,PDVER,CALCVER,FLINEVER;
99 
100 int PTYPE,GRIDTYPE,DTYPE,FLTYPE,NPTYPE,AVG2DTYPE,AVG1DTYPE,ENERTYPE,LOSSTYPE,SPTYPE,TSTYPE,LOGDTTYPE,STEPTYPE,PERFTYPE,ADTYPE,PDTYPE,CALCTYPE,FLINETYPE,MODETYPE,EXPANDTYPE,NPCOMPUTETYPE;
101 
102 char dumpnamelist[NUMDUMPTYPES][MAXFILENAME];
103 SFTYPE DTdumpgen[NUMDUMPTYPES];
104 long dumpcntgen[NUMDUMPTYPES];
105 //SFTYPE DTd;
106 //SFTYPE DTavg;
107 //SFTYPE DTener;
108 //SFTYPE DTi;
109 //SFTYPE DTdebug;
110 long DTr;
111 long DTfake;
112 //long dump_cnt;
113 //long avg_cnt;
114 //long debug_cnt;
115 //long image_cnt;
116 long rdump_cnt;
117 //long fieldline_cnt; // assumed to keep track with images (as in diag.c), so no need to include in restart()
118 
119 int nstroke; // OPENMPMARK: Bad in some inversion codes
120 
121 
123 FTYPE t_transition_in,t_transition_out;
124 int global_enerregiondef[NUMENERREGIONS][NUMUPDOWN][NDIM];
125 
126 
127 
128 
129 /* global flags */
130 int failed;
131 int lim[NDIM],fluxmethod[NPR],FLUXB,UTOPRIMVERSION,TIMEORDER,TIMETYPE,DOENOFLUX,avgscheme[NDIM];
132 int dofluxreconevolvepointfield,emffixedstencil,extrazones4emf,splitmaem,unewisavg;
133 int do_transverse_flux_integration[NPR],do_conserved_integration[NPR],do_source_integration[NPR];
134 int useghostplusactive;
135 FTYPE defcon;
136 int global_preinterpolate;
137 
138 /* diagnostics */
140 SFTYPE frdot[N1][NPR];
141 SFTYPE pdottermsjet2[COMPDIM*2][NUMFLUXTERMS][NPR];
142 CTYPE failfloorcountlocal[2][NUMTSCALES][NUMFAILFLOORFLAGS]; // don't track this separately in jet
143 CTYPE failfloorcountlocal_tot[2][NUMTSCALES][NUMFAILFLOORFLAGS]; // don't track this separately in jet
144 
146 int dothisenerreg[NUMENERREGIONS];
147 int dofluxreg[NUMENERREGIONS][COMPDIM*2];
148 int enerposreg[NUMENERREGIONS][COMPDIM*2];
152 SFTYPE fladdreg[NUMENERREGIONS][NPR];
153 SFTYPE fladdreg_tot[NUMENERREGIONS][NPR];
154 SFTYPE fladdtermsreg[NUMENERREGIONS][NUMFAILFLOORFLAGS][NPR];
155 SFTYPE fladdtermsreg_tot[NUMENERREGIONS][NUMFAILFLOORFLAGS][NPR];
156 SFTYPE Ureg_init[NUMENERREGIONS][NPR];
157 SFTYPE Ureg_init_tot[NUMENERREGIONS][NPR];
158 SFTYPE pcumreg[NUMENERREGIONS][COMPDIM*2][NPR];
159 SFTYPE pcumreg_tot[NUMENERREGIONS][COMPDIM*2][NPR];
160 SFTYPE pdotreg[NUMENERREGIONS][COMPDIM*2][NPR];
161 SFTYPE pdottermsreg[NUMENERREGIONS][COMPDIM*2][NUMFLUXTERMS][NPR];
162 SFTYPE sourceaddreg[NUMENERREGIONS][NPR];
163 SFTYPE sourceaddreg_tot[NUMENERREGIONS][NPR];
164 SFTYPE sourceaddtermsreg[NUMENERREGIONS][NUMSOURCES][NPR];
165 SFTYPE sourceaddtermsreg_tot[NUMENERREGIONS][NUMSOURCES][NPR];
166 SFTYPE dissreg[NUMENERREGIONS][NUMDISSVERSIONS];
167 SFTYPE dissreg_tot[NUMENERREGIONS][NUMDISSVERSIONS];
168 //SFTYPE horizonflux[NPR];
169 //SFTYPE horizoncum[NPR];
170 //SFTYPE horizonflux_tot[NPR];
171 //SFTYPE horizoncum_tot[NPR];
172 
175 SFTYPE pdotreg_tot[NUMENERREGIONS][COMPDIM*2][NPR];
176 SFTYPE pdottermsreg_tot[NUMENERREGIONS][COMPDIM*2][NUMFLUXTERMS][NPR];
177 
181 int *doflux;
182 int *enerpos;
183 SFTYPE *fladd;
184 SFTYPE *fladd_tot;
185 SFTYPE (*fladdterms)[NPR];
186 SFTYPE (*fladdterms_tot)[NPR];
187 SFTYPE *U_init;
188 SFTYPE *U_init_tot;
189 SFTYPE (*pcum)[NPR];
190 SFTYPE (*pcum_tot)[NPR];
191 SFTYPE (*pdot)[NPR];
192 SFTYPE (*pdotterms)[NUMFLUXTERMS][NPR];
193 SFTYPE *sourceadd;
194 SFTYPE *sourceadd_tot;
195 SFTYPE (*sourceaddterms)[NPR];
196 SFTYPE (*sourceaddterms_tot)[NPR];
197 SFTYPE *diss;
198 SFTYPE *diss_tot;
199 
201 SFTYPE (*pdot_tot)[NPR];
202 SFTYPE (*pdotterms_tot)[NUMFLUXTERMS][NPR];
203 
204 
206 
207 
208 /* Jon's addition */
209 int horizoni,horizoncpupos1;
210 long realnstep;
211 int partialstep;
212 int mpicombine;
213 int mpicombinetype;
214 int truempicombinetype;
215 int halftimep;
216 int appendold;
217 int whocalleducon; // OPENMPNOTE: Ensure those are set as threadprivate [noted only called outside parallel regions]
219 int whichrestart;
220 int whichfake;
221 long restartsteps[2];
222 long fakesteps[2];
223 int binaryoutput,sortedoutput;
224 int CHECKCONT,DOTSTEPDIAG,DOLOGSTEP,DOLOGPERF;
225 int NDTCCHECK,NZCCHECK,NDTDOTCCHECK,NGOCHECK,NTIMECHECK,NDTPERFDUMPCHECK;
226 SFTYPE PERFWALLTIME,ZCPSESTIMATE;
227 
228 FTYPE zerouuperbaryon;
229 
230 long steptofaildump,steptofailmap;
231 int ifail,jfail,kfail; // OPENMPNOTE: Ensure those are set as private [noted only for diagnostics outside parallel regions]
232 int dofailmap,dofaildump,restartonfail;
234 FTYPE h_over_r;
236 FTYPE h_over_r_jet;
237 int BCtype[COMPDIM*2];
238 int rescaletype;
239 int cooling;
240 int DOENERDIAG,DOGDUMPDIAG,DORDUMPDIAG,DODUMPDIAG,DOAVGDIAG, DOIMAGEDIAG,DOAREAMAPDIAG;
241 int GAMMIEDUMP,GAMMIEIMAGE,GAMMIEENER,DODIAGS,RESTARTMODE,WHICHFILE,POSDEFMETRIC,DOENODEBUGEVERYSUBSTEP,DODIAGEVERYSUBSTEP;
242 int INVERTFROMAVERAGEIFFAILED,LIMIT_AC_PRIM_FRAC_CHANGE,LIMIT_AC_FRAC_CHANGE; //atch
243 int PARAMODWENO;
244 FTYPE MAX_AC_FRAC_CHANGE,MAX_AC_PRIM_FRAC_CHANGE;
245 FTYPE RHOMIN,UUMIN,SPECIFICENTROPYMIN,RHOMINLIMIT,UUMINLIMIT,SPECIFICENTROPYMINLIMIT;
246 FTYPE prMAX[NPR];
247 FTYPE prfloorcoef[NPR];
248 FTYPE BSQORHOLIMIT,BSQOULIMIT,UORHOLIMIT,GAMMAMAX,GAMMADAMP,GAMMAFAIL;
249 FTYPE ERADLIMIT,GAMMAMAXRAD,GAMMAMAXRADFAIL;
250 FTYPE UULIMIT;
251 FTYPE SAFE;
252 int debugfail;
253 FTYPE uttdiscr; // OPENMPNOTE: Ensure those are set as threadprivate [noted only for WHICHVEL==VEL3]
254 int jonchecks;
255 int dnumcolumns[NUMDUMPTYPES];
256 int dnumversion[NUMDUMPTYPES];
257 struct blink * blinkptr0[NUMDUMPTYPES];
258 struct blink * cpulinkptr0[NUMDUMPTYPES];
259 int DOCOLSPLIT[NUMDUMPTYPES];
260 int docolsplit; // global var for now CHANGINGMARK
261 int nextcol;
262 int doevolvemetricsubsteps, gravityskipstep;
263 FTYPE gravitydtglobal, sourcedtglobal, wavedtglobal;
264 int waveglobaldti[NDIM],waveglobaldtj[NDIM],waveglobaldtk[NDIM];
265 int didstorepositiondata,didstoremetricdata;
266 
267 /* physical consts */
268 FTYPE msun,lsun,rsun,G,H,C,qe,Na,malpha,mn,me,kb,arad,sigmasb,sigmamat,mevocsq,ergPmev,mp,Q,R,Re,hpl,hbar,K,K2;
269 SFTYPE a,MBH,QBH,EP3,THETAROT,THETAROTMETRIC,THETAROTPRIMITIVES;
270 FTYPE Mfactor,Jfactor,rhofactor;
271 SFTYPE dabh,dE,dJ,dEold,dJold;
272 FTYPE mb,mbcsq,mbwithrhounit,amu,a0,MBH0,QBH0,EP30,THETAROT0,Mdot,Mdotc,Mcgs,Ccode;
273 FTYPE Lunit,Tunit,Vunit,rhounit,rhomassunit,Munit,mdotunit,energyunit,edotunit,Pressureunit,Tempunit,Bunit,massunitPmsun;
274 int rho0unittype;
275 FTYPE ledd,leddcode;
276 
277 int NUMBUFFERS;
278 
279 int advancepassnumber;
280 
284 int nprstart,nprend; // normally 0 and NPR-1
285 int nprlist[MAXNPR]; // maximum is NPR elements
286 
288 int npr2interpstart,npr2interpend; // normally 0 and NPR2INTERP-1
289 int npr2interplist[MAXNPR]; // maximum is NPR2INTERP elements
290 
292 int npr2notinterpstart,npr2notinterpend; // normally 0 and -1
293 int npr2notinterplist[MAXNPR]; // maximum is NPR2INTERP elements
294 
296 int nprboundstart,nprboundend; // normally 0 and NPRBOUND-1
297 int nprboundlist[MAXNPR]; // maximum is NPRBOUND elements
298 
300 int nprfluxboundstart,nprfluxboundend; // normally 0 and NPRFLUXBOUND-1
301 int nprfluxboundlist[MAXNPR]; // maximum is NPRFLUXBOUND elements
302 
304 int nprdumpstart,nprdumpend; // normally 0 and NPRDUMP-1
305 int nprdumplist[MAXNPR]; // maximum is NPRDUMP elements
306 
308 int nprinvertstart,nprinvertend; // normally 0 and NPRINVERT-1
309 int nprinvertlist[MAXNPR]; // maximum is NPRINVERT elements
310 
311 
312 
313 int nogdetlist[MAXNPR]; // holds result for WHICHEOM stuff
314 
315 
316 int fluxloop[NDIM][NUMFLUXLOOPNUMBERS];
317 int emffluxloop[NDIM][NUMFLUXLOOPNUMBERS];
318 int Uconsloop[NUMFLUXLOOPNUMBERS];
319 int emfUconsloop[NUMFLUXLOOPNUMBERS];
320 int Uconsevolveloop[NUMFLUXLOOPNUMBERS];
321 int a_interporder[NUMINTERPS];
322 int *interporder;
323 int a_usedqarray[NUMINTERPS];
324 int *usedqarray;
325 
326 
328 //int dirglobal,locglobal,iglobal,jglobal,kglobal,iterglobal,interporfluxglobal;
329 
331 int special3dspc;
332 int dofull2pi;
333 
338 FTYPE global_npow,global_npow2,global_cpow2,global_rbr;
342 FTYPE global_x10,global_x20,global_fracdisk,global_fracjet,global_jetnu,global_rsjet,global_r0grid,global_r0jet,global_rjetend,global_r0disk,global_rdiskend;
346 FTYPE global_fracphi; //phi-extent measured in units of 2*PI, i.e. 0.25 means PI/2; only used if dofull2pi == 0
347 
348 
349 
350 int numbercpu[ 3+1 ];
351 
353 FTYPE nu,ss,ucrit,Ttpow,jetalpha;
354 
355 
357 FTYPE (*ptr_pressure_rho0_u[NUMEOSS])(FTYPE *EOSextra, FTYPE rho0, FTYPE u);
358 FTYPE (*ptr_compute_u_from_entropy[NUMEOSS])(FTYPE *EOSextra, FTYPE rho0, FTYPE entropy);
359 FTYPE (*ptr_u_rho0_p[NUMEOSS])(FTYPE *EOSextra, FTYPE rho0, FTYPE p);
360 FTYPE (*ptr_u_rho0_T[NUMEOSS])(FTYPE *EOSextra, FTYPE rho0, FTYPE T);
361 FTYPE (*ptr_dpdu_rho0_u[NUMEOSS])(FTYPE *EOSextra, FTYPE rho0, FTYPE u);
362 FTYPE (*ptr_dpdrho0_rho0_u[NUMEOSS])(FTYPE *EOSextra, FTYPE rho0, FTYPE u);
363 FTYPE (*ptr_cs2_compute[NUMEOSS])(FTYPE *EOSextra, FTYPE rho0, FTYPE u);
364 
365 FTYPE (*ptr_compute_entropy[NUMEOSS])(FTYPE *EOSextra, FTYPE rho0, FTYPE u);
366 FTYPE (*ptr_compute_dSdrho[NUMEOSS])(FTYPE *EOSextra, FTYPE rho0, FTYPE u);
367 FTYPE (*ptr_compute_dSdu[NUMEOSS])(FTYPE *EOSextra, FTYPE rho0, FTYPE u);
368 
369 FTYPE (*ptr_compute_specificentropy_wmrho0[NUMEOSS])(FTYPE *EOSextra, FTYPE rho0, FTYPE wmrho0);
370 FTYPE (*ptr_compute_dspecificSdrho_wmrho0[NUMEOSS])(FTYPE *EOSextra, FTYPE rho0, FTYPE wmrho0);
371 FTYPE (*ptr_compute_dspecificSdwmrho0_wmrho0[NUMEOSS])(FTYPE *EOSextra, FTYPE rho0, FTYPE wmrho0);
372 
373 FTYPE (*ptr_pressure_wmrho0[NUMEOSS])(FTYPE *EOSextra, FTYPE rho0, FTYPE wmrho0);
374 FTYPE (*ptr_compute_idwmrho0dp[NUMEOSS])(FTYPE *EOSextra, FTYPE rho0, FTYPE wmrho0);
375 FTYPE (*ptr_compute_idrho0dp[NUMEOSS])(FTYPE *EOSextra, FTYPE rho0, FTYPE wmrho0);
376 FTYPE (*ptr_compute_qdot[NUMEOSS])(FTYPE *EOSextra, FTYPE rho0, FTYPE u);
377 int (*ptr_compute_sources_EOS[NUMEOSS])(FTYPE *EOSextra, FTYPE *pr, struct of_geom *geom, struct of_state *q, FTYPE *Ui, FTYPE *dUother, FTYPE(*dUcomp)[NPR]);
378 void (*ptr_compute_allextras[NUMEOSS])(int justnum, FTYPE *EOSextra, FTYPE rho0, FTYPE u, int *numextrasreturn, FTYPE *extras);
379 int (*ptr_get_extrasprocessed[NUMEOSS])(int doall, FTYPE *EOSextra, FTYPE *pr, FTYPE *extras, FTYPE *processed);
380 FTYPE (*ptr_compute_temp[NUMEOSS])(FTYPE *EOSextra, FTYPE rho0, FTYPE u);
381 void (*ptr_compute_EOS_parms[NUMEOSS])(FTYPE (*EOSextra)[NSTORE2][NSTORE3][NUMEOSGLOBALS], FTYPE (*prim)[NSTORE2][NSTORE3][NPR]);
382 void (*ptr_compute_EOS_parms_full[NUMEOSS])(FTYPE (*EOSextra)[NSTORE2][NSTORE3][NUMEOSGLOBALS], FTYPE (*prim)[NSTORE2][NSTORE3][NPR]);
383 void (*ptr_store_EOS_parms[NUMEOSS])(int numparms, FTYPE *EOSextra, FTYPE *parlist);
384 void (*ptr_get_EOS_parms[NUMEOSS])(int*numparms, FTYPE *EOSextra, FTYPE *parlist);
385 void (*ptr_fix_primitive_eos_scalars[NUMEOSS])(FTYPE *EOSextra, FTYPE *pr);
386 void (*ptr_getall_forinversion[NUMEOSS])(int eomtype, int whichd, FTYPE *EOSextra, FTYPE quant1, FTYPE quant2, FTYPE *fun, FTYPE *dfunofrho, FTYPE *dfunofu);
387 
388 FTYPE SQRTMINNUMREPRESENT;
389 
390 FTYPE NUMEPSILONPOW23;
391 
393 int AVOIDADVANCESHIFTX1DN,AVOIDADVANCESHIFTX1UP,AVOIDADVANCESHIFTX2DN,AVOIDADVANCESHIFTX2UP,AVOIDADVANCESHIFTX3DN,AVOIDADVANCESHIFTX3UP,GLOBALBCMOVEDWITHACTIVESECTION;
394 
395 
396 #include "defs.user.h"
397 
398 
399 
400 int crapdebug;
401 
402 
403 
405 int ENDOPENMPSHAREDLIST;
406 
407 
408 #pragma omp threadprivate(OPENMPGLOBALPRIVATELIST)
Generated on Fri May 20 2016 15:52:32 for HARM by   doxygen 1.8.3.1