SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
bcx1,bcx2= a s s a2 a2 a s s s : a s s a2 a2 a s s s
bcy1,bcy2= s a s a2 a2 s a s s : s a s a2 a2 s a s s
bcz1,bcz2= s s a a2 a2 s s a s : s s a a2 a2 s s a s
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= 5.0000000000000000 5.0000000000000000 5.0000000000000000
Vbox= 125.00000000000000
rsnap: read snapshot var.dat in 9.9999998928979039E-004 seconds
rprint_list: nname= 18
rprint_list: ix,iy,iz,iz2= -1 -1 -1 -1
rprint_list: nnamev= 0
rprint_list: nnamez= -1
rprint_list: nnamez= -1
rprint_list: nnamey= -1
rprint_list: nnamex= -1
rprint_list: nnamer= 0
rprint_list: nnamexz= -1
rprint_list: nnamexy= -1
rprint_list: nnamerz= 7
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_shock: run through parse list
rprint_viscosity: run through parse list
units_general: unit_velocity= 1.0000000000000000
units_general: unit_density= 1.0000000000000000
units_general: unit_length= 1.0000000000000000
units_general: unit_magnetic= 3.5449077018110318
units_eos: unit_temperature= 4.8108883747190188E-009
units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 0.40546510810816422 1.0000000000000000 0.59999999999999998 0.40000000000000008
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
timestep based on CFL cond; cdt= 0.40000000000000002
initialize_gravity: smoothed 1/r potential
xmask_den= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000
zmask_den= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000
diffusion: nothing (i.e. no mass diffusion)
initialize_density: no need to read initial stratification for lanti_shockdiffusion=F.
initialize_hydro: outer radius rdampext= 3.0000000000000000
xmask_hyd= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000
zmask_hyd= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000
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.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000
zmask_mag= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000
resistivity: shell
initialize_shock: prenormalised shock_factor sum= 1.0000000000000000
viscous force: nu*(del2u+graddivu/3+2S.glnrho)
viscous force: nu_shock*(XXXXXXXXXXX)
xmask_vis= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000
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
pencil_consistency_check: checking pencil case
pencil_consistency_check: checking requested pencils
pencil_consistency_check: performing small pencil check
pencil_consistency_check: the small pencil check reported no problems
pencil_consistency_check: checking dependence on pencil initialization
pencil_consistency_check: results are independent of pencil initialization
pencil_consistency_check: checking diagnostics pencils
pencil_consistency_check: performing small pencil check
pencil_consistency_check: the small pencil check reported no problems
pencil_consistency_check: set lpencil_check_diagnos_opti=T in run.in to report
pencil_consistency_check: optimisation possibilities for diagnostics pencils
pencil_consistency_check: checking dependence on pencil initialization
pencil_consistency_check: diagnostics are independent of pencil initialization
pencil_consistency_check: all tests passed
pde: ENTER
SVN: ------- v. ( ) $Id$
calc_pencils_energy: max(advec_cs2) = 97.165773730143172
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>
duu_dt: max(advec_uu) = 0.0000000000000000
dlnrho_dt: SOLVE
Bcs for lnrho, x: < a2>, y: < a2>, z: < a2>
dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000
dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000
denergy_dt: SOLVE denergy_dt
Bcs for ss, x: < a2>, y: < a2>, z: < a2>
denergy_dt: lnTT,cs2,cp1= 0.23417276267693732 0.84257521444799810 1.0000000000000000
calc_heatcond_constchi: chi= 2.0000000000000000E-003
calc_heatcond_constchi: added thdiff
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=shell
daa_dt: max(diffus_eta) = 0.23064000000000004
daa_dt: max(diffus_eta2) = 0.0000000000000000
daa_dt: max(diffus_eta3) = 0.0000000000000000
pde: freezing variables for r < 1.0000000000000000 : T T T T T T T T F
----it-------------t------------------dt----------------urms---------------umax----------------brms-----------------bmax-----------------epsK-----------------epsM---------------rhom----------------ssm---------------Marms--------------Mamax---------------vArms----------------vAmax---------------beta1m-------------beta1max------------shockmax------
0 0.00000000000000 3.4497939497199E-02 0.00000000000000 0.00000000000000 1.50660298999331E-01 8.20930484218245E-01 0.00000000000000E+00 1.4294321038045E-03 0.42939736357543 0.26020575672839 0.00000000000000 0.00000000000000 1.73334582581480E-01 8.87864862046937E-01 2.8165751012175E-02 7.2928923995986E-01 0.0000000000000E+00
10 0.29855533670282 2.8257550665232E-02 0.04764564133289 0.40577670415882 1.39183892100697E-01 6.22213189947349E-01 3.15306739590431E-04 6.4535215365133E-04 0.42884797589611 0.26046415600019 0.05101782770303 0.44907090468680 1.65282712316009E-01 7.93195199643802E-01 2.6431654430068E-02 7.0646434777711E-01 1.5673505902845E-02
20 0.57345337016128 2.7283730342637E-02 0.06427752670822 0.47724661391863 1.32621148662506E-01 6.52753636881266E-01 4.00025775709997E-04 5.0067504020646E-04 0.42912870431812 0.26040376021498 0.06955837137219 0.51587043234429 1.57758867539400E-01 7.66829511404838E-01 2.3786291625314E-02 5.8696870057860E-01 2.0640067151489E-02
30 0.85146537533974 2.8531335315585E-02 0.07009142680047 0.40185082925566 1.27698449121874E-01 6.27389573392630E-01 3.47822177363789E-04 4.9360252277450E-04 0.43009770049199 0.26019708656788 0.07581379112222 0.43668790624682 1.52742608432759E-01 7.05307353832752E-01 2.2262302711566E-02 4.8627595728718E-01 1.2286599585020E-02
40 1.14120661364245 2.9435740297643E-02 0.07363472282460 0.35922877242974 1.24145150929859E-01 6.32719199940289E-01 3.25982651050151E-04 4.7744315714213E-04 0.43111509841149 0.25982664053757 0.07971595832275 0.39665690490745 1.49356058324093E-01 6.80523449336110E-01 2.1273237144203E-02 4.3386353791040E-01 1.4291159956566E-02
Simulation finished after 50 time-steps
Writing final snapshot at time t = 1.4375952214241658
Wall clock time [hours] = 2.726E-03 (+/- 2.7778E-07)
Wall clock time/timestep/meshpoint [microsec] = 5.990601
Wall clock time/timestep/local meshpoint [microsec] = 5.990601
Rhs wall clock time/timestep/local meshpoint [microsec] = 6.726685
Maximum used memory per cpu [MBytes] = 32.375
Maximum used memory [MBytes] = 32.375
real 11.48
user 11.46
sys 0.01