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 <<a href="mailto:jfg@ucolick.org">jfg@ucolick.org</a>> 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've attached two flash.par files, the first is from the initial run, the second is for the restarted run. I'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't want to attach all of them)<div class="Ih2E3d"><br>
<br>
-- <br>
James Guillochon<br>
Department of Astronomy & 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'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 <<a href="mailto:jfg@ucolick.org" target="_blank">jfg@ucolick.org</a> <mailto:<a href="mailto:jfg@ucolick.org" target="_blank">jfg@ucolick.org</a>>> wrote:<br>
<br>
Hey guys,<br>
<br>
It seems that the velz variable is getting shifted relative to the<br>
other variables when I restart a simulation. This shift appears to<br>
be equal to the nzb variable, or 8 blocks, at the highest refinement<br>
level. I have attached two screenshots: The first is prior to the<br>
restart, the second is immediately after the restart.<br>
<br>
The only variable that changes between restarts is the max<br>
refinement level; however the shift in velz is visible before the<br>
grid is refined. No other variables are shifted. This is a real<br>
problem for me since my problem assumes symmetry about the z axis!<br>
<br>
Any ideas guys?<br>
<br>
-- James Guillochon<br>
Department of Astronomy & Astrophysics<br>
University of California, Santa Cruz<br></div>
<a href="mailto:jfg@ucolick.org" target="_blank">jfg@ucolick.org</a> <mailto:<a href="mailto:jfg@ucolick.org" target="_blank">jfg@ucolick.org</a>><br>
<br>
<br>
!DSPAM:10135,4829ef7c2382245126300!<br>
</blockquote>
<br>lrefine_min = 1<br>
lrefine_max = 10<br>
<br>
<br>
basenm = "tidal_"<br>
restart = .true.<br>
checkpointFileIntervalTime = 200.0<br>
checkpointFileNumber = 9<br>
plotFileIntervalTime = 0.5<br>
<br>
xmax = 6.e11<br>
ymax = 6.e11<br>
zmax = 6.e11<br>
<br>
ptmass = 2.e39<br>
ptxpos = 3.e11<br>
ptypos = 3.e11<br>
ptzpos = 3.e11<br>
<br>
small = 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 = 1000000<br>
tmax = 2200.<br>
<br>
gamma = 1.6667<br>
<br>
eintSwitch = 1.e-4<br>
#if beta is < 0, the code uses periDist instead.<br>
sim_periBeta = 1.<br>
sim_periDist = 1.e12<br>
sim_periTime = 2000.<br>
sim_rhoAmbient = 1.0e-5<br>
sim_pAmbient = 1.0e+5<br>
sim_xctr = 3.e11<br>
sim_yctr = 3.e11<br>
sim_zctr = 3.e11<br>
sim_nsubzones = 10<br>
<br>
<br>
xl_boundary_type = "outflow"<br>
xr_boundary_type = "outflow"<br>
yl_boundary_type = "outflow"<br>
yr_boundary_type = "outflow"<br>
zl_boundary_type = "outflow"<br>
zr_boundary_type = "outflow"<br>
grav_boundary_type = "isolated"<br>
mpole_lmax = 8<br>
#useGravity = .true.<br>
<br>
plot_var_1 = "dens"<br>
plot_var_2 = "temp"<br>
plot_var_3 = "velx"<br>
plot_var_4 = "vely"<br>
plot_var_5 = "velz"<br>
plot_var_6 = "gpot"<br>
<br>
refine_var_1 = "dens"<br>
#refine_var_2 = "pres"<br>
refine_filter_1 = 5.e-2<br>
#refine_filter_2 = 5.e-2<br>
refine_val_cutoff_1 = 1.e0<br>
<br>lrefine_min = 1<br>
lrefine_max = 8<br>
<br>
<br>
basenm = "tidal_"<br>
restart = .false.<br>
checkpointFileIntervalTime = 200.0<br>
plotFileIntervalTime = 10.0<br>
<br>
xmax = 6.e11<br>
ymax = 6.e11<br>
zmax = 6.e11<br>
<br>
ptmass = 2.e39<br>
ptxpos = 3.e11<br>
ptypos = 3.e11<br>
ptzpos = 3.e11<br>
<br>
small = 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 = 1000000<br>
tmax = 4000.<br>
<br>
gamma = 1.6667<br>
<br>
eintSwitch = 1.e-4<br>
#if beta is < 0, the code uses periDist instead.<br>
sim_periBeta = 1.<br>
sim_periDist = 1.e12<br>
sim_periTime = 2000.<br>
sim_rhoAmbient = 1.0e-5<br>
sim_pAmbient = 1.0e+5<br>
sim_xctr = 3.e11<br>
sim_yctr = 3.e11<br>
sim_zctr = 3.e11<br>
sim_nsubzones = 10<br>
<br>
<br>
xl_boundary_type = "outflow"<br>
xr_boundary_type = "outflow"<br>
yl_boundary_type = "outflow"<br>
yr_boundary_type = "outflow"<br>
zl_boundary_type = "outflow"<br>
zr_boundary_type = "outflow"<br>
grav_boundary_type = "isolated"<br>
mpole_lmax = 8<br>
#useGravity = .true.<br>
<br>
plot_var_1 = "dens"<br>
plot_var_2 = "temp"<br>
plot_var_3 = "velx"<br>
plot_var_4 = "vely"<br>
plot_var_5 = "velz"<br>
plot_var_6 = "gpot"<br>
<br>
refine_var_1 = "dens"<br>
#refine_var_2 = "pres"<br>
refine_filter_1 = 5.e-2<br>
#refine_filter_2 = 5.e-2<br>
refine_val_cutoff_1 = 1.e0<br>
<br></blockquote></div><br>