===================================================================== WARNING: using SINGLE PRECISION, you'd better know what you're doing! ===================================================================== SVN: ------- v. ( ) $Id$ WARNING: find_namelist: namelist "_pars" is missing!! SVN: ------- v. ( ) $Id$ WARNING: find_namelist: namelist "_pars" is missing!! SVN: ------- v. ( ) $Id$ bcx1,bcx2= a s s a2 : a s s a2 bcy1,bcy2= s a s a2 : s a s a2 bcz1,bcz2= p p p p : p p p p 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: ------- v. ( ) $Id$ The verbose level is ip= 13 (ldebug= F ) This is a 3-D run nxgrid, nygrid, nzgrid= 32 32 32 Lx, Ly, Lz= 4.19999981 4.19999981 4.00000000 Vbox= 70.5599899 rsnap: read snapshot var.dat in 5.8282199999999998E-004 seconds rprint_list: nname= 10 rprint_list: ix,iy,iz,iz2= -1 -1 -1 -1 rprint_list: nnamev= 1 rprint_list: nnamez= -1 rprint_list: nnamez= -1 rprint_list: nnamey= -1 rprint_list: nnamex= -1 rprint_list: nnamer= 0 rprint_list: nnamexz= -1 rprint_list: nnamexy= -1 rprint_list: nnamerz= 6 rprint_general: run through parse list rprint_hydro: run through parse list rprint_hydro: run through parse list rprint_density: run through parse list rprint_viscosity: run through parse list setup_slices: slice_position = p setup_slices: ix_loc,iy_loc,iz_loc,iz2_loc, (video files) = 4 4 4 35 units_general: unit_velocity= 1.0000000000000000 units_general: unit_density= 1.0000000000000000 units_general: unit_length= 1.0000000000000000 units_general: unit_magnetic= 3.5449078083038330 units_eos: unit_temperature= 1.2027220775801925E-008 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.00000000 2.19722462 3.00000000 9.00000000 1.00000000 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 timestep based on CFL cond; cdt= 0.400000006 xmask_den= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 zmask_den= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 diffusion: nothing (i.e. no mass diffusion) initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. initialize_hydro: set outer radius rdampext=r_ext xmask_hyd= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 zmask_hyd= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 select_eos_variable: Using lnrho and cs2 viscous force: nu*(del2u+graddivu/3+2S.glnrho) xmask_vis= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 choose_pencils: finding out which pencils are needed for the pencil case. write_pencil_info: pencil information written to the file pencils.list. pencil_consistency_check: checking pencil case 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: performing full pencil check (takes a while) pencil_consistency_check: possible overcalculation... pencil sij2 ( 114) is requested, but does not appear to be required! 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: performing full pencil check (takes a while) 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: max(advec_uu) = 0.356653839 calc_pencils_energy: max(advec_cs2) = 1556.61243 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: < a2>, y: < a2>, z: < p> dlnrho_dt: max(diffus_diffrho ) = 0.00000000 dlnrho_dt: max(diffus_diffrho3) = 0.00000000 ----it--------t--------dt-------urms------umax------epsK------rhom-----dtc-----dtu-----dtnu-- 0 0.000 9.98E-03 0.02455 0.04249 1.353E-03 1.00000 0.9844 0.0139 0.0575 2 0.019 8.32E-03 0.21181 0.82631 1.379E-02 1.00000 0.8202 0.1786 0.0479 4 0.035 8.14E-03 0.23773 0.92588 1.747E-02 1.00000 0.8031 0.1958 0.0469 6 0.051 8.11E-03 0.24346 0.94320 1.810E-02 1.00000 0.7999 0.1990 0.0468 8 0.068 8.10E-03 0.24562 0.94798 1.812E-02 1.00000 0.7989 0.2000 0.0467 Simulation finished after 10 time-steps Writing final snapshot at time t = 8.3788161631673574E-002 Wall clock time [hours] = 2.564E-04 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 2.816401 Wall clock time/timestep/local meshpoint [microsec] = 2.816401 Rhs wall clock time/timestep/local meshpoint [microsec] = 53.87893 After substep wall clock time/timestep/local meshpoint [microsec] = 2.9103830E-05 Before boundary wall clock time/timestep/local meshpoint [microsec] = 5.4249540E-03 Diagnostics wall clock time/timestep/local meshpoint [microsec] = 0.000000 Timestep wall clock time/timestep/local meshpoint [microsec] = 2.810155 Timestep (without waiting) wall clock time/timestep/local meshpoint [microsec] = 2.810155 Timestep wall clock time/timestep/meshpoint [microsec] = 2.810155 Timestep (without waiting) wall clock time/timestep/meshpoint [microsec] = 2.810155 Diagnostics wall clock time/timestep/meshpoint [microsec] = 0.000000 Rhs wall clock time/timestep/meshpoint [microsec] = 53.87893 After substep wall clock time/timestep/meshpoint [microsec] = 2.9103830E-05 Before boundary wall clock time/timestep/meshpoint [microsec] = 5.4249540E-03 Maximum used memory per cpu [MBytes] = 30.129 Maximum used memory [MBytes] = 30.129 real 19.23 user 18.85 sys 0.12