SVN: ------- v. ( ) $Id$ WARNING: find_namelist: namelist "special_init_pars" is missing!! bcx1,bcx2= p p p p p : p p p p p bcy1,bcy2= p p p p p : p p p p p bcz1,bcz2= p p p p p : p p p p p lperi= T T T 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: ------- 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= 320 320 1 Lx, Ly, Lz= 5.2000000000000002 5.2000000000000002 0.52000000000000002 Vbox= 27.040000000000003 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 2.5000000000000005E-003 1.0000000000000000 diffusion: nothing (i.e. no mass diffusion) 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 initialize_shock: prenormalised shock_factor sum= 1.0000000000000000 viscous force: nu*(del2u+graddivu/3+2S.glnrho) viscous force: nu_shock*(XXXXXXXXXXX) rt_int,rt_ext 0.40000000000000002 2.5000000000000000 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_energy: max(advec_cs2) = 7.3054519246106828 duu_dt: SOLVE Bcs for ux, x: < p>, y: < p>, z: < p> Bcs for uy, x: < p>, y: < p>, z: < p> Bcs for uz, x: < p>, y: < p>, z: < p> duu_dt: max(advec_uu) = 48.537096255183982 dlnrho_dt: SOLVE Bcs for lnrho, x: < p>, y: < p>, z: < p> dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000E+000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000E+000 dspecial_dt: SOLVE dSPECIAL_dt pde: freezing variables for r < 0.40000000000000002 : T T T T F pde: freezing variables for r > 2.5000000000000000 : T T T T F ----it--------t---------dt------totmass----totenergy---totangmom------TTp-------torqint_1----torqext_1-------xq2----------yq2----------xq1----------yq1---------vxq2---------vyq2---------vxq1---------vyq1----------urm----------upm----- 0 0.000 2.78E-03 2.704E+01 -3.457E-01 1.283E+00 1.725E-03 0.000000000 0.000000000 1.000E-04 0.000E+00 -9.999E-01 0.000E+00 0.000E+00 1.000E-04 0.000E+00 -9.999E-01 0.000000000 4.421E-05 5 0.014 2.78E-03 2.704E+01 -3.457E-01 1.283E+00 1.725E-03 -0.000003360 0.000001803 9.999E-05 1.390E-06 -9.998E-01 -1.389E-02 -1.390E-06 9.999E-05 1.389E-02 -9.998E-01 -0.000000003 4.421E-05 10 0.028 2.78E-03 2.704E+01 -3.457E-01 1.283E+00 1.725E-03 -0.000007881 0.000005840 9.996E-05 2.779E-06 -9.995E-01 -2.779E-02 -2.779E-06 9.996E-05 2.779E-02 -9.995E-01 -0.000000005 4.421E-05 15 0.042 2.78E-03 2.704E+01 -3.457E-01 1.283E+00 1.725E-03 -0.000004136 0.000003298 9.991E-05 4.168E-06 -9.990E-01 -4.167E-02 -4.168E-06 9.991E-05 4.167E-02 -9.990E-01 -0.000000006 4.421E-05 20 0.056 2.78E-03 2.704E+01 -3.457E-01 1.283E+00 1.725E-03 -0.000001121 0.000002133 9.985E-05 5.556E-06 -9.984E-01 -5.555E-02 -5.556E-06 9.985E-05 5.555E-02 -9.984E-01 -0.000000007 4.421E-05 25 0.069 2.78E-03 2.704E+01 -3.457E-01 1.283E+00 1.725E-03 0.000000218 0.000001666 9.976E-05 6.942E-06 -9.975E-01 -6.942E-02 -6.942E-06 9.976E-05 6.942E-02 -9.975E-01 -0.000000007 4.421E-05 30 0.083 2.78E-03 2.704E+01 -3.457E-01 1.283E+00 1.725E-03 0.000000973 0.000000124 9.965E-05 8.328E-06 -9.964E-01 -8.327E-02 -8.328E-06 9.965E-05 8.327E-02 -9.964E-01 -0.000000007 4.421E-05 35 0.097 2.78E-03 2.704E+01 -3.457E-01 1.283E+00 1.725E-03 0.000001394 -0.000001840 9.953E-05 9.712E-06 -9.952E-01 -9.711E-02 -9.712E-06 9.953E-05 9.711E-02 -9.952E-01 -0.000000006 4.421E-05 40 0.111 2.78E-03 2.704E+01 -3.457E-01 1.283E+00 1.725E-03 0.000001571 -0.000003050 9.938E-05 1.109E-05 -9.937E-01 -1.109E-01 -1.109E-05 9.938E-05 1.109E-01 -9.937E-01 -0.000000005 4.421E-05 Simulation finished after 41 time-steps Writing final snapshot at time t = 0.11394614606703260 Wall clock time [hours] = 5.867E-03 (+/- 2.7778E-07) Wall clock time/timestep/meshpoint [microsec] = 5.030488 Maximum used memory per cpu [MBytes] = 188.453 Maximum used memory [MBytes] = 188.453 real 23.74 user 23.47 sys 0.19