SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) 8cc7b1488 SVN: ------- v. ( ) 8cc7b1488 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$ The verbose level is ip= 14 (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_general_by_label: no persistent value of dt found! WARNING: input_persist_general_by_label: no persistent value of it 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 = 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 1.4425635000000020E-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! initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. 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 gammacr1= 0.33333333333333326 initialize_interstellar: t_next_SNI = 0.35137243E-01 initialize_interstellar: eampl_SN, kampl_SN = 0.2063E+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_next_SNI, t_next_SNII= 0.35137243E-01 0.00000000E+00 viscous force: nu_shock*(XXXXXXXXXXX) viscous force: nu*(del2u+graddivu/3+2S.glnrho) viscous force: nu*(del6u+S.glnrho) farraymin 3.7083367875272092E-005 4.5858332163458895E-005 4.1952600786889477E-005 2.5019037577498498E-005 7.4119482392732625E-004 5.0000000000000004E-006 5.0000000000000004E-006 5.0000000000000004E-006 5.0000000000000004E-006 5.0000000000000004E-006 5.0000000000000004E-006 5.0000000000000004E-006 pde: ENTER SVN: ------- v. ( ) $Id$ check_SN: ENTER explode_SN: step, time = 1 0.00000E+00 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: < ism> dlnrho_dt: diffrho_shock= 1.0000000000000000 denergy_dt: SOLVE denergy_dt Bcs for ss, x: < she>, y: < p>, z: < ism> denergy_dt: lnTT,cs2,cp1= 8.7494711002314851 6307.3512691139531 0.0000000000000000 calc_heatcond_shock: use shock diffusion calc_heatcond_shock: chi_shock= 1.0000000000000000 calc_heatcond_shock: added thdiff calc_heatcond_hyper3: chi_hyper3= 8.5000000000000001E-013 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 SOLVE decr_dt Bcs for ecr, x: < she>, y: < p>, z: < s> dfcr_dt: SOLVE Bcs for Fecx, x: < she>, y: < p>, z: < s> Bcs for Fecy, x: < she>, y: < p>, z: < s> Bcs for Fecz, x: < she>, y: < p>, z: < s> shearing: Sshear,Sshear1= -25.000000000000000 -25.000000000000000 shearing: qshear,qshear0= 1.0000000000000000 0.0000000000000000 ---it--------t----------dt-------urms-------umax-------rhomin------rhomax------ekin-------ethm-------ecrm------ecrmax------TTmin------TTmax----fviscmax--mesh3Remax- 0 0.000000000 1.000E-09 1.6911E-03 5.1874E-03 2.34889E+00 8.26857E+01 4.9045E-05 2.8735E+06 3.8987E+01 3.2524E+04 4.1047E+03 3.1287E+07 8.905E+01 2.042E-05 2 0.000000002 1.000E-09 2.4092E-01 8.3609E+00 2.34889E+00 8.26857E+01 2.3101E+00 2.8713E+06 3.8987E+01 3.2524E+04 3.9885E+03 3.1287E+07 3.496E+06 3.291E-02 4 0.000000004 1.015E-09 4.8010E-01 1.6628E+01 2.34889E+00 8.26857E+01 9.1723E+00 2.8712E+06 3.8987E+01 3.2524E+04 3.8702E+03 3.1279E+07 5.375E+07 6.546E-02 6 0.000000006 1.077E-09 7.4782E-01 2.5649E+01 2.34889E+00 8.26857E+01 2.2239E+01 2.8797E+06 3.8987E+01 3.2524E+04 3.7253E+03 3.1240E+07 2.922E+08 1.010E-01 8 0.000000009 1.431E-09 9.8561E-01 3.2976E+01 2.34889E+00 8.26857E+01 3.8552E+01 2.9009E+06 3.8987E+01 3.2524E+04 3.5686E+03 3.1147E+07 7.153E+08 1.298E-01 Simulation finished after 10 time-steps Writing final snapshot at time t = 1.1490119256527313E-008 Wall clock time [hours] = 7.254E-04 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 5.100343 Wall clock time/timestep/local meshpoint [microsec] = 20.40137 Rhs wall clock time/timestep/local meshpoint [microsec] = 12.38742 Maximum used memory per cpu [MBytes] = 41.863 Maximum used memory [MBytes] = 160.867 real 2.99 user 9.98 sys 0.59