SVN: ------- v. ( ) $Id$
bcx1,bcx2= a s s s s : a s s s s
bcy1,bcy2= s a s s s : s a s s s
bcz1,bcz2= p p p p p : p p p p p
lperi= F F T
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
initialize_mpicomm: enabled MPI
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$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: pointmasses.f90,v v. 1.1 (2019/02/02 03:54:41) wlyra
The verbose level is ip= 14 (ldebug= F )
This is a 3-D run
nxgrid, nygrid, nzgrid= 8 4 8
Lx, Ly, Lz= 2.1000000000000001 0.17453000000000007 6.2831840000000003
Vbox= 5.6808689339176519
rsnap: read snapshot var.dat in 1.9137144088745117E-002 seconds
WARNING: setup_slice: slice_position=w may be wrong for nprocx>1!
setup_slices: slice_position = w
setup_slices: iz_loc, (video files) = 7
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 -4.6051701859880909 0.10000000000000001 1.0000000000000004E-002 1.0000000000000000
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
diffusion: mesh hyperdiffusion
initialize_density: no need to read initial stratification for lanti_shockdiffusion=F.
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)
sh: 1: Syntax error: Bad fd number
pencil_consistency_check: checking pencil case
At line 1154 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'xxq' of procedure 'get_evr'
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_hydro: call gij_etc
calc_pencils_energy: max(advec_cs2) = 46.699031468702778
duu_dt: SOLVE
Bcs for ux, x: < a>, y: < s>, z: < p>
Bcs for uy, x: < s>, y: < a>, z: < p>
Bcs for uz, x: < s>, y: < s>, z: < p>
duu_dt: max(advec_uu) = 4.9637728963999841
dlnrho_dt: SOLVE
Bcs for lnrho, x: < s>, y: < s>, z: < p>
dlnrho_dt: diffrho_hyper3_mesh= 5.0000000000000000
dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000
dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000
----it--------t---------dt-------dtc------dtv------dtu--------xq1----------yq1---------zq1----------xq2----------yq2----------zq2---------vzq1---------vzq2-----
0 0.000 3.13E-02 5.3E-01 1.0E+00 3.9E-01 9.990E-01 1.571E+00 0.00E+00 1.000E-03 1.571E+00 -3.142E+00 9.990E-01 1.000E-03
5 0.154 3.04E-02 5.2E-01 1.0E+00 4.0E-01 9.990E-01 1.571E+00 1.54E-01 1.000E-03 1.571E+00 -2.988E+00 9.990E-01 1.000E-03
10 0.307 3.11E-02 5.3E-01 1.0E+00 3.9E-01 9.990E-01 1.571E+00 3.07E-01 1.000E-03 1.571E+00 -2.835E+00 9.990E-01 1.000E-03
15 0.462 3.06E-02 5.2E-01 1.0E+00 4.0E-01 9.990E-01 1.571E+00 4.62E-01 1.000E-03 1.571E+00 -2.680E+00 9.990E-01 1.000E-03
20 0.615 3.09E-02 5.3E-01 1.0E+00 3.9E-01 9.990E-01 1.571E+00 6.15E-01 1.000E-03 1.571E+00 -2.526E+00 9.990E-01 1.000E-03
25 0.770 3.09E-02 5.3E-01 1.0E+00 3.9E-01 9.990E-01 1.571E+00 7.70E-01 1.000E-03 1.571E+00 -2.371E+00 9.990E-01 1.000E-03
30 0.925 3.10E-02 5.3E-01 1.0E+00 3.9E-01 9.990E-01 1.571E+00 9.25E-01 1.000E-03 1.571E+00 -2.217E+00 9.990E-01 1.000E-03
35 1.080 3.10E-02 5.3E-01 1.0E+00 3.9E-01 9.990E-01 1.571E+00 1.08E+00 1.000E-03 1.571E+00 -2.061E+00 9.990E-01 1.000E-03
40 1.235 3.09E-02 5.3E-01 1.0E+00 4.0E-01 9.990E-01 1.571E+00 1.23E+00 1.000E-03 1.571E+00 -1.907E+00 9.990E-01 1.000E-03
45 1.390 3.12E-02 5.3E-01 1.0E+00 3.9E-01 9.990E-01 1.571E+00 1.39E+00 1.000E-03 1.571E+00 -1.752E+00 9.990E-01 1.000E-03
Simulation finished after 50 time-steps
Writing final snapshot at time t = 1.5453963072360808
Wall clock time [hours] = 4.021E-04 (+/- 2.7778E-10)
Wall clock time/timestep/meshpoint [microsec] = 113.0941
Wall clock time/timestep/local meshpoint [microsec] = 452.3763
Rhs wall clock time/timestep/local meshpoint [microsec] = 266.1775
Maximum used memory per cpu [MBytes] = 9.578
Maximum used memory [MBytes] = 36.387
real 4.65
user 1.08
sys 0.82