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 6.7268239999999924E-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 1.4667216102325049 1.7000000000000000 1.7340000000000000 0.40000000000000008 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 initialize_gravity: RGB-1Msun-11.7 gravity 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 K-profile nothing nothing nothing heat conduction: K-profile resistivity: constant eta viscous force: nu*(del2u+graddivu/3+2S.glnrho) 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 cool_prof ( 52) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil ee ( 55) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cp1 ( 60) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil TT ( 63) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cv ( 72) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cv1 ( 73) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil r_mn1 ( 97) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil phix ( 98) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil phiy ( 99) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil pomx ( 100) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil pomy ( 101) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil rcyl_mn ( 102) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil rcyl_mn1 ( 103) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil evr ( 105) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil evth ( 107) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil aa ( 163) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil uxb ( 196) 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) = 2342.0964849193329 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.35629982835442653 0.95202376644071307 1.0000000000000000 calc_heatcond_arrays: hcond0= 0.0000000000000000 calc_heatcond_arrays: lgravz= F calc_heatcond: lhcond_global= T calc_heatcond: 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) = 4.9202479338842977 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----inertiaxx_car-inertiayy_car-inertiazz_car- 0 0.000 1.06E-02 0.00006197 1.280E-09 1.289E-09 1.271E-09 0.00053 2.720E-03 -6.136E-10 4.287E-08 8.256E-05 0.12894 0.24142 1.398840 0.13446550 5.013E-10 6.4236E-09 1.36075E-01 1.36075E-01 1.36075E-01 1 0.011 1.06E-02 0.00007728 2.002E-09 1.999E-09 1.971E-09 0.00148 3.343E-03 -2.403E-09 6.588E-08 1.247E-04 0.12894 0.24129 1.398840 0.13445010 8.420E-10 9.4447E-09 1.36075E-01 1.36075E-01 1.36075E-01 2 0.021 1.05E-02 0.00009854 3.266E-09 3.237E-09 3.208E-09 0.00212 4.289E-03 -3.514E-09 1.034E-07 2.052E-04 0.12894 0.24117 1.398840 0.13443517 1.460E-09 1.3816E-08 1.36075E-01 1.36075E-01 1.36075E-01 3 0.032 1.05E-02 0.00011280 4.274E-09 4.232E-09 4.217E-09 0.00214 4.840E-03 -3.822E-09 1.257E-07 2.613E-04 0.12894 0.24105 1.398840 0.13442082 1.922E-09 1.5733E-08 1.36075E-01 1.36075E-01 1.36075E-01 4 0.042 1.05E-02 0.00012631 5.343E-09 5.307E-09 5.303E-09 0.00193 5.051E-03 -3.598E-09 1.399E-07 2.846E-04 0.12894 0.24094 1.398840 0.13440709 2.277E-09 1.7042E-08 1.36074E-01 1.36074E-01 1.36074E-01 5 0.053 1.05E-02 0.00014131 6.675E-09 6.651E-09 6.644E-09 0.00178 5.035E-03 -3.120E-09 1.512E-07 2.826E-04 0.12894 0.24084 1.398840 0.13439401 2.586E-09 1.7975E-08 1.36073E-01 1.36073E-01 1.36073E-01 6 0.063 1.05E-02 0.00015711 8.243E-09 8.228E-09 8.211E-09 0.00151 4.879E-03 -2.583E-09 1.593E-07 2.652E-04 0.12894 0.24074 1.398840 0.13438157 2.850E-09 1.7831E-08 1.36072E-01 1.36072E-01 1.36072E-01 7 0.074 1.05E-02 0.00017374 1.008E-08 1.006E-08 1.004E-08 0.00162 4.645E-03 -2.093E-09 1.671E-07 2.402E-04 0.12894 0.24064 1.398840 0.13436975 3.097E-09 1.6835E-08 1.36070E-01 1.36070E-01 1.36070E-01 8 0.084 1.05E-02 0.00019187 1.229E-08 1.227E-08 1.225E-08 0.00180 4.375E-03 -1.696E-09 1.802E-07 2.128E-04 0.12894 0.24056 1.398840 0.13435853 3.381E-09 1.5614E-08 1.36068E-01 1.36068E-01 1.36068E-01 9 0.095 1.05E-02 0.00021181 1.498E-08 1.495E-08 1.494E-08 0.00197 4.096E-03 -1.395E-09 2.034E-07 1.863E-04 0.12894 0.24047 1.398840 0.13434788 3.752E-09 1.4511E-08 1.36066E-01 1.36066E-01 1.36066E-01 Simulation finished after 10 time-steps Writing final snapshot at time t = 0.10545807208801390 Wall clock time [hours] = 5.246E-02 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 72.04129 Wall clock time/timestep/local meshpoint [microsec] = 576.3303 Rhs wall clock time/timestep/local meshpoint [microsec] = 513.0668 Maximum used memory per cpu [MBytes] = 72.656 Maximum used memory [MBytes] = 565.676 real 189.28 user 716.94 sys 35.54