SVN: ------- v. ( ) $Id$
SVN: start.in,v v. 1.1 (2014/04/24 15:15:40) wlyra
SVN: run.in,v v. 1.1 (2014/04/24 15:15:40) wlyra
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= 14 (ldebug= F )
This is a 3-D run
nxgrid, nygrid, nzgrid= 32 16 16
Lx, Ly, Lz= 3.8750000000000000 0.60000000000000009 0.59999999999999998
Vbox= 13.577075619635460
rsnap: read snapshot var.dat in 0.0000000000000000 seconds
setup_slices: slice_position = d
setup_slices: ix_loc,iy_loc,iz_loc,iy2_loc, (video files) = 14 14 4 11
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= 1.2027220936797545E-008
units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 -5.3185200738655558 7.0000000000000007E-002 4.9000000000000016E-003 1.0000000000000000
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
WARNING: initialize_diagnostics: volume averages are calculated wrongly for nonequidistant grids unless lproper_averages=T.!
initialize_gravity: non-smoothed newtonian gravity
diffusion: mesh hyperdiffusion
diffusion: shock diffusion
WARNING: initialize_density: For diffusion energy/momentum correction should use lmassdiff_fix=T!
initialize_density: no need to read initial stratification for lanti_shockdiffusion=F.
WARNING: initialize_energy: llocal_iso=T. Make sure you have the appropriate INITIAL_CONDITION in Makefile.local!
select_eos_variable: Using rho and cs2
resistivity: hyper3 resolution-invariant
resistivity: shock
resistivity: x-dependent
viscous force: nu_hyper3_mesh/pi^5 *(Deltav)^6/Deltaq
viscous force: nu_shock*(XXXXXXXXXXX)
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_hydro: call gij_etc
calc_pencils_energy: max(advec_cs2) = 9.1631313484622883
duu_dt: SOLVE
Bcs for ux, x: < a>, y: < a>, z: < p>
Bcs for uy, x: < s>, y: < out>, z: < p>
Bcs for uz, x: < e3>, y: < s>, z: < p>
duu_dt: max(advec_uu) = 27.718808614853046
dlnrho_dt: SOLVE
Bcs for lnrho, x: < e3>, y: < str>, z: < p>
dlnrho_dt: diffrho_shock= 10.000000000000000
dlnrho_dt: diffrho_hyper3_mesh= 5.0000000000000000
dlnrho_dt: max(diffus_diffrho ) = 0.10171453368340980
dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000
daa_dt: SOLVE
Bcs for Ax, x: < s>, y: < a>, z: < p>
Bcs for Ay, x: < a2>, y: < s>, z: < p>
Bcs for Az, x: < a>, y: < a>, z: < p>
daa_dt: iresistivity=hyper3-mesh shock xdep
daa_dt: max(diffus_eta) = 1.0662901980223569
daa_dt: max(diffus_eta2) = 0.0000000000000000
daa_dt: max(diffus_eta3) = 0.0000000000000000
---it-------t---------dt--------umax------rhomin------rhom------rhomax------orms-------oz2m-------urms--------ux2m-------uy2m-------uz2m----------u2m----------b2m---------brms---------bmax--------bx2m-------by2m--------bz2m-------uxuzm-----uxuym----uyuzm-----bxbzm-----bxbym-----bybzm------axm--------amax--------vA2m-----vArms-----vAmax---
0 0.0000 1.251E-02 9.935E-01 1.554E-05 4.683E-02 9.618E-01 1.187E-01 2.82E-07 4.94317E-01 1.095E-09 1.18E-09 2.44350E-01 2.44350E-01 5.26290E-06 2.29410E-03 9.28843E-03 0.000E+00 0.000E+00 5.26290E-06 2.89E-08 2.1E-11 3.35E-07 0.00E+00 0.00E+00 0.00E+00 0.00E+00 1.01448E-02 6.92E-05 8.32E-03 1.36E-02
1 0.0125 1.251E-02 9.936E-01 1.555E-05 4.683E-02 9.617E-01 1.187E-01 2.81E-07 4.94317E-01 1.095E-09 1.44E-09 2.44350E-01 2.44350E-01 5.26252E-06 2.29402E-03 9.28720E-03 5.234E-16 4.641E-16 5.26252E-06 -6.83E-07 2.2E-11 3.35E-07 -1.93E-17 1.46E-18 -1.50E-17 2.29E-11 1.01446E-02 6.92E-05 8.32E-03 1.36E-02
2 0.0250 1.251E-02 9.936E-01 1.556E-05 4.683E-02 9.616E-01 1.187E-01 2.81E-07 4.94317E-01 1.104E-09 2.15E-09 2.44350E-01 2.44350E-01 5.26214E-06 2.29394E-03 9.28599E-03 2.062E-15 1.820E-15 5.26214E-06 -1.39E-06 2.3E-11 3.34E-07 -7.95E-17 6.29E-18 -5.97E-17 4.21E-11 1.01445E-02 6.92E-05 8.32E-03 1.36E-02
3 0.0375 1.251E-02 9.936E-01 1.557E-05 4.683E-02 9.615E-01 1.187E-01 2.82E-07 4.94317E-01 1.122E-09 3.32E-09 2.44350E-01 2.44350E-01 5.26175E-06 2.29385E-03 9.28473E-03 4.568E-15 4.014E-15 5.26175E-06 -2.10E-06 2.3E-11 3.34E-07 -1.74E-16 1.33E-17 -1.33E-16 6.27E-11 1.01444E-02 6.92E-05 8.32E-03 1.36E-02
4 0.0500 1.251E-02 9.936E-01 1.558E-05 4.683E-02 9.614E-01 1.187E-01 2.84E-07 4.94317E-01 1.149E-09 4.92E-09 2.44350E-01 2.44350E-01 5.26137E-06 2.29377E-03 9.28335E-03 7.994E-15 6.983E-15 5.26137E-06 -2.80E-06 2.4E-11 3.33E-07 -3.03E-16 2.39E-17 -2.40E-16 8.67E-11 1.01442E-02 6.92E-05 8.32E-03 1.36E-02
5 0.0625 1.251E-02 9.935E-01 1.559E-05 4.683E-02 9.613E-01 1.187E-01 2.87E-07 4.94317E-01 1.184E-09 6.96E-09 2.44350E-01 2.44350E-01 5.26098E-06 2.29368E-03 9.28195E-03 1.230E-14 1.068E-14 5.26098E-06 -3.50E-06 2.4E-11 3.33E-07 -4.69E-16 4.16E-17 -3.86E-16 1.08E-10 1.01441E-02 6.92E-05 8.32E-03 1.36E-02
6 0.0751 1.251E-02 9.936E-01 1.560E-05 4.683E-02 9.612E-01 1.187E-01 2.90E-07 4.94317E-01 1.227E-09 9.42E-09 2.44350E-01 2.44350E-01 5.26058E-06 2.29360E-03 9.28084E-03 1.743E-14 1.510E-14 5.26058E-06 -4.19E-06 2.5E-11 3.32E-07 -6.74E-16 6.93E-17 -5.62E-16 1.34E-10 1.01439E-02 6.92E-05 8.32E-03 1.36E-02
7 0.0876 1.250E-02 9.936E-01 1.561E-05 4.683E-02 9.610E-01 1.187E-01 2.94E-07 4.94317E-01 1.278E-09 1.23E-08 2.44350E-01 2.44350E-01 5.26018E-06 2.29351E-03 9.27984E-03 2.335E-14 2.018E-14 5.26018E-06 -4.88E-06 2.5E-11 3.32E-07 -9.09E-16 1.07E-16 -7.61E-16 1.72E-10 1.01438E-02 6.92E-05 8.32E-03 1.36E-02
8 0.1001 1.250E-02 9.935E-01 1.563E-05 4.683E-02 9.608E-01 1.187E-01 2.99E-07 4.94317E-01 1.336E-09 1.56E-08 2.44349E-01 2.44350E-01 5.25978E-06 2.29342E-03 9.27873E-03 3.002E-14 2.588E-14 5.25978E-06 -5.56E-06 2.6E-11 3.31E-07 -1.18E-15 1.56E-16 -9.88E-16 2.18E-10 1.01437E-02 6.92E-05 8.32E-03 1.36E-02
9 0.1126 1.251E-02 9.935E-01 1.564E-05 4.683E-02 9.607E-01 1.187E-01 3.04E-07 4.94317E-01 1.401E-09 1.94E-08 2.44349E-01 2.44349E-01 5.25937E-06 2.29333E-03 9.27745E-03 3.740E-14 3.212E-14 5.25937E-06 -6.23E-06 2.6E-11 3.31E-07 -1.49E-15 2.21E-16 -1.24E-15 2.60E-10 1.01435E-02 6.92E-05 8.32E-03 1.36E-02
Simulation finished after 10 time-steps
Writing final snapshot at time t = 0.12507011389665129
Wall clock time [hours] = 1.731E-04 (+/- 2.7778E-07)
Wall clock time/timestep/meshpoint [microsec] = 7.604980
Wall clock time/timestep/local meshpoint [microsec] = 7.604980
Rhs wall clock time/timestep/local meshpoint [microsec] = 10.30273
Maximum used memory per cpu [MBytes] = 21.500
Maximum used memory [MBytes] = 21.500
real 1.24
user 1.11
sys 0.01