Pencil1: iapp, nprocs, ncpus= 0 4 4 SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) 20a96485a SVN: ------- v. ( ) 5d67d8a3f initialize_mpicomm: enabled MPI 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= 1963 (ldebug= F ) This is a 3-D run nxgrid, nygrid, nzgrid= 16 16 200 Lx, Ly, Lz= 6.4000000000000001E-002 6.4000000000000001E-002 0.79600000000000004 Vbox= 3.2604159999999999E-003 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 = T 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 = 3.5137242741072820E-002 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.0000000000000000E+000 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! setup_slices: slice_position = m setup_slices: ix_loc,iy_loc, (video files) = 11 11 lcourant_dt T dt 1.0000000000000000E-008 units_general: unit_velocity= 97781.229948950713 units_general: unit_density= 1.6737236000000002E-024 units_general: unit_length= 3.0856775800000002E+021 units_general: unit_magnetic= 4.4843788716874900E-007 units_eos: unit_temperature= 36.637181156923013 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.5000000000000000 0.0000000000000000E+000 1.0000000000000000 0.31859999999999994 0.31859999999999994 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 initialize_gravity: no x-gravity initialize_gravity: no y-gravity diffusion: shock diffusion select_eos_variable: Using rho and ss initialize_energy: cs2top,cs2cool,cs2cool2= 1.0000000000000000 0.0000000000000000E+000 0.0000000000000000E+000 initialize_energy: now set cs2cool=cs2top initialize_energy: nheatc_max,iheatcond= 4 shock hyper3 nothing nothing heat conduction: shock heat conduction: hyperdiffusivity of ss resistivity: hyper3 resistivity: constant eta initialize_interstellar: t_next_SNI = 0.35137243E-01 initialize_interstellar: solar_mass = 0.4045E-07 initialize_interstellar: eampl_SN, kampl_SN = 0.2127E+01 0.0000E+00 initialize_interstellar: Lxyzmin = 0.1149E-01 select_cooling: WSW cooling fct initialize_interstellar: unit_Lambda 0.1810E+18 initialize_interstellar: unit_Gamma 0.3030E-06 initialize_interstellar: t_interval_SNI, SNI rate = 0.6105E-01 0.3999E+04 initialize_interstellar: average_SNI_heating = 0 initialize_interstellar: average_SNII_heating =0 initialize_interstellar: nseed,seed 1 1812 initialize_interstellar: finished initialize_interstellar: t_next_SNI, t_next_SNII= 0.35137243E-01 0.00000000E+00 viscous force: nu_shock*(XXXXXXXXXXX) viscous force: nu*(del6u+S.glnrho) viscous force: nu*(del2u+graddivu/3+2S.glnrho) 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 rho ( 13) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil grho ( 16) 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 farraymin 1.0000000000000001E-009 1.0000000000000001E-009 1.0000000000000001E-009 2.5019037577498496E-008 7.4119482392732623E-007 1.0000000000000001E-009 1.0000000000000001E-009 1.0000000000000001E-009 pde: ENTER SVN: ------- v. ( ) $Id$ check_SN: ENTER check_SNI: ENTER check_SNIIb: ENTER position_SN_gaussianz: ENTER position_SN_gaussianz: ENTER position_SN_gaussianz: ENTER position_SN_gaussianz: ENTER position_SN_gaussianz: ENTER position_SN_gaussianz: ENTER position_SN_gaussianz: ENTER position_SN_gaussianz: ENTER position_SN_gaussianz: ENTER position_SN_gaussianz: ENTER WARNING: calc_pencils_viscosity: viscous heating not implemented for lvisc_hyper3_nu_const! 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: < s> coriolis_cartesian: add Coriolis force; Omega= 25.000000000000000 dlnrho_dt: SOLVE Bcs for lnrho, x: < she>, y: < p>, z: < s> dlnrho_dt: diffrho_shock= 1.0000000000000000 denergy_dt: SOLVE denergy_dt Bcs for ss, x: < she>, y: < p>, z: < s> denergy_dt: lnTT,cs2,cp1= 5.6757182395299903 291.69777219535422 1.1728390000000000 calc_heatcond_shock: use shock diffusion calc_heatcond_shock: chi_shock= 1.0000000000000000 calc_heatcond_shock: added thdiff calc_heatcond_hyper3: chi_hyper3= 7.5000000000000000E-012 calc_heatcond_hyper3: added thdiff calc_heat_cool_interstellar: ENTER daa_dt: SOLVE Bcs for Ax, x: < she>, y: < p>, z: < s> Bcs for Ay, x: < she>, y: < p>, z: < s> Bcs for Az, x: < she>, y: < p>, z: < s> daa_dt: iresistivity=hyper3 eta-const shearing: Sshear,Sshear1= -25.000000000000000 -25.000000000000000 shearing: qshear,qshear0= 1.0000000000000000 0.0000000000000000E+000 ---it--------t----------dt--------urms-------umax-------rhomin------rhomax------TTmin-------TTm-------TTmax------brms-------bmax---- 0 0.000000000 1.0000E-08 0.0000E+00 0.0000E+00 2.33610E-02 8.22356E-01 1.8983E+02 1.1303E+03 2.0256E+06 5.6913E-04 2.1548E-03 1 0.000000010 1.5394E-09 6.8747E-02 3.4102E+00 2.33610E-02 8.22356E-01 1.8973E+02 1.1306E+03 2.0249E+06 5.6850E-04 2.1521E-03 WARNING: powerhel: computation of wavevector is wrong for non-cubical domains! WARNING: powerhel: computation of wavevector is wrong for non-cubical domains! 2 0.000000012 1.1603E-09 7.9331E-02 3.9352E+00 2.33610E-02 8.22356E-01 1.8971E+02 1.1307E+03 2.0248E+06 5.6841E-04 2.1517E-03 WARNING: powerhel: computation of wavevector is wrong for non-cubical domains! WARNING: powerhel: computation of wavevector is wrong for non-cubical domains! 3 0.000000013 1.0882E-09 8.7308E-02 4.3308E+00 2.33610E-02 8.22356E-01 1.8970E+02 1.1307E+03 2.0247E+06 5.6833E-04 2.1514E-03 WARNING: powerhel: computation of wavevector is wrong for non-cubical domains! WARNING: powerhel: computation of wavevector is wrong for non-cubical domains! 4 0.000000014 1.1381E-09 9.4789E-02 4.7019E+00 2.33610E-02 8.22356E-01 1.8969E+02 1.1308E+03 2.0246E+06 5.6827E-04 2.1511E-03 WARNING: powerhel: computation of wavevector is wrong for non-cubical domains! WARNING: powerhel: computation of wavevector is wrong for non-cubical domains! Simulation finished after 5 time-steps Writing final snapshot at time t = 1.4925933468854715E-008 Wall clock time [hours] = 1.048E-03 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 14.73097 Maximum used memory per cpu [MBytes] = 46.066 Maximum used memory [MBytes] = 175.297 real 73.07 user 280.25 sys 2.73