Pencil1: iapp, nprocs, ncpus= 0 4 4 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 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 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) 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.0000000000000000E+000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000E+000 ----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] = 3.416E-05 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 9.608739 Maximum used memory per cpu [MBytes] = 27.105 Maximum used memory [MBytes] = 105.879 real 0.63 user 0.94 sys 0.16