James,<br><br>We will definitely need all the Simulation_ and Grid_ files to have any chance of finding out<br>the problem. However, without having all the files it is not clear that we can solve your problem, because<br>we run many problems with restart in our nightly test suite, and all of them start transparently. It seems <br>
very likely that your restart issue is specific to your problem.<br><br>If you are really uncomfortable about sharing the files, then you can try the following steps, at each step <br>making sure that you have transparent restart :<br>
1. Start with a simple problem like Sedov.<br>2. Add self gravity to it.<br>3. Add your custom gravity routines.<br>4. Start your own problem with gravity turned off, and default mark_refine_derefine<br>5. Add each custom component one by one.<br>
<br>That is basically the strategy I would use for identifying the problem.<br><br>Let me also assure you that if you were to give the files to us, only one person will have access to them,<br>and he/she will delete them as soon as they are done with them. <br>
<br>Anshu<br><br><div class="gmail_quote">On Tue, May 13, 2008 at 4:01 PM, James Guillochon &lt;<a href="mailto:jfg@ucolick.org">jfg@ucolick.org</a>&gt; wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
I just checked and the problem occurs even if the max refinement level is the same before and after the restart.<br>
<br>
I&#39;ve attached two flash.par files, the first is from the initial run, the second is for the restarted run. I&#39;ve modified the following built-in files:<br>
<br>
Gravity_accelAtCoords.F90<br>
Gravity_accelListOfBlocks.F90<br>
Gravity_accelOneRow.F90<br>
Gravity_data.F90<br>
Gravity_init.F90<br>
Gravity_potentialListOfBlocks.F90<br>
Grid_data.F90<br>
Grid_init.F90<br>
Grid_markRefineDerefine.F90<br>
Simulation_data.F90<br>
Simulation_initBlock.F90<br>
Simulation_init.F90<br>
<br>
Which in particular would you like me to attach? (I don&#39;t want to attach all of them)<div class="Ih2E3d"><br>
<br>
-- <br>
James Guillochon<br>
Department of Astronomy &amp; Astrophysics<br>
University of California, Santa Cruz<br>
<a href="mailto:jfg@ucolick.org" target="_blank">jfg@ucolick.org</a><br>
<br>
<br></div>
Anshu Dubey wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="Ih2E3d">
Hi James,<br>
<br>
Can you share your setup, flash.par and any other customized file with us ? We have never encountered the problem<br>
that you are seeing, and we&#39;d really like to be able to reproduce it ourselves to figure out what is going on.<br>
<br>
Anshu<br>
<br></div><div class="Ih2E3d">
On Tue, May 13, 2008 at 2:35 PM, James Guillochon &lt;<a href="mailto:jfg@ucolick.org" target="_blank">jfg@ucolick.org</a> &lt;mailto:<a href="mailto:jfg@ucolick.org" target="_blank">jfg@ucolick.org</a>&gt;&gt; wrote:<br>

<br>
 &nbsp; &nbsp;Hey guys,<br>
<br>
 &nbsp; &nbsp;It seems that the velz variable is getting shifted relative to the<br>
 &nbsp; &nbsp;other variables when I restart a simulation. This shift appears to<br>
 &nbsp; &nbsp;be equal to the nzb variable, or 8 blocks, at the highest refinement<br>
 &nbsp; &nbsp;level. I have attached two screenshots: The first is prior to the<br>
 &nbsp; &nbsp;restart, the second is immediately after the restart.<br>
<br>
 &nbsp; &nbsp;The only variable that changes between restarts is the max<br>
 &nbsp; &nbsp;refinement level; however the shift in velz is visible before the<br>
 &nbsp; &nbsp;grid is refined. No other variables are shifted. This is a real<br>
 &nbsp; &nbsp;problem for me since my problem assumes symmetry about the z axis!<br>
<br>
 &nbsp; &nbsp;Any ideas guys?<br>
<br>
 &nbsp; &nbsp;--  &nbsp; &nbsp;James Guillochon<br>
 &nbsp; &nbsp;Department of Astronomy &amp; Astrophysics<br>
 &nbsp; &nbsp;University of California, Santa Cruz<br></div>
 &nbsp; &nbsp;<a href="mailto:jfg@ucolick.org" target="_blank">jfg@ucolick.org</a> &lt;mailto:<a href="mailto:jfg@ucolick.org" target="_blank">jfg@ucolick.org</a>&gt;<br>
<br>
<br>
!DSPAM:10135,4829ef7c2382245126300!<br>
</blockquote>
<br>lrefine_min &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 1<br>
lrefine_max &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 10<br>
<br>
<br>
basenm &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;tidal_&quot;<br>
restart &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = .true.<br>
checkpointFileIntervalTime = 200.0<br>
checkpointFileNumber = 9<br>
plotFileIntervalTime &nbsp; &nbsp;= 0.5<br>
<br>
xmax &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 6.e11<br>
ymax &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 6.e11<br>
zmax &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 6.e11<br>
<br>
ptmass &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 2.e39<br>
ptxpos &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 3.e11<br>
ptypos &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 3.e11<br>
ptzpos &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 3.e11<br>
<br>
small = &nbsp;1.e-30<br>
smalle = 1.e-30<br>
smallt = 1.e-30<br>
smallu = 1.e-30<br>
smallx = 1.e-30<br>
smally = 1.e-30<br>
smlrho = 1.e-30<br>
smallp = 1.e-30<br>
#<br>
#cfl = 0.8<br>
#cvisc = 0.e0<br>
#iplm = 1<br>
nriem = 1000<br>
#eos_maxNewton = 1000<br>
#convertToConsvdForMeshCalls = .true.<br>
dtinit = 1e-2<br>
#dtmin = 1e-1<br>
#dtmax = 1e-2<br>
#memory_stat_freq = 1<br>
<br>
nend &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1000000<br>
tmax &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 2200.<br>
<br>
gamma &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 1.6667<br>
<br>
eintSwitch &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1.e-4<br>
#if beta is &lt; 0, the code uses periDist instead.<br>
sim_periBeta &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1.<br>
sim_periDist &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1.e12<br>
sim_periTime &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 2000.<br>
sim_rhoAmbient &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1.0e-5<br>
sim_pAmbient &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1.0e+5<br>
sim_xctr &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 3.e11<br>
sim_yctr &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 3.e11<br>
sim_zctr &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 3.e11<br>
sim_nsubzones &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 10<br>
<br>
<br>
xl_boundary_type &nbsp; &nbsp; &nbsp; &nbsp;= &quot;outflow&quot;<br>
xr_boundary_type &nbsp; &nbsp; &nbsp; &nbsp;= &quot;outflow&quot;<br>
yl_boundary_type &nbsp; &nbsp; &nbsp; &nbsp;= &quot;outflow&quot;<br>
yr_boundary_type &nbsp; &nbsp; &nbsp; &nbsp;= &quot;outflow&quot;<br>
zl_boundary_type &nbsp; &nbsp; &nbsp; &nbsp;= &quot;outflow&quot;<br>
zr_boundary_type &nbsp; &nbsp; &nbsp; &nbsp;= &quot;outflow&quot;<br>
grav_boundary_type &nbsp; &nbsp; &nbsp;= &quot;isolated&quot;<br>
mpole_lmax &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 8<br>
#useGravity &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= .true.<br>
<br>
plot_var_1 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;dens&quot;<br>
plot_var_2 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;temp&quot;<br>
plot_var_3 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;velx&quot;<br>
plot_var_4 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;vely&quot;<br>
plot_var_5 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;velz&quot;<br>
plot_var_6 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;gpot&quot;<br>
<br>
refine_var_1 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;dens&quot;<br>
#refine_var_2 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;pres&quot;<br>
refine_filter_1 &nbsp; &nbsp; &nbsp; &nbsp; = 5.e-2<br>
#refine_filter_2 &nbsp; &nbsp; &nbsp; &nbsp; = 5.e-2<br>
refine_val_cutoff_1 &nbsp; &nbsp; = 1.e0<br>
<br>lrefine_min &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 1<br>
lrefine_max &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 8<br>
<br>
<br>
basenm &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;tidal_&quot;<br>
restart &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = .false.<br>
checkpointFileIntervalTime = 200.0<br>
plotFileIntervalTime &nbsp; &nbsp;= 10.0<br>
<br>
xmax &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 6.e11<br>
ymax &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 6.e11<br>
zmax &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 6.e11<br>
<br>
ptmass &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 2.e39<br>
ptxpos &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 3.e11<br>
ptypos &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 3.e11<br>
ptzpos &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 3.e11<br>
<br>
small = &nbsp;1.e-30<br>
smalle = 1.e-30<br>
smallt = 1.e-30<br>
smallu = 1.e-30<br>
smallx = 1.e-30<br>
smally = 1.e-30<br>
smlrho = 1.e-30<br>
smallp = 1.e-30<br>
#<br>
#cfl = 0.8<br>
#cvisc = 0.e0<br>
#iplm = 1<br>
nriem = 1000<br>
#eos_maxNewton = 1000<br>
#convertToConsvdForMeshCalls = .true.<br>
dtinit = 1e-2<br>
#dtmin = 1e-1<br>
#dtmax = 1e-2<br>
#memory_stat_freq = 1<br>
<br>
nend &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1000000<br>
tmax &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 4000.<br>
<br>
gamma &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 1.6667<br>
<br>
eintSwitch &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1.e-4<br>
#if beta is &lt; 0, the code uses periDist instead.<br>
sim_periBeta &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1.<br>
sim_periDist &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1.e12<br>
sim_periTime &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 2000.<br>
sim_rhoAmbient &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1.0e-5<br>
sim_pAmbient &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1.0e+5<br>
sim_xctr &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 3.e11<br>
sim_yctr &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 3.e11<br>
sim_zctr &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 3.e11<br>
sim_nsubzones &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 10<br>
<br>
<br>
xl_boundary_type &nbsp; &nbsp; &nbsp; &nbsp;= &quot;outflow&quot;<br>
xr_boundary_type &nbsp; &nbsp; &nbsp; &nbsp;= &quot;outflow&quot;<br>
yl_boundary_type &nbsp; &nbsp; &nbsp; &nbsp;= &quot;outflow&quot;<br>
yr_boundary_type &nbsp; &nbsp; &nbsp; &nbsp;= &quot;outflow&quot;<br>
zl_boundary_type &nbsp; &nbsp; &nbsp; &nbsp;= &quot;outflow&quot;<br>
zr_boundary_type &nbsp; &nbsp; &nbsp; &nbsp;= &quot;outflow&quot;<br>
grav_boundary_type &nbsp; &nbsp; &nbsp;= &quot;isolated&quot;<br>
mpole_lmax &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 8<br>
#useGravity &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= .true.<br>
<br>
plot_var_1 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;dens&quot;<br>
plot_var_2 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;temp&quot;<br>
plot_var_3 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;velx&quot;<br>
plot_var_4 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;vely&quot;<br>
plot_var_5 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;velz&quot;<br>
plot_var_6 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;gpot&quot;<br>
<br>
refine_var_1 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;dens&quot;<br>
#refine_var_2 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;pres&quot;<br>
refine_filter_1 &nbsp; &nbsp; &nbsp; &nbsp; = 5.e-2<br>
#refine_filter_2 &nbsp; &nbsp; &nbsp; &nbsp; = 5.e-2<br>
refine_val_cutoff_1 &nbsp; &nbsp; = 1.e0<br>
<br></blockquote></div><br>