=====================================================================
WARNING: using SINGLE PRECISION, you'd better know what you're doing!
=====================================================================
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
bcx1,bcx2= a s s a2 a2 s s a : a s s a2 a2 s s a
bcy1,bcy2= s a s a2 a2 s s a : s a s a2 a2 s s a
bcz1,bcz2= s s a s s s s s : s s a a2 a2 s s a
initialize_mpicomm: enabled MPI
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= 64 64 32
Lx, Ly, Lz= 2.20000005 2.20000005 2.20000005
Vbox= 10.6480007
rsnap: read snapshot var.dat in 4.0471499999999999E-004 seconds
rprint_list: nname= 18
rprint_list: ix,iy,iz,iz2= -1 -1 -1 -1
rprint_list: nnamev= 5
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_viscosity: run through parse list
setup_slices: slice_position = e
setup_slices: ix_loc,iz_loc,iz2_loc, (video files) = 35 4 19
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.19209221E-06 0.816497028 0.400000453 0.399999976
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
timestep based on CFL cond; cdt= 0.400000006
initialize_gravity: smoothed 1/r potential in spherical shell
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 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
diffusion: nothing (i.e. no mass diffusion)
initialize_density: no need to read initial stratification for lanti_shockdiffusion=F.
initialize_hydro: set outer radius rdampext=r_ext
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 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
WARNING: initialize_hydro: need to set orms in print.in to get othresh. Disabled othresh!
select_eos_variable: Using lnrho and ss
initialize_energy: set boundary temperatures for spherical shell problem
initialize_energy: g0,mpoly,beta1 1.20000005 1.00000000 1.50000012
initialize_energy: TT_int, TT_ext= 2.50000000 1.00000000
initialize_energy: cs2_ext, cs2_int= 0.666666627 1.66666651
initialize_energy: nheatc_max,iheatcond= 4 K-const nothing nothing nothing
heat conduction: K=cte
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 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
resistivity: shell
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 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
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) = 1225.75745
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.632456005
duu_dt: max(advec_uu) = 0.999772191
dlnrho_dt: SOLVE
Bcs for lnrho, x: < a2>, y: < a2>, z: < s:a2>
dlnrho_dt: max(diffus_diffrho ) = 0.00000000
dlnrho_dt: max(diffus_diffrho3) = 0.00000000
denergy_dt: SOLVE denergy_dt
Bcs for ss, x: < a2>, y: < a2>, z: < s:a2>
denergy_dt: lnTT,cs2,cp1= 0.00000000 0.666666627 1.00000000
calc_heatcond_constK: hcond= 0.00000000
calc_heatcond_constK: added thdiff
calc_heat_cool: lgravz, lgravr, lgravx, lspherical_coords= F T F F
daa_dt: SOLVE
Bcs for Ax, x: < s>, y: < s>, z: < s>
Bcs for Ay, x: < s>, y: < s>, z: < s>
Bcs for Az, x: < a>, y: < a>, z: < s:a>
daa_dt: iresistivity=shell
daa_dt: max(diffus_eta) = 0.813999176
daa_dt: max(diffus_eta2) = 0.00000000
daa_dt: max(diffus_eta3) = 0.00000000
----it--------t--------dt-------urms------umax------brms------bmax------epsK------epsM-----rhom------ssm------ethm-----dtc-----dtu-----dtb-----dtnu---dtchi---dteta--
0 0.000 7.08E-03 0.01547 0.02717 4.064E-05 1.331E-04 4.402E-05 5.488E-09 1.2077 0.02953 0.9739 0.9802 0.0198 0.0001 0.0027 0.0037 0.0192
2 0.014 7.09E-03 0.01172 0.04736 3.978E-05 1.329E-04 1.981E-05 5.061E-09 1.2076 0.02954 0.9738 0.9851 0.0355 0.0001 0.0027 0.0037 0.0192
4 0.028 7.09E-03 0.01257 0.06671 3.897E-05 1.327E-04 1.563E-05 4.672E-09 1.2075 0.02958 0.9736 0.9865 0.0450 0.0001 0.0027 0.0038 0.0192
6 0.043 7.10E-03 0.01394 0.07307 3.819E-05 1.325E-04 1.715E-05 4.317E-09 1.2073 0.02966 0.9734 0.9863 0.0487 0.0000 0.0028 0.0038 0.0193
8 0.057 7.11E-03 0.01505 0.07590 3.746E-05 1.323E-04 1.844E-05 3.992E-09 1.2072 0.02977 0.9732 0.9890 0.0568 0.0000 0.0028 0.0038 0.0193
10 0.071 7.12E-03 0.01573 0.07496 3.676E-05 1.320E-04 1.738E-05 3.696E-09 1.2070 0.02990 0.9729 0.9905 0.0561 0.0000 0.0028 0.0038 0.0193
12 0.085 7.13E-03 0.01616 0.07487 3.610E-05 1.317E-04 1.615E-05 3.424E-09 1.2069 0.03005 0.9727 0.9917 0.0556 0.0000 0.0028 0.0038 0.0194
14 0.099 7.14E-03 0.01655 0.07907 3.547E-05 1.314E-04 1.585E-05 3.176E-09 1.2067 0.03021 0.9725 0.9930 0.0583 0.0000 0.0028 0.0039 0.0194
16 0.114 7.15E-03 0.01693 0.07939 3.487E-05 1.311E-04 1.588E-05 2.949E-09 1.2066 0.03037 0.9723 0.9942 0.0590 0.0000 0.0028 0.0039 0.0194
Simulation finished after 17 time-steps
Writing final snapshot at time t = 0.12093318905681372
Wall clock time [hours] = 7.710E-04 (+/- 2.7778E-13)
Wall clock time/timestep/meshpoint [microsec] = 1.245670
Wall clock time/timestep/local meshpoint [microsec] = 4.982681
Rhs wall clock time/timestep/local meshpoint [microsec] = 7.060812
Maximum used memory per cpu [MBytes] = 39.453
Maximum used memory [MBytes] = 156.270
real 4.36
user 16.23
sys 0.14