Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified 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 4.2785699999999954E-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) = 1835.8157533965004 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.11274230134979168 0.74622896165608876 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.290E-09 1.271E-09 1.275E-09 0.00046 2.718E-03 6.958E-10 0.22182 0.22182 10 0.056 5.55E-03 0.00085 2.431E-07 2.430E-07 2.429E-07 0.04339 2.872E-02 -5.869E-09 0.22172 0.22181 20 0.111 5.47E-03 0.00141 6.592E-07 6.587E-07 6.574E-07 0.07583 5.099E-02 -9.994E-09 0.22162 0.22181 30 0.165 5.47E-03 0.00180 1.076E-06 1.076E-06 1.073E-06 0.10035 6.350E-02 -1.237E-08 0.22153 0.22181 40 0.220 5.44E-03 0.00202 1.358E-06 1.358E-06 1.349E-06 0.09908 6.838E-02 -1.030E-08 0.22144 0.22180 50 0.274 5.43E-03 0.00212 1.506E-06 1.506E-06 1.491E-06 0.10142 6.727E-02 -1.040E-08 0.22136 0.22180 60 0.328 5.41E-03 0.00223 1.672E-06 1.671E-06 1.643E-06 0.11917 6.471E-02 -8.359E-09 0.22127 0.22179 70 0.383 5.43E-03 0.00246 2.023E-06 2.023E-06 1.985E-06 0.10900 6.572E-02 -6.172E-09 0.22119 0.22178 80 0.437 5.48E-03 0.00264 2.335E-06 2.336E-06 2.298E-06 0.08446 7.037E-02 -8.227E-09 0.22111 0.22177 90 0.492 5.46E-03 0.00295 2.927E-06 2.927E-06 2.877E-06 0.10542 7.591E-02 -7.613E-09 0.22104 0.22176 Simulation finished after 100 time-steps Writing final snapshot at time t = 0.54663755173572848 Wall clock time [hours] = 4.121E-02 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 5.659636 Wall clock time/timestep/local meshpoint [microsec] = 5.659636 Rhs wall clock time/timestep/local meshpoint [microsec] = 5.492318 Maximum used memory per cpu [MBytes] = 107.543 Maximum used memory [MBytes] = 107.543 real 148.70 user 146.65 sys 0.75