gandalfcode / gandalf

GANDALF (Graphical Astrophysics code for N-body Dynamics And Lagrangian Fluids)
GNU General Public License v2.0
44 stars 12 forks source link

[run_id].diag what are the headers for this data? #185

Open speed-wagon opened 4 years ago

speed-wagon commented 4 years ago

The diagnostic file contains some good data regarding the simulation. However the file contains some 30 columns. Whilst I am able to guess some of the values like the time step, time and dt...it would be nice to now exactly what all the headers were.

I havent been able to locate any info on this in the docs and looking into the source code has only lead me to find ~20 of the 30 column headers.

Can anyone help with this?

rbooth200 commented 4 years ago

Is what you are looking for in the function: void Simulation::RecordDiagnostics(void)

in https://github.com/gandalfcode/gandalf/blob/master/src/Common/SimAnalysis.hpp

speed-wagon commented 4 years ago

That is the function were I obtained the first 20 or so headers, but my .diag file contains 30 columns.

The simlation is the accretion disk. Unsure if this adds addional simulation specific columns

rbooth200 commented 4 years ago

The lines

for (k=0; k<3; k++) outfile << diag.angmom[k]*simunits.angmom.outscale << "     ";
  for (k=0; k<ndim; k++) outfile << diag.rcom[k]*simunits.r.outscale << "     ";
  for (k=0; k<ndim; k++) outfile << diag.vcom[k]*simunits.v.outscale << "     ";
  for (k=0; k<ndim; k++) outfile << diag.mom[k]*simunits.mom.outscale << "     ";
  for (k=0; k<ndim; k++) outfile << diag.force[k] << "     "; //*simunits.f.o

Output 3 items each for a 3D simulation.