I am working on a case that involves temperature and extra scalar fields. These fields are present in the "case0.f00001" file automatically. However, when include the following code on usercheck:
if (istep.eq.0) call outpost(vx, vy, vz, pr, t, ' ')
to save the initial condition, the extra scalar fields are missing from the saved file, just temperature is present.
If I print the control flags, they are set correctly:
ifto= T
ifpsco= T T F F
What I am doing wrong?
Thanks,
R