<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<div class="moz-cite-prefix">Hi,<br>
I believe so. But I am working on laser and cylindrical target. So
I specified a cylindrical target in cylindrical coordinate. There
shouldn't be a problem to specify a spherical target. The problem
is, the ray tracing is only good for Cartesian (2D and maybe 3D)
and 2D cylindrical (r-z) coordinate if I remember correctly.
Therefore, you may not be able to do the full 3D simulation. I may
try using cylindrical coordinate so that you can have at least 2D
spherical simulation with azimuthal symmetric. <br>
I know a guy who is also doing similar problem. He may have
experience on it.<br>
<br>
Po-Yu<br>
<br>
<br>
<pre class="moz-signature" 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 class="moz-txt-link-abbreviated" href="mailto:pchang@lle.rochester.edu">pchang@lle.rochester.edu</a> </pre>
<br>
<br>
On 2/6/2015 12:53 PM, Majid Masnavi wrote:<br>
</div>
<blockquote
cite="mid:D1CBEF83-D79F-469E-AD9D-6F7CE2576FB5@creol.ucf.edu"
type="cite">
<meta http-equiv="Content-Type" content="text/html;
charset=windows-1252">
<div>Hello,</div>
<div><br>
</div>
<div>I would like to use the code to solve interaction of laser
and spherical target. Is there any way to define spherical
target?</div>
<div><br>
</div>
<div>Cheers</div>
<div>M<br>
<br>
Sent from my iPhone</div>
<div><br>
On Feb 6, 2015, at 11:38, Di Sheng <<a moz-do-not-send="true"
href="mailto:disheng222@gmail.com">disheng222@gmail.com</a>>
wrote:<br>
<br>
</div>
<blockquote type="cite">
<div>
<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 moz-do-not-send="true"
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 moz-do-not-send="true" 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 moz-do-not-send="true"
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
moz-do-not-send="true"
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 moz-do-not-send="true"
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>
</div>
</blockquote>
</blockquote>
<br>
</body>
</html>