===================================================================== WARNING: using SINGLE PRECISION, you'd better know what you're doing! ===================================================================== SVN: ------- v. ( ) $Id$ WARNING: find_namelist: namelist "_pars" is missing!! SVN: ------- v. ( ) $Id$ WARNING: find_namelist: namelist "_pars" is missing!! SVN: ------- v. ( ) $Id$ initialize_mpicomm: enabled MPI WARNING: run: quad precision suppressed! 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= 14 (ldebug= F ) This is a 3-D run nxgrid, nygrid, nzgrid= 32 32 32 Lx, Ly, Lz= 3.30882072 3.30882072 3.20542002 Vbox= 35.0938835 rsnap: read snapshot var.dat in 5.1414099999999997E-004 seconds units_general: unit_velocity= 1000000.0000000000 units_general: unit_density= 1.0000000000000001E-007 units_general: unit_length= 100000000.00000000 units_general: unit_magnetic= 3.5449078083038330 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 initialize_eos: ENTER initialize_eos: chiH= 2.1789600494381816E-040 13.600000381469725 initialize_eos: chiH_= 1.2080410653356137E-041 0.75400000810623169 initialize_gravity: no x-gravity initialize_gravity: no y-gravity initialize_gravity: constant gravz= -2.75000000 diffusion: nothing (i.e. no mass diffusion) initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. heat conduction: nothing heat conduction: nothing heat conduction: nothing WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_temperature_ionization. The values provided are for one-atomic ideal gas. Use at own risk! resistivity: hyper3 resistivity: shock calc_angle_weights: correction_factor = 1.00000000 initialize_radiation: ndir = 6 WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_temperature_ionization. The values provided are for one-atomic ideal gas. Use at own risk! initialize_shock: prenormalised shock_factor sum= 1.00000000 viscous force: nu*(del6u+S.glnrho) viscous force: nu_shock*(XXXXXXXXXXX) WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_temperature_ionization. The values provided are for one-atomic ideal gas. Use at own risk! 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 aa ( 154) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil b2 ( 164) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil va2 ( 180) 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 pde: ENTER SVN: ------- v. ( ) $Id$ calc_pencils_hydro: max(advec_uu) = 4.48127270 calc_pencils_energy: max(p%advec_cs2) = 623.263977 WARNING: calc_pencils_viscosity: viscous heating not implemented for lvisc_hyper3_nu_const! calc_pencils_magnetic_pencpar: B_ext = 0.00000000 0.00000000 0.892060995 calc_pencils_magnetic_pencpar: logic = F duu_dt: SOLVE Bcs for ux, x: < p>, y: < p>, z: < s> Bcs for uy, x: < p>, y: < p>, z: < s> Bcs for uz, x: < p>, y: < p>, z: < a> dlnrho_dt: SOLVE Bcs for lnrho, x: < p>, y: < p>, z: < StS> dlnrho_dt: max(diffus_diffrho ) = 0.00000000 dlnrho_dt: max(diffus_diffrho3) = 0.00000000 denergy_dt: SOLVE denergy_dt Bcs for lnTT, x: < p>, y: < p>, z: < StS> denergy_dt: cs2 = 2.22124910 daa_dt: SOLVE Bcs for Ax, x: < p>, y: < p>, z: < s:pwd> Bcs for Ay, x: < p>, y: < p>, z: < s:nil> Bcs for Az, x: < p>, y: < p>, z: < a:nil> daa_dt: iresistivity=hyper3 shock daa_dt: max(diffus_eta) = 45.3062172 daa_dt: max(diffus_eta2) = 0.00000000 daa_dt: max(diffus_eta3) = 0.490921199 ---it------t--------dt--------urms------umax-----vArms------vAmax----rhom---rhomin---rhomax-----ssm-------TTm-----TTmin----TTmax---Qradrms---Qradmax---dtb---dtc---dtu-- 0 0.0000 2.781E-03 0.17130 0.49030 0.01345 0.05838 62.57 2.335E-03 266.23 0.0016598 12726.6 4888.7 21398.6 1.847E-01 9.619E-01 0.00 0.08 0.02 10 0.0606 9.368E-03 0.07887 0.22699 0.01345 0.06360 62.57 1.934E-03 286.06 0.0016600 12733.9 4669.3 21780.2 2.422E-01 4.414E+00 0.01 0.28 0.04 20 0.1506 9.315E-03 0.05976 0.16652 0.01340 0.06677 62.56 1.914E-03 278.13 0.0016599 12729.7 4564.8 21589.9 2.107E-01 1.745E+00 0.01 0.28 0.03 30 0.2451 9.506E-03 0.05109 0.13362 0.01330 0.06459 62.56 2.094E-03 274.24 0.0016596 12725.6 4480.2 21554.0 2.120E-01 1.611E+00 0.01 0.28 0.02 40 0.3481 1.058E-02 0.04607 0.12537 0.01314 0.06292 62.56 2.213E-03 271.61 0.0016593 12720.3 4301.2 21495.1 2.178E-01 2.977E+00 0.01 0.32 0.02 50 0.4461 9.681E-03 0.04327 0.12663 0.01298 0.06184 62.56 2.281E-03 271.72 0.0016591 12714.3 4063.0 21491.2 2.284E-01 1.105E+01 0.01 0.29 0.02 60 0.5376 8.960E-03 0.04173 0.13332 0.01283 0.06138 62.55 2.286E-03 270.89 0.0016589 12707.9 4048.4 21475.8 2.248E-01 7.573E+00 0.01 0.27 0.02 70 0.6252 8.594E-03 0.04133 0.14867 0.01270 0.05976 62.56 2.155E-03 270.68 0.0016589 12701.1 4063.7 21479.9 2.273E-01 6.031E+00 0.01 0.26 0.02 80 0.7031 7.523E-03 0.04193 0.16742 0.01261 0.05972 62.56 2.127E-03 269.52 0.0016589 12694.7 4080.0 21467.8 2.344E-01 9.031E+00 0.00 0.22 0.02 90 0.7758 7.475E-03 0.04337 0.18172 0.01257 0.05963 62.56 2.196E-03 269.58 0.0016590 12688.5 4114.8 21458.6 2.315E-01 4.106E+00 0.00 0.22 0.02 100 0.8487 7.515E-03 0.04569 0.18893 0.01258 0.05985 62.56 2.285E-03 269.83 0.0016593 12682.2 4148.9 21463.0 2.362E-01 4.821E+00 0.00 0.22 0.02 110 0.9259 7.744E-03 0.04914 0.19275 0.01269 0.06124 62.57 2.242E-03 270.18 0.0016596 12675.3 4109.2 21466.1 2.463E-01 7.100E+00 0.01 0.23 0.02 120 1.0030 7.640E-03 0.05332 0.20210 0.01292 0.06352 62.57 2.195E-03 269.45 0.0016602 12667.9 4157.7 21450.9 2.464E-01 4.567E+00 0.01 0.23 0.02 130 1.0781 7.704E-03 0.05754 0.21797 0.01328 0.06521 62.57 1.967E-03 270.37 0.0016608 12659.6 4132.3 21469.5 2.484E-01 2.284E+00 0.01 0.23 0.03 140 1.1544 7.610E-03 0.06126 0.22476 0.01381 0.06790 62.58 1.744E-03 269.71 0.0016615 12649.9 4063.2 21459.9 2.530E-01 2.497E+00 0.01 0.23 0.03 Simulation finished after 150 time-steps Writing final snapshot at time t = 1.2259751888341270 Wall clock time [hours] = 4.921E-03 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 3.604048 Wall clock time/timestep/local meshpoint [microsec] = 14.41619 Rhs wall clock time/timestep/local meshpoint [microsec] = 21.98326 Maximum used memory per cpu [MBytes] = 33.957 Maximum used memory [MBytes] = 134.672 real 44.67 user 176.26 sys 0.96