forked from pfhopkins/gizmo-public
-
Notifications
You must be signed in to change notification settings - Fork 0
/
begrun.c
2701 lines (2280 loc) · 112 KB
/
begrun.c
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
#include <mpi.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <sys/types.h>
#include <sys/stat.h>
#include <unistd.h>
#include <gsl/gsl_rng.h>
#include <ctype.h>
#include "allvars.h"
#include "proto.h"
#include "kernel.h"
/*! \file begrun.c
* \brief initial set-up of a simulation run
*
* This file contains various functions to initialize a simulation run. In
* particular, the parameterfile is read in and parsed, the initial
* conditions or restart files are read, and global variables are initialized
* to their proper values.
*/
/*! This function performs the initial set-up of the simulation. First, the
* parameterfile is set, then routines for setting units, reading
* ICs/restart-files are called, auxialiary memory is allocated, etc.
*/
/*
* This file was originally part of the GADGET3 code developed by
* Volker Springel. The code has been modified
* in part by Phil Hopkins (phopkins@caltech.edu) for GIZMO. The modifications
* mostly center on added functionality for new modules, elimination of unnecessary
* variables, implementing the DEVELOPER_MODE options, and re-organizing the read order
* to allow easier manipulation on restarts.
*/
void begrun(void)
{
struct global_data_all_processes all;
#ifdef _OPENMP
int tid;
#endif
if(ThisTask == 0)
{
printf("Running on %d MPI tasks.\n", NTask);
#ifdef _OPENMP
#pragma omp parallel private(tid)
{
#pragma omp master
printf("Using %d OpenMP threads\n", omp_get_num_threads());
tid = omp_get_thread_num();
}
#endif
printf("\nSize of particle structure %d [bytes]\n", (int) sizeof(struct particle_data));
printf("Size of hydro-cell structure %d [bytes]\n\n", (int) sizeof(struct gas_cell_data));
}
#ifdef CHIMES_TURB_DIFF_IONS
// Check that TURB_DIFF_METALS and TURB_DIFF_METALS_LOWORDER
// have also been switched on.
#ifndef TURB_DIFF_METALS
if (ThisTask == 0)
{
printf("ERROR: CHIMES_TURB_DIFF_IONS requires TURB_DIFF_METALS, but this is missing. Aborting.\n");
endrun(6572);
}
#endif // !(TURB_DIFF_METALS)
#ifndef TURB_DIFF_METALS_LOWORDER
if (ThisTask == 0)
{
printf("ERROR: CHIMES_TURB_DIFF_IONS requires TURB_DIFF_METALS_LOWORDER, but this is missing. Aborting.\n");
endrun(6573);
}
#endif // !(TURB_DIFF_METALS_LOWORDER)
#endif // CHIMES_TURB_DIFF_IONS
read_parameter_file(ParameterFile); /* ... read in parameters for this run */
mymalloc_init();
#ifdef DEBUG
write_pid_file();
enable_core_dumps_and_fpu_exceptions();
#endif
#ifdef GR_TABULATED_COSMOLOGY
#ifdef GR_TABULATED_COSMOLOGY_W
fwa_init();
#endif
#endif
set_units();
set_cosmo_factors_for_current_time();
All.Time = All.TimeBegin;
#ifdef COOLING
InitCool();
#endif
#ifdef BOX_PERIODIC
ewald_init();
#endif
#ifdef BOX_PERIODIC
boxSize = All.BoxSize;
boxHalf = 0.5 * All.BoxSize;
#endif
#ifdef BOX_LONG_X
boxSize_X = All.BoxSize * BOX_LONG_X;
boxHalf_X = 0.5 * boxSize_X;
#endif
#ifdef BOX_LONG_Y
boxSize_Y = All.BoxSize * BOX_LONG_Y;
boxHalf_Y = 0.5 * boxSize_Y;
#endif
#ifdef BOX_LONG_Z
boxSize_Z = All.BoxSize * BOX_LONG_Z;
boxHalf_Z = 0.5 * boxSize_Z;
#endif
#ifdef BOX_SHEARING
#ifdef BOX_LONG_X
Shearing_Box_Vel_Offset = BOX_SHEARING_Q * BOX_SHEARING_OMEGA_BOX_CENTER * All.BoxSize * BOX_LONG_X;
#else
Shearing_Box_Vel_Offset = BOX_SHEARING_Q * BOX_SHEARING_OMEGA_BOX_CENTER * All.BoxSize;
#endif
calc_shearing_box_pos_offset();
#endif
/* begin pre-definitions for special boundaries */
#if defined(BOX_REFLECT_X) || defined(BOX_REFLECT_Y) || defined(BOX_REFLECT_Z) || defined(BOX_OUTFLOW_X) || defined(BOX_OUTFLOW_Y) || defined(BOX_OUTFLOW_Z)
special_boundary_condition_xyz_def_reflect[0]=special_boundary_condition_xyz_def_reflect[1]=special_boundary_condition_xyz_def_reflect[2]=BOX_VALUE_FOR_NOTHING_SPECIAL_BOUNDARY_; /* sets arbitrary value code for 'nothing special' */
special_boundary_condition_xyz_def_outflow[0]=special_boundary_condition_xyz_def_outflow[1]=special_boundary_condition_xyz_def_outflow[2]=BOX_VALUE_FOR_NOTHING_SPECIAL_BOUNDARY_; /* sets arbitrary value code for 'nothing special' */
#if defined(BOX_REFLECT_X)
#if CHECK_IF_PREPROCESSOR_HAS_NUMERICAL_VALUE_(BOX_REFLECT_X)
special_boundary_condition_xyz_def_reflect[0] = BOX_REFLECT_X; /* set to user definition */
#else
special_boundary_condition_xyz_def_reflect[0] = 0; /* assume special boundary applies to both 'ends' of box, if not specified by user */
#endif
#endif
#if defined(BOX_REFLECT_Y)
#if CHECK_IF_PREPROCESSOR_HAS_NUMERICAL_VALUE_(BOX_REFLECT_Y)
special_boundary_condition_xyz_def_reflect[1] = BOX_REFLECT_Y; /* set to user definition */
#else
special_boundary_condition_xyz_def_reflect[1] = 0; /* assume special boundary applies to both 'ends' of box, if not specified by user */
#endif
#endif
#if defined(BOX_REFLECT_Z)
#if CHECK_IF_PREPROCESSOR_HAS_NUMERICAL_VALUE_(BOX_REFLECT_Z)
special_boundary_condition_xyz_def_reflect[2] = BOX_REFLECT_Z; /* set to user definition */
#else
special_boundary_condition_xyz_def_reflect[2] = 0; /* assume special boundary applies to both 'ends' of box, if not specified by user */
#endif
#endif
#if defined(BOX_OUTFLOW_X)
#if CHECK_IF_PREPROCESSOR_HAS_NUMERICAL_VALUE_(BOX_OUTFLOW_X)
special_boundary_condition_xyz_def_outflow[0] = BOX_OUTFLOW_X; /* set to user definition */
#else
special_boundary_condition_xyz_def_outflow[0] = 0; /* assume special boundary applies to both 'ends' of box, if not specified by user */
#endif
#endif
#if defined(BOX_OUTFLOW_Y)
#if CHECK_IF_PREPROCESSOR_HAS_NUMERICAL_VALUE_(BOX_OUTFLOW_Y)
special_boundary_condition_xyz_def_outflow[1] = BOX_OUTFLOW_Y; /* set to user definition */
#else
special_boundary_condition_xyz_def_outflow[1] = 0; /* assume special boundary applies to both 'ends' of box, if not specified by user */
#endif
#endif
#if defined(BOX_OUTFLOW_Z)
#if CHECK_IF_PREPROCESSOR_HAS_NUMERICAL_VALUE_(BOX_OUTFLOW_Z)
special_boundary_condition_xyz_def_outflow[2] = BOX_OUTFLOW_Z; /* set to user definition */
#else
special_boundary_condition_xyz_def_outflow[2] = 0; /* assume special boundary applies to both 'ends' of box, if not specified by user */
#endif
#endif
#endif /* end set of clauses to deal with causal flags for special boundary conditions */
random_generator = gsl_rng_alloc(gsl_rng_ranlxd1);
gsl_rng_set(random_generator, 42 + ThisTask); /* start-up seed */
set_random_numbers();
#ifdef PMGRID
if(RestartFlag != 3 && RestartFlag != 4) {long_range_init();}
#endif
#ifdef SUBFIND
GrNr = -1;
#endif
#ifdef EOS_TABULATED
int ierr = eos_init(All.EosTable);
if(ierr) {printf("error initializing the eos"); endrun(1);}
#endif
#ifdef EOS_TILLOTSON
tillotson_eos_init();
#endif
#ifdef NUCLEAR_NETWORK
network_init(All.EosSpecies, All.NetworkRates, All.NetworkPartFunc, All.NetworkMasses, All.NetworkWeakrates, &All.nd);
network_workspace_init(&All.nd, &All.nw);
#endif
#ifdef TURB_DRIVING
init_turb();
#endif
#if defined(DM_SIDM)
init_geofactor_table();
#endif
All.TimeLastRestartFile = CPUThisRun;
if(RestartFlag == 0 || RestartFlag == 2 || RestartFlag == 3 || RestartFlag == 4 || RestartFlag == 5 || RestartFlag == 6)
{
init(); /* ... read in initial model */
}
else
{
all = All; /* save global variables. (will be read from restart file) */
restart(RestartFlag); /* ... read restart file. Note: This also resets
all variables in the struct `All'.
However, during the run, some variables in the parameter
file are allowed to be changed, if desired. These need to
copied in the way below.
Note: All.PartAllocFactor is treated in restart() separately.
*/
set_random_numbers();
All.MinSizeTimestep = all.MinSizeTimestep;
All.MaxSizeTimestep = all.MaxSizeTimestep;
All.BufferSize = all.BufferSize;
All.TimeLimitCPU = all.TimeLimitCPU;
All.ResubmitOn = all.ResubmitOn;
All.SnapFormat = all.SnapFormat;
All.TimeBetSnapshot = all.TimeBetSnapshot;
All.TimeBetStatistics = all.TimeBetStatistics;
All.CpuTimeBetRestartFile = all.CpuTimeBetRestartFile;
All.ErrTolIntAccuracy = all.ErrTolIntAccuracy;
All.MinGasHsmlFractional = all.MinGasHsmlFractional;
All.MinGasTemp = all.MinGasTemp;
#ifdef CHIMES
All.ChimesThermEvolOn = all.ChimesThermEvolOn;
#endif
/* allow softenings to be modified during the run */
if(All.ComovingIntegrationOn)
{
All.SofteningGasMaxPhys = all.SofteningGasMaxPhys;
All.SofteningHaloMaxPhys = all.SofteningHaloMaxPhys;
All.SofteningDiskMaxPhys = all.SofteningDiskMaxPhys;
All.SofteningBulgeMaxPhys = all.SofteningBulgeMaxPhys;
All.SofteningStarsMaxPhys = all.SofteningStarsMaxPhys;
All.SofteningBndryMaxPhys = all.SofteningBndryMaxPhys;
}
All.SofteningGas = all.SofteningGas;
All.SofteningHalo = all.SofteningHalo;
All.SofteningDisk = all.SofteningDisk;
All.SofteningBulge = all.SofteningBulge;
All.SofteningStars = all.SofteningStars;
All.SofteningBndry = all.SofteningBndry;
All.MaxHsml = all.MaxHsml;
All.MaxRMSDisplacementFac = all.MaxRMSDisplacementFac;
All.ErrTolForceAcc = all.ErrTolForceAcc;
All.NumFilesPerSnapshot = all.NumFilesPerSnapshot;
All.NumFilesWrittenInParallel = all.NumFilesWrittenInParallel;
All.TreeDomainUpdateFrequency = all.TreeDomainUpdateFrequency;
All.OutputListOn = all.OutputListOn;
All.CourantFac = all.CourantFac;
All.OutputListLength = all.OutputListLength;
memcpy(All.OutputListTimes, all.OutputListTimes, sizeof(double) * All.OutputListLength);
memcpy(All.OutputListFlag, all.OutputListFlag, sizeof(char) * All.OutputListLength);
#ifdef GALSF
All.CritPhysDensity = all.CritPhysDensity;
All.MaxSfrTimescale = all.MaxSfrTimescale;
#endif
#ifdef SPHAV_CD10_VISCOSITY_SWITCH
All.ArtBulkViscConst = all.ArtBulkViscConst;
All.ViscosityAMin = all.ViscosityAMin;
All.ViscosityAMax = all.ViscosityAMax;
#endif
#ifdef TURB_DIFFUSION
All.TurbDiffusion_Coefficient = all.TurbDiffusion_Coefficient;
#endif
#ifdef SPHAV_ARTIFICIAL_CONDUCTIVITY
All.ArtCondConstant = all.ArtCondConstant;
#endif
#if defined(SPH_TP12_ARTIFICIAL_RESISTIVITY)
All.ArtMagDispConst = all.ArtMagDispConst;
#endif
#ifdef DIVBCLEANING_DEDNER
All.DivBcleanParabolicSigma = all.DivBcleanParabolicSigma;
All.DivBcleanHyperbolicSigma = all.DivBcleanHyperbolicSigma;
All.FastestWaveSpeed = 0.0;
All.FastestWaveDecay = 0.0;
#endif
#ifdef BLACK_HOLES
All.BlackHoleEddingtonFactor = all.BlackHoleEddingtonFactor;
All.SeedBlackHoleMass = all.SeedBlackHoleMass;
All.BlackHoleNgbFactor = all.BlackHoleNgbFactor;
All.BlackHoleMaxAccretionRadius = all.BlackHoleMaxAccretionRadius;
All.BlackHoleRadiativeEfficiency = all.BlackHoleRadiativeEfficiency;
All.BlackHoleFeedbackFactor = all.BlackHoleFeedbackFactor;
#if defined(BH_SEED_FROM_FOF) || defined(BH_SEED_FROM_LOCALGAS)
All.SeedBlackHoleMassSigma = all.SeedBlackHoleMassSigma;
All.SeedBlackHoleMinRedshift = all.SeedBlackHoleMinRedshift;
#ifdef BH_SEED_FROM_LOCALGAS
All.SeedBlackHolePerUnitMass = all.SeedBlackHolePerUnitMass;
#endif
#endif
#ifdef BH_ALPHADISK_ACCRETION
All.SeedAlphaDiskMass = all.SeedAlphaDiskMass;
#endif
#ifdef BH_SEED_FROM_FOF
All.MinFoFMassForNewSeed = all.MinFoFMassForNewSeed;
#endif
#if defined(BH_WIND_CONTINUOUS) || defined(BH_WIND_KICK) || defined(BH_WIND_SPAWN)
All.BAL_f_accretion = all.BAL_f_accretion;
All.BAL_v_outflow = all.BAL_v_outflow;
#endif
#if defined(SINGLE_STAR_FB_JETS)
All.BAL_f_launch_v = all.BAL_f_launch_v;
#endif
#if defined(BH_COSMIC_RAYS)
All.BH_CosmicRay_Injection_Efficiency = all.BH_CosmicRay_Injection_Efficiency;
#endif
#ifdef BH_WIND_SPAWN
All.BAL_internal_temperature = all.BAL_internal_temperature;
All.BAL_wind_particle_mass = all.BAL_wind_particle_mass; // dangeous to change this, as it is also part of the merger criterion!
#endif
#ifdef BH_PHOTONMOMENTUM
All.BH_Rad_MomentumFactor = all.BH_Rad_MomentumFactor;
#endif
#endif // blackholes
#ifdef GALSF_FB_FIRE_RT_LOCALRP
All.RP_Local_Momentum_Renormalization = all.RP_Local_Momentum_Renormalization;
#endif
#ifdef GALSF_FB_FIRE_RT_HIIHEATING
All.HIIRegion_fLum_Coupled = all.HIIRegion_fLum_Coupled;
#endif
#ifdef RT_LEBRON
All.PhotonMomentum_Coupled_Fraction = all.PhotonMomentum_Coupled_Fraction;
#endif
#ifdef GALSF_FB_FIRE_RT_LONGRANGE
All.PhotonMomentum_fUV = all.PhotonMomentum_fUV;
All.PhotonMomentum_fOPT = all.PhotonMomentum_fOPT;
#endif
#ifdef COSMIC_RAY_FLUID
#if (CRFLUID_DIFFUSION_MODEL == 0)
All.CosmicRayDiffusionCoeff = all.CosmicRayDiffusionCoeff;
#endif
#endif
#ifdef GALSF_FB_FIRE_AGE_TRACERS
All.AgeTracerRateNormalization = all.AgeTracerRateNormalization;
#ifdef GALSF_FB_FIRE_AGE_TRACERS_CUSTOM
strcpy(All.AgeTracerListFilename, all.AgeTracerListFilename);
#else
All.AgeTracerBinStart = all.AgeTracerBinStart;
All.AgeTracerBinEnd = all.AgeTracerBinEnd;
#endif
#endif
#ifdef CR_DYNAMICAL_INJECTION_IN_SNE
All.CosmicRay_SNeFraction = all.CosmicRay_SNeFraction;
#endif
#ifdef GR_TABULATED_COSMOLOGY
All.DarkEnergyConstantW = all.DarkEnergyConstantW;
#endif
All.MaxNumNgbDeviation = all.MaxNumNgbDeviation;
#ifdef AGS_HSML_CALCULATION_IS_ACTIVE
/* Allow the tolerance over the number of neighbours to vary during the run:
If it was initially set to a very strict value, convergence in ngb-iteration may at some point fail */
All.AGS_MaxNumNgbDeviation = all.AGS_MaxNumNgbDeviation;
#endif
strcpy(All.ResubmitCommand, all.ResubmitCommand);
strcpy(All.OutputListFilename, all.OutputListFilename);
strcpy(All.OutputDir, all.OutputDir);
strcpy(All.RestartFile, all.RestartFile);
/*
strcpy(All.EnergyFile, all.EnergyFile);
strcpy(All.InfoFile, all.InfoFile);
strcpy(All.CpuFile, all.CpuFile);
strcpy(All.TimingsFile, all.TimingsFile);
strcpy(All.TimebinFile, all.TimebinFile);
*/
strcpy(All.SnapshotFileBase, all.SnapshotFileBase);
#ifdef COOL_GRACKLE
strcpy(All.GrackleDataFile, all.GrackleDataFile);
#endif
#ifdef EOS_TABULATED
strcpy(All.EosTable, all.EosTable);
#endif
#ifdef NUCLEAR_NETWORK
strcpy(All.EosSpecies, all.EosSpecies);
strcpy(All.NetworkRates, all.NetworkRates);
strcpy(All.NetworkPartFunc, all.NetworkPartFunc);
strcpy(All.NetworkMasses, all.NetworkMasses);
strcpy(All.NetworkWeakrates, all.NetworkWeakrates);
All.nd = all.nd;
All.nw = all.nw;
All.NetworkTempThreshold = all.NetworkTempThreshold;
#endif
if(All.TimeMax != all.TimeMax) {readjust_timebase(All.TimeMax, all.TimeMax);}
}
#ifdef GALSF_EFFECTIVE_EQS
init_clouds();
#endif
char contfname[1000];
sprintf(contfname, "%scont", All.OutputDir);
unlink(contfname);
open_outputfiles();
#ifdef PMGRID
long_range_init_regionsize();
#endif
reconstruct_timebins();
#ifndef BOX_SHEARING
#if (NUMDIMS==2)
int i;
for(i = 0; i < NumPart; i++)
{
P[i].Pos[2] = P[i].Vel[2] = P[i].GravAccel[2] = 0;
if(P[i].Type == 0) {SphP[i].VelPred[2] = SphP[i].HydroAccel[2] = 0;}
}
#endif
#endif
#ifdef RT_RAD_PRESSURE_OUTPUT
{int i; for(i=0;i<NumPart;i++) {int k; for(k=0;k<3;k++) {SphP[i].Rad_Accel[k]=0;}}}
#endif
#ifdef RADTRANSFER
#if defined(RT_EVOLVE_INTENSITIES)
rt_init_intensity_directions();
#endif
#if defined(RT_DIFFUSION_CG)
All.Radiation_Ti_begstep = 0;
#endif
#ifdef RT_CHEM_PHOTOION
rt_get_sigma();
#endif
#endif
#if defined(RADTRANSFER) || defined(RT_USE_GRAVTREE)
rt_define_effective_frequencies_in_bands();
#endif
#ifdef RADTRANSFER
rt_set_simple_inits(RestartFlag);
#endif
if(All.ComovingIntegrationOn) {init_drift_table();}
if(RestartFlag == 2)
{All.Ti_nextoutput = find_next_outputtime(All.Ti_Current + 100);}
else if(RestartFlag == 1)
{All.Ti_nextoutput = find_next_outputtime(All.Ti_Current + 1);}
else
{All.Ti_nextoutput = find_next_outputtime(All.Ti_Current);}
All.TimeLastRestartFile = CPUThisRun;
}
/*! Computes conversion factors between internal code units and the cgs-system
*/
void set_units(void)
{
/* convert some physical input parameters to internal units */
if(All.G <= 0) {All.G = GRAVITY_G_CGS * UNIT_MASS_IN_CGS / (UNIT_LENGTH_IN_CGS * UNIT_VEL_IN_CGS*UNIT_VEL_IN_CGS);}
#ifdef GR_TABULATED_COSMOLOGY_G
All.Gini = All.G;
All.G = All.Gini * dGfak(All.TimeBegin);
#endif
All.Hubble_H0_CodeUnits = H0_CGS * UNIT_TIME_IN_CGS;
if(ThisTask == 0)
{
printf("\nCode units to be used: make sure you check these are correct! \n");
printf(" Hubble H0 (internal units) = %g \n", All.Hubble_H0_CodeUnits);
printf(" Gravity G (internal units) = %g \n", All.G);
printf(" unit Mass in g = %g \n", UNIT_MASS_IN_CGS);
printf(" unit Length in cm = %g \n", UNIT_LENGTH_IN_CGS);
printf(" unit Time in s = %g \n", UNIT_TIME_IN_CGS);
printf(" unit Velocity in cm/s = %g \n", UNIT_VEL_IN_CGS);
printf(" unit Energy in erg = %g \n", UNIT_ENERGY_IN_CGS);
printf(" unit Density in g/cm^3 = %g \n", UNIT_DENSITY_IN_CGS);
printf(" unit Pressure in erg/cm^3 = %g \n", UNIT_PRESSURE_IN_CGS);
printf(" unit Luminosity in erg/s = %g \n", UNIT_LUM_IN_CGS);
printf(" unit Flux in erg/s/cm^2 = %g \n", UNIT_FLUX_IN_CGS);
printf(" unit B[internal] in gauss = %g \n", UNIT_B_IN_GAUSS);
printf("\n");
}
#if !defined(SELFGRAVITY_OFF)
if(ThisTask==0)
{
printf("\nThe chosen interaction and softening kernel function is =%d \n",(int)KERNEL_FUNCTION);
printf(" Gravity will be exactly Newtonian beyond =%g times the given Plummer-equivalent softenings \n",(1./KERNEL_FAC_FROM_FORCESOFT_TO_PLUMMER));
printf("\n");
}
#endif
double meanweight = 4.0 / (1 + 3 * HYDROGEN_MASSFRAC); /* assumes fully-atomic otherwise */
#ifdef COOL_LOW_TEMPERATURES
meanweight = 1. / ( HYDROGEN_MASSFRAC*0.5 + (1-HYDROGEN_MASSFRAC)/4. + 1./(16.+12.)); /* assumes fully-molecular if low-temp cooling enabled */
#endif
All.MinEgySpec = All.MinGasTemp / (meanweight * (GAMMA_DEFAULT-1) * U_TO_TEMP_UNITS);
#if defined(GALSF)
/* for historical reasons, we need to convert to "All.MaxSfrTimescale", defined as the SF timescale in code units at the critical physical
density given above. use the dimensionless SfEffPerFreeFall (which has been read in) to calculate this. This must be done -BEFORE- calling set_units_sfr) */
#ifndef GALSF_EFFECTIVE_EQS
All.MaxSfrTimescale = (1/All.MaxSfrTimescale) * sqrt(3.*M_PI / (32. * All.G * (All.CritPhysDensity / UNIT_DENSITY_IN_NHCGS)));
#endif
set_units_sfr();
#endif
#ifdef DM_FUZZY
/* For Schroedinger equation: this encodes the coefficient with the mass of the particle: units vel*L = hbar / particle_mass. This is the key variable used throughout */
All.ScalarField_hbar_over_mass = 591569.0 / ((double)All.ScalarField_hbar_over_mass * UNIT_VEL_IN_CGS * UNIT_LENGTH_IN_CGS);
#endif
#if defined(CONDUCTION_SPITZER) || defined(VISCOSITY_BRAGINSKII)
/* Note: Because we replace \nabla(T) in the conduction equation with \nabla(u), our conduction coefficient is not the usual kappa, but
* rather kappa*(gamma-1)*mu/kB. We therefore need to multiply with another factor of (meanweight_ion / k_B * (gamma-1)) */
double meanweight_ion = 4.0 / (8 - 5 * (1 - HYDROGEN_MASSFRAC)); /* mean weight in code units, assuming full ionization */
double u_to_temp = meanweight_ion * (GAMMA_DEFAULT-1.) * U_TO_TEMP_UNITS; /* for full ionization, assume gas has a monatomic ideal eos gamma=5/3 */
/* Kappa_Spitzer definition taken from Zakamska & Narayan 2003 ( ApJ 582:162-169, Eq. (5) ) */
double coulomb_log = 37.8; // Sarazin value (recommendation from PIC calculations) //
double coefficient = (1.84e-5/coulomb_log) * pow(u_to_temp,3.5) * ((UNIT_TIME_IN_CGS*UNIT_TIME_IN_CGS*UNIT_TIME_IN_CGS) / (UNIT_LENGTH_IN_CGS*UNIT_MASS_IN_CGS)); // ok, this multiplied by the specific energy (u_code)^(3/2) gives the diffusity of u_code, as needed (density term is included in said diffusivity)
#ifdef CONDUCTION_SPITZER
All.ConductionCoeff *= coefficient;
#endif
#ifdef VISCOSITY_BRAGINSKII
All.ShearViscosityCoeff *= coefficient * 0.636396*sqrt(ELECTRONMASS_CGS/(PROTONMASS_CGS*meanweight_ion)); // the viscosity coefficient eta is identical in these units up to the order-unity constant, and multiplied by sqrt[m_electron/m_ion] //
All.BulkViscosityCoeff = 0; // no bulk viscosity in the Braginskii-Spitzer formulation //
#endif
/* factor used for determining saturation */
All.ElectronFreePathFactor = 8 * pow(3.0, 1.5) * pow((GAMMA_DEFAULT-1), 2) / pow(3 + 5 * HYDROGEN_MASSFRAC, 2)
/ (1 + HYDROGEN_MASSFRAC) / sqrt(M_PI) / coulomb_log * pow(PROTONMASS_CGS, 3) / pow(ELECTRONCHARGE_CGS, 4) / (UNIT_DENSITY_IN_CGS) * pow(UNIT_SPECEGY_IN_CGS, 2);
/* If the above value is multiplied with u^2/rho in code units (with rho being the physical density), then
* one gets the electron mean free path in centimeters. Since we want to compare this with another length
* scale in code units, we now add an additional factor to convert back to code units. */
All.ElectronFreePathFactor /= UNIT_LENGTH_IN_CGS;
#endif
}
/*! This function opens various log-files that report on the status and performance of the simulation.
On restart from restart-files, (start-option 1), the code will append to these files. */
void open_outputfiles(void)
{
char mode[2], buf[200];
if(RestartFlag == 0) {strcpy(mode, "w");} else {strcpy(mode, "a");}
if(ThisTask == 0) {mkdir(All.OutputDir, 02755);}
MPI_Barrier(MPI_COMM_WORLD);
#ifdef BLACK_HOLES /* Note: This is done by everyone [all tasks can write to these log-files], even if it might be empty */
if(ThisTask == 0) {sprintf(buf, "%sblackhole_details", All.OutputDir); mkdir(buf, 02755);}
MPI_Barrier(MPI_COMM_WORLD);
#if !defined(IO_REDUCED_MODE) || defined(BH_OUTPUT_MOREINFO)
sprintf(buf, "%sblackhole_details/blackhole_details_%d.txt", All.OutputDir, ThisTask);
if(!(FdBlackHolesDetails = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
#endif
#ifdef OUTPUT_SINK_ACCRETION_HIST
sprintf(buf, "%sblackhole_details/bhswallow_%d.txt", All.OutputDir, ThisTask);
if(!(FdBhSwallowDetails = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
#endif
#ifdef OUTPUT_SINK_FORMATION_PROPS
sprintf(buf, "%sblackhole_details/bhformation_%d.txt", All.OutputDir, ThisTask);
if(!(FdBhFormationDetails = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
#endif
#ifdef BH_OUTPUT_MOREINFO
sprintf(buf, "%sblackhole_details/bhmergers_%d.txt", All.OutputDir, ThisTask);
if(!(FdBhMergerDetails = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
#ifdef BH_WIND_KICK
sprintf(buf, "%sblackhole_details/bhwinds_%d.txt", All.OutputDir, ThisTask);
if(!(FdBhWindDetails = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
#endif
#endif // bh-output-more-info if
#endif // black-holes if
if(ThisTask != 0) {return;} /* only the root processors writes to the log files listed below */
sprintf(buf, "%s%s", All.OutputDir, "cpu.txt");
if(!(FdCPU = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
#ifndef IO_REDUCED_MODE
sprintf(buf, "%s%s", All.OutputDir, "timebin.txt");
if(!(FdTimebin = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
sprintf(buf, "%s%s", All.OutputDir, "info.txt");
if(!(FdInfo = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
sprintf(buf, "%s%s", All.OutputDir, "energy.txt");
if(!(FdEnergy = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
sprintf(buf, "%s%s", All.OutputDir, "timings.txt");
if(!(FdTimings = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
sprintf(buf, "%s%s", All.OutputDir, "balance.txt");
if(!(FdBalance = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
fprintf(FdBalance, "\n");
fprintf(FdBalance, "Treewalk1 = '%c' / '%c'\n", CPU_Symbol[CPU_TREEWALK1], CPU_SymbolImbalance[CPU_TREEWALK1]);
fprintf(FdBalance, "Treewalk2 = '%c' / '%c'\n", CPU_Symbol[CPU_TREEWALK2], CPU_SymbolImbalance[CPU_TREEWALK2]);
fprintf(FdBalance, "Treewait1 = '%c' / '%c'\n", CPU_Symbol[CPU_TREEWAIT1], CPU_SymbolImbalance[CPU_TREEWAIT1]);
fprintf(FdBalance, "Treewait2 = '%c' / '%c'\n", CPU_Symbol[CPU_TREEWAIT2], CPU_SymbolImbalance[CPU_TREEWAIT2]);
fprintf(FdBalance, "Treesend = '%c' / '%c'\n", CPU_Symbol[CPU_TREESEND], CPU_SymbolImbalance[CPU_TREESEND]);
fprintf(FdBalance, "Treerecv = '%c' / '%c'\n", CPU_Symbol[CPU_TREERECV], CPU_SymbolImbalance[CPU_TREERECV]);
fprintf(FdBalance, "Treebuild = '%c' / '%c'\n", CPU_Symbol[CPU_TREEBUILD], CPU_SymbolImbalance[CPU_TREEBUILD]);
fprintf(FdBalance, "Treehmaxupdate = '%c' / '%c'\n", CPU_Symbol[CPU_TREEHMAXUPDATE], CPU_SymbolImbalance[CPU_TREEHMAXUPDATE]);
fprintf(FdBalance, "Treemisc = '%c' / '%c'\n", CPU_Symbol[CPU_TREEMISC], CPU_SymbolImbalance[CPU_TREEMISC]);
fprintf(FdBalance, "Domain decomp = '%c' / '%c'\n", CPU_Symbol[CPU_DOMAIN], CPU_SymbolImbalance[CPU_DOMAIN]);
fprintf(FdBalance, "Peano-Hilbert = '%c' / '%c'\n", CPU_Symbol[CPU_PEANO], CPU_SymbolImbalance[CPU_PEANO]);
fprintf(FdBalance, "Density compute= '%c' / '%c'\n", CPU_Symbol[CPU_DENSCOMPUTE], CPU_SymbolImbalance[CPU_DENSCOMPUTE]);
fprintf(FdBalance, "Density imbal = '%c' / '%c'\n", CPU_Symbol[CPU_DENSWAIT], CPU_SymbolImbalance[CPU_DENSWAIT]);
fprintf(FdBalance, "Density commu = '%c' / '%c'\n", CPU_Symbol[CPU_DENSCOMM], CPU_SymbolImbalance[CPU_DENSCOMM]);
fprintf(FdBalance, "Density misc = '%c' / '%c'\n", CPU_Symbol[CPU_DENSMISC], CPU_SymbolImbalance[CPU_DENSMISC]);
fprintf(FdBalance, "Hydro compute = '%c' / '%c'\n", CPU_Symbol[CPU_HYDCOMPUTE], CPU_SymbolImbalance[CPU_HYDCOMPUTE]);
fprintf(FdBalance, "Hydro imbalance= '%c' / '%c'\n", CPU_Symbol[CPU_HYDWAIT], CPU_SymbolImbalance[CPU_HYDWAIT]);
fprintf(FdBalance, "Hydro comm = '%c' / '%c'\n", CPU_Symbol[CPU_HYDCOMM], CPU_SymbolImbalance[CPU_HYDCOMM]);
fprintf(FdBalance, "Hydro misc = '%c' / '%c'\n", CPU_Symbol[CPU_HYDMISC], CPU_SymbolImbalance[CPU_HYDMISC]);
fprintf(FdBalance, "Drifts = '%c' / '%c'\n", CPU_Symbol[CPU_DRIFT], CPU_SymbolImbalance[CPU_DRIFT]);
fprintf(FdBalance, "Kicks = '%c' / '%c'\n", CPU_Symbol[CPU_TIMELINE], CPU_SymbolImbalance[CPU_TIMELINE]);
fprintf(FdBalance, "Potential = '%c' / '%c'\n", CPU_Symbol[CPU_POTENTIAL], CPU_SymbolImbalance[CPU_POTENTIAL]);
fprintf(FdBalance, "PM-gravity = '%c' / '%c'\n", CPU_Symbol[CPU_MESH], CPU_SymbolImbalance[CPU_MESH]);
fprintf(FdBalance, "Snapshot dump = '%c' / '%c'\n", CPU_Symbol[CPU_SNAPSHOT], CPU_SymbolImbalance[CPU_SNAPSHOT]);
fprintf(FdBalance, "Blackhole = '%c' / '%c'\n", CPU_Symbol[CPU_BLACKHOLES], CPU_SymbolImbalance[CPU_BLACKHOLES]);
fprintf(FdBalance, "Cooling & SFR = '%c' / '%c'\n", CPU_Symbol[CPU_COOLINGSFR], CPU_SymbolImbalance[CPU_COOLINGSFR]);
fprintf(FdBalance, "Coolimbal check= '%c' / '%c'\n", CPU_Symbol[CPU_COOLSFRIMBAL], CPU_SymbolImbalance[CPU_COOLSFRIMBAL]);
fprintf(FdBalance, "FoF & subfind = '%c' / '%c'\n", CPU_Symbol[CPU_FOF], CPU_SymbolImbalance[CPU_FOF]);
fprintf(FdBalance, "Grain/PIC part = '%c' / '%c'\n", CPU_Symbol[CPU_DRAGFORCE], CPU_SymbolImbalance[CPU_DRAGFORCE]);
fprintf(FdBalance, "Mech/Thermal FB= '%c' / '%c'\n", CPU_Symbol[CPU_SNIIHEATING], CPU_SymbolImbalance[CPU_SNIIHEATING]);
fprintf(FdBalance, "HII-module = '%c' / '%c'\n", CPU_Symbol[CPU_HIIHEATING], CPU_SymbolImbalance[CPU_HIIHEATING]);
fprintf(FdBalance, "Local wind kick= '%c' / '%c'\n", CPU_Symbol[CPU_LOCALWIND], CPU_SymbolImbalance[CPU_LOCALWIND]);
fprintf(FdBalance, "RHD-nonfluxops = '%c' / '%c'\n", CPU_Symbol[CPU_RTNONFLUXOPS], CPU_SymbolImbalance[CPU_RTNONFLUXOPS]);
fprintf(FdBalance, "AGS-nongas-comp= '%c' / '%c'\n", CPU_Symbol[CPU_AGSDENSCOMPUTE], CPU_SymbolImbalance[CPU_AGSDENSCOMPUTE]);
fprintf(FdBalance, "AGS-imbal = '%c' / '%c'\n", CPU_Symbol[CPU_AGSDENSWAIT], CPU_SymbolImbalance[CPU_AGSDENSWAIT]);
fprintf(FdBalance, "AGS-comm = '%c' / '%c'\n", CPU_Symbol[CPU_AGSDENSCOMM], CPU_SymbolImbalance[CPU_AGSDENSCOMM]);
fprintf(FdBalance, "AGS-misc = '%c' / '%c'\n", CPU_Symbol[CPU_AGSDENSMISC], CPU_SymbolImbalance[CPU_AGSDENSMISC]);
fprintf(FdBalance, "DynDiffusn-comp= '%c' / '%c'\n", CPU_Symbol[CPU_DYNDIFFCOMPUTE], CPU_SymbolImbalance[CPU_DYNDIFFCOMPUTE]);
fprintf(FdBalance, "DynDiffusn-imbl= '%c' / '%c'\n", CPU_Symbol[CPU_DYNDIFFWAIT], CPU_SymbolImbalance[CPU_DYNDIFFWAIT]);
fprintf(FdBalance, "DynDiffusn-comm= '%c' / '%c'\n", CPU_Symbol[CPU_DYNDIFFCOMM], CPU_SymbolImbalance[CPU_DYNDIFFCOMM]);
fprintf(FdBalance, "DynDiffusn-misc= '%c' / '%c'\n", CPU_Symbol[CPU_DYNDIFFMISC], CPU_SymbolImbalance[CPU_DYNDIFFMISC]);
fprintf(FdBalance, "MultiDiff-comp = '%c' / '%c'\n", CPU_Symbol[CPU_IMPROVDIFFCOMPUTE], CPU_SymbolImbalance[CPU_IMPROVDIFFCOMPUTE]);
fprintf(FdBalance, "MultiDiff-imbl = '%c' / '%c'\n", CPU_Symbol[CPU_IMPROVDIFFWAIT], CPU_SymbolImbalance[CPU_IMPROVDIFFWAIT]);
fprintf(FdBalance, "MultiDiff-comm = '%c' / '%c'\n", CPU_Symbol[CPU_IMPROVDIFFCOMM], CPU_SymbolImbalance[CPU_IMPROVDIFFCOMM]);
fprintf(FdBalance, "MultiDiff-misc = '%c' / '%c'\n", CPU_Symbol[CPU_IMPROVDIFFMISC], CPU_SymbolImbalance[CPU_IMPROVDIFFMISC]);
fprintf(FdBalance, "Miscellaneous = '%c' / '%c'\n", CPU_Symbol[CPU_MISC], CPU_SymbolImbalance[CPU_MISC]);
fprintf(FdBalance, "\n");
#endif
#ifdef GALSF
sprintf(buf, "%s%s", All.OutputDir, "sfr.txt");
if(!(FdSfr = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
#endif
#ifdef GALSF_FB_FIRE_RT_LOCALRP
sprintf(buf, "%s%s", All.OutputDir, "MomWinds.txt");
if(!(FdMomWinds = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
#endif
#ifdef GALSF_FB_FIRE_RT_HIIHEATING
sprintf(buf, "%s%s", All.OutputDir, "HIIheating.txt");
if(!(FdHIIHeating = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
#endif
#ifdef GALSF_FB_MECHANICAL
sprintf(buf, "%s%s", All.OutputDir, "SNeIIheating.txt");
if(!(FdSneIIHeating = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
#endif
#if defined(RT_CHEM_PHOTOION) && !defined(IO_REDUCED_MODE)
sprintf(buf, "%s%s", All.OutputDir, "rt_photoion_chem.txt");
if(!(FdRad = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
#endif
#ifdef BLACK_HOLES
sprintf(buf, "%s%s", All.OutputDir, "blackholes.txt");
if(!(FdBlackHoles = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
#endif
#if defined(TURB_DRIVING) && !defined(IO_REDUCED_MODE)
sprintf(buf, "%s%s", All.OutputDir, "turb.txt");
if(!(FdTurb = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
#endif
#if defined(GR_TABULATED_COSMOLOGY) && !defined(IO_REDUCED_MODE)
sprintf(buf, "%s%s", All.OutputDir, "darkenergy.txt");
if(!(FdDE = fopen(buf, mode))) {printf("error in opening file '%s'\n", buf); endrun(1);}
else if(RestartFlag == 0)
{
fprintf(FdDE, "nstep time H(a) ");
#ifndef GR_TABULATED_COSMOLOGY_W
fprintf(FdDE, "w0 Omega_L ");
#else
fprintf(FdDE, "w(a) Omega_L ");
#endif
#ifdef GR_TABULATED_COSMOLOGY_G
fprintf(FdDE, "dH dG ");
#endif
fprintf(FdDE, "\n"); fflush(FdDE);
}
#endif
}
/*! This function parses the parameterfile in a simple way. Each paramater is
* defined by a keyword (`tag'), and can be either of type douple, int, or
* character string. The routine makes sure that each parameter appears
* exactly once in the parameterfile, otherwise error messages are
* produced that complain about the missing parameters.
*/
void read_parameter_file(char *fname)
{
#define REAL 1
#define STRING 2
#define INT 3
#define MAXTAGS 300
FILE *fd, *fdout;
char buf[200], buf1[200], buf2[200], buf3[400], tag[MAXTAGS][50], alternate_tag[MAXTAGS][50];
int i, j, nt, id[MAXTAGS], pnum, errorFlag = 0;
void *addr[MAXTAGS];
double safe_memorypertask = mpi_report_comittable_memory(0,0); /* used for some checks below */
#ifdef CHIMES
double Tdust_buf, Tmol_buf, relTol_buf, absTol_buf, expTol_buf, z_reion_buf;
#endif
if(sizeof(long long) != 8)
{
if(ThisTask == 0) {printf("\nType `long long' is not 64 bit on this platform. Stopping.\n\n");}
endrun(0);
}
if(sizeof(int) != 4)
{
if(ThisTask == 0) {printf("\nType `int' is not 32 bit on this platform. Stopping.\n\n");}
endrun(0);
}
if(sizeof(float) != 4)
{
if(ThisTask == 0) {printf("\nType `float' is not 32 bit on this platform. Stopping.\n\n");}
endrun(0);
}
if(sizeof(double) != 8)
{
if(ThisTask == 0) {printf("\nType `double' is not 64 bit on this platform. Stopping.\n\n");}
endrun(0);
}
if(ThisTask == 0) /* read parameter file on process 0 */
{
nt = 0;
for(j=0;j<MAXTAGS;j++) {strcpy(alternate_tag[j], "-");}
strcpy(tag[nt], "InitCondFile");
strcpy(alternate_tag[nt], "Initial_Conditions_File");
addr[nt] = All.InitCondFile;
id[nt++] = STRING;
strcpy(tag[nt], "OutputDir");
strcpy(alternate_tag[nt], "Output_Directory");
addr[nt] = All.OutputDir;
id[nt++] = STRING;
strcpy(tag[nt], "SnapshotFileBase");
strcpy(alternate_tag[nt], "Snapshot_Filename_Base");
addr[nt] = All.SnapshotFileBase;
id[nt++] = STRING;
strcpy(tag[nt], "RestartFile");
strcpy(alternate_tag[nt], "Restart_Filename_Base");
addr[nt] = All.RestartFile;
id[nt++] = STRING;
#ifdef DEVELOPER_MODE
strcpy(tag[nt], "ResubmitOn");
strcpy(alternate_tag[nt], "Use_Automatic_Shell_Resubmission");
addr[nt] = &All.ResubmitOn;
id[nt++] = INT;
strcpy(tag[nt], "ResubmitCommand");
strcpy(alternate_tag[nt], "Shell_Resubmission_Command");
addr[nt] = All.ResubmitCommand;
id[nt++] = STRING;
#endif
strcpy(tag[nt], "OutputListFilename");
strcpy(alternate_tag[nt], "Snapshot_Times_Table_Filename");
addr[nt] = All.OutputListFilename;
id[nt++] = STRING;
strcpy(tag[nt], "OutputListOn");
strcpy(alternate_tag[nt], "Use_Tabulated_Snapshot_Times");
addr[nt] = &All.OutputListOn;
id[nt++] = INT;
strcpy(tag[nt], "Omega0");
strcpy(alternate_tag[nt], "Omega_Matter");
addr[nt] = &All.OmegaMatter;
id[nt++] = REAL;
strcpy(tag[nt], "OmegaBaryon");
strcpy(alternate_tag[nt], "Omega_Baryon");
addr[nt] = &All.OmegaBaryon;
id[nt++] = REAL;
strcpy(tag[nt], "OmegaLambda");
strcpy(alternate_tag[nt], "Omega_Lambda");
addr[nt] = &All.OmegaLambda;
id[nt++] = REAL;
strcpy(tag[nt], "OmegaRadiation");
strcpy(alternate_tag[nt], "Omega_Radiation");
addr[nt] = &All.OmegaRadiation;
id[nt++] = REAL;
strcpy(tag[nt], "HubbleParam");
strcpy(alternate_tag[nt], "Hubble_Param_Little_h");
addr[nt] = &All.HubbleParam;
id[nt++] = REAL;
strcpy(tag[nt], "BoxSize");
strcpy(alternate_tag[nt], "Box_Size_In_Code_Units");
addr[nt] = &All.BoxSize;
id[nt++] = REAL;
strcpy(tag[nt], "MaxMemSize");
strcpy(alternate_tag[nt], "Max_Memory_Per_MPI_Task_in_MB");
addr[nt] = &All.MaxMemSize;
id[nt++] = INT;
strcpy(tag[nt], "TimeOfFirstSnapshot");
strcpy(alternate_tag[nt], "Simulation_Time_of_First_Snapshot");
addr[nt] = &All.TimeOfFirstSnapshot;
id[nt++] = REAL;
strcpy(tag[nt], "CpuTimeBetRestartFile");
strcpy(alternate_tag[nt], "Walltime_in_Seconds_Between_Restartfiles");
addr[nt] = &All.CpuTimeBetRestartFile;
id[nt++] = REAL;
#ifdef DEVELOPER_MODE
strcpy(tag[nt], "TimeBetStatistics");
strcpy(alternate_tag[nt], "Time_Between_Internal_Diagnostic_Statistics");
addr[nt] = &All.TimeBetStatistics;
id[nt++] = REAL;
#endif
strcpy(tag[nt], "TimeBegin");
strcpy(alternate_tag[nt], "Initial_Simulation_Time");
addr[nt] = &All.TimeBegin;
id[nt++] = REAL;
strcpy(tag[nt], "TimeMax");
strcpy(alternate_tag[nt], "Final_Simulation_Time");
addr[nt] = &All.TimeMax;
id[nt++] = REAL;
strcpy(tag[nt], "TimeBetSnapshot");
strcpy(alternate_tag[nt], "ScaleFac_Between_Snapshots");
addr[nt] = &All.TimeBetSnapshot;
id[nt++] = REAL;
strcpy(tag[nt], "UnitVelocity_in_cm_per_s");
strcpy(alternate_tag[nt], "UnitVelocity_in_centimeterspersecond");
addr[nt] = &All.UnitVelocity_in_cm_per_s;
id[nt++] = REAL;
strcpy(tag[nt], "UnitLength_in_cm");
strcpy(alternate_tag[nt], "UnitLength_in_centimeters");
addr[nt] = &All.UnitLength_in_cm;
id[nt++] = REAL;
strcpy(tag[nt], "UnitMass_in_g");
strcpy(alternate_tag[nt], "UnitMass_in_grams");
addr[nt] = &All.UnitMass_in_g;
id[nt++] = REAL;
#ifdef MAGNETIC
strcpy(tag[nt], "UnitMagneticField_in_gauss");
strcpy(alternate_tag[nt], "UnitMagneticField_in_Gauss_for_ICSnapshotIO");
addr[nt] = &All.UnitMagneticField_in_gauss;
id[nt++] = REAL;
#endif
strcpy(tag[nt], "TreeDomainUpdateFrequency");
strcpy(alternate_tag[nt], "TreeRebuild_ActiveFraction");
addr[nt] = &All.TreeDomainUpdateFrequency;
id[nt++] = REAL;
#ifdef DEVELOPER_MODE
strcpy(tag[nt], "ErrTolIntAccuracy");
addr[nt] = &All.ErrTolIntAccuracy;
id[nt++] = REAL;
strcpy(tag[nt], "ErrTolTheta");
addr[nt] = &All.ErrTolTheta;
id[nt++] = REAL;
strcpy(tag[nt], "CourantFac");
addr[nt] = &All.CourantFac;
id[nt++] = REAL;
strcpy(tag[nt], "ErrTolForceAcc");
addr[nt] = &All.ErrTolForceAcc;
id[nt++] = REAL;
strcpy(tag[nt], "MaxRMSDisplacementFac");
addr[nt] = &All.MaxRMSDisplacementFac;
id[nt++] = REAL;
#ifdef HYDRO_SPH
strcpy(tag[nt], "ArtBulkViscConst");
addr[nt] = &All.ArtBulkViscConst;
id[nt++] = REAL;
#ifdef SPHAV_ARTIFICIAL_CONDUCTIVITY
strcpy(tag[nt], "ArtCondConstant");
addr[nt] = &All.ArtCondConstant;
id[nt++] = REAL;
#endif
#ifdef SPHAV_CD10_VISCOSITY_SWITCH
strcpy(tag[nt], "ViscosityAMin");
addr[nt] = &All.ViscosityAMin;
id[nt++] = REAL;
strcpy(tag[nt], "ViscosityAMax");
addr[nt] = &All.ViscosityAMax;
id[nt++] = REAL;
#endif
#ifdef SPH_TP12_ARTIFICIAL_RESISTIVITY
strcpy(tag[nt], "ArtificialResistivityMax");
addr[nt] = &All.ArtMagDispConst;
id[nt++] = REAL;
#endif
#endif
#ifdef DIVBCLEANING_DEDNER
strcpy(tag[nt], "DivBcleaningParabolicSigma");
addr[nt] = &All.DivBcleanParabolicSigma;
id[nt++] = REAL;
strcpy(tag[nt], "DivBcleaningHyperbolicSigma");
addr[nt] = &All.DivBcleanHyperbolicSigma;
id[nt++] = REAL;
#endif
#endif // closes DEVELOPER_MODE check
#ifdef GRAIN_FLUID
#ifdef GRAIN_RDI_TESTPROBLEM