the check for it in connection with generating the first output. This was a bug introduced by checking for the restart variable when calling output.
doing a restart. This behaviour let to overwriting of the original results, not necessarily with the correct values (e.g. IonRates files).
No change, just a local conflict on laptop. Conflicts: mat_ini_cubep3m.F90
These changes may be out of date compared to changes made on sfinx.
in a bug in make_LLS so I am not sure the code has changed fundamentally but this version has been tested and validated.
is not reached. Useful for the reporting to the log file, making it consistent with the behaviour for homogeneous LLS case.
case of LLS. If the mean free path is less than a certain value, LLS opacity is not applied.
position dependent LLSs. Previous versions of position dependent LLS were never tested or used and were wrong.
Do not use for other cases!!!!
…clock. Possible improvement to consider: make timestep_wallclock function return time in specified units (seconds, minutes, hours, etc..)
cause no change in the output being produced by the code
and mass weighted mean ionization fractions for writing to PhotonCounts2.dat (unit 95) were wrongly defined as arrays running from 0 to 2. This mistake was due to synchronization with the Helium version of the code (where both the HeII and HeIII fractions are calculated and stored in these arrays). volfrac and massfrac are now scalars again.
photonstatistics and output modules to a separate directory, called compressed. The modules use a compressed version of the ionized/neutral fractions but have not been kept up to date. They are to retained for future use but are for now moved out of the way to create a cleaner source directory.
was missing, leading to all nodes writing out the data!!!!
Added more comments to the write_photonstatistics subroutine Increased the length of various filename variables to 512.
it just checks that the intermediate redshift is indeed intermediate between the current and the next redshift (intermediate redshift when restarting).
to accommodate full path names in file_admin.f90. Likewise in file_admin.f90 to 400 characters
the mpi_bcast of the phih_grid array was erroneously removed. This has now been restored. Additionally the two quantities used in checking convergence of the global ionization fractions (prev_sum_xh1_int, prev_sum_xh0_int) are being mpi_bcast.
…on ALLFRAC case. It should not be 1.0 - sum_xh1_int, but real(mesh(1)*mesh(2)*mesh(3)) - sum_xh1_int. The global average fractions are not normalized during this calculation hence this factor instead of 1.0 This is now fixed.
This module now uses the new temperature structures.
…ates. The code was still using the temper variable which has become undefined. This can lead to divbyzero errors. SHould use get_temperature_point. This has now been added to this subroutine.