Saving initial condition with multiple scalar fields.

42 views
Skip to first unread message

R Frantz

unread,
Apr 11, 2023, 4:30:16 AM4/11/23
to Nek5000
Dear Neks,

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

YuHsiang Lan

unread,
Apr 11, 2023, 5:17:27 AM4/11/23
to Nek5000
Hi,

The outpost only dump single temperature field
Use outpost2 instead so you can specified the nfldt.
  nfldt = 3  ! temperature + 2 scalars
  call outpost2(vx,vy,vz,pr,t,nfldt,'   ')

Hope this helps,
Yu-Hsiang
--

R Frantz

unread,
Apr 11, 2023, 11:16:20 AM4/11/23
to Nek5000
Yes, it solves the problem. Thank you!

tuanliu wenxin

unread,
Feb 5, 2026, 6:51:47 AM (17 hours ago) Feb 5
to Nek5000
Thank you very much! This help me a lot! 
Reply all
Reply to author
Forward
0 new messages