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. ( ) 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 8.1611529999999766E-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 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] = 6.865E-04 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 4.826697 Wall clock time/timestep/local meshpoint [microsec] = 19.30679 Rhs wall clock time/timestep/local meshpoint [microsec] = 9.657308 Maximum used memory per cpu [MBytes] = 50.172 Maximum used memory [MBytes] = 192.934 real 3.13 user 8.12 sys 0.53