Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified 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 2.2933979000000049E-002 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 ( 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 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 ( 162) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil uxb ( 195) 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) = 1822.8096178714193 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.10175662536339702 0.73807599698072957 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.284E-09 1.292E-09 1.275E-09 0.00053 2.723E-03 -6.152E-10 4.298E-08 8.277E-05 0.01945 0.21632 0.212759 0.033968 1.308E-10 2.5147E-10 1 0.006 5.56E-03 0.00123 5.038E-07 5.043E-07 5.040E-07 0.01941 5.359E-03 7.016E-10 9.201E-06 2.633E-04 0.01945 0.21633 0.212759 0.033968 7.275E-08 1.3409E-07 2 0.011 5.55E-03 0.00241 1.934E-06 1.935E-06 1.934E-06 0.03882 9.431E-03 1.853E-09 3.583E-05 7.763E-04 0.01945 0.21634 0.212759 0.033969 2.882E-07 5.2800E-07 3 0.017 5.58E-03 0.00336 3.754E-06 3.754E-06 3.754E-06 0.04977 1.309E-02 2.611E-09 6.749E-05 1.461E-03 0.01945 0.21634 0.212759 0.033973 4.644E-07 7.6874E-07 4 0.022 5.60E-03 0.00398 5.277E-06 5.277E-06 5.276E-06 0.04388 1.593E-02 2.826E-09 8.897E-05 2.109E-03 0.01945 0.21635 0.212759 0.033980 4.675E-07 6.1490E-07 5 0.028 5.46E-03 0.00440 6.449E-06 6.449E-06 6.448E-06 0.11341 1.835E-02 2.538E-09 1.069E-04 2.755E-03 0.01945 0.21635 0.212759 0.033987 8.017E-07 1.4453E-06 6 0.033 5.26E-03 0.00473 7.464E-06 7.464E-06 7.462E-06 0.20383 2.074E-02 2.136E-09 1.341E-04 3.561E-03 0.01945 0.21636 0.212759 0.033994 1.844E-06 3.9427E-06 7 0.039 5.14E-03 0.00472 7.437E-06 7.436E-06 7.432E-06 0.21076 2.278E-02 1.781E-09 1.408E-04 4.382E-03 0.01945 0.21636 0.212759 0.034002 2.207E-06 4.2943E-06 8 0.044 5.24E-03 0.00449 6.731E-06 6.730E-06 6.724E-06 0.07945 2.589E-02 7.625E-10 1.352E-04 5.963E-03 0.01945 0.21637 0.212759 0.034011 1.814E-06 4.1616E-06 9 0.049 5.00E-03 0.00540 9.708E-06 9.708E-06 9.699E-06 0.22048 3.153E-02 -4.957E-10 2.577E-04 9.531E-03 0.01945 0.21637 0.212759 0.034017 5.861E-06 1.8193E-05 Simulation finished after 10 time-steps Writing final snapshot at time t = 5.3976341561671880E-002 Wall clock time [hours] = 4.211E-02 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 57.82981 Wall clock time/timestep/local meshpoint [microsec] = 462.6385 Rhs wall clock time/timestep/local meshpoint [microsec] = 348.9684 Maximum used memory per cpu [MBytes] = 73.891 Maximum used memory [MBytes] = 588.871 real 152.51 user 1043.60 sys 2.02