Pencil1: iapp, nprocs, ncpus= 0 4 4 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 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 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) 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.0000000000000000E+000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000E+000 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] = 3.031E-05 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 5.328572 Maximum used memory per cpu [MBytes] = 25.070 Maximum used memory [MBytes] = 98.445 real 0.54 user 0.71 sys 0.14