SVN: ------- v. ( ) $Id$ bcx1,bcx2= s a2 s a2 a2 s s s : s a2 s a2 a2 s s s bcy1,bcy2= p p p p p p p p : p p p p p p p p bcz1,bcz2= p p p p p p p p : p p p p p p p p lperi= F T T SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ bcx1,bcx2= s a2 s a2 a2 s s s : s a2 s a2 a2 s s s bcy1,bcy2= p p p p p p p p : p p p p p p p p bcz1,bcz2= p p p p p p p p : p p p p p p p p 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$ SVN: ------- v. ( ) $Id$ The verbose level is ip= 10 (ldebug= F ) This is a 2-D run nxgrid, nygrid, nzgrid= 32 32 1 Lx, Ly, Lz= 2.1000000000000001 6.2831853071795862 0.52000000000000002 Vbox= 19.132299260361840 rsnap: read snapshot var.dat in 0.0000000000000000E+000 seconds rprint_list: nname= 16 rprint_list: ix,iy,iz,iz2= -1 -1 -1 -1 rprint_list: nnamev= 2 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= -1 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 = p setup_slices: ix_loc,iy_loc,iz_loc, (video files) = 4 4 4 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= 1.2027220936797545E-008 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 -5.9914645471079817 5.0000000000000003E-002 2.5000000000000005E-003 1.0000000000000000 initialize_grid, r_int,r_ext= 0.40000000000000002 2.5000000000000000 timestep based on CFL cond; cdt= 0.40000000000000002 WARNING: initialize_diagnostics: volume averages are calculated wrongly for nonequidistant grids unless lproper_averages=T.! dVol_rel1= 6.9521690767519468E-004 dA_xy_rel1= 6.9521690767519468E-004 dA_xz_rel1= 3.1250000000000000E-002 dA_yz_rel1= 3.1250000000000000E-002 initialize_gravity: non-smoothed newtonian gravity 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 diffusion: Dhyper/pi^4 *(Delta(rho))^6/Deltaq^2 initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. 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 WARNING: initialize_energy: llocal_iso=T. Make sure you have the appropriate INITIAL_CONDITION in Makefile.local! select_eos_variable: Using rho and cs2 viscous force: nu_hyper3/pi^4 *(Deltav)^6/Deltaq^2 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 Lz= 0.52000000000000002 z= -0.26000000000000001 -0.26000000000000001 -0.26000000000000001 -0.26000000000000001 -0.26000000000000001 -0.26000000000000001 -0.26000000000000001 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 life_support_on: needed for pencil consistency check. 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 rho ( 7) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil grho ( 10) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil potself ( 175) 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 life_support_off: end of pencil consistency check. pencil_consistency_check: all tests passed pde: ENTER SVN: ------- v. ( ) $Id$ calc_pencils_energy: max(advec_cs2) = 4.8764099727013486 duu_dt: SOLVE Bcs for ux, x: < s>, y: < p>, z: < p> Bcs for uy, x: < a2>, y: < p>, z: < p> Bcs for uz, x: < s>, y: < p>, z: < p> duu_dt: max(advec_uu) = 19.754211595431695 dlnrho_dt: SOLVE Bcs for lnrho, x: < a2>, y: < p>, z: < p> dlnrho_dt: diffrho_hyper3= 1.0000000000000000E-003 dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000E+000 dlnrho_dt: max(diffus_diffrho3) = 1.8373871318870829E-002 ----it--------t---------dt-------dtc-----rhomin-rhomax---rhom------ux2m-------uy2m---potselfm---potself2m-gpotselfxm-gpotselfx2m--qtoomre-qtoomremin-qtoomremax- 0 0.000 1.82E-02 1.00E-01 0.2530 3.953 0.7141 2.002E-07 0.7313 0.000E+00 0.000E+00 3.147E-02 1.829E-03 1.25816 1.59E+00 1.59E+00 wvid: written slices in 9.9999993108212948E-004 seconds 5 0.091 1.82E-02 1.00E-01 0.2529 3.958 0.7141 2.009E-07 0.7313 0.000E+00 0.000E+00 3.147E-02 1.829E-03 1.25816 1.58E+00 1.60E+00 wvid: written slices in 0.0000000000000000E+000 seconds 10 0.182 1.81E-02 1.00E-01 0.2529 3.962 0.7141 2.110E-07 0.7313 0.000E+00 0.000E+00 3.147E-02 1.829E-03 1.25817 1.58E+00 1.61E+00 wvid: written slices in 0.0000000000000000E+000 seconds 15 0.272 1.82E-02 1.00E-01 0.2529 3.964 0.7141 2.258E-07 0.7312 0.000E+00 0.000E+00 3.147E-02 1.829E-03 1.25818 1.57E+00 1.61E+00 wvid: written slices in 0.0000000000000000E+000 seconds 20 0.363 1.82E-02 1.00E-01 0.2529 3.960 0.7141 2.422E-07 0.7312 0.000E+00 0.000E+00 3.147E-02 1.828E-03 1.25820 1.56E+00 1.62E+00 wvid: written slices in 0.0000000000000000E+000 seconds 25 0.454 1.82E-02 1.00E-01 0.2528 3.964 0.7140 2.590E-07 0.7312 0.000E+00 0.000E+00 3.146E-02 1.828E-03 1.25822 1.56E+00 1.63E+00 wvid: written slices in 1.0000001639127731E-003 seconds 30 0.545 1.82E-02 1.00E-01 0.2528 3.961 0.7140 2.755E-07 0.7312 0.000E+00 0.000E+00 3.146E-02 1.828E-03 1.25825 1.55E+00 1.63E+00 wvid: written slices in 0.0000000000000000E+000 seconds 35 0.636 1.82E-02 1.00E-01 0.2528 3.962 0.7140 2.917E-07 0.7312 0.000E+00 0.000E+00 3.146E-02 1.827E-03 1.25828 1.54E+00 1.64E+00 wvid: written slices in 0.0000000000000000E+000 seconds 40 0.727 1.82E-02 1.00E-01 0.2528 3.959 0.7139 3.072E-07 0.7312 0.000E+00 0.000E+00 3.146E-02 1.827E-03 1.25831 1.53E+00 1.65E+00 45 0.818 1.82E-02 1.00E-01 0.2528 3.958 0.7139 3.219E-07 0.7312 0.000E+00 0.000E+00 3.145E-02 1.826E-03 1.25835 1.53E+00 1.67E+00 wvid: written slices in 0.0000000000000000E+000 seconds Simulation finished after 46 time-steps Writing final snapshot at time t = 0.83583441206134002 Wall clock time [hours] = 2.000E-04 (+/- 2.7778E-07) Wall clock time/timestep/meshpoint [microsec] = 15.28533 Maximum used memory per cpu [MBytes] = 8.195 Maximum used memory [MBytes] = 8.195 real 3.99 user 3.82 sys 0.02