SVN: ------- v. ( ) $Id$
SVN: start.in,v v. 1.2 (2020/05/19 07:33:50) wlyra
SVN: run.in,v v. 1.2 (2020/05/19 07:33:50) wlyra
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 2-D run
nxgrid, nygrid, nzgrid= 32 64 1
Lx, Ly, Lz= 2.1000000000000001 6.2831853071795862 0.52000000000000002
Vbox= 19.132299260361840
rsnap: read snapshot var.dat in 1.1284300000000000E-004 seconds
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 5.0000000000000013E-006 1.0000000000000000
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
WARNING: initialize_diagnostics: volume averages are calculated wrongly for nonequidistant grids unless lproper_averages=T.!
diffusion: mesh hyperdiffusion
Note: rho_const or lnrho_const are not used
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)
viscous force: nu*(del2u+graddivu/3+2S.glnrho)
sh: 1: Syntax error: Bad fd number
Adding gas+dust gravity to the massive particles
pde: ENTER
SVN: ------- v. ( ) $Id$
calc_pencils_hydro: call gij_etc
calc_pencils_energy: max(advec_cs2) = 15.230660441023506
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) = 40.137349438062927
dlnrho_dt: SOLVE
Bcs for lnrho, x: < a2>, y: < p>, z: < p>
dlnrho_dt: diffrho_hyper3_mesh= 30.000000000000000
dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000
dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000
At line 1154 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'xxq' of procedure 'get_evr'
----it--------t-----------dt-----------rhomax---------rhomin----------urms-----------uy2m------------xq1------------yq1------------xq2------------yq2-----------vxq1-----------vyq1-----------vxq2-----------vyq2---------torqint_2------torqext_2---
0 0.000 7.86103E-03 7.90569E-03 5.05964E-04 8.27856E-01 6.85345E-01 3.90000E-05 3.14159E+00 1.29996E+00 0.00000E+00 0.00000E+00 2.20140E-05 0.00000E+00 7.33777E-01 1.32349E-22 1.58819E-22
1 0.008 7.86100E-03 7.90570E-03 5.05964E-04 8.27856E-01 6.85345E-01 3.89998E-05 3.14603E+00 1.29996E+00 4.43711E-03 -4.18323E-08 2.20141E-05 -1.26082E-03 7.33710E-01 -4.13893E-09 -6.29650E-07
2 0.016 7.86098E-03 7.90573E-03 5.05964E-04 8.27856E-01 6.85346E-01 3.89993E-05 3.15047E+00 1.29994E+00 8.87355E-03 -8.36590E-08 2.20144E-05 -2.55600E-03 7.33567E-01 -7.87999E-09 -7.93691E-07
3 0.024 7.86096E-03 7.90577E-03 5.05964E-04 8.27856E-01 6.85346E-01 3.89984E-05 3.15491E+00 1.29992E+00 1.33092E-02 -1.25475E-07 2.20148E-05 -3.88624E-03 7.33429E-01 -1.08822E-08 -6.84094E-07
4 0.031 7.86095E-03 7.90582E-03 5.05964E-04 8.27857E-01 6.85347E-01 3.89971E-05 3.15934E+00 1.29988E+00 1.77442E-02 -1.67281E-07 2.20155E-05 -5.23886E-03 7.33325E-01 -1.29027E-08 -5.26176E-07
5 0.039 7.86093E-03 7.90588E-03 5.05964E-04 8.27857E-01 6.85347E-01 3.89955E-05 3.16378E+00 1.29983E+00 2.21788E-02 -2.09079E-07 2.20163E-05 -6.60403E-03 7.33257E-01 -1.38146E-08 -3.91620E-07
6 0.047 7.86091E-03 7.90595E-03 5.05964E-04 8.27857E-01 6.85348E-01 3.89935E-05 3.16822E+00 1.29978E+00 2.66132E-02 -2.50872E-07 2.20173E-05 -7.97617E-03 7.33219E-01 -1.36030E-08 -2.87931E-07
7 0.055 7.86090E-03 7.90602E-03 5.05964E-04 8.27857E-01 6.85348E-01 3.89912E-05 3.17267E+00 1.29971E+00 3.10477E-02 -2.92660E-07 2.20185E-05 -9.35227E-03 7.33205E-01 -1.23454E-08 -2.08252E-07
8 0.063 7.86088E-03 7.90610E-03 5.05964E-04 8.27858E-01 6.85348E-01 3.89885E-05 3.17711E+00 1.29963E+00 3.54824E-02 -3.34444E-07 2.20199E-05 -1.07306E-02 7.33212E-01 -1.01857E-08 -1.45016E-07
9 0.071 7.86087E-03 7.90619E-03 5.05964E-04 8.27858E-01 6.85349E-01 3.89855E-05 3.18155E+00 1.29954E+00 3.99175E-02 -3.76225E-07 2.20215E-05 -1.21103E-02 7.33237E-01 -7.30967E-09 -9.23766E-08
Simulation finished after 10 time-steps
Writing final snapshot at time t = 7.8609409231685606E-002
Wall clock time [hours] = 1.193E-05 (+/- 2.7778E-13)
Wall clock time/timestep/meshpoint [microsec] = 2.097327
Wall clock time/timestep/local meshpoint [microsec] = 8.389310
Rhs wall clock time/timestep/local meshpoint [microsec] = 6.784216
Maximum used memory per cpu [MBytes] = 25.098
Maximum used memory [MBytes] = 99.812
real 0.39
user 0.37
sys 0.09