[FLASH-USERS] small numbers in guard cellsod
张高原
zgy0106 at yeah.net
Mon Sep 22 13:51:44 EDT 2014
Hi,
I am trying to setup a simulation using FLASH code. It runs well in a test problem, Sod Shock-Tube, using the original initial parameters. But when I run it using physical parameters in units of "cgs", it raises a warning:
"after gc filling: min. unk(DENS_VAR)=0.1019999999999999981868E-24 PE=0 block=21 type=1 ",
without any output file except the "0000" one. If I change the "plotfileIntervalStep" to 1, it produces some 'plt' files, but the result is obviously wrong. Following is what the screen output looks like. The "flash.par" file and a result density image are in the attachments.
16 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
15 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
14 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
13 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
12 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
11 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
10 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
9 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
8 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
7 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
6 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
5 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
4 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
3 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
2 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
1 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
WARNING after gc filling: min. unk(DENS_VAR)=0.1019999999999999981868E-24 PE=0 block=21 type=1
...
1 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
WARNING after gc filling: min. unk(DENS_VAR)=0.9999999999999999071852E-10 PE=0 block=37 type=1
16 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
15 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
14 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
13 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
12 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
11 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
10 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
9 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
8 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
7 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
6 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
5 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
4 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 ^C0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
3 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
2 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
1 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
WARNING after gc filling: min. unk(DENS_VAR)=0.9999999999999999071852E-10 PE=0 block=41 type=1
From the warning message, I think it's because of the value of density. Is it too small? Or any thoughts on how to make it better?
Thanks in advance.
Cheers,
Gao-Yuan
--
ZHANG Gao-Yuan
Department of Astronomy, Nanjing University
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20140923/48198715/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: flash.par
Type: application/octet-stream
Size: 4657 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20140923/48198715/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: sodclump_hdf5_plt_cnt_0001_Slice_z_dens.png
Type: image/png
Size: 64190 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20140923/48198715/attachment.png>
More information about the flash-users
mailing list