<div dir="ltr">Hi Po-Yu, <div>Many thanks!</div><div>It's very helpful.</div><div><br></div><div>Best</div><div>Sheng</div></div><div class="gmail_extra"><br><div class="gmail_quote">2015-02-06 10:27 GMT-06:00 pchang <span dir="ltr"><<a href="mailto:pchang@lle.rochester.edu" target="_blank">pchang@lle.rochester.edu</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
<div>Hi Sheng,<br>
For your first question, since I never got the error "Assertion
`dataset >= 0' failed", I can't say too much about it.<br>
<br>
About the group frequency, you may try the example "LaserSlab"
(section 25.7.5 in the flash user's guide) or look at the
flash.par in that example. Following are the section in my
flash.par . The group bounds (rt_mgdBounds) needs to match what is
given in the table. You can find the format of the table in
section 21.4.3 in the flash user's guide.<br>
<br>
Cheers,<br>
<br>
Po-Yu<br>
<br>
########################################<br>
# #<br>
# RADIATION/OPACITY PARAMETERS #<br>
# #<br>
########################################<br>
rt_useMGD = .true.<br>
rt_mgdNumGroups = 12<br>
rt_mgdBounds_1 = 1.0e-02<br>
rt_mgdBounds_2 = 3.831e-02<br>
rt_mgdBounds_3 = 1.468e-01<br>
rt_mgdBounds_4 = 5.623e-01<br>
rt_mgdBounds_5 = 2.154e+00<br>
rt_mgdBounds_6 = 8.254e+00<br>
rt_mgdBounds_7 = 3.162e+01<br>
rt_mgdBounds_8 = 1.212e+02<br>
rt_mgdBounds_9 = 4.642e+02<br>
rt_mgdBounds_10 = 1.778e+03<br>
rt_mgdBounds_11 = 6.813e+03<br>
rt_mgdBounds_12 = 2.610e+04<br>
rt_mgdBounds_13 = 1.0e+05<br>
rt_mgdFlMode = "fl_harmonic"<br>
rt_mgdFlCoef = 1.0<br>
<br>
rt_mgdXlBoundaryType = "reflecting"<br>
rt_mgdXrBoundaryType = "vacuum"<br>
rt_mgdYlBoundaryType = "vacuum"<br>
rt_mgdYrBoundaryType = "vacuum"<br>
rt_mgdZlBoundaryType = "reflecting"<br>
rt_mgdZrBoundaryType = "reflecting"<br>
<br>
useOpacity = .true.<br>
<br>
### SET Gas (DD) OPACITY OPTIONS ###<br>
op_gasAbsorb = "op_tabpa"<br>
op_gasEmiss = "op_tabpe"<br>
op_gasTrans = "op_tabro"<br>
op_gasFileType = "ionmix4"<br>
op_gasFileName =
"DD_0p01Ar_T0p01-1e5_Rho1e16-1e26_vu12_0p01-1e5_Zmin0p01.cn4"<br>
<br>
### SET CHAMBER (DD) OPACITY OPTIONS ###<br>
op_chamAbsorb = "op_tabpa"<br>
op_chamEmiss = "op_tabpe"<br>
op_chamTrans = "op_tabro"<br>
op_chamFileType = "ionmix4"<br>
op_chamFileName = "DD_T0p01-1e5_Rho1e16-1e26_vu12_0p01-1e5.cn4"<br>
<br>
### SET TARGET (CH) OPACITY OPTIONS ###<br>
op_targAbsorb = "op_tabpa"<br>
op_targEmiss = "op_tabpe"<br>
op_targTrans = "op_tabro"<br>
op_targFileType = "ionmix4"<br>
op_targFileName = "CH_T0p01-1e5_Rho1e16-1e26_vu12_0p01-1e5.cn4"<span class=""><br>
<br>
<br>
<br>
<pre cols="72">--
Po-Yu Chang
Postdoctoral Associate
Department of Mechanical Engineering,
Fusion Science Center for Extreme States of Matter,
Laboratory for Laser Energetics
University of Rochester
250 East River Road
Rochester, New York 14623
Phone: (585) 273-5179
FAX: (585) 275-5960
<a href="mailto:pchang@lle.rochester.edu" target="_blank">pchang@lle.rochester.edu</a> </pre>
<br>
<br></span><div><div class="h5">
On 2/6/2015 11:04 AM, Di Sheng wrote:<br>
</div></div></div><div><div class="h5">
<blockquote type="cite">
<div dir="ltr">Hi Po-Yu,
<div>Thanks a lot for your help!</div>
<div>As for the first question, after adding "+hdf5typeio"
(i.e., ./setup -auto Cellular +hdf5typeio parallelIO=False), I
still got the same errors (Assertion `dataset >= 0'
failed). In fact, when I use "./setup -auto Cellular +ug", the
errors will disappear and the execution goes well. There must
be some reason for it. I mean, should I adopt some
"-with-unit=...." for Cellular or some more settings for setup
? </div>
<div><br>
</div>
<div>As for the second question, I didn't set the group in the
flash.par. Could you tell me how to specify the group
frequency? </div>
<div>Did I miss any other parameters? </div>
<div><br>
</div>
<div>ps. My objective is just to run this application normally,
so as to evaluate my own MPI middleware. That's it. I'll be
happy as long as it can be run, no matter what kind of
settings/parameters are used. </div>
<div><br>
</div>
<div>Many thanks again!</div>
<div><br>
</div>
<div>Best</div>
<div>Sheng</div>
<div><a href="http://www.mcs.anl.gov/%7Eshdi" target="_blank">http://www.mcs.anl.gov/~shdi</a></div>
<div class="gmail_extra"><br>
<div class="gmail_quote">2015-02-06 9:45 GMT-06:00 pchang <span dir="ltr"><<a href="mailto:pchang@lle.rochester.edu" target="_blank">pchang@lle.rochester.edu</a>></span>:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Hi
Sheng,<br>
For your first question, you may try adding "+hdf5typeio"
in your setup.<br>
For your second question, did you specify the
radiation/opacity group frequency in the flash.par and
xe-100grp-lte.imx? Are the specified the same?<br>
<br>
Po-Yu<br>
<br>
<br>
-- <br>
Po-Yu Chang<br>
<br>
Postdoctoral Associate<br>
Department of Mechanical Engineering,<br>
Fusion Science Center for Extreme States of Matter,<br>
Laboratory for Laser Energetics<br>
University of Rochester<br>
250 East River Road<br>
Rochester, New York 14623<br>
Phone: (585) 273-5179<br>
FAX: (585) 275-5960<br>
<a href="mailto:pchang@lle.rochester.edu" target="_blank">pchang@lle.rochester.edu</a>
<div>
<div><br>
<br>
<br>
<br>
On 2/5/2015 6:08 PM, Di Sheng wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Hi
everyone,<br>
I am a newbie of FLASH. I ran into the following
problems which already bothered me for a long while.
Could any one answer any one of them? many thanks!<br>
<br>
1. How to deal with "Assertion `dataset >= 0'
failed"?<br>
<br>
"./setup -auto Cellular" success, and "make" is also
successful. However, when I run "mpirun -np 1
flash4", it shows the following errors:<br>
<br>
HDF5-DIAG: Error detected in HDF5 (1.8.14) thread 0:<br>
#000: H5Ddeprec.c line 178 in H5Dcreate1(): not a
location ID<br>
major: Invalid arguments to routine<br>
minor: Inappropriate type<br>
#001: H5Gloc.c line 253 in H5G_loc(): invalid
object ID<br>
major: Invalid arguments to routine<br>
minor: Bad value<br>
flash4: io_h5create_dataset.c:57:
io_h5create_dataset: Assertion `dataset >= 0'
failed.<br>
<br>
Program received signal SIGABRT: Process abort
signal.<br>
.......<br>
<br>
The same situation goes to the "Jeans" application
(SimulationMain/Jeans).<br>
<br>
2. As for MGDInfinite, I setup it as follows:<br>
./setup -auto MGDInfinite -1d +hdf5typeio +mtmmmt
+mgd +uhd3t -with-unit=physics/materialProperties/Opacity/OpacityMain/Constant
species=xe mgd_meshgroups=4 parallelIO=FALSE<br>
<br>
Then, I add the following lines to the flash.par:<br>
sim_xeMassFrac = 1<br>
<br>
eos_xeEosType = "eos_tab"<br>
eos_xeSubType = "ionmix4"<br>
eos_xeTableFile = "xe-100grp-lte.imx"<br>
However, when I run "mpirun -np 1 flash4", I got
errors as follows:<br>
<br>
[eos_tabBrowseIonmix4Tables] IONMIX4 file found:
xe-100grp-lte.imx<br>
in eos_inittabulated, tableName = xe-100grp-lte.imx<br>
in eos_inittabulated, groupName = -none-<br>
At line 158 of file eos_tabReadIonmix4Tables.F90
(unit = 7, file = 'xe-100grp-lte.imx')<br>
Fortran runtime error: Bad value during integer read<br>
<br>
------------<br>
Best<br>
Sheng<br>
<br>
</blockquote>
<br>
<br>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>
</blockquote>
<br>
<br>
</div></div></div>
</blockquote></div><br></div>