Pencil1: iapp, nprocs, ncpus= 0 4 4 SVN: ------- v. ( ) $Id$ SVN: start.in,v v. 1.1 (2008-11-30 03:16:25) wlyra SVN: run.in,v v. 1.1 (2008-11-30 03:16:24) 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= 16 (ldebug= F ) This is a 2-D run nxgrid, nygrid, nzgrid= 32 64 1 Lx, Ly, Lz= 1.0000000000000000 1.0000000000000000 0.20000000000000001 Vbox= 1.5000000000000000 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 -9.2103403719761818 1.0000000000000000E-002 1.0000000000000009E-004 1.0000000000000000 initialize_gravity: shear with Omega proto r^-q, q= 0.0000000000000000E+000 diffusion: nothing (i.e. no mass diffusion) initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. select_eos_variable: Using rho and cs2 resistivity: constant eta WARNING: initialize_magnetic: Resistivity coefficient eta is zero! 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: performing full pencil check (takes a while) pencil_consistency_check: possible overcalculation... pencil rho ( 7) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil grho ( 10) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil ugrho ( 12) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil ugu ( 109) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil uuadvec_guu ( 134) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil del2a ( 159) is requested, but does not appear to be required! 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: performing full pencil check (takes a while) 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) = 0.50570000000000004 duu_dt: SOLVE Bcs for ux, x: < f>, y: < p>, z: < p> Bcs for uy, x: < f>, y: < p>, z: < p> Bcs for uz, x: < f>, y: < p>, z: < p> duu_dt: max(advec_uu) = 5.3398363172274199E-014 dlnrho_dt: SOLVE Bcs for lnrho, x: < f>, y: < p>, z: < p> dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000E+000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000E+000 daa_dt: SOLVE Bcs for Ax, x: < f>, y: < p>, z: < p> Bcs for Ay, x: < f>, y: < p>, z: < p> Bcs for Az, x: < f>, y: < p>, z: < p> daa_dt: iresistivity=eta-const daa_dt: max(diffus_eta) = 0.0000000000000000E+000 daa_dt: max(diffus_eta2) = 0.0000000000000000E+000 daa_dt: max(diffus_eta3) = 0.0000000000000000E+000 ---it-------t----------dt----------umax--------rhomin--------rhom--------rhomax-------nshift--------amax----------a2m----------azm---------arms---------bmax---------bx2m---------by2m----------b2m---------brms----- 0 0.0000 5.625E-01 2.000E+00 1.000E+00 1.000E+00 1.000E+00 3.600E+01 2.801E-04 2.827E-09 1.885E-05 5.317E-05 1.102E-03 9.202E-08 9.102E-08 1.830E-07 4.278E-04 5 2.8124 5.625E-01 2.000E+00 9.990E-01 1.000E+00 1.001E+00 3.600E+01 2.801E-04 2.827E-09 1.885E-05 5.317E-05 1.102E-03 9.202E-08 9.102E-08 1.830E-07 4.278E-04 10 5.6249 5.625E-01 2.000E+00 9.981E-01 1.000E+00 1.002E+00 3.600E+01 2.802E-04 2.827E-09 1.885E-05 5.317E-05 1.103E-03 9.202E-08 9.102E-08 1.830E-07 4.278E-04 15 8.4373 5.625E-01 2.000E+00 9.971E-01 1.000E+00 1.003E+00 3.600E+01 2.802E-04 2.827E-09 1.885E-05 5.317E-05 1.103E-03 9.202E-08 9.101E-08 1.830E-07 4.278E-04 20 11.2498 5.625E-01 2.000E+00 9.961E-01 1.000E+00 1.004E+00 3.600E+01 2.802E-04 2.827E-09 1.885E-05 5.317E-05 1.103E-03 9.201E-08 9.101E-08 1.830E-07 4.278E-04 Simulation finished after 21 time-steps Writing final snapshot at time t = 11.812254875183900 Wall clock time [hours] = 4.609E-04 (+/- 2.7778E-10) Wall clock time/timestep/meshpoint [microsec] = 38.57615 Maximum used memory per cpu [MBytes] = 9.723 Maximum used memory [MBytes] = 37.395 real 10.43 user 10.34 sys 2.98