Pencil1: iapp, nprocs, ncpus= 0 2 2 SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ initialize_mpicomm: enabled MPI 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.0000000000000000 1.0000000000000000 2.0000000000000000 Vbox= 2.0000000000000000 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.5449077018110318 units_eos: unit_temperature= 9.6217761721314082E-009 units_eos: cp, lnTT0, cs0, pp0, Rgas= 2.0000000000000000 -1.3862951936203205 0.57735000000000003 0.19999982149999285 0.79999995199999807 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= 0.33333302250000002 0.33333299999999999 0.0000000000000000E+000 initialize_energy: now set cs2top=(cool*cs2cool + cool2*cs2cool2)/(cool+cool2) initialize_energy: Calculated Fbot = 1.0000000600000059E-002 initialize_energy: nheatc_max,iheatcond= 4 K-profile nothing nothing nothing heat conduction: K-profile cooling_profile,z2,wcool,cs2cool=gaussian 1.0000000000000000 0.20000000000000001 0.33333299999999999 viscous force: nu*(del2u+graddivu/3+2S.glnrho) At line 152 of file eos_dummies.inc Fortran runtime warning: An array temporary was created for argument 'f' of procedure 'get_gamma_etc' At line 176 of file eos_dummies.inc Fortran runtime warning: An array temporary was created for argument 'f' of procedure 'eoscalc_point_f' pencil_consistency_check: checking pencil case At line 152 of file eos_dummies.inc Fortran runtime warning: An array temporary was created for argument 'f' of procedure 'get_gamma_etc' At line 176 of file eos_dummies.inc Fortran runtime warning: An array temporary was created for argument 'f' of procedure 'eoscalc_point_f' 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 ( 12) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil rho ( 13) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil ee ( 54) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cp1 ( 59) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil TT ( 62) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cv ( 71) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cv1 ( 72) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil z_mn ( 94) 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) = 3071.9133269640456 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.10000000000000001 duu_dt: max(advec_uu) = 5.9250845432077044E-005 dlnrho_dt: SOLVE Bcs for lnrho, x: < p>, y: < p>, z: < a2> dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000E+000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000E+000 denergy_dt: SOLVE denergy_dt Bcs for ss, x: < p>, y: < p>, z: < c1:cT> denergy_dt: lnTT,cs2,cp1= 6.8312952739877485E-003 1.3424727748122127 0.50000000000000000 calc_heatcond: hcond0= 1.6000000000000000E-002 calc_heatcond: lgravz= T calc_heatcond: Fbot,Ftop= 1.0000000600000059E-002 0.0000000000000000E+000 calc_heatcond: added thdiff calc_heat_cool: lgravz, lgravr, lgravx, lspherical_coords= T F F F cooling_profile,z2,wcool,cs2cool=gaussian 1.0000000000000000 0.20000000000000001 0.33333299999999999 --it-----------t-------------------dt----------------urms--------------umax-------------rhom------------ssbycpm------------dtc-------------dtu-------------dtnu-----------dtchi------ 0 0.0000000000000E+00 6.7934342668984E-03 0.00633494973307 0.09559838043254 14.4708006727069 -0.44600582855207 0.9782842680130 0.0251658508978 0.2072676794831 0.3454470813524 10 6.7934015851104E-02 6.7933327514499E-03 0.00555072815237 0.07230049197946 14.4708248641281 -0.44635726637377 0.9782849189789 0.0192098244316 0.2072645822467 0.3454423045229 20 1.3586674889854E-01 6.7932219627419E-03 0.00529708235339 0.04707023056431 14.4708460967685 -0.44673948348581 0.9782853113495 0.0191772568469 0.2072612020833 0.3454530326315 30 2.0379898731185E-01 6.7932487784274E-03 0.00561242542619 0.04134314153769 14.4708291937098 -0.44712008079588 0.9782834660600 0.0171112925353 0.2072620202298 0.3455191423222 40 2.7173167069761E-01 6.7932232805153E-03 0.00584389886042 0.04486324704298 14.4708024985952 -0.44747109877967 0.9782834560116 0.0134086078292 0.2072612422885 0.3456218264845 Simulation finished after 41 time-steps Writing final snapshot at time t = 0.27852489397812757 Wall clock time [hours] = 2.392E-03 (+/- 2.7778E-10) Wall clock time/timestep/meshpoint [microsec] = 6.408755 Maximum used memory per cpu [MBytes] = 17.629 Maximum used memory [MBytes] = 32.895 real 48.95 user 69.71 sys 3.79