Pencil1: iapp, nprocs, ncpus= 0 4 4 SVN: ------- v. ( ) $Id$ SVN: start.in,v v. 1.2 (2016/10/14 23:26:55) wlyra SVN: run.in,v v. 1.2 (2016/10/14 23:26:55) 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= 256 768 1 Lx, Ly, Lz= 1.2000000000000002 6.2831853071795862 0.20000000000000001 Vbox= 7.5398223686155044 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= 4.8108883747190188E-009 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 -5.5859994389998171 5.0000000000000003E-002 1.5000000000000011E-003 0.40000000000000008 diffusion: mesh hyperdiffusion initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. select_eos_variable: Using rho and ss initialize_energy: nheatc_max,iheatcond= 4 hyper3-mesh nothing nothing nothing heat conduction: hyperdiffusivity of ss viscous force: nu_hyper3_mesh/pi^5 *(Deltav)^6/Deltaq viscous force: nu_shock*(XXXXXXXXXXX) pde: ENTER SVN: ------- v. ( ) $Id$ calc_pencils_hydro: call gij_etc calc_pencils_energy: max(advec_cs2) = 865.83658027986576 WARNING: calc_pencils_viscosity: viscous heating not implemented for lvisc_hyper3_mesh! 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) = 482.79792988503829 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.0000000000000000E+000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000E+000 denergy_dt: SOLVE denergy_dt Bcs for ss, x: < a2>, y: < p>, z: < p> denergy_dt: lnTT,cs2,cp1= -4.6697087071256620 6.2500000000000003E-003 0.0000000000000000E+000 calc_heatcond_hyper3_mesh: chi_hyper3= 0.0000000000000000E+000 calc_heatcond_hyper3: added thdiff 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-----------ux2m-----------uy2m------------xq1------------yq1------------xq2------------yq2-----------vxq1-----------vyq1-----------vyq2---------torqext_2------torqint_2--- 0 0.000 7.21558E-04 1.00000E+00 1.00000E+00 9.99250E-01 0.00000E+00 9.98500E-01 1.26000E-05 -3.14159E+00 9.99987E-01 0.00000E+00 0.00000E+00 1.26000E-05 9.99987E-01 -3.38813E-20 2.71051E-20 5 0.004 7.21558E-04 1.00002E+00 9.99999E-01 9.99250E-01 5.81655E-13 9.98500E-01 1.26000E-05 -3.13798E+00 9.99987E-01 3.60779E-03 2.73034E-11 1.26000E-05 9.99987E-01 6.10713E-12 6.00591E-12 10 0.007 7.21558E-04 1.00007E+00 9.99997E-01 9.99250E-01 2.32465E-12 9.98500E-01 1.26000E-05 -3.13438E+00 9.99987E-01 7.21558E-03 5.46914E-11 1.26000E-05 9.99987E-01 -1.38165E-11 -1.02103E-11 15 0.011 7.21558E-04 1.00017E+00 9.99993E-01 9.99250E-01 5.22572E-12 9.98500E-01 1.26000E-05 -3.13077E+00 9.99987E-01 1.08234E-02 8.21628E-11 1.26000E-05 9.99987E-01 9.14435E-12 1.73769E-11 20 0.014 7.21557E-04 1.00030E+00 9.99987E-01 9.99250E-01 9.28119E-12 9.98500E-01 1.26000E-05 -3.12716E+00 9.99987E-01 1.44312E-02 1.09717E-10 1.26000E-05 9.99987E-01 -3.30354E-11 -9.68500E-12 25 0.018 7.21557E-04 1.00044E+00 9.99980E-01 9.99250E-01 1.44870E-11 9.98500E-01 1.26000E-05 -3.12355E+00 9.99987E-01 1.80389E-02 1.37352E-10 1.26000E-05 9.99987E-01 -1.32815E-11 2.86737E-11 30 0.022 7.21557E-04 1.00067E+00 9.99972E-01 9.99250E-01 2.08387E-11 9.98500E-01 1.26000E-05 -3.11995E+00 9.99987E-01 2.16467E-02 1.65068E-10 1.26000E-05 9.99987E-01 -6.56673E-11 9.75896E-12 35 0.025 7.21557E-04 1.00084E+00 9.99962E-01 9.99250E-01 2.83314E-11 9.98500E-01 1.26000E-05 -3.11634E+00 9.99987E-01 2.52545E-02 1.92864E-10 1.26000E-05 9.99987E-01 -7.24493E-11 4.49498E-11 40 0.029 7.21556E-04 1.00117E+00 9.99952E-01 9.99250E-01 3.69599E-11 9.98500E-01 1.26000E-05 -3.11273E+00 9.99987E-01 2.88623E-02 2.20739E-10 1.26000E-05 9.99987E-01 -1.23903E-10 5.22231E-11 45 0.032 7.21556E-04 1.00143E+00 9.99940E-01 9.99250E-01 4.67186E-11 9.98500E-01 1.26000E-05 -3.10912E+00 9.99987E-01 3.24701E-02 2.48692E-10 1.26000E-05 9.99987E-01 -1.74388E-10 7.59112E-11 Simulation finished after 46 time-steps Writing final snapshot at time t = 3.3191633335019483E-002 Wall clock time [hours] = 1.800E-02 (+/- 2.7778E-10) Wall clock time/timestep/meshpoint [microsec] = 7.166665 Maximum used memory per cpu [MBytes] = 66.203 Maximum used memory [MBytes] = 263.766 real 66.62 user 104.48 sys 5.62