Pencil1: iapp, nprocs, ncpus= 0 4 4 SVN: ------- v. ( ) $Id$ bcx1,bcx2= a a2 s a2 s : a a2 s a2 s bcy1,bcy2= p p p p p : p p p p p bcz1,bcz2= p p p p p : p p p p p lperi= F T T SVN: start.in,v v. 1.1 (2017/04/26 17:31:11) wlyra SVN: run.in,v v. 1.1 (2017/04/26 17:31:11) wlyra bcx1,bcx2= a a2 s a2 s : a a2 s a2 s bcy1,bcy2= p p p p p : p p p p p bcz1,bcz2= p p p p p : p p p p p 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= 10 (ldebug= F ) This is a 2-D run nxgrid, nygrid, nzgrid= 16 48 1 Lx, Ly, Lz= 2.1000000000000001 6.2831853071795862 2.0000000000000000 Vbox= 19.132299260361840 rsnap: read snapshot var.dat in 6.6572790000000003E-003 seconds rprint_list: nname= 7 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= -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_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= 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 dVol_rel1= 9.5785440613026837E-004 dA_xy_rel1= 9.5785440613026837E-004 dA_xz_rel1= 6.2500000000000000E-002 dA_yz_rel1= 2.0833333333333332E-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 zmask_den= 1.0000000000000000 diffusion: mesh hyperdiffusion 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 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_mesh/pi^5 *(Deltav)^6/Deltaq viscous force: nu_shock*(XXXXXXXXXXX) xmask_vis= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 Lz= 2.0000000000000000 z= 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 Lz= 2.0000000000000000 Lz= 2.0000000000000000 z= 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 z= 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 Lz= 2.0000000000000000 z= 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 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_hydro: call gij_etc calc_pencils_energy: max(advec_cs2) = 2.5986041829730091 duu_dt: SOLVE Bcs for ux, x: < a>, 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) = 8.4814791505693786E-015 dlnrho_dt: SOLVE Bcs for lnrho, x: < a2>, y: < p>, z: < p> dlnrho_dt: diffrho_hyper3_mesh= 20.000000000000000 dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000E+000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000E+000 ---it-------t----------dt----------umax-----rhomax-rhomin--nshift-- 0 0.000 1.340E-01 1.868E+00 1.0 1.0 3.017 5 0.670 1.340E-01 1.868E+00 1.0 1.0 3.017 10 1.340 1.339E-01 1.868E+00 1.1 1.0 3.016 15 2.009 1.339E-01 1.868E+00 1.2 0.9 3.015 20 2.678 1.338E-01 1.868E+00 1.2 0.9 3.013 25 3.347 1.338E-01 1.869E+00 1.2 0.9 3.013 30 4.016 1.337E-01 1.868E+00 1.2 0.9 3.012 35 4.685 1.337E-01 1.868E+00 1.3 0.8 3.011 40 5.353 1.337E-01 1.868E+00 1.3 0.8 3.010 45 6.021 1.337E-01 1.868E+00 1.2 0.8 3.010 Simulation finished after 50 time-steps Writing final snapshot at time t = 6.6897194356860243 Wall clock time [hours] = 6.814E-05 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 6.388097 Maximum used memory per cpu [MBytes] = 24.105 Maximum used memory [MBytes] = 95.422 real 0.62 user 1.22 sys 0.09