SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ 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$ SVN: ------- v. ( ) $Id$ The verbose level is ip= 15 (ldebug= F ) This is a 2-D run nxgrid, nygrid, nzgrid= 20 1 201 Lx, Ly, Lz= 2.5960000000000001 1.0000000000000000 1.0000000000000000 Vbox= 2.5960000000000001 rsnap: read snapshot var.dat in 1.0000000474974513E-003 seconds setup_slices: slice_position = p setup_slices: ix_loc,iy_loc,iz_loc,iz2_loc, (video files) = 4 4 4 204 units_general: unit_velocity= 1.0000000000000000 units_general: unit_density= 1.0000000000000000 units_general: unit_length= 1.0000000000000000 units_general: unit_magnetic= 3.5449077018110318 units_eos: unit_temperature= 4.8108883747190188E-009 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 -2.0792692765183309 0.28870000000000001 5.0008614000000007E-002 0.40000000000000008 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 initialize_gravity: no x-gravity initialize_gravity: no y-gravity initialize_gravity: constant gravz= -1.0000000000000000 diffusion: nothing (i.e. no mass diffusion) initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. select_eos_variable: Using lnrho and ss initialize_energy: cs2top,cs2cool,cs2cool2= 8.3333333333333412E-002 0.0000000000000000 0.0000000000000000 initialize_energy: now set cs2cool=cs2top initialize_energy: nheatc_max,iheatcond= 4 K-const nothing nothing nothing heat conduction: K=cte viscous force: mu/rho*(del2u+graddivu/3) 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: 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: 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_energy: max(advec_cs2) = 36721.074578012238 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> duu_dt: max(advec_uu) = 4.6642984322891115E-003 dlnrho_dt: SOLVE Bcs for lnrho, x: < p>, y: < p>, z: < a2> dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000 denergy_dt: SOLVE denergy_dt Bcs for ss, x: < p>, y: < p>, z: < cT> denergy_dt: lnTT,cs2,cp1= 0.31845373111853470 0.91666666666666663 1.0000000000000000 calc_heatcond_constK: hcond= 0.0000000000000000 calc_heatcond_constK: added thdiff ----it--------t---------dt----------urms----------umax--------uzrmaxs--------uzmax--------ssm---------ekin----------ethm--------rhom------dtc---- 0 0.000 2.982E-05 5.23651E-04 1.4604E-03 1.29749E-03 1.2527E-03 0.32706 8.2481E-07 3.3307E+00 5.99897 0.01429 100 0.003 2.982E-05 1.86336E-04 5.1876E-04 3.91339E-04 3.9134E-04 0.32706 1.1853E-07 3.3307E+00 5.99897 0.01429 200 0.006 2.982E-05 1.54259E-04 3.8661E-04 3.18620E-04 3.1862E-04 0.32706 8.1323E-08 3.3307E+00 5.99897 0.01429 300 0.009 2.982E-05 1.37293E-04 3.3508E-04 3.04707E-04 3.0471E-04 0.32706 6.4604E-08 3.3307E+00 5.99896 0.01429 400 0.012 2.982E-05 1.26004E-04 3.0450E-04 2.90043E-04 2.9004E-04 0.32706 5.4627E-08 3.3307E+00 5.99896 0.01429 500 0.015 2.982E-05 1.17663E-04 2.9474E-04 2.74442E-04 2.7444E-04 0.32706 4.7831E-08 3.3307E+00 5.99896 0.01429 600 0.018 2.982E-05 1.11100E-04 2.8889E-04 2.58670E-04 2.5867E-04 0.32706 4.2819E-08 3.3307E+00 5.99896 0.01429 700 0.021 2.982E-05 1.05715E-04 2.8561E-04 2.43478E-04 2.4348E-04 0.32706 3.8917E-08 3.3307E+00 5.99896 0.01429 800 0.024 2.982E-05 1.01160E-04 2.8411E-04 2.29127E-04 2.2913E-04 0.32706 3.5758E-08 3.3307E+00 5.99896 0.01429 900 0.027 2.982E-05 9.72177E-05 2.8242E-04 2.15672E-04 2.1567E-04 0.32706 3.3127E-08 3.3307E+00 5.99897 0.01429 Simulation finished after 1000 time-steps Writing final snapshot at time t = 2.9819166890472513E-002 Wall clock time [hours] = 4.142E-03 (+/- 2.7778E-07) Wall clock time/timestep/meshpoint [microsec] = 3.708955 Wall clock time/timestep/local meshpoint [microsec] = 3.708955 Rhs wall clock time/timestep/local meshpoint [microsec] = 3.373881 Maximum used memory per cpu [MBytes] = 17.262 Maximum used memory [MBytes] = 17.262 real 15.04 user 15.00 sys 0.03