SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) no cvsid is given i SVN: ------- v. ( ) no cvsid is given i initialize_mpicomm: enabled MPI WARNING: run: quad precision suppressed! 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= 64 64 64 Lx, Ly, Lz= 2.2000000000000002 2.2000000000000002 2.2000000000000002 Vbox= 10.648000000000003 rsnap: read snapshot var.dat in 7.6318719999999840E-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= 4.8108883747190188E-009 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 2.7627751007914565 3.2500000000000000 6.3375000000000004 0.40000000000000008 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 initialize_grid: Set up half grid x12, y12, z12 initialize_gravity: smoothed and flattened 1/r potential mass diffusion: slope limited diffusion mass diffusion: using 2nd order WARNING: initialize_density: set lgamma_is_1=F if gamma/=1 with entropy! initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. initialize_hydro: outer radius rdampext= 1.0300000000000000 select_eos_variable: Using lnrho and ss initialize_energy: nheatc_max,iheatcond= 4 kramers chit nothing nothing heat conduction: kramers heat conduction: chit resistivity: constant eta viscous force: nu*(del2u+graddivu/3+2S.glnrho) Changed iss_run_aver bc to nil: Will time average the boundary of ss Changed iss_run_aver bc to nil: Will time average the boundary of ss Changed iss_run_aver bc to nil: Will time average the boundary of ss Changed iss_run_aver bc to nil: Will time average the boundary of ss Changed iss_run_aver bc to nil: Will time average the boundary of ss Changed iss_run_aver bc to nil: Will time average the boundary of ss 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 lnrho ( 12) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil rho ( 13) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil ee ( 56) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cp1 ( 61) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cv ( 73) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cv1 ( 74) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil r_mn1 ( 98) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil phix ( 99) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil phiy ( 100) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil pomx ( 101) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil pomy ( 102) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil rcyl_mn ( 103) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil rcyl_mn1 ( 104) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil evr ( 106) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil evth ( 108) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil aa ( 165) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil uxb ( 198) 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_hydro: upwinding advection term calc_pencils_hydro: max(advec_uu) = 0.0000000000000000 calc_pencils_energy: max(advec_cs2) = 1835.8157533964998 duu_dt: SOLVE Bcs for ux, x: < a>, y: < s>, z: < s> Bcs for uy, x: < s>, y: < a>, z: < s> Bcs for uz, x: < s>, y: < s>, z: < a> Bcs for sld_char, x: < cop>, y: < cop>, z: < cop> coriolis_cartesian: add Coriolis force; Omega= 1.0000000000000000 dlnrho_dt: SOLVE Bcs for lnrho, x: < a2>, y: < a2>, z: < a2> dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000 denergy_dt: SOLVE denergy_dt Bcs for ss, x: < sT>, y: < sT>, z: < sT> denergy_dt: lnTT,cs2,cp1= 0.11274230134979213 0.74622896165608921 1.0000000000000000 calc_heatcond_kramers: added thdiff calc_heatcond_chit: chi_t0= 0.0000000000000000 calc_heatcond_chit: chi_t1= 2.9999999999999997E-004 calc_heatcond_chit: added thdiff calc_heat_cool: lgravz, lgravr, lgravx, lspherical_coords= F T F F daa_dt: SOLVE Bcs for Ax, x: < a>, y: < s>, z: < s> Bcs for Ay, x: < s>, y: < a>, z: < s> Bcs for Az, x: < s>, y: < s>, z: < a> Bcs for sld_char, x: < cop>, y: < cop>, z: < cop> daa_dt: iresistivity=eta-const daa_dt: use upwinding in advection term daa_dt: max(diffus_eta) = 0.73803719008264457 daa_dt: max(diffus_eta2) = 0.0000000000000000 daa_dt: max(diffus_eta3) = 0.0000000000000000 ----it--------t--------dt-------urms------ux2m-------uy2m-------uz2m-------umax------orms--------oum------u2sphm-----o2sphm------rhom------ssm-------mass-------ethm-------ekin--------epsK---- 0 0.000 5.59E-03 0.00006 1.280E-09 1.289E-09 1.271E-09 0.00053 2.720E-03 -6.136E-10 4.287E-08 8.256E-05 0.01943 0.22182 0.212586 0.034054 1.306E-10 2.5109E-10 1 0.006 5.56E-03 0.00032 3.493E-08 3.520E-08 3.498E-08 0.01690 4.293E-03 3.211E-10 1.172E-06 2.057E-04 0.01943 0.22181 0.212586 0.034053 3.777E-08 6.3468E-08 2 0.011 5.54E-03 0.00064 1.375E-07 1.379E-07 1.375E-07 0.03396 7.006E-03 1.098E-09 4.606E-06 5.477E-04 0.01943 0.22180 0.212586 0.034054 1.521E-07 2.5681E-07 3 0.017 5.58E-03 0.00080 2.108E-07 2.110E-07 2.106E-07 0.04328 9.640E-03 1.577E-09 7.051E-06 1.037E-03 0.01943 0.22179 0.212586 0.034056 2.227E-07 3.1729E-07 4 0.022 5.60E-03 0.00076 1.935E-07 1.934E-07 1.931E-07 0.03739 1.186E-02 1.672E-09 6.459E-06 1.569E-03 0.01943 0.22178 0.212586 0.034060 1.787E-07 1.2927E-07 5 0.028 5.47E-03 0.00113 4.253E-07 4.251E-07 4.248E-07 0.10831 1.369E-02 1.474E-09 1.422E-05 2.092E-03 0.01943 0.22177 0.212586 0.034065 4.532E-07 8.1118E-07 6 0.033 5.27E-03 0.00184 1.124E-06 1.123E-06 1.122E-06 0.19524 1.537E-02 1.496E-09 3.758E-05 2.636E-03 0.01943 0.22176 0.212586 0.034069 1.332E-06 2.8488E-06 7 0.039 5.16E-03 0.00202 1.363E-06 1.361E-06 1.361E-06 0.20147 1.712E-02 2.070E-09 4.558E-05 3.269E-03 0.01943 0.22175 0.212586 0.034074 1.592E-06 3.0315E-06 8 0.044 5.26E-03 0.00192 1.227E-06 1.225E-06 1.226E-06 0.07166 1.993E-02 2.330E-09 4.103E-05 4.431E-03 0.01943 0.22174 0.212586 0.034081 1.131E-06 2.6179E-06 9 0.049 5.01E-03 0.00353 4.142E-06 4.141E-06 4.143E-06 0.22044 2.325E-02 1.999E-09 1.387E-04 6.033E-03 0.01943 0.22173 0.212586 0.034085 4.512E-06 1.3998E-05 Simulation finished after 10 time-steps Writing final snapshot at time t = 5.4035487148736398E-002 Wall clock time [hours] = 5.684E-02 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 78.06281 Wall clock time/timestep/local meshpoint [microsec] = 624.5025 Rhs wall clock time/timestep/local meshpoint [microsec] = 551.5570 Maximum used memory per cpu [MBytes] = 69.691 Maximum used memory [MBytes] = 552.566 real 205.10 user 785.58 sys 29.66