Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified SVN: ------- v. ( ) $Id$ WARNING: find_namelist: namelist "_pars" is missing!! SVN: ------- v. ( ) no cvsid is given i WARNING: find_namelist: namelist "_pars" is missing!! 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 4.3638770000000021E-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_gravity: smoothed and flattened 1/r potential diffusion: nothing (i.e. no mass diffusion) 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) 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) = 1822.6663507387650 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> 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.10150335253343412 0.73788908605497561 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> 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-------ssm----ss_run_averm- 0 0.000 5.59E-03 0.00006 1.294E-09 1.275E-09 1.278E-09 0.00046 2.721E-03 6.976E-10 0.21632 0.21632 10 0.056 5.55E-03 0.00246 2.013E-06 2.013E-06 2.002E-06 0.04388 3.125E-02 -6.115E-09 0.21638 0.21633 20 0.111 5.47E-03 0.00269 2.423E-06 2.423E-06 2.402E-06 0.07535 5.145E-02 -1.039E-08 0.21644 0.21633 30 0.165 5.48E-03 0.00287 2.745E-06 2.744E-06 2.734E-06 0.10139 6.401E-02 -1.358E-08 0.21650 0.21633 40 0.220 5.46E-03 0.00300 3.001E-06 3.000E-06 2.981E-06 0.10067 6.866E-02 -1.170E-08 0.21655 0.21634 50 0.274 5.42E-03 0.00311 3.232E-06 3.232E-06 3.213E-06 0.10239 6.782E-02 -1.263E-08 0.21661 0.21634 60 0.328 5.41E-03 0.00325 3.543E-06 3.542E-06 3.504E-06 0.11956 6.541E-02 -9.247E-09 0.21666 0.21634 70 0.383 5.46E-03 0.00353 4.165E-06 4.166E-06 4.121E-06 0.09438 6.641E-02 -9.463E-09 0.21671 0.21635 80 0.437 5.46E-03 0.00371 4.607E-06 4.609E-06 4.562E-06 0.09604 7.028E-02 -8.999E-09 0.21676 0.21636 90 0.492 5.47E-03 0.00406 5.501E-06 5.502E-06 5.446E-06 0.09814 7.635E-02 -9.471E-09 0.21681 0.21636 Simulation finished after 100 time-steps Writing final snapshot at time t = 0.54652932930709552 Wall clock time [hours] = 4.080E-02 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 5.603175 Wall clock time/timestep/local meshpoint [microsec] = 5.603175 Rhs wall clock time/timestep/local meshpoint [microsec] = 5.444678 Maximum used memory per cpu [MBytes] = 107.785 Maximum used memory [MBytes] = 107.785 real 147.23 user 146.05 sys 0.76