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 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$ SVN: ------- v. ( ) 20a96485a SVN: ------- v. ( ) 5d67d8a3f initialize_mpicomm: enabled MPI 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= 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: 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.2622625000000000E-002 seconds setup_slices: slice_position = m setup_slices: ix_loc,iy_loc, (video files) = 11 11 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.0000000000000000 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 WARNING: initialize_density: set lgamma_is_1=F if gamma/=1 with entropy! select_eos_variable: Using rho and ss initialize_energy: cs2top,cs2cool,cs2cool2= 1.0000000000000000 0.0000000000000000 0.0000000000000000 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 WARNING: initialize_power_spectrum: computation of wavevector wrong for non-cubical domains! 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.13000000E-07 viscous force: nu_shock*(XXXXXXXXXXX) viscous force: nu*(del6u+S.glnrho) viscous force: nu*(del2u+graddivu/3+2S.glnrho) sh: 1: Syntax error: Bad fd number 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: possible overcalculation... pencil ee ( 55) 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-005 1.0000000000000001E-005 1.0000000000000001E-005 2.5019037577498498E-004 7.4119482392732625E-003 1.0000000000000001E-005 1.0000000000000001E-005 1.0000000000000001E-005 pde: ENTER SVN: ------- v. ( ) $Id$ check_SN: ENTER check_SNI: ENTER check_SNIIb: 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.0000000000000000 ---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.1326E+03 1.8483E+06 5.6913E-04 2.1548E-03 1 0.000000010 3.0877E-09 6.7720E-02 3.0082E+00 2.33610E-02 8.22356E-01 1.8973E+02 1.1328E+03 1.8478E+06 5.6850E-04 2.1521E-03 WARNING: powerhel: computation of wavevector wrong for non-cubical domains! WARNING: powerhel: computation of wavevector wrong for non-cubical domains! explode_SN: step, time = 3 0.13088E-07 2 0.000000013 5.0606E-09 8.8615E-02 3.9369E+00 2.33610E-02 8.22356E-01 1.8970E+02 2.0377E+03 1.8476E+06 5.6831E-04 2.1513E-03 WARNING: powerhel: computation of wavevector wrong for non-cubical domains! WARNING: powerhel: computation of wavevector wrong for non-cubical domains! 3 0.000000018 2.1880E-09 1.2740E-01 5.4591E+00 2.33610E-02 8.22356E-01 1.8965E+02 2.0379E+03 1.8474E+06 5.6799E-04 2.1499E-03 WARNING: powerhel: computation of wavevector wrong for non-cubical domains! WARNING: powerhel: computation of wavevector wrong for non-cubical domains! 4 0.000000020 3.7717E-09 1.4588E-01 6.1170E+00 2.33610E-02 8.22356E-01 1.8962E+02 2.0380E+03 1.8473E+06 5.6785E-04 2.1493E-03 WARNING: powerhel: computation of wavevector wrong for non-cubical domains! WARNING: powerhel: computation of wavevector wrong for non-cubical domains! Simulation finished after 5 time-steps Writing final snapshot at time t = 2.4107948475976662E-008 Wall clock time [hours] = 6.919E-04 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 9.729686 Wall clock time/timestep/local meshpoint [microsec] = 38.91874 Rhs wall clock time/timestep/local meshpoint [microsec] = 690.3028 Maximum used memory per cpu [MBytes] = 45.410 Maximum used memory [MBytes] = 175.434 real 56.03 user 219.70 sys 2.81