SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ bcx1,bcx2= a s s a2 a2 a s s s : a s s a2 a2 a s s s bcy1,bcy2= s a s a2 a2 s a s s : s a s a2 a2 s a s s bcz1,bcz2= s s a a2 a2 s s a s : s s a a2 a2 s s a s 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= 5.0000000000000000 5.0000000000000000 5.0000000000000000 Vbox= 125.00000000000000 rprint_list: nname= 18 rprint_list: ix,iy,iz,iz2= -1 -1 -1 -1 rprint_list: nnamev= 0 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_shock: run through parse list rprint_viscosity: run through parse list 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 0.40546510810816422 1.0000000000000000 0.59999999999999998 0.40000000000000008 timestep based on CFL cond; cdt= 0.40000000000000002 initialize_gravity: smoothed 1/r potential xmask_den= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 zmask_den= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 diffusion: nothing (i.e. no mass diffusion) initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. initialize_hydro: outer radius rdampext= 3.0000000000000000 xmask_hyd= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 zmask_hyd= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 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.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 zmask_mag= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 resistivity: shell initialize_shock: prenormalised shock_factor sum= 1.0000000000000000 viscous force: nu*(del2u+graddivu/3+2S.glnrho) viscous force: nu_shock*(XXXXXXXXXXX) xmask_vis= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 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) = 97.165773730143172 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> duu_dt: max(advec_uu) = 0.0000000000000000E+000 dlnrho_dt: SOLVE Bcs for lnrho, x: < a2>, y: < a2>, z: < a2> dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000E+000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000E+000 denergy_dt: SOLVE denergy_dt Bcs for ss, x: < a2>, y: < a2>, z: < a2> denergy_dt: lnTT,cs2,cp1= 0.23417276267693740 0.84257521444799821 1.0000000000000000 calc_heatcond_constchi: chi= 2.0000000000000000E-003 calc_heatcond_constchi: added thdiff 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=shell daa_dt: max(diffus_eta) = 0.23064000000000001 daa_dt: max(diffus_eta2) = 0.0000000000000000E+000 daa_dt: max(diffus_eta3) = 0.0000000000000000E+000 pde: freezing variables for r < 1.0000000000000000 : T T T T T T T T F ----it-------------t------------------dt----------------urms---------------umax----------------brms-----------------bmax-----------------epsK-----------------epsM---------------rhom----------------ssm---------------Marms--------------Mamax---------------vArms----------------vAmax---------------beta1m-------------beta1max------------shockmax------ 0 0.00000000000000 3.4497939497199E-02 0.00000000000000 0.00000000000000 1.50660298999331E-01 8.20930484218245E-01 0.00000000000000E+00 1.4294321038045E-03 0.42939736357543 0.26020575672839 0.00000000000000 0.00000000000000 1.73334582581480E-01 8.87864862046937E-01 2.8165751012175E-02 7.2928923995987E-01 0.0000000000000E+00 10 0.29855533670282 2.8257550665232E-02 0.04764564133289 0.40577670415882 1.39183892100697E-01 6.22213189947349E-01 3.15306739590431E-04 6.4535215365133E-04 0.42884797589611 0.26046415600019 0.05101782770303 0.44907090468680 1.65282712316009E-01 7.93195199643802E-01 2.6431654430068E-02 7.0646434777711E-01 1.5673505902845E-02 20 0.57345337016128 2.7283730342637E-02 0.06427752670822 0.47724661391863 1.32621148662506E-01 6.52753636881266E-01 4.00025775709997E-04 5.0067504020646E-04 0.42912870431812 0.26040376021498 0.06955837137219 0.51587043234429 1.57758867539400E-01 7.66829511404838E-01 2.3786291625314E-02 5.8696870057860E-01 2.0640067151489E-02 30 0.85146537533974 2.8531335315585E-02 0.07009142680047 0.40185082925566 1.27698449121874E-01 6.27389573392631E-01 3.47822177363789E-04 4.9360252277450E-04 0.43009770049199 0.26019708656788 0.07581379112222 0.43668790624682 1.52742608432759E-01 7.05307353832751E-01 2.2262302711566E-02 4.8627595728718E-01 1.2286599585020E-02 40 1.14120661364245 2.9435740297644E-02 0.07363472282460 0.35922877242974 1.24145150929859E-01 6.32719199940288E-01 3.25982651050151E-04 4.7744315714213E-04 0.43111509841149 0.25982664053757 0.07971595832275 0.39665690490745 1.49356058324093E-01 6.80523449336109E-01 2.1273237144203E-02 4.3386353791040E-01 1.4291159956565E-02 Simulation finished after 50 time-steps Writing final snapshot at time t = 1.4375952214241661 Wall clock time [hours] = 9.924E-03 (+/- 2.7778E-07) Wall clock time/timestep/meshpoint [microsec] = 21.80603 Maximum used memory per cpu [MBytes] = 25.434 Maximum used memory [MBytes] = 25.434 real 40.81 user 39.22 sys 0.04