Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified 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 WARNING: run: quad precision suppressed! 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 5.0795300000000001E-004 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) pencil_consistency_check: checking pencil case At line 1153 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_hydro: max(advec_uu) = 4.9637728963999841 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> 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.5453963072360806 Wall clock time [hours] = 1.533E-05 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 4.310576 Wall clock time/timestep/local meshpoint [microsec] = 17.24230 Rhs wall clock time/timestep/local meshpoint [microsec] = 61.63880 Maximum used memory per cpu [MBytes] = 35.754 Maximum used memory [MBytes] = 140.766 real 1.53 user 0.71 sys 0.22