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.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 6.7932219999999998E-003 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
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.00000000E+00
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-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.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.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 wrong for non-cubical domains!
WARNING: powerhel: computation of wavevector 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 wrong for non-cubical domains!
WARNING: powerhel: computation of wavevector 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 wrong for non-cubical domains!
WARNING: powerhel: computation of wavevector 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 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 = 1.4925933468854715E-008
Wall clock time [hours] = 5.349E-04 (+/- 2.7778E-13)
Wall clock time/timestep/meshpoint [microsec] = 7.521577
Wall clock time/timestep/local meshpoint [microsec] = 30.08631
Rhs wall clock time/timestep/local meshpoint [microsec] = 462.8937
Maximum used memory per cpu [MBytes] = 44.984
Maximum used memory [MBytes] = 172.941
real 37.98
user 149.48
sys 1.26