Pencil1: iapp, nprocs, ncpus= 0 8 8 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 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= 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.577075619635462 setup_slices: slice_position = d setup_slices: ix_loc,iy_loc,iz_loc, (video files) = 14 11 4 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 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 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) 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.1631313484628354 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.721359820810992 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.12941331054539285 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000E+000 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.0939889748843974 daa_dt: max(diffus_eta2) = 0.0000000000000000E+000 daa_dt: max(diffus_eta3) = 0.0000000000000000E+000 ---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.936E-01 1.554E-05 4.683E-02 9.618E-01 1.187E-01 2.86E-07 4.94317E-01 1.103E-09 1.20E-09 2.44349E-01 2.44349E-01 5.26290E-06 2.29410E-03 9.28843E-03 0.000E+00 0.000E+00 5.26290E-06 -1.40E-07 2.1E-11 -5.66E-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.254E-02 9.937E-01 1.555E-05 4.683E-02 9.617E-01 1.191E-01 2.85E-07 4.94271E-01 1.108E-09 1.45E-09 2.44304E-01 2.44304E-01 5.26252E-06 2.29402E-03 9.28702E-03 5.524E-16 5.291E-16 5.26252E-06 -9.80E-07 2.1E-11 -5.66E-07 -3.03E-17 1.19E-19 1.06E-17 1.42E-11 1.01446E-02 6.92E-05 8.32E-03 1.36E-02 2 0.0250 1.257E-02 9.939E-01 1.556E-05 4.683E-02 9.616E-01 1.196E-01 2.95E-07 4.94217E-01 1.157E-09 2.11E-09 2.44250E-01 2.44250E-01 5.26214E-06 2.29394E-03 9.28554E-03 2.175E-15 2.072E-15 5.26214E-06 -2.11E-06 2.3E-11 -5.65E-07 -1.13E-16 -1.10E-18 4.24E-17 3.02E-11 1.01445E-02 6.92E-05 8.32E-03 1.36E-02 3 0.0376 1.262E-02 9.941E-01 1.557E-05 4.683E-02 9.614E-01 1.202E-01 3.43E-07 4.94151E-01 1.338E-09 3.14E-09 2.44185E-01 2.44185E-01 5.26175E-06 2.29385E-03 9.28417E-03 4.816E-15 4.577E-15 5.26175E-06 -3.59E-06 2.5E-11 -5.64E-07 -2.45E-16 -9.79E-18 1.04E-16 4.72E-11 1.01444E-02 6.92E-05 8.32E-03 1.36E-02 4 0.0502 1.269E-02 9.944E-01 1.558E-05 4.683E-02 9.613E-01 1.212E-01 5.00E-07 4.94072E-01 1.847E-09 4.47E-09 2.44108E-01 2.44108E-01 5.26136E-06 2.29377E-03 9.28284E-03 8.434E-15 8.008E-15 5.26136E-06 -5.50E-06 2.8E-11 -5.63E-07 -4.27E-16 -2.71E-17 1.99E-16 6.22E-11 1.01442E-02 6.92E-05 8.32E-03 1.36E-02 5 0.0629 1.279E-02 9.946E-01 1.559E-05 4.683E-02 9.612E-01 1.225E-01 9.15E-07 4.93980E-01 3.075E-09 6.00E-09 2.44016E-01 2.44016E-01 5.26096E-06 2.29368E-03 9.28150E-03 1.300E-14 1.233E-14 5.26096E-06 -7.92E-06 3.1E-11 -5.64E-07 -6.54E-16 -4.76E-17 3.20E-16 7.32E-11 1.01441E-02 6.92E-05 8.32E-03 1.36E-02 6 0.0757 1.292E-02 9.949E-01 1.561E-05 4.683E-02 9.611E-01 1.244E-01 1.88E-06 4.93870E-01 5.743E-09 7.67E-09 2.43908E-01 2.43908E-01 5.26055E-06 2.29359E-03 9.28021E-03 1.849E-14 1.751E-14 5.26055E-06 -1.09E-05 3.5E-11 -5.65E-07 -9.31E-16 -6.96E-17 4.88E-16 7.98E-11 1.01439E-02 6.92E-05 8.32E-03 1.36E-02 7 0.0886 1.300E-02 9.952E-01 1.562E-05 4.683E-02 9.609E-01 1.271E-01 3.92E-06 4.93741E-01 1.113E-08 9.40E-09 2.43780E-01 2.43780E-01 5.26013E-06 2.29350E-03 9.27885E-03 2.490E-14 2.354E-14 5.26013E-06 -1.47E-05 4.2E-11 -5.65E-07 -1.27E-15 -9.99E-17 7.13E-16 8.42E-11 1.01438E-02 6.92E-05 8.32E-03 1.36E-02 8 0.1016 1.300E-02 9.955E-01 1.563E-05 4.683E-02 9.608E-01 1.310E-01 7.96E-06 4.93591E-01 2.131E-08 1.12E-08 2.43632E-01 2.43632E-01 5.25970E-06 2.29340E-03 9.27735E-03 3.217E-14 3.037E-14 5.25970E-06 -1.92E-05 5.2E-11 -5.63E-07 -1.66E-15 -1.43E-16 9.61E-16 8.75E-11 1.01437E-02 6.92E-05 8.32E-03 1.36E-02 9 0.1146 1.300E-02 9.958E-01 1.564E-05 4.683E-02 9.607E-01 1.365E-01 1.54E-05 4.93417E-01 3.946E-08 1.31E-08 2.43460E-01 2.43460E-01 5.25926E-06 2.29331E-03 9.27571E-03 4.022E-14 3.792E-14 5.25926E-06 -2.46E-05 6.4E-11 -5.58E-07 -2.09E-15 -1.95E-16 1.22E-15 9.31E-11 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.12762353338483248 Wall clock time [hours] = 5.326E-04 (+/- 2.7778E-10) Wall clock time/timestep/meshpoint [microsec] = 23.40576 Maximum used memory per cpu [MBytes] = 13.109 Maximum used memory [MBytes] = 104.332 real 4.65 user 5.28 sys 0.92