Pencil1: iapp, nprocs, ncpus= 0 4 4 ===================================================================== WARNING: using SINGLE PRECISION, you'd better know what you're doing! ===================================================================== SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ bcx1,bcx2= a s s a2 a2 s s a : a s s a2 a2 s s a bcy1,bcy2= s a s a2 a2 s s a : s a s a2 a2 s s a bcz1,bcz2= s s a s s s s s : s s a a2 a2 s s a initialize_mpicomm: enabled MPI 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= 13 (ldebug= F ) This is a 3-D run nxgrid, nygrid, nzgrid= 64 64 32 Lx, Ly, Lz= 2.20000005 2.20000005 2.20000005 Vbox= 10.6480007 rprint_list: nname= 18 rprint_list: ix,iy,iz,iz2= -1 -1 -1 -1 rprint_list: nnamev= 5 rprint_list: nnamez= -1 rprint_list: nnamey= -1 rprint_list: nnamex= -1 rprint_list: nnamer= 0 rprint_list: nnamexz= -1 rprint_list: nnamexy= -1 rprint_list: nnamerz= 7 rprint_general: run through parse list rprint_hydro: run through parse list rprint_hydro: run through parse list rprint_density: run through parse list rprint_viscosity: run through parse list setup_slices: slice_position = e setup_slices: ix_loc,iz_loc,iz2_loc, (video files) = 35 4 19 units_general: unit_velocity= 1.0000000000000000 units_general: unit_density= 1.0000000000000000 units_general: unit_length= 1.0000000000000000 units_general: unit_magnetic= 3.5449078083038330 units_eos: unit_temperature= 4.8108880235694808E-009 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.00000000 1.19209221E-06 0.816497028 0.400000453 0.399999976 timestep based on CFL cond; cdt= 0.400000006 initialize_gravity: smoothed 1/r potential in spherical shell xmask_den= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 zmask_den= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 diffusion: nothing (i.e. no mass diffusion) initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. initialize_hydro: set outer radius rdampext=r_ext xmask_hyd= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 zmask_hyd= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 WARNING: initialize_hydro: need to set orms in print.in to get othresh. Disabled othresh! select_eos_variable: Using lnrho and ss initialize_energy: set boundary temperatures for spherical shell problem initialize_energy: g0,mpoly,beta1 1.20000005 1.00000000 1.50000012 initialize_energy: TT_int, TT_ext= 2.50000000 1.00000000 initialize_energy: cs2_ext, cs2_int= 0.666666627 1.66666651 initialize_energy: nheatc_max,iheatcond= 4 K-const nothing nothing nothing heat conduction: K=cte xmask_mag= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 zmask_mag= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 resistivity: shell viscous force: nu*(del2u+graddivu/3+2S.glnrho) xmask_vis= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 choose_pencils: finding out which pencils are needed for the pencil case. write_pencil_info: pencil information written to the file pencils.list. 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: 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: 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_energy: max(advec_cs2) = 1225.75745 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= 0.632456005 duu_dt: max(advec_uu) = 0.999772191 dlnrho_dt: SOLVE Bcs for lnrho, x: < a2>, y: < a2>, z: < s:a2> dlnrho_dt: max(diffus_diffrho ) = 0.00000000E+00 dlnrho_dt: max(diffus_diffrho3) = 0.00000000E+00 denergy_dt: SOLVE denergy_dt Bcs for ss, x: < a2>, y: < a2>, z: < s:a2> denergy_dt: lnTT,cs2,cp1= 9.47388733E-15 0.666666627 1.00000000 calc_heatcond_constK: hcond= 0.00000000E+00 calc_heatcond_constK: added thdiff calc_heat_cool: lgravz, lgravr, lgravx, lspherical_coords= F T F F daa_dt: SOLVE Bcs for Ax, x: < s>, y: < s>, z: < s> Bcs for Ay, x: < s>, y: < s>, z: < s> Bcs for Az, x: < a>, y: < a>, z: < s:a> daa_dt: iresistivity=shell daa_dt: max(diffus_eta) = 0.813999176 daa_dt: max(diffus_eta2) = 0.00000000E+00 daa_dt: max(diffus_eta3) = 0.00000000E+00 ----it--------t--------dt-------urms------umax------brms------bmax------epsK------epsM-----rhom------ssm------ethm-----dtc-----dtu-----dtb-----dtnu---dtchi---dteta-- 0 0.000 7.08E-03 0.01547 0.02717 4.064E-05 1.331E-04 4.402E-05 5.488E-09 1.2077 0.02953 0.9739 0.9802 0.0198 0.0001 0.0027 0.0037 0.0192 2 0.014 7.09E-03 0.01172 0.04736 3.978E-05 1.329E-04 1.981E-05 5.061E-09 1.2076 0.02954 0.9738 0.9851 0.0355 0.0001 0.0027 0.0037 0.0192 4 0.028 7.09E-03 0.01257 0.06671 3.897E-05 1.327E-04 1.563E-05 4.672E-09 1.2075 0.02958 0.9736 0.9865 0.0450 0.0001 0.0027 0.0038 0.0192 6 0.043 7.10E-03 0.01394 0.07307 3.819E-05 1.325E-04 1.715E-05 4.317E-09 1.2073 0.02966 0.9734 0.9863 0.0487 0.0000 0.0028 0.0038 0.0193 8 0.057 7.11E-03 0.01505 0.07590 3.746E-05 1.323E-04 1.844E-05 3.992E-09 1.2072 0.02977 0.9732 0.9890 0.0568 0.0000 0.0028 0.0038 0.0193 10 0.071 7.12E-03 0.01573 0.07496 3.676E-05 1.320E-04 1.738E-05 3.696E-09 1.2070 0.02990 0.9729 0.9905 0.0561 0.0000 0.0028 0.0038 0.0193 12 0.085 7.13E-03 0.01616 0.07487 3.610E-05 1.317E-04 1.615E-05 3.424E-09 1.2069 0.03005 0.9727 0.9917 0.0556 0.0000 0.0028 0.0038 0.0194 14 0.099 7.14E-03 0.01655 0.07907 3.547E-05 1.314E-04 1.585E-05 3.176E-09 1.2067 0.03021 0.9725 0.9930 0.0583 0.0000 0.0028 0.0039 0.0194 16 0.114 7.15E-03 0.01693 0.07939 3.487E-05 1.311E-04 1.588E-05 2.949E-09 1.2066 0.03037 0.9723 0.9942 0.0590 0.0000 0.0028 0.0039 0.0194 Simulation finished after 17 time-steps Writing final snapshot at time t = 0.12093319022096694 Wall clock time [hours] = 6.555E-03 (+/- 2.7778E-10) Wall clock time/timestep/meshpoint [microsec] = 10.59102 Maximum used memory per cpu [MBytes] = 22.832 Maximum used memory [MBytes] = 90.090 real 35.11 user 53.38 sys 0.95