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 8.0000000307336450E-003 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: 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
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
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
dL_y_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: nothing (i.e. no mass diffusion)
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
choose_pencils: finding out which pencils are needed for the pencil case.
write_pencil_info: pencil information written to the file pencils.list.
Lz= 0.52000000000000002
z= -0.26000000000000001 -0.26000000000000001 -0.26000000000000001 -0.26000000000000001 -0.26000000000000001 -0.26000000000000001 -0.26000000000000001
sh: 1: Syntax error: Bad fd number
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 ( 13) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil grho ( 16) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil potself ( 183) 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: max(diffus_diffrho ) = 0.0000000000000000
dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000
----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 0.0000000000000000 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.0000000000000000 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.0000000000000000 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 9.9999993108212948E-004 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.0000000000000000 seconds
25 0.454 1.82E-02 1.00E-01 0.2528 3.964 0.7140 2.589E-07 0.7312 0.000E+00 0.000E+00 3.146E-02 1.828E-03 1.25822 1.56E+00 1.64E+00
wvid: written slices in 9.9999993108212948E-004 seconds
30 0.545 1.82E-02 1.00E-01 0.2528 3.961 0.7140 2.754E-07 0.7312 0.000E+00 0.000E+00 3.146E-02 1.828E-03 1.25825 1.54E+00 1.63E+00
wvid: written slices in 0.0000000000000000 seconds
35 0.636 1.82E-02 1.00E-01 0.2528 3.962 0.7140 2.915E-07 0.7312 0.000E+00 0.000E+00 3.146E-02 1.827E-03 1.25828 1.53E+00 1.64E+00
wvid: written slices in 0.0000000000000000 seconds
40 0.727 1.82E-02 1.00E-01 0.2528 3.959 0.7139 3.071E-07 0.7312 0.000E+00 0.000E+00 3.146E-02 1.827E-03 1.25832 1.52E+00 1.66E+00
45 0.818 1.82E-02 1.00E-01 0.2528 3.958 0.7139 3.218E-07 0.7312 0.000E+00 0.000E+00 3.145E-02 1.826E-03 1.25836 1.52E+00 1.67E+00
wvid: written slices in 0.0000000000000000 seconds
Simulation finished after 46 time-steps
Writing final snapshot at time t = 0.83582879615265770
Wall clock time [hours] = 5.722E-05 (+/- 2.7778E-07)
Wall clock time/timestep/meshpoint [microsec] = 4.373302
Wall clock time/timestep/local meshpoint [microsec] = 4.373302
Rhs wall clock time/timestep/local meshpoint [microsec] = 8.130944
Maximum used memory per cpu [MBytes] = 14.250
Maximum used memory [MBytes] = 14.250
real 1.86
user 1.22
sys 0.02