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$ 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= 1.00000000 1.00000000 2.00000000 Vbox= 2.00000000 setup_slices: slice_position = p setup_slices: ix_loc,iy_loc,iz_loc,iz2_loc, (video files) = 4 4 4 35 units_general: unit_velocity= 1.0000000000000000 units_general: unit_density= 1.0000000000000000 units_general: unit_length= 1.0000000000000000 units_general: unit_magnetic= 3.5449078083038330 units_eos: unit_temperature= 4.8108880235694808E-009 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.00000000 -0.693148017 0.577350020 0.199999824 0.399999976 WARNING: initialize_diagnostics: volume averages are calculated wrongly for nonequidistant grids unless lproper_averages=T.! initialize_gravity: no x-gravity initialize_gravity: no y-gravity initialize_gravity: constant gravz= -1.00000000 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= 0.333333045 0.333332986 0.00000000E+00 initialize_energy: now set cs2top=(cool*cs2cool + cool2*cs2cool2)/(cool+cool2) initialize_energy: Calculated Fbot = 1.00000007E-02 initialize_energy: nheatc_max,iheatcond= 4 K-profile nothing nothing nothing heat conduction: K-profile viscous force: nu*(del2u+graddivu/3+2S.glnrho) 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 lnrho ( 6) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil rho ( 7) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil ee ( 48) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cp1 ( 53) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil TT ( 56) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cv ( 65) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cv1 ( 66) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil z_mn ( 88) 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_energy: max(advec_cs2) = 3050.53662 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> coriolis_cartesian: add Coriolis force; Omega= 0.100000001 duu_dt: max(advec_uu) = 0.00000000E+00 dlnrho_dt: SOLVE Bcs for lnrho, x: < p>, y: < p>, z: < a2> dlnrho_dt: max(diffus_diffrho ) = 0.00000000E+00 dlnrho_dt: max(diffus_diffrho3) = 0.00000000E+00 denergy_dt: SOLVE denergy_dt Bcs for ss, x: < p>, y: < p>, z: < c1:cT> denergy_dt: lnTT,cs2,cp1= 0.777028382 1.44999945 1.00000000 calc_heatcond: hcond0= 8.00000038E-03 calc_heatcond: lgravz= T calc_heatcond: Fbot,Ftop= 1.00000007E-02 0.00000000E+00 calc_heatcond: added thdiff calc_heat_cool: lgravz, lgravr, lgravx, lspherical_coords= T F F F cooling_profile,z2,wcool,cs2cool=gaussian 1.00000000 0.200000003 0.333332986 --it-----t-------dt------urms----umax----rhom----ssm----dtc---dtu---dtnu-dtchi- 0 0.00 7.097E-03 0.0053 0.0995 11.0770 -0.4380 0.980 0.019 0.269 0.448 10 0.07 7.097E-03 0.0047 0.0731 11.0771 -0.4385 0.980 0.015 0.269 0.448 20 0.14 7.098E-03 0.0046 0.0458 11.0771 -0.4390 0.980 0.017 0.269 0.448 30 0.21 7.098E-03 0.0051 0.0399 11.0771 -0.4396 0.980 0.016 0.269 0.448 40 0.28 7.098E-03 0.0056 0.0455 11.0770 -0.4400 0.980 0.013 0.269 0.448 Simulation finished after 41 time-steps Writing final snapshot at time t = 0.29100462200585753 Wall clock time [hours] = 3.586E-03 (+/- 2.7778E-07) Wall clock time/timestep/meshpoint [microsec] = 9.609315 Maximum used memory per cpu [MBytes] = 10.699 Maximum used memory [MBytes] = 10.699 real 75.32 user 70.89 sys 0.02