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 rsnap: read snapshot var.dat in 6.3571929931640625E-003 seconds 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 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.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.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.0939889748843974 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.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.251E-02 9.936E-01 1.555E-05 4.683E-02 9.617E-01 1.187E-01 2.85E-07 4.94317E-01 1.105E-09 1.45E-09 2.44349E-01 2.44349E-01 5.26252E-06 2.29402E-03 9.28702E-03 5.522E-16 5.293E-16 5.26252E-06 -8.53E-07 2.1E-11 -5.67E-07 -3.12E-17 1.40E-20 1.09E-17 1.46E-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.85E-07 4.94317E-01 1.116E-09 2.16E-09 2.44349E-01 2.44349E-01 5.26214E-06 2.29394E-03 9.28555E-03 2.170E-15 2.074E-15 5.26214E-06 -1.56E-06 2.1E-11 -5.67E-07 -1.18E-16 -8.40E-19 4.33E-17 3.23E-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.614E-01 1.187E-01 2.86E-07 4.94317E-01 1.136E-09 3.31E-09 2.44349E-01 2.44349E-01 5.26176E-06 2.29385E-03 9.28418E-03 4.795E-15 4.572E-15 5.26176E-06 -2.27E-06 2.0E-11 -5.68E-07 -2.55E-16 -7.24E-18 1.01E-16 5.08E-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.613E-01 1.187E-01 2.88E-07 4.94317E-01 1.165E-09 4.91E-09 2.44349E-01 2.44349E-01 5.26137E-06 2.29377E-03 9.28286E-03 8.372E-15 7.967E-15 5.26137E-06 -2.97E-06 2.0E-11 -5.68E-07 -4.41E-16 -2.02E-17 1.86E-16 6.20E-11 1.01442E-02 6.92E-05 8.32E-03 1.36E-02 5 0.0625 1.251E-02 9.936E-01 1.559E-05 4.683E-02 9.612E-01 1.187E-01 2.91E-07 4.94317E-01 1.202E-09 6.94E-09 2.44349E-01 2.44349E-01 5.26098E-06 2.29368E-03 9.28154E-03 1.285E-14 1.221E-14 5.26098E-06 -3.67E-06 1.9E-11 -5.69E-07 -6.77E-16 -3.87E-17 2.96E-16 6.38E-11 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.611E-01 1.187E-01 2.94E-07 4.94317E-01 1.247E-09 9.39E-09 2.44349E-01 2.44349E-01 5.26058E-06 2.29360E-03 9.28028E-03 1.819E-14 1.726E-14 5.26058E-06 -4.37E-06 1.7E-11 -5.70E-07 -9.61E-16 -6.26E-17 4.38E-16 5.76E-11 1.01439E-02 6.92E-05 8.32E-03 1.36E-02 7 0.0876 1.250E-02 9.935E-01 1.562E-05 4.683E-02 9.610E-01 1.187E-01 2.99E-07 4.94317E-01 1.299E-09 1.23E-08 2.44349E-01 2.44349E-01 5.26018E-06 2.29351E-03 9.27895E-03 2.435E-14 2.308E-14 5.26018E-06 -5.05E-06 1.5E-11 -5.71E-07 -1.30E-15 -9.41E-17 6.10E-16 4.88E-11 1.01438E-02 6.92E-05 8.32E-03 1.36E-02 8 0.1001 1.250E-02 9.936E-01 1.563E-05 4.683E-02 9.608E-01 1.187E-01 3.03E-07 4.94317E-01 1.359E-09 1.56E-08 2.44349E-01 2.44349E-01 5.25978E-06 2.29342E-03 9.27752E-03 3.127E-14 2.960E-14 5.25978E-06 -5.74E-06 1.3E-11 -5.71E-07 -1.68E-15 -1.33E-16 7.94E-16 3.98E-11 1.01437E-02 6.92E-05 8.32E-03 1.36E-02 9 0.1126 1.250E-02 9.936E-01 1.564E-05 4.683E-02 9.608E-01 1.187E-01 3.09E-07 4.94317E-01 1.426E-09 1.93E-08 2.44349E-01 2.44349E-01 5.25938E-06 2.29333E-03 9.27595E-03 3.893E-14 3.679E-14 5.25938E-06 -6.41E-06 9.9E-12 -5.72E-07 -2.09E-15 -1.81E-16 9.94E-16 3.91E-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.12507230821533089 Wall clock time [hours] = 5.081E-04 (+/- 2.7778E-10) Wall clock time/timestep/meshpoint [microsec] = 22.32813 Wall clock time/timestep/local meshpoint [microsec] = 178.6250 Rhs wall clock time/timestep/local meshpoint [microsec] = 194.9572 Maximum used memory per cpu [MBytes] = 13.238 Maximum used memory [MBytes] = 104.656 real 5.43 user 5.30 sys 0.67