SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
bcx1,bcx2= a s s a2 sT a s s : a s s a2 sT a s s
bcy1,bcy2= s a s a2 sT s a s : s a s a2 sT s a s
bcz1,bcz2= s s a a2 sT s s a : s s a a2 sT s s a
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: gravity_r.f90,v v. 1.1 (2018/08/24 15:48:10) wlyra
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
The verbose level is ip= 13 (ldebug= F )
This is a 3-D run
nxgrid, nygrid, nzgrid= 32 32 32
Lx, Ly, Lz= 3.00000000 3.00000000 3.00000000
Vbox= 27.0000000
rsnap: read snapshot var.dat in 0.0000000000000000 seconds
rprint_list: nname= 13
rprint_list: ix,iy,iz,iz2= -1 -1 11 -1
rprint_list: nnamev= 2
rprint_list: nnamez= 4
rprint_list: nnamez= 4
rprint_list: nnamey= 2
rprint_list: nnamex= 2
rprint_list: nnamer= 2
rprint_list: nnamexz= 2
rprint_list: nnamexy= 2
rprint_list: nnamerz= 10
rprint_general: run through parse list
rprint_hydro: run through parse list
rprint_hydro: run through parse list
rprint_density: run through parse list
rprint_viscosity: run through parse list
setup_slices: slice_position = m
setup_slices: ix_loc,iy_loc,iz_loc,iz2_loc, (video files) = 19 19 14 35
units_general: unit_velocity= 1.0000000000000000
units_general: unit_density= 1.0000000000000000
units_general: unit_length= 1.0000000000000000
units_general: unit_magnetic= 3.5449078083038330
units_eos: unit_temperature= 4.8108880235694808E-009
units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.00000000 1.21639538 1.50000000 1.35000002 0.399999976
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
timestep based on CFL cond; cdt= 0.400000006
TIMEAVG: defaulting idx no. 1 to 1
TIMEAVG: defaulting idx no. 2 to 2
TIMEAVG: defaulting idx no. 3 to 3
TIMEAVG: defaulting idx no. 4 to 4
TIMEAVG: defaulting idx no. 5 to 5
TIMEAVG: defaulting idx no. 6 to 6
TIMEAVG: defaulting idx no. 7 to 7
TIMEAVG: defaulting idx no. 8 to 8
initialize_gravity: M5 dwarf gravity potential
xmask_den= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000
zmask_den= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000
diffusion: nothing (i.e. no mass diffusion)
initialize_density: no need to read initial stratification for lanti_shockdiffusion=F.
initialize_hydro: outer radius rdampext= 1.10000002
xmask_hyd= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000
zmask_hyd= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000
select_eos_variable: Using lnrho and ss
initialize_energy: nheatc_max,iheatcond= 4 chi-const nothing nothing nothing
heat conduction: constant chi
xmask_mag= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000
zmask_mag= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000
resistivity: constant eta
viscous force: nu*(del2u+graddivu/3+2S.glnrho)
xmask_vis= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000
choose_pencils: finding out which pencils are needed for the pencil case.
write_pencil_info: pencil information written to the file pencils.list.
sh: 1: Syntax error: Bad fd number
pde: ENTER
SVN: ------- v. ( ) $Id$
calc_pencils_energy: max(advec_cs2) = 433.583130
duu_dt: SOLVE
Bcs for ux, x: < a>, y: < s>, z: < s>
Bcs for uy, x: < s>, y: < a>, z: < s>
Bcs for uz, x: < s>, y: < s>, z: < a>
coriolis_cartesian: add Coriolis force; Omega= 0.500000000
duu_dt: max(advec_uu) = 0.00000000
dlnrho_dt: SOLVE
Bcs for lnrho, x: < a2>, y: < a2>, z: < a2>
dlnrho_dt: max(diffus_diffrho ) = 0.00000000
dlnrho_dt: max(diffus_diffrho3) = 0.00000000
denergy_dt: SOLVE denergy_dt
Bcs for ss, x: < sT>, y: < sT>, z: < sT>
denergy_dt: lnTT,cs2,cp1= 0.708186388 1.35353708 1.00000000
calc_heatcond_constchi: chi= 1.00000005E-03
calc_heatcond_constchi: added thdiff
calc_heat_cool: lgravz, lgravr, lgravx, lspherical_coords= F T F F
daa_dt: SOLVE
Bcs for Ax, x: < a>, y: < s>, z: < s>
Bcs for Ay, x: < s>, y: < a>, z: < s>
Bcs for Az, x: < s>, y: < s>, z: < a>
daa_dt: iresistivity=eta-const
daa_dt: max(diffus_eta) = 0.192200035
daa_dt: max(diffus_eta2) = 0.00000000
daa_dt: max(diffus_eta3) = 0.00000000
----it--------t--------dt-------urms------umax------brms-------bmax-------rhom------ssm------dtc-----dtu-----dtnu---dtchi--
0 0.000 1.41E-02 0.01140 0.1249 5.999E-06 6.287E-05 0.34611 0.15062 0.9412 0.0779 0.0150 0.0250
10 0.143 1.43E-02 0.00916 0.1024 5.560E-06 5.650E-05 0.34610 0.15063 0.9600 0.0547 0.0152 0.0254
20 0.286 1.44E-02 0.00810 0.0890 5.176E-06 5.133E-05 0.34610 0.15065 0.9659 0.0500 0.0153 0.0256
30 0.430 1.44E-02 0.00726 0.0730 4.832E-06 4.915E-05 0.34610 0.15067 0.9738 0.0413 0.0154 0.0257
40 0.574 1.44E-02 0.00650 0.0682 4.520E-06 4.710E-05 0.34610 0.15068 0.9757 0.0391 0.0154 0.0257
Simulation finished after 41 time-steps
Writing final snapshot at time t = 0.58838927699252963
Wall clock time [hours] = 1.667E-03 (+/- 2.7778E-07)
Wall clock time/timestep/meshpoint [microsec] = 4.465987
Wall clock time/timestep/local meshpoint [microsec] = 4.465987
Rhs wall clock time/timestep/local meshpoint [microsec] = 4.326425
Maximum used memory per cpu [MBytes] = 18.250
Maximum used memory [MBytes] = 18.250
real 6.02
user 6.01
sys 0.01