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= shesheshesheshesheshe : shesheshesheshesheshe bcy1,bcy2= p p p p p p p : p p p p p p p bcz1,bcz2= s s ubsismisms s : s s ubsismisms s WARNING: run: quad precision suppressed! SVN: ------- v. ( ) $Id$ WARNING: register_io: HDF5 files are possibly locked; writing may fail.! 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= 10 (ldebug= F ) This is a 3-D run nxgrid, nygrid, nzgrid= 32 32 64 Lx, Ly, Lz= 0.12500000000000000 0.12500000000000000 0.25000000000000000 Vbox= 3.9062500000000000E-003 WARNING: input_persist_general_by_label: no persistent value of dt found! WARNING: input_persist_general_by_label: no persistent value of eps_rkf found! WARNING: input_persist_interstellar: lSNI from run.in overwritten. Set l_persist_overwrite_lSNI=T to update! input_persist_interstellar: lSNI = T WARNING: input_persist_interstellar: lSNII from run.in overwritten. Set l_persist_overwrite_lSNII=T to update! input_persist_interstellar: lSNII = F WARNING: input_persist_interstellar: t_next_SNI from run.in overwritten. Set l_persist_overwrite_tSNI=T to update! input_persist_interstellar: t_next_SNI = 0.0000000000000000 WARNING: input_persist_interstellar: t_next_SNII from run.in overwritten. Set l_persist_overwrite_tSNII=T to update! input_persist_interstellar: t_next_SNII = 0.0000000000000000 WARNING: input_persist_interstellar: x_cluster from run.in overwritten. Set l_persist_overwrite_xcluster=T to update! WARNING: input_persist_interstellar: y_cluster from run.in overwritten. Set l_persist_overwrite_ycluster=T to update! WARNING: input_persist_interstellar: z_cluster from run.in overwritten. Set l_persist_overwrite_zcluster=T to update! WARNING: input_persist_interstellar: t_cluster from run.in overwritten. Set l_persist_overwrite_tcluster=T to update! rsnap: read snapshot var.dat in 1.9999998621642590E-003 seconds rprint_list: nname= 15 rprint_list: ix,iy,iz,iz2= -1 -1 -1 -1 rprint_list: nnamev= 0 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= -1 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_shock: 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 67 units_general: unit_velocity= 97781.229948950713 units_general: unit_density= 1.6737236000000002E-024 units_general: unit_length= 3.0856775800000002E+021 units_general: unit_magnetic= 3.5449077018110318 unit temperature 1.0000000000000000 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 initialize_modules: unit_velocity, unit_density, etc, are in cgs register: k_B,m_p,m_e,eV= 2.936550E-67 3.401436E-65 1.852480E-68 3.407721E-63 timestep based on CFL cond; cdt= 0.40000000000000002 dVol_rel1= 1.5258789062500000E-005 dA_xy_rel1= 9.7656250000000000E-004 dA_xz_rel1= 4.8828125000000000E-004 dA_yz_rel1= 4.8828125000000000E-004 dL_y_rel1= 3.1250000000000000E-002 initialize_eos: ENTER initialize_eos: reference values for ionization initialize_eos: yHmin,yHmax,yMetals= 2.2250738585072014E-308 0.99999999999999978 0.0000000000000000 initialize_eos: TT_ion,ss_ion,kappa0= 157821.25804611665 6.1758201460340879E-003 22087.988206997077 initialize_eos: lnrho_e,lnrho_H,lnrho_He,lnrho_e_= 53.708496399739161 64.982454679406558 67.051181795731779 49.369847344164853 initialize_gravity: no x-gravity initialize_gravity: no y-gravity xmask_den= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 zmask_den= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 diffusion: nothing (i.e. no mass diffusion) initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. xmask_hyd= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 zmask_hyd= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 initialize_energy: cs2top,cs2cool,cs2cool2= 3.9084999999999999E+037 0.0000000000000000 0.0000000000000000 initialize_energy: now set cs2cool=cs2top initialize_energy: nheatc_max,iheatcond= 4 nothing nothing nothing nothing heat conduction: nothing initialize_interstellar: t_next_SNI = 0.00000000E+00 initialize_interstellar: mu not set by eos: 0.62000000000000000 initialize_interstellar: eampl_SN, kampl_SN = 0.2127E+01 0.0000E+00 select_cooling: default RBN cooling fct initialize_interstellar: unit_Lambda 0.1810E+18 initialize_interstellar: unit_Gamma 0.3030E-06 WARNING: initialize_interstellar: using T/e instead of cv1 for diagnostics. Not yet implemented cv1 for ionization! initialize_interstellar: t_next_SNI, t_next_SNII= 0.00000000E+00 0.00000000E+00 initialize_shear: Sshear,Sshear1= -25.000000000000000 -25.000000000000000 initialize_shear: qshear,qshear0= 1.0000000000000000 0.0000000000000000 viscous force: nu*(del2u+graddivu/3+2S.glnrho) viscous force: nu_shock*(XXXXXXXXXXX) xmask_vis= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_ionization. The values provided are for one-atomic ideal gas. Use at own risk! choose_pencils: finding out which pencils are needed for the pencil case. write_pencil_info: pencil information written to the file pencils.list. Lz= 0.25000000000000000 z= -0.13690476190476189 -0.13293650793650794 -0.12896825396825395 -0.12500000000000000 -0.12103174603174602 -0.11706349206349206 -0.11309523809523808 -0.10912698412698413 -0.10515873015873015 -0.10119047619047619 -9.7222222222222210E-002 -9.3253968253968256E-002 -8.9285714285714274E-002 -8.5317460317460320E-002 -8.1349206349206338E-002 -7.7380952380952384E-002 -7.3412698412698402E-002 -6.9444444444444448E-002 -6.5476190476190466E-002 -6.1507936507936511E-002 -5.7539682539682530E-002 -5.3571428571428575E-002 -4.9603174603174593E-002 -4.5634920634920639E-002 -4.1666666666666657E-002 -3.7698412698412703E-002 -3.3730158730158721E-002 -2.9761904761904767E-002 -2.5793650793650785E-002 -2.1825396825396831E-002 -1.7857142857142849E-002 -1.3888888888888895E-002 -9.9206349206349131E-003 -5.9523809523809590E-003 -1.9841269841269771E-003 1.9841269841269771E-003 5.9523809523809590E-003 9.9206349206349131E-003 1.3888888888888895E-002 1.7857142857142849E-002 2.1825396825396831E-002 2.5793650793650785E-002 2.9761904761904767E-002 3.3730158730158721E-002 3.7698412698412703E-002 4.1666666666666657E-002 4.5634920634920639E-002 4.9603174603174593E-002 5.3571428571428575E-002 5.7539682539682530E-002 6.1507936507936511E-002 6.5476190476190466E-002 6.9444444444444448E-002 7.3412698412698402E-002 7.7380952380952384E-002 8.1349206349206338E-002 8.5317460317460320E-002 8.9285714285714274E-002 9.3253968253968256E-002 9.7222222222222210E-002 0.10119047619047619 0.10515873015873015 0.10912698412698413 0.11309523809523808 0.11706349206349206 0.12103174603174602 0.12500000000000000 0.12896825396825395 0.13293650793650791 0.13690476190476186 boundcond_shear: use shearing sheet boundary condition pencil_consistency_check: checking pencil case life_support_on: needed for pencil consistency check. 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 lnrho ( 12) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil ee ( 55) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cs2 ( 58) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cp1tilde ( 61) 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 life_support_off: end of pencil consistency check. pencil_consistency_check: all tests passed pde: ENTER SVN: ------- v. ( ) $Id$ check_SN: ENTER explode_SN: step, time = 1 0.00000E+00 boundcond_shear: use shearing sheet boundary condition boundcond_shear: use shearing sheet boundary condition boundcond_shear: use shearing sheet boundary condition boundcond_shear: use shearing sheet boundary condition calc_pencils_hydro: max(advec_uu) = 0.0000000000000000 calc_pencils_energy: max(advec_cs2) = 17661506.211830802 duu_dt: SOLVE Bcs for ux, x: < she>, y: < p>, z: < s> Bcs for uy, x: < she>, y: < p>, z: < s> Bcs for uz, x: < she>, y: < p>, z: < ubs> coriolis_cartesian: add Coriolis force; Omega= 25.000000000000000 dlnrho_dt: SOLVE Bcs for lnrho, x: < she>, y: < p>, z: < ism> dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000 denergy_dt: SOLVE denergy_dt Bcs for ss, x: < she>, y: < p>, z: < ism> denergy_dt: lnTT,cs2,cp1= 8.3442678054305901 90.769188378982889 1.1728390000000000 calc_heat_cool_interstellar: ENTER shearing: Sshear,Sshear1= -25.000000000000000 -25.000000000000000 shearing: qshear,qshear0= 1.0000000000000000 0.0000000000000000 --it-------t--------dt-------urms----ekintot----ethtot-----rhom------ssm------TTmin-----TTmax------TTm------dtc-----dtu-----dtnu---dtchi-- 0 0.0000E+00 6.74E-06 0.000E+00 0.000E+00 4.741E+00 6.457E-01 7.604E-01 3.595E+03 8.255E+05 3.929E+04 0.9930 0.0000 0.0874 0.0000 2 1.4126E-05 8.58E-06 8.966E-01 1.223E-03 2.988E+00 6.457E-01 7.509E-01 3.501E+03 4.984E+05 9.576E+03 0.9817 0.0684 0.3599 0.0000 4 3.4408E-05 1.11E-05 1.153E+00 2.001E-03 2.780E+00 6.457E-01 7.475E-01 3.372E+03 3.069E+04 6.071E+03 0.3158 0.4605 0.7195 0.0000 6 6.0063E-05 1.87E-05 1.076E+00 1.707E-03 2.729E+00 6.457E-01 7.450E-01 3.366E+03 1.443E+04 5.262E+03 0.3639 0.4927 0.5586 0.0000 8 1.0038E-04 2.21E-05 1.038E+00 1.559E-03 2.680E+00 6.457E-01 7.419E-01 3.364E+03 1.142E+04 4.560E+03 0.3826 0.3928 0.7298 0.0000 10 1.4437E-04 2.22E-05 1.005E+00 1.433E-03 2.640E+00 6.457E-01 7.387E-01 3.368E+03 1.025E+04 4.070E+03 0.3642 0.3472 0.7103 0.0000 explode_SN: step, time = 12 0.16657E-03 12 1.7095E-04 4.42E-06 1.236E+00 2.162E-03 3.663E+00 6.457E-01 7.505E-01 3.419E+03 1.758E+06 2.249E+04 0.9492 0.0778 0.1826 0.0000 14 1.7975E-04 4.11E-06 1.839E+00 4.618E-03 3.054E+00 6.458E-01 7.468E-01 3.417E+03 1.684E+06 1.209E+04 0.8641 0.1413 0.5053 0.0000 16 1.8738E-04 3.38E-06 2.392E+00 7.084E-03 2.897E+00 6.458E-01 7.451E-01 3.415E+03 1.626E+06 9.372E+03 0.6991 0.2184 0.6929 0.0000 18 1.9375E-04 2.70E-06 2.851E+00 8.915E-03 2.838E+00 6.458E-01 7.441E-01 3.414E+03 1.704E+06 8.381E+03 0.5712 0.3104 0.5387 0.0000 20 1.9879E-04 1.99E-06 3.277E+00 1.019E-02 2.805E+00 6.458E-01 7.434E-01 3.413E+03 2.725E+06 7.951E+03 0.5331 0.3533 0.4629 0.0000 22 2.0247E-04 1.44E-06 3.688E+00 1.111E-02 2.789E+00 6.458E-01 7.430E-01 3.412E+03 5.266E+06 7.884E+03 0.5359 0.3626 0.4372 0.0000 24 2.0514E-04 1.05E-06 4.121E+00 1.184E-02 2.779E+00 6.458E-01 7.427E-01 3.412E+03 1.043E+07 8.088E+03 0.5488 0.3585 0.4192 0.0000 26 2.0708E-04 7.65E-07 4.620E+00 1.248E-02 2.774E+00 6.458E-01 7.424E-01 3.411E+03 2.094E+07 8.564E+03 0.5676 0.3462 0.4050 0.0000 28 2.0851E-04 6.01E-07 5.090E+00 1.303E-02 2.771E+00 6.458E-01 7.423E-01 3.411E+03 3.926E+07 9.354E+03 0.6106 0.3229 0.3575 0.0000 Simulation finished after 30 time-steps Writing final snapshot at time t = 2.0966551948464212E-004 Wall clock time [hours] = 6.149E-03 (+/- 2.7778E-07) Wall clock time/timestep/meshpoint [microsec] = 11.25844 Wall clock time/timestep/local meshpoint [microsec] = 11.25844 Rhs wall clock time/timestep/local meshpoint [microsec] = 26.69830 After substep wall clock time/timestep/local meshpoint [microsec] = 5.0862639E-04 Before boundary wall clock time/timestep/local meshpoint [microsec] = 15.66772 Diagnostics wall clock time/timestep/local meshpoint [microsec] = 0.000000 Timestep wall clock time/timestep/local meshpoint [microsec] = 11.24878 Timestep (without waiting) wall clock time/timestep/local meshpoint [microsec] = 11.24878 Timestep wall clock time/timestep/meshpoint [microsec] = 11.24878 Timestep (without waiting) wall clock time/timestep/meshpoint [microsec] = 11.24878 Diagnostics wall clock time/timestep/meshpoint [microsec] = 0.000000 Rhs wall clock time/timestep/meshpoint [microsec] = 26.69830 After substep wall clock time/timestep/meshpoint [microsec] = 5.0862639E-04 Before boundary wall clock time/timestep/meshpoint [microsec] = 15.66772 Maximum used memory per cpu [MBytes] = 45.422 Maximum used memory [MBytes] = 45.422 real 105.59 user 104.67 sys 0.85