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 1.0231899999999999E-004 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: 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
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
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
dL_y_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
choose_pencils: finding out which pencils are needed for the pencil case.
write_pencil_info: pencil information written to the file pencils.list.
Lz= 2.0000000000000000
z= 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000
sh: 1: Syntax error: Bad fd number
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) = 1.2722218725854067E-014
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.0000000000000000
dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000
---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.014
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.010
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.6897194356860323
Wall clock time [hours] = 2.836E-05 (+/- 2.7778E-13)
Wall clock time/timestep/meshpoint [microsec] = 2.658313
Wall clock time/timestep/local meshpoint [microsec] = 10.63325
Rhs wall clock time/timestep/local meshpoint [microsec] = 4.163019
Maximum used memory per cpu [MBytes] = 24.992
Maximum used memory [MBytes] = 97.637
real 0.49
user 0.59
sys 0.10