Pencil1: iapp, nprocs, ncpus= 0 4 4 SVN: ------- v. ( ) $Id$ SVN: start.in,v v. 1.1 (2012-10-04 16:04:16) wlyra SVN: run.in,v v. 1.1 (2012-10-04 16:04:16) 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 32 8 Lx, Ly, Lz= 1.6000000000000001 6.2831853999999998 0.20000000000000001 Vbox= 2.4127431935999999 setup_slices: slice_position = m setup_slices: iz_loc,iz2_loc, (video files) = 7 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 -4.6051701859880909 0.10000000000000001 1.0000000000000004E-002 1.0000000000000000 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) pde: ENTER SVN: ------- v. ( ) $Id$ calc_pencils_hydro: call gij_etc calc_pencils_energy: max(advec_cs2) = 53.437612850949485 duu_dt: SOLVE Bcs for ux, x: < a>, y: < p>, z: < p> Bcs for uy, x: < a2>, y: < p>, z: < p> Bcs for uz, x: < s>, y: < p>, z: < p> duu_dt: max(advec_uu) = 19.878441331084431 dlnrho_dt: SOLVE Bcs for lnrho, x: < s>, y: < p>, z: < p> dlnrho_dt: diffrho_shock= 1.0000000000000000 dlnrho_dt: diffrho_hyper3_mesh= 5.0000000000000000 dlnrho_dt: max(diffus_diffrho ) = 9.9917901892293704E-003 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000E+000 daa_dt: SOLVE Bcs for Ax, x: < a2>, y: < p>, z: < p> Bcs for Ay, x: < a2>, y: < p>, z: < p> Bcs for Az, x: < a>, y: < p>, z: < p> daa_dt: iresistivity=hyper3-mesh shock xdep daa_dt: max(diffus_eta) = 2.8184037178372451E-002 daa_dt: max(diffus_eta2) = 0.0000000000000000E+000 daa_dt: max(diffus_eta3) = 0.0000000000000000E+000 ----it--------t---------dt-------rhom---rhomax--rhomin----u2m-------um2-------urms------umax------uy2m--------b2m---------bm2--------brms---------bmax---------axm---------amax---- 0 0.000 1.40E-02 0.814539 3.95 0.3536 0.812499 2.437499 0.901387 1.561249 0.812499 3.42214E-05 1.625E-04 5.84991E-03 1.2746E-02 0.00000E+00 5.1257E-03 5 0.070 1.40E-02 0.814543 3.95 0.3536 0.812496 2.437693 0.901386 1.561311 0.812496 3.42162E-05 1.624E-04 5.84946E-03 1.2743E-02 3.11538E-04 5.1304E-03 10 0.140 1.40E-02 0.814552 3.94 0.3537 0.812492 2.438001 0.901383 1.561410 0.812492 3.42112E-05 1.623E-04 5.84903E-03 1.2738E-02 6.08749E-04 5.1461E-03 15 0.211 1.40E-02 0.814563 3.93 0.3537 0.812487 2.438238 0.901380 1.561486 0.812487 3.42063E-05 1.622E-04 5.84862E-03 1.2734E-02 8.96859E-04 5.1740E-03 20 0.281 1.40E-02 0.814572 3.93 0.3537 0.812482 2.438376 0.901378 1.561530 0.812482 3.42014E-05 1.621E-04 5.84820E-03 1.2732E-02 1.17895E-03 5.2237E-03 25 0.351 1.40E-02 0.814578 3.93 0.3537 0.812478 2.438437 0.901376 1.561550 0.812478 3.41965E-05 1.621E-04 5.84778E-03 1.2731E-02 1.45687E-03 6.1669E-03 30 0.421 1.40E-02 0.814581 3.92 0.3537 0.812475 2.438452 0.901374 1.561554 0.812475 3.41916E-05 1.621E-04 5.84736E-03 1.2733E-02 1.73178E-03 7.2040E-03 35 0.491 1.40E-02 0.814581 3.92 0.3537 0.812471 2.438444 0.901372 1.561552 0.812471 3.41867E-05 1.622E-04 5.84694E-03 1.2736E-02 2.00443E-03 8.2623E-03 40 0.561 1.40E-02 0.814579 3.92 0.3537 0.812469 2.438429 0.901370 1.561547 0.812469 3.41818E-05 1.623E-04 5.84652E-03 1.2740E-02 2.27532E-03 9.3343E-03 45 0.631 1.40E-02 0.814576 3.92 0.3537 0.812466 2.438415 0.901369 1.561542 0.812466 3.41770E-05 1.625E-04 5.84611E-03 1.2746E-02 2.54478E-03 1.0415E-02 Simulation finished after 46 time-steps Writing final snapshot at time t = 0.64552808391132555 Wall clock time [hours] = 4.551E-04 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 4.347794 Maximum used memory per cpu [MBytes] = 26.430 Maximum used memory [MBytes] = 104.641 real 2.08 user 6.71 sys 0.25