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 rsnap: read snapshot var.dat in 1.0890960693359375E-003 seconds 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 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 initialize_gravity: shear with Omega proto r^-q, q= 0.0000000000000000 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! 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: performing full pencil check (takes a while) pencil_consistency_check: possible overcalculation... pencil rho ( 13) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil grho ( 16) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil ugrho ( 18) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil ugu ( 115) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil uuadvec_guu ( 140) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil del2a ( 165) 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.0000000000000000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000 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.0000000000000000 daa_dt: max(diffus_eta2) = 0.0000000000000000 daa_dt: max(diffus_eta3) = 0.0000000000000000 ---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] = 3.698E-04 (+/- 2.7778E-10) Wall clock time/timestep/meshpoint [microsec] = 30.95289 Wall clock time/timestep/local meshpoint [microsec] = 123.8116 Rhs wall clock time/timestep/local meshpoint [microsec] = 276.8426 Maximum used memory per cpu [MBytes] = 9.574 Maximum used memory [MBytes] = 37.582 real 9.05 user 10.70 sys 2.84