SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ bcx1,bcx2= a s s a2 sT a s s : a s s a2 sT a s s bcy1,bcy2= s a s a2 sT s a s : s a s a2 sT s a s bcz1,bcz2= s s a a2 sT s s a : s s a a2 sT s s a 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= 32 32 32 Lx, Ly, Lz= 3.00000000 3.00000000 3.00000000 Vbox= 27.0000000 rprint_list: nname= 13 rprint_list: ix,iy,iz,iz2= -1 -1 11 -1 rprint_list: nnamev= 2 rprint_list: nnamez= 4 rprint_list: nnamey= 2 rprint_list: nnamex= 2 rprint_list: nnamer= 2 rprint_list: nnamexz= 2 rprint_list: nnamexy= 2 rprint_list: nnamerz= 10 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 = m setup_slices: ix_loc,iy_loc,iz_loc,iz2_loc, (video files) = 19 19 14 35 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.21639538 1.50000000 1.35000002 0.399999976 timestep based on CFL cond; cdt= 0.400000006 TIMEAVG: defaulting idx no. 1 to 1 TIMEAVG: defaulting idx no. 2 to 2 TIMEAVG: defaulting idx no. 3 to 3 TIMEAVG: defaulting idx no. 4 to 4 TIMEAVG: defaulting idx no. 5 to 5 TIMEAVG: defaulting idx no. 6 to 6 TIMEAVG: defaulting idx no. 7 to 7 TIMEAVG: defaulting idx no. 8 to 8 initialize_gravity: M5 dwarf gravity potential 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 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 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: outer radius rdampext= 1.10000002 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 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 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 select_eos_variable: Using lnrho and ss initialize_energy: nheatc_max,iheatcond= 4 chi-const nothing nothing nothing heat conduction: constant chi 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 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 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: constant eta 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 choose_pencils: finding out which pencils are needed for the pencil case. write_pencil_info: pencil information written to the file pencils.list. pde: ENTER SVN: ------- v. ( ) $Id$ calc_pencils_energy: max(advec_cs2) = 433.583099 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.500000000 duu_dt: max(advec_uu) = 0.00000000E+00 dlnrho_dt: SOLVE Bcs for lnrho, x: < a2>, y: < a2>, z: < 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: < sT>, y: < sT>, z: < sT> denergy_dt: lnTT,cs2,cp1= 0.708186328 1.35353696 1.00000000 calc_heatcond_constchi: chi= 1.00000005E-03 calc_heatcond_constchi: 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: max(diffus_eta) = 0.192200035 daa_dt: max(diffus_eta2) = 0.00000000E+00 daa_dt: max(diffus_eta3) = 0.00000000E+00 ----it--------t--------dt-------urms------umax------brms-------bmax-------rhom------ssm------dtc-----dtu-----dtnu---dtchi-- 0 0.000 1.41E-02 0.01140 0.1249 5.999E-06 6.287E-05 0.34611 0.15062 0.9412 0.0779 0.0150 0.0250 10 0.143 1.43E-02 0.00916 0.1024 5.560E-06 5.650E-05 0.34610 0.15063 0.9600 0.0547 0.0152 0.0254 20 0.286 1.44E-02 0.00810 0.0890 5.176E-06 5.133E-05 0.34610 0.15065 0.9659 0.0500 0.0153 0.0256 30 0.430 1.44E-02 0.00726 0.0730 4.832E-06 4.915E-05 0.34610 0.15067 0.9738 0.0413 0.0154 0.0257 40 0.574 1.44E-02 0.00650 0.0682 4.520E-06 4.710E-05 0.34610 0.15068 0.9757 0.0391 0.0154 0.0257 Simulation finished after 41 time-steps Writing final snapshot at time t = 0.58838926372118294 Wall clock time [hours] = 6.283E-03 (+/- 2.7778E-07) Wall clock time/timestep/meshpoint [microsec] = 16.83528 Maximum used memory per cpu [MBytes] = 11.500 Maximum used memory [MBytes] = 11.500 real 23.01 user 21.54 sys 0.05