<html><head><meta http-equiv="Content-Type" content="text/html charset=windows-1252"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">One quick suggestion I have is to add +pm4dev to you setup line. That brings in the most current grid package.<div><br></div><div>Sean<br><div apple-content-edited="true">
<div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><span class="Apple-style-span" style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; border-collapse: separate; border-spacing: 0px; -webkit-text-decorations-in-effect: none;"><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; border-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-stroke-width: 0px;"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div><br class="Apple-interchange-newline">--------------------------------------------------------</div><div>Sean M. Couch</div><div>Hubble Fellow</div><div>Flash Center for Computational Science</div><div>Department of Astronomy & Astrophysics</div><div>The University of Chicago</div><div>5747 S Ellis Ave, Jo 315</div><div>Chicago, IL 60637</div><div>(773) 702-3899</div><div><a href="http://www.flash.uchicago.edu/~smc">www.flash.uchicago.edu/~smc</a></div><div><br></div></div></span></span></div><br class="Apple-interchange-newline"><br class="Apple-interchange-newline">
</div>
<br><div><div>On Jul 7, 2014, at 9:30 AM, Luke Zoltan Kelley <<a href="mailto:lkelley@cfa.harvard.edu">lkelley@cfa.harvard.edu</a>> wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><meta http-equiv="Content-Type" content="text/html charset=windows-1252"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div><div>Sure, I wasn't sure If flash-users allowed attachments so I've added the log contents:</div><div><br></div><div><br></div><div> FLASH log file: 07-07-2014 10:26:23.961 Run number: 1</div><div> ==============================================================================</div><div> Number of MPI tasks: 1</div><div> MPI version: 3</div><div> MPI subversion: 0</div><div> Dimensionality: 3</div><div> Max Number of Blocks/Proc: 200</div><div> Number x zones: 8</div><div> Number y zones: 8</div><div> Number z zones: 8</div><div> Setup stamp: Mon Jul 7 10:25:52 2014</div><div> Build stamp: Mon Jul 7 10:25:56 2014</div><div> System info: Darwin daedalus.local 13.3.0 Darwin Kernel Version 13.3.0: Tue Jun 3 21:27:35 P</div><div> Version: FLASH 4.2.2_release </div><div> Build directory: /Users/lzkelley/Applications/flash/flash4.2.2/object_starwind</div><div> Setup syntax: /Users/lzkelley/Applications/flash/flash4.2.2/bin/setup.py starwind -auto -objdir=object_starwind -3d +spherical -noclobber -debug</div><div> f compiler flags: /opt/local//bin/mpif90 -ggdb -c -O0 -fdefault-real-8 -fdefault-double-8 -pedantic -Wall -Wextra -Waliasing -Wsurprising -Wconversion -Wunderflow -ffpe-trap=invalid,zero,overflow -fbounds-check -fimplicit-none -fstack-protector-all -ffree-line-length-0 -Wno-unused -Wno-unused-parameter -DMAXBLOCKS=200 -DNXB=8 -DNYB=8 -DNZB=8 -DN_DIM=3</div><div> c compiler flags: /opt/local//bin/mpicc -I /opt/local//include -DH5_USE_16_API -ggdb -c -O0 -Wno-div-by-zero -Wundef -Wconversion -Wstrict-prototypes -Wunreachable-code -pedantic -Wall -Wextra -Winit-self -Wfloat-equal -Wunsafe-loop-optimizations -Wpadded -fstack-protector-all -DMAXBLOCKS=200 -DNXB=8 -DNYB=8 -DNZB=8 -DN_DIM=3 -DDarwin</div><div> ==============================================================================</div><div> Comment: FLASH 3 run</div><div> ==============================================================================</div><div> FLASH Units used:</div><div> Driver/DriverMain/Split</div><div> Grid/GridBoundaryConditions/OneRow</div><div> Grid/GridMain/paramesh/interpolation/Paramesh4/prolong</div><div> Grid/GridMain/paramesh/interpolation/prolong</div><div> Grid/GridMain/paramesh/paramesh4/Paramesh4dev/PM4_package/headers</div><div> Grid/GridMain/paramesh/paramesh4/Paramesh4dev/PM4_package/mpi_source</div><div> Grid/GridMain/paramesh/paramesh4/Paramesh4dev/PM4_package/source</div><div> Grid/GridMain/paramesh/paramesh4/Paramesh4dev/PM4_package/utilities/multigrid</div><div> Grid/GridSolvers/Multipole_new</div><div> Grid/localAPI</div><div> IO/IOMain/hdf5/serial/PM</div><div> IO/localAPI</div><div> Multispecies/MultispeciesMain</div><div> PhysicalConstants/PhysicalConstantsMain</div><div> RuntimeParameters/RuntimeParametersMain</div><div> Simulation/SimulationComposition/Burn</div><div> Simulation/SimulationMain/starwind</div><div> flashUtilities/contiguousConversion</div><div> flashUtilities/general</div><div> flashUtilities/interpolation/oneDim</div><div> flashUtilities/nameValueLL</div><div> flashUtilities/sorting/quicksort</div><div> flashUtilities/system/memoryUsage/legacy</div><div> monitors/Logfile/LogfileMain</div><div> monitors/Timers/TimersMain/MPINative</div><div> physics/Eos/EosMain/Gamma</div><div> physics/Eos/localAPI</div><div> physics/Hydro/HydroMain/split/PPM/PPMKernel</div><div> ==============================================================================</div><div> RuntimeParameters:</div><div> </div><div> ==============================================================================</div><div>bndpriorityone = 1</div><div>bndprioritythree = 3</div><div>bndprioritytwo = 2</div><div>checkpointfileintervalstep = 0</div><div>checkpointfilenumber = 0</div><div>dr_abortpause = 2</div><div>drift_break_inst = 0</div><div>drift_trunc_mantissa = 2</div><div>drift_verbose_inst = 0</div><div>eos_entrelescalechoice = 6</div><div>fileformatversion = 9</div><div>forcedplotfilenumber = 0</div><div>gr_lrefmaxtimevalue_1 = -1</div><div>gr_lrefmaxtimevalue_10 = -1</div><div>gr_lrefmaxtimevalue_11 = -1</div><div>gr_lrefmaxtimevalue_12 = -1</div><div>gr_lrefmaxtimevalue_13 = -1</div><div>gr_lrefmaxtimevalue_14 = -1</div><div>gr_lrefmaxtimevalue_15 = -1</div><div>gr_lrefmaxtimevalue_16 = -1</div><div>gr_lrefmaxtimevalue_17 = -1</div><div>gr_lrefmaxtimevalue_18 = -1</div><div>gr_lrefmaxtimevalue_19 = -1</div><div>gr_lrefmaxtimevalue_2 = -1</div><div>gr_lrefmaxtimevalue_20 = -1</div><div>gr_lrefmaxtimevalue_3 = -1</div><div>gr_lrefmaxtimevalue_4 = -1</div><div>gr_lrefmaxtimevalue_5 = -1</div><div>gr_lrefmaxtimevalue_6 = -1</div><div>gr_lrefmaxtimevalue_7 = -1</div><div>gr_lrefmaxtimevalue_8 = -1</div><div>gr_lrefmaxtimevalue_9 = -1</div><div>gr_pmrpdivergencefree = 1</div><div>gr_pmrpifaceoff = 0</div><div>gr_pmrpl2p5d = 0</div><div>gr_pmrpmaxblocks = -1</div><div>gr_pmrpmflags = 1</div><div>gr_pmrpnboundaries = 6</div><div>gr_pmrpndim = 3</div><div>gr_pmrpnedgevar1 = -1</div><div>gr_pmrpnfacevar = -1</div><div>gr_pmrpnfielddivf = -1</div><div>gr_pmrpnfluxvar = -1</div><div>gr_pmrpnguard = -1</div><div>gr_pmrpnguardwork = -1</div><div>gr_pmrpnvar = -1</div><div>gr_pmrpnvarwork = 1</div><div>gr_pmrpnvarcorn = 0</div><div>gr_pmrpnvaredge = 0</div><div>gr_pmrpnxb = -1</div><div>gr_pmrpnyb = -1</div><div>gr_pmrpnzb = -1</div><div>gr_restrictallmethod = 3</div><div>gr_sanitizedatamode = 1</div><div>gr_sanitizeverbosity = 5</div><div>iprocs = 1</div><div>igodu = 0</div><div>interpol_order = 2</div><div>iplm = 0</div><div>irenorm = 0</div><div>jprocs = 1</div><div>kprocs = 1</div><div>lrefine_del = 0</div><div>lrefine_max = 4 [CHANGED]</div><div>lrefine_min = 3 [CHANGED]</div><div>lrefine_min_init = 1</div><div>max_particles_per_blk = 100</div><div>memory_stat_freq = 100000</div><div>meshcopycount = 1</div><div>min_particles_per_blk = 1</div><div>mpole_innerzonesize = 16</div><div>mpole_lmax = 8 [CHANGED]</div><div>mpole_maxradialzones = 1</div><div>nbegin = 1</div><div>nblockx = 1</div><div>nblocky = 1</div><div>nblockz = 1</div><div>nend = 1000000000 [CHANGED]</div><div>nrefs = 2</div><div>nriem = 10000 [CHANGED]</div><div>nsteptotalsts = 5</div><div>outputsplitnum = 1</div><div>plotfileintervalstep = 0</div><div>plotfilenumber = 0</div><div>ppmeintcompfluxconstructionmeth = 0</div><div>ppmeintfluxconstructionmeth = 0</div><div>ppmenercompfluxconstructionmeth = 0</div><div>ppmenerfluxconstructionmeth = 0</div><div>refine_var_count = 4</div><div>rolling_checkpoint = 10000</div><div>sim_maxblocks = 40000</div><div>sim_nsubzones = 4 [CHANGED]</div><div>sweeporder = 123</div><div>wr_integrals_freq = 1</div><div>cfl = 0.800E+00</div><div>checkpointfileintervaltime = 0.100E+06 [CHANGED]</div><div>checkpointfileintervalz = 0.180+309</div><div>cvisc = 0.100E+00</div><div>derefine_cutoff_1 = 0.200E+00</div><div>derefine_cutoff_2 = 0.200E+00</div><div>derefine_cutoff_3 = 0.200E+00</div><div>derefine_cutoff_4 = 0.200E+00</div><div>dp_sh = 0.330E+00</div><div>dp_sh_md = 0.330E+00</div><div>dr_tstepslowstartfactor = 0.100E+00</div><div>dtinit = 0.100E-03 [CHANGED]</div><div>dtmax = 0.100E+09 [CHANGED]</div><div>dtmin = 0.100E-09</div><div>eintswitch = 0.100E-03 [CHANGED]</div><div>eos_singlespeciesa = 0.100E+01</div><div>eos_singlespeciesz = 0.100E+01</div><div>epsiln = 0.330E+00</div><div>gamma = 0.167E+01</div><div>gr_lrefinemaxredlogbase = 0.100E+02</div><div>gr_lrefinemaxredradiusfact = 0.000E+00</div><div>gr_lrefinemaxredtref = 0.000E+00</div><div>gr_lrefinemaxredtimescale = 0.100E+01</div><div>gr_lrefmaxtime_1 = -0.100E+01</div><div>gr_lrefmaxtime_10 = -0.100E+01</div><div>gr_lrefmaxtime_11 = -0.100E+01</div><div>gr_lrefmaxtime_12 = -0.100E+01</div><div>gr_lrefmaxtime_13 = -0.100E+01</div><div>gr_lrefmaxtime_14 = -0.100E+01</div><div>gr_lrefmaxtime_15 = -0.100E+01</div><div>gr_lrefmaxtime_16 = -0.100E+01</div><div>gr_lrefmaxtime_17 = -0.100E+01</div><div>gr_lrefmaxtime_18 = -0.100E+01</div><div>gr_lrefmaxtime_19 = -0.100E+01</div><div>gr_lrefmaxtime_2 = -0.100E+01</div><div>gr_lrefmaxtime_20 = -0.100E+01</div><div>gr_lrefmaxtime_3 = -0.100E+01</div><div>gr_lrefmaxtime_4 = -0.100E+01</div><div>gr_lrefmaxtime_5 = -0.100E+01</div><div>gr_lrefmaxtime_6 = -0.100E+01</div><div>gr_lrefmaxtime_7 = -0.100E+01</div><div>gr_lrefmaxtime_8 = -0.100E+01</div><div>gr_lrefmaxtime_9 = -0.100E+01</div><div>larget = 0.100E+11</div><div>mpole_innerzoneresolution = 0.100E+00</div><div>mpole_zoneexponent_1 = 0.100E+01</div><div>mpole_zoneexponent_2 = 0.100E+01</div><div>mpole_zoneexponent_3 = 0.100E+01</div><div>mpole_zoneexponent_4 = 0.100E+01</div><div>mpole_zoneradiusfraction_1 = 0.100E+01</div><div>mpole_zoneradiusfraction_2 = 0.100E+01</div><div>mpole_zoneradiusfraction_3 = 0.100E+01</div><div>mpole_zoneradiusfraction_4 = 0.100E+01</div><div>mpole_zonescalar_1 = 0.100E+01</div><div>mpole_zonescalar_2 = 0.100E+01</div><div>mpole_zonescalar_3 = 0.100E+01</div><div>mpole_zonescalar_4 = 0.100E+01</div><div>nusts = 0.100E+00</div><div>omg1 = 0.750E+00</div><div>omg2 = 0.100E+02</div><div>plotfileintervaltime = 0.100E-02 [CHANGED]</div><div>plotfileintervalz = 0.180+309</div><div>refine_cutoff_1 = 0.800E+00</div><div>refine_cutoff_2 = 0.800E+00</div><div>refine_cutoff_3 = 0.800E+00</div><div>refine_cutoff_4 = 0.800E+00</div><div>refine_filter_1 = 0.100E-01</div><div>refine_filter_2 = 0.100E-01</div><div>refine_filter_3 = 0.100E-01</div><div>refine_filter_4 = 0.100E-01</div><div>rieman_tol = 0.100E-13 [CHANGED]</div><div>rss_limit = -0.100E+01</div><div>sim_fluffdampcoeff = 0.950E+00</div><div>sim_fluffdampcutoff = 0.100E-17 [CHANGED]</div><div>sim_fluffrefinecutoff = 0.100E+02</div><div>sim_fluidgamma = 0.167E+01 [CHANGED]</div><div>sim_objcentdens = 0.165E+03 [CHANGED]</div><div>sim_objcentpres = 0.246E+18 [CHANGED]</div><div>sim_objmass = 0.999E+00 [CHANGED]</div><div>sim_objpolyn = 0.300E+01</div><div>sim_objradius = 0.102E+01 [CHANGED]</div><div>sim_pambient = 0.100E-03</div><div>sim_relaxrate = 0.900E+00 [CHANGED]</div><div>sim_rhoambient = 0.100E-18 [CHANGED]</div><div>sim_rotangle = 0.000E+00</div><div>sim_rotfac = 0.000E+00</div><div>sim_smallt = 0.100E+05</div><div>sim_starptmass = 0.200E+34</div><div>sim_tambient = 0.100E+01</div><div>sim_trelax = 0.100E+05 [CHANGED]</div><div>sim_tspinup = 0.000E+00</div><div>sim_winddens = 0.100E-02 [CHANGED]</div><div>sim_windduration = 0.250E+04 [CHANGED]</div><div>sim_windpeak = 0.300E+05 [CHANGED]</div><div>sim_windrampup = 0.300E+01</div><div>sim_windtemp = 0.100E+05 [CHANGED]</div><div>sim_windvel = 0.300E+08</div><div>sim_xcenter = 0.100E+12 [CHANGED]</div><div>sim_ycenter = 0.100E+12 [CHANGED]</div><div>sim_zcenter = 0.100E+12 [CHANGED]</div><div>small = 0.100E-14 [CHANGED]</div><div>smalle = 0.100E+11 [CHANGED]</div><div>smallp = 0.100E-03 [CHANGED]</div><div>smallt = 0.100E-14 [CHANGED]</div><div>smallu = 0.100E-14 [CHANGED]</div><div>smallx = 0.100E-29 [CHANGED]</div><div>smlrho = 0.100E-19 [CHANGED]</div><div>tinitial = 0.000E+00</div><div>tmax = 0.100E+01 [CHANGED]</div><div>tstep_change_factor = 0.200E+01</div><div>vgrid = 0.000E+00</div><div>wall_clock_checkpoint = 0.432E+05</div><div>wall_clock_time_limit = 0.605E+06</div><div>x_refine_center = 0.000E+00</div><div>xmax = 0.200E+12 [CHANGED]</div><div>xmin = 0.000E+00</div><div>y_refine_center = 0.000E+00</div><div>ymax = 0.900E+02 [CHANGED]</div><div>ymin = 0.000E+00</div><div>zfinal = 0.000E+00</div><div>zinitial = -0.100E+01</div><div>z_refine_center = 0.000E+00</div><div>zmax = 0.180E+03 [CHANGED]</div><div>zmin = 0.000E+00</div><div>unitsystem = none </div><div>basenm = starwind_ [CHANGED]</div><div>eosmode = dens_ie </div><div>eosmodeinit = dens_pres [CHANGED]</div><div>geometry = spherical [CHANGED]</div><div>gr_pmrpoutputdir = ./ </div><div>grav_boundary_type = isolated </div><div>hy_eosmodeafter = dens_ie </div><div>hy_fluxrepresentation = fluxes </div><div>log_file = starwind.log [CHANGED]</div><div>mpole_zonetype_1 = exponential </div><div>mpole_zonetype_2 = exponential </div><div>mpole_zonetype_3 = exponential </div><div>mpole_zonetype_4 = logarithmic </div><div>output_directory = </div><div>pc_unitsbase = CGS </div><div>plot_grid_var_1 = none </div><div>plot_grid_var_10 = none </div><div>plot_grid_var_11 = none </div><div>plot_grid_var_12 = none </div><div>plot_grid_var_2 = none </div><div>plot_grid_var_3 = none </div><div>plot_grid_var_4 = none </div><div>plot_grid_var_5 = none </div><div>plot_grid_var_6 = none </div><div>plot_grid_var_7 = none </div><div>plot_grid_var_8 = none </div><div>plot_grid_var_9 = none </div><div>plot_var_1 = dens [CHANGED]</div><div>plot_var_10 = none </div><div>plot_var_11 = none </div><div>plot_var_12 = none </div><div>plot_var_2 = temp [CHANGED]</div><div>plot_var_3 = velx [CHANGED]</div><div>plot_var_4 = vely [CHANGED]</div><div>plot_var_5 = velz [CHANGED]</div><div>plot_var_6 = none </div><div>plot_var_7 = none </div><div>plot_var_8 = none </div><div>plot_var_9 = none </div><div>prof_file = profile.dat </div><div>refine_var_1 = dens [CHANGED]</div><div>refine_var_2 = pres [CHANGED]</div><div>refine_var_3 = none </div><div>refine_var_4 = none </div><div>run_comment = FLASH 3 run </div><div>run_number = 1 </div><div>stats_file = flash.dat </div><div>xl_boundary_type = reflecting [CHANGED]</div><div>xr_boundary_type = diode [CHANGED]</div><div>yl_boundary_type = periodic </div><div>yr_boundary_type = periodic </div><div>zl_boundary_type = outflow [CHANGED]</div><div>zr_boundary_type = outflow [CHANGED]</div><div>allowdtstsdominate = F</div><div>alwayscomputeuservars = T</div><div>alwaysrestrictcheckpoint = T</div><div>bytepack = F</div><div>charlimiting = T</div><div>chkguardcellsinput = F</div><div>chkguardcellsoutput = F</div><div>converttoconsvdformeshcalls = F</div><div>converttoconsvdinmeshinterp = T</div><div>corners = F</div><div>dr_printtsteploc = T</div><div>dr_shortenlaststepbeforetmax = T [CHANGED]</div><div>drift_tuples = F</div><div>eachprocwritesownabortlog = F</div><div>eachprocwritessummary = F</div><div>earlyblockdistadjustment = T</div><div>enablemaskedgcfill = F</div><div>flux_correct = T</div><div>geometryoverride = F</div><div>gr_bcenableapplymixedgds = T</div><div>gr_lrefinemaxbytime = F</div><div>gr_lrefinemaxreddobylogr = F</div><div>gr_lrefinemaxreddobytime = F</div><div>gr_pmrpadvancealllevels = F</div><div>gr_pmrpamrerrorchecking = F</div><div>gr_pmrpcartesianpm = F</div><div>gr_pmrpconserve = F</div><div>gr_pmrpconsvfluxdensities = T</div><div>gr_pmrpconsvfluxes = F</div><div>gr_pmrpcurvilinear = F</div><div>gr_pmrpcurvilinearconserve = F</div><div>gr_pmrpcylindricalpm = F</div><div>gr_pmrpdiagonals = T</div><div>gr_pmrpedgevalue = T</div><div>gr_pmrpedgevalueinteg = F</div><div>gr_pmrpemptycells = F</div><div>gr_pmrpforceconsistency = T</div><div>gr_pmrplsingularline = F</div><div>gr_pmrpnopermanentguardcells = F</div><div>gr_pmrppolarpm = F</div><div>gr_pmrppredcorr = F</div><div>gr_pmrpsphericalpm = F</div><div>gr_pmrptimingmpi = F</div><div>gr_pmrptimingmpix = F</div><div>gr_pmrpvardt = F</div><div>hybrid_riemann = F</div><div>ignoreforcedplot = F</div><div>io_writemscalarintegrals = F</div><div>leveque = F</div><div>mpole_2dsymmetryplane = F</div><div>mpole_3daxisymmetry = F</div><div>mpole_dumpmoments = F</div><div>mpole_ignoreinnerzone = F</div><div>mpole_multithreading = F</div><div>mpole_printradialinfo = F</div><div>plotfilegridquantitydp = F</div><div>plotfilemetadatadp = F</div><div>ppm_modifystates = F</div><div>reducegcellfills = F</div><div>refine_on_particle_count = F</div><div>restart = F</div><div>summaryoutputonly = F</div><div>threadblocklistbuild = F</div><div>threaddelepblocklist = F</div><div>threaddelepwithinblock = F</div><div>threaddriverblocklist = F</div><div>threaddriverwithinblock = F</div><div>threadeoswithinblock = F</div><div>threadhydroblocklist = F</div><div>threadhydrowithinblock = F</div><div>threadraytracebuild = F</div><div>threadwithinblockbuild = F</div><div>typematchedxfer = T</div><div>unbiased_geometry = F</div><div>updatehydrofluxes = T</div><div>useburn = F</div><div>usecollectivehdf5 = T</div><div>useconductivity = F</div><div>usecool = F</div><div>usecosmology = F</div><div>usedeleptonize = F</div><div>usediffuse = F</div><div>useenergydeposition = F</div><div>useflame = F</div><div>usegravity = F</div><div>useheat = F</div><div>useheatexchange = F</div><div>usehydro = T</div><div>useionize = F</div><div>uselegacylabels = T</div><div>usemagneticresistivity = F</div><div>usemassdiffusivity = F</div><div>useopacity = F</div><div>useparticles = F</div><div>usepolytrope = F</div><div>useprimordialchemistry = F</div><div>useradtrans = F</div><div>usests = F</div><div>usestsfordiffusion = F</div><div>usestir = F</div><div>usetreecol = F</div><div>useturb = T</div><div>useviscosity = F</div><div>use_cma_advection = F</div><div>use_cma_flattening = F</div><div>use_cma_steepening = F</div><div>use_steepening = T</div><div>writestatsummary = T</div><div> </div><div> ==============================================================================</div><div> </div><div> Known units of measurement:</div><div> </div><div> Unit CGS Value Base Unit</div><div> 1 cm 1.0000 cm</div><div> 2 s 1.0000 s</div><div> 3 g 1.0000 g</div><div> 4 K 1.0000 K</div><div> 5 esu 1.0000 esu</div><div> 6 mol 1.0000 mol</div><div> 7 m 100.00 cm</div><div> 8 km 1.00000E+05 cm</div><div> 9 pc 3.08568E+18 cm</div><div> 10 kpc 3.08568E+21 cm</div><div> 11 Mpc 3.08568E+24 cm</div><div> 12 Gpc 3.08568E+27 cm</div><div> 13 Rsun 6.96000E+10 cm</div><div> 14 AU 1.49598E+13 cm</div><div> 15 yr 3.15569E+07 s</div><div> 16 Myr 3.15569E+13 s</div><div> 17 Gyr 3.15569E+16 s</div><div> 18 kg 1000.0 g</div><div> 19 Msun 1.98892E+33 g</div><div> 20 amu 1.66054E-24 g</div><div> 21 eV 11605. K</div><div> 22 C 2.99792E+09 esu</div><div> 23 LFLY 3.08568E+24 cm</div><div> 24 TFLY 2.05759E+17 s</div><div> 25 MFLY 9.88470E+45 g</div><div> 26 clLength 3.08568E+24 cm</div><div> 27 clTime 3.15569E+16 s</div><div> 28 clMass 1.98892E+48 g</div><div> 29 clTemp 1.16044E+07 K</div><div>-----------End of Units--------------------</div><div> </div><div> Known physical constants:</div><div> </div><div> Constant Name Constant Value cm s g K esu mol</div><div> 1 Newton 6.67428E-08 3.0 -2.0 -1.0 0.0 0.0 0.0 </div><div> 2 speed of light 2.99792E+10 1.0 -1.0 0.0 0.0 0.0 0.0 </div><div> 3 Planck 6.62607E-27 2.0 -1.0 1.0 0.0 0.0 0.0 </div><div> 4 electron charge 4.80320E-10 0.0 0.0 0.0 0.0 1.0 0.0 </div><div> 5 electron mass 9.10938E-28 0.0 0.0 1.0 0.0 0.0 0.0 </div><div> 6 proton mass 1.67262E-24 0.0 0.0 1.0 0.0 0.0 0.0 </div><div> 7 fine-structure 7.29735E-03 0.0 0.0 0.0 0.0 0.0 0.0 </div><div> 8 Avogadro 6.02214E+23 0.0 0.0 0.0 0.0 0.0 -1.0 </div><div> 9 Boltzmann 1.38065E-16 2.0 -2.0 1.0 -1.0 0.0 0.0 </div><div> 10 ideal gas constant 8.31447E+07 2.0 -2.0 1.0 -1.0 0.0 -1.0 </div><div> 11 Wien 0.28978 1.0 0.0 0.0 1.0 0.0 0.0 </div><div> 12 Stefan-Boltzmann 5.67040E-05 0.0 -3.0 1.0 -4.0 0.0 0.0 </div><div> 13 Radiation Constant 7.56577E-15 -1.0 -2.0 1.0 -4.0 0.0 0.0 </div><div> 14 pi 3.1416 0.0 0.0 0.0 0.0 0.0 0.0 </div><div> 15 e 2.7183 0.0 0.0 0.0 0.0 0.0 0.0 </div><div> 16 Euler 0.57722 0.0 0.0 0.0 0.0 0.0 0.0 </div><div> ==============================================================================</div><div> </div><div> Multifluid database contents:</div><div> </div><div>Initially defined values of species:</div><div>Name Index Total Positive Neutral Negative bind Ener Gamma eosType</div><div><br></div><div>Name Index OpacityLowTemp Zmin Subtype ZFreeTableFile EnerTableFile PresTableFile</div><div><br></div><div>Species Constituents</div><div> </div><div> ==============================================================================</div><div> [ 07-07-2014 10:26:23.962 ] [amr_set_runtime_parameters] NOTE: Enabling curvilinear suppport because of the selected geometry</div><div> [ 07-07-2014 10:26:23.999 ] [gr_initGeometry] checking BCs for idir: 1</div><div> [ 07-07-2014 10:26:23.999 ] [gr_initGeometry] checking BCs for idir: 2</div><div> [ 07-07-2014 10:26:23.999 ] [gr_initGeometry] checking BCs for idir: 3</div><div> Object mass (m_sun): 9.98681760E-01</div><div> Object radius: 7.08274944E+10</div><div> Object central density: 1.65245400E+02</div><div> Ambient CFL timestep: 6.12372436E+01</div><div> Fluff CFL timestep: 1.93649167E+01</div><div> Obj. radius: 7.08274944E+10</div><div> [ 07-07-2014 10:26:24.007 ] [GRID amr_refine_derefine]: initiating refinement</div><div> [ 07-07-2014 10:26:24.007 ] [GRID amr_refine_derefine]: redist. phase. tot blks requested: 9</div><div> [GRID amr_refine_derefine] min blks 9 max blks 9 tot blks 9</div><div> [GRID amr_refine_derefine] min leaf blks 8 max leaf blks 8 tot leaf blks 8</div><div> [ 07-07-2014 10:26:24.008 ] [GRID amr_refine_derefine]: refinement complete</div><div> [ 07-07-2014 10:26:24.010 ] [GRID gr_expandDomain]: iteration=1, create level=3</div><div> INFO: Grid_fillGuardCells is ignoring masking.</div><div> [ 07-07-2014 10:26:24.097 ] [mpi_amr_comm_setup]: buffer_dim_send=1, buffer_dim_recv=1</div><div> [ 07-07-2014 10:26:24.139 ] WARNING after gc filling: min. unk(DENS_VAR)=0.000000000000000000000 PE=0 block=1 type=2</div><div> [ 07-07-2014 10:26:24.140 ] WARNING after gc filling: min. unk(ENER_VAR)=0.000000000000000 PE=0 block=1 type=2</div><div> [ 07-07-2014 10:26:24.140 ] WARNING after gc filling: min. unk(EINT_VAR)=0.000000000000000 PE=0 block=1 type=2</div><div> [ 07-07-2014 10:26:24.141 ] WARNING after gc filling: min. unk(DENS_VAR)=0.000000000000000000000 PE=0 block=2 type=1</div><div> [ 07-07-2014 10:26:24.141 ] WARNING after gc filling: min. unk(ENER_VAR)=0.000000000000000 PE=0 block=2 type=1</div><div> [ 07-07-2014 10:26:24.142 ] WARNING after gc filling: min. unk(EINT_VAR)=0.000000000000000 PE=0 block=2 type=1</div><div> [ 07-07-2014 10:26:24.142 ] WARNING after gc filling: min. unk(DENS_VAR)=0.000000000000000000000 PE=0 block=3 type=1</div><div> [ 07-07-2014 10:26:24.142 ] WARNING after gc filling: min. unk(ENER_VAR)=0.000000000000000 PE=0 block=3 type=1</div><div> [ 07-07-2014 10:26:24.143 ] WARNING after gc filling: min. unk(EINT_VAR)=0.000000000000000 PE=0 block=3 type=1</div><div> [ 07-07-2014 10:26:24.143 ] WARNING after gc filling: min. unk(DENS_VAR)=0.000000000000000000000 PE=0 block=4 type=1</div><div> [ 07-07-2014 10:26:24.143 ] WARNING after gc filling: min. unk(ENER_VAR)=0.000000000000000 PE=0 block=4 type=1</div><div> [ 07-07-2014 10:26:24.144 ] WARNING after gc filling: min. unk(EINT_VAR)=0.000000000000000 PE=0 block=4 type=1</div><div> [ 07-07-2014 10:26:24.144 ] WARNING after gc filling: min. unk(DENS_VAR)=0.000000000000000000000 PE=0 block=5 type=1</div><div> [ 07-07-2014 10:26:24.145 ] WARNING after gc filling: min. unk(ENER_VAR)=0.000000000000000 PE=0 block=5 type=1</div><div> [ 07-07-2014 10:26:24.145 ] WARNING after gc filling: min. unk(EINT_VAR)=0.000000000000000 PE=0 block=5 type=1</div><div> [ 07-07-2014 10:26:24.145 ] WARNING after gc filling: min. unk(DENS_VAR)=0.000000000000000000000 PE=0 block=6 type=1</div><div> [ 07-07-2014 10:26:24.146 ] WARNING after gc filling: min. unk(ENER_VAR)=0.000000000000000 PE=0 block=6 type=1</div><div> [ 07-07-2014 10:26:24.146 ] WARNING after gc filling: min. unk(EINT_VAR)=0.000000000000000 PE=0 block=6 type=1</div><div> [ 07-07-2014 10:26:24.146 ] WARNING after gc filling: min. unk(DENS_VAR)=0.000000000000000000000 PE=0 block=7 type=1</div><div> [ 07-07-2014 10:26:24.146 ] WARNING after gc filling: min. unk(ENER_VAR)=0.000000000000000 PE=0 block=7 type=1</div><div> [ 07-07-2014 10:26:24.147 ] WARNING after gc filling: min. unk(EINT_VAR)=0.000000000000000 PE=0 block=7 type=1</div><div> [ 07-07-2014 10:26:24.147 ] WARNING after gc filling: min. unk(DENS_VAR)=0.000000000000000000000 PE=0 block=8 type=1</div><div> [ 07-07-2014 10:26:24.147 ] WARNING after gc filling: min. unk(ENER_VAR)=0.000000000000000 PE=0 block=8 type=1</div><div> [ 07-07-2014 10:26:24.148 ] WARNING after gc filling: min. unk(EINT_VAR)=0.000000000000000 PE=0 block=8 type=1</div><div> [ 07-07-2014 10:26:24.148 ] WARNING after gc filling: min. unk(DENS_VAR)=0.000000000000000000000 PE=0 block=9 type=1</div><div> [ 07-07-2014 10:26:24.148 ] WARNING after gc filling: min. unk(ENER_VAR)=0.000000000000000 PE=0 block=9 type=1</div><div> [ 07-07-2014 10:26:24.148 ] WARNING after gc filling: min. unk(EINT_VAR)=0.000000000000000 PE=0 block=9 type=1</div></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><br><div><div>On Jul 7, 2014, at 10:17 AM, Sean Couch <<a href="mailto:smc@flash.uchicago.edu">smc@flash.uchicago.edu</a>> wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><meta http-equiv="Content-Type" content="text/html charset=windows-1252"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">Hi Luke, <div><br></div><div>Could you send along your log file so that we can see exactly what setup and runtime parameters you’re using?</div><div><br></div><div>Cheers,</div><div>Sean<br><div apple-content-edited="true">
<div style="font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><span class="Apple-style-span" style="font-family: Helvetica; border-collapse: separate; border-spacing: 0px;"><span class="Apple-style-span" style="border-collapse: separate; font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; border-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-stroke-width: 0px;"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div><br class="Apple-interchange-newline">--------------------------------------------------------</div><div>Sean M. Couch</div><div>Hubble Fellow</div><div>Flash Center for Computational Science</div><div>Department of Astronomy & Astrophysics</div><div>The University of Chicago</div><div>5747 S Ellis Ave, Jo 315</div><div>Chicago, IL 60637</div><div>(773) 702-3899</div><div><a href="http://www.flash.uchicago.edu/~smc">www.flash.uchicago.edu/~smc</a></div><div><br></div></div></span></span></div><br class="Apple-interchange-newline"><br class="Apple-interchange-newline">
</div>
<br><div><div>On Jul 6, 2014, at 11:06 AM, Luke Zoltan Kelley <<a href="mailto:lkelley@cfa.harvard.edu">lkelley@cfa.harvard.edu</a>> wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">Hello flash-users!<br><br>I'm new to FLASH and trying to convert a simulation from cartesian to (3D) spherical, but I'm running into some errors. If anyone has had recent success using 3D spherical, I would be very appreciative of any tips they had, or issues they overcame. In particular, I'm getting a warning from source/Grid/GridMain/paramesh/paramesh4/gr_sanitizeDataAfterInterp.F90 that my guard-cells aren't being set correctly:<br><br><br>WARNING after gc filling: min. unk(DENS_VAR)=0.000000000000000000000 PE=0 block=1 type=2 <br>16 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 <br>15 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 <br>14 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 <br>13 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 <br>12 0.10E-18 0.63E-01 0.63E-01 0.62E-01 0.62E-01 0.63E-01 0.63E-01 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18<br>11 0.10E-18 0.63E-01 0.62E-01 0.63E-01 0.63E-01 0.62E-01 0.63E-01 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18<br>10 0.10E-18 0.62E-01 0.62E-01 0.63E-01 0.63E-01 0.62E-01 0.62E-01 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18<br> 9 0.10E-18 0.63E-01 0.63E-01 0.62E-01 0.62E-01 0.63E-01 0.63E-01 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18<br> 8 0.10E-18 0.62E-01 0.63E-01 0.63E-01 0.63E-01 0.63E-01 0.62E-01 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18<br> 7 0.10E-18 0.63E-01 0.62E-01 0.62E-01 0.62E-01 0.62E-01 0.63E-01 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18<br> 6 0.10E-18 0.62E-01 0.63E-01 0.62E-01 0.62E-01 0.63E-01 0.62E-01 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18<br> 5 0.10E-18 0.62E-01 0.63E-01 0.63E-01 0.63E-01 0.63E-01 0.62E-01 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18 0.10E-18<br> 4 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 <br> 3 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 <br> 2 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 <br> 1 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 <br><br><br>These warnings are given for DENS_VAR, ENER_VAR, EINT_VAR and all blocks.<br>(presumably) for this reason I'm getting arithmetic exceptions when source/physics/Eos/EosMain/Gamma/eos_idealGamma.F90 tries to divide by the (zero) densities, the backtrace I get is:<br><br><br>Program received signal EXC_ARITHMETIC, Arithmetic exception.<br>0x000000010024cf6e in eos_idealgamma_ (mode=103, veclen=8, eosdata=(), vecbegin=Cannot access memory at address 0x0<br>) at eos_idealGamma.F90:341<br>341<span class="Apple-tab-span" style="white-space:pre"> </span> eosData(dens+ilo:dens+ihi)<br>(gdb) bt<br>#0 0x000000010024cf6e in eos_idealgamma_ (mode=103, veclen=8, eosdata=(), vecbegin=Cannot access memory at address 0x0<br>) at eos_idealGamma.F90:341<br>#1 0x000000010001149b in eos_ (mode=103, veclen=8, eosdata=(), massfrac=(), mask=(.FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE.), vecbegin=Cannot access memory at address 0x0<br>) at Eos.F90:265<br>#2 0x000000010001edea in eos_wrapped_ (mode=103, range=(( 5, 12) ( 1, 4) ( 1, 1) ), blockid=1, griddatastruct=Cannot access memory at address 0x0<br>) at Eos_wrapped.F90:220<br>#3 0x0000000100019b15 in complexskipping.1886 () at Eos_guardCells.F90:43<br>#4 0x000000010001afd8 in eos_guardcells_ (corners=.TRUE., layers=(4, 4, 4), skipsrl=.TRUE.) at Eos_guardCells.F90:120<br>#5 0x00000001000322d7 in grid_fillguardcells_ (griddatastruct=380, idir=-1, minlayers=Cannot access memory at address 0x0<br>) at Grid_fillGuardCells.F90:482<br>#6 0x000000010006340e in grid_markrefinederefine_ () at Grid_markRefineDerefine.F90:98<br>#7 0x000000010029eec2 in gr_expanddomain_ (particlesinitialized=.FALSE.) at gr_expandDomain.F90:212<br>#8 0x0000000100061385 in grid_initdomain_ (restart=.FALSE., particlesinitialized=.FALSE.) at Grid_initDomain.F90:98<br>#9 0x000000010000d405 in driver_initflash_ () at Driver_initFlash.F90:156<br>#10 0x000000010001f19c in flash () at Flash.F90:49<br>#11 0x000000010001f20e in main (argc=1, argv=0x7fff5fbff870 '/Users/lzkelley/Applications/flash/flash4.2.2/object_starwind/flash4\000') at Flash.F90:43<br><br><br>In my Simulation_initBlock.F90, I am looping over the guard-cells and setting all densities to non-zero values --- if I print out the densities during the initialization loop I get the appropriate densities from i,j,k = {1,1,1} up to {16,16,16} --- so it seems like the guard-cells are being overwritten to zero...<br>This does not occur in 3D cartesian, and also does *not* occur when I run SodSpherical (which seems to work properly, strangely enough).<br><br>I believe I've also found a bug in source/physics/Hydro/HydroMain/split/PPM/PPMKernel/avisco.F90:416, when sweeping along "y" (theta) in 3D spherical the denominator [ 'sin(xl(i))' ] is allowed to be zero and causes arithmetic exceptions.<br><br>Any help on this particular problem, or on getting spherical running in general would be greatly appreciated!<br>Thanks,<br>Luke</blockquote></div><br></div></div></blockquote></div><br></div></blockquote></div><br></div></body></html>