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 rsnap: read snapshot var.dat in 9.9999993108212948E-004 seconds 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 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.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.00000000 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 cooling_profile,z2,wcool,cs2cool=gaussian 1.00000000 0.200000003 0.333332986 viscous force: nu*(del2u+graddivu/3+2S.glnrho) sh: 1: Syntax error: Bad fd number At line 156 of file eos_dummies.inc Fortran runtime warning: An array temporary was created for argument 'f' of procedure 'get_gamma_etc' At line 180 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 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 ( 55) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cp1 ( 60) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil TT ( 63) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cv ( 72) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cv1 ( 73) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil z_mn ( 95) 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.00000000 dlnrho_dt: SOLVE Bcs for lnrho, x: < p>, y: < p>, z: < a2> dlnrho_dt: max(diffus_diffrho ) = 0.00000000 dlnrho_dt: max(diffus_diffrho3) = 0.00000000 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.00000000 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.0064 0.0995 14.4127 -0.4525 0.980 0.019 0.269 0.448 10 0.07 7.097E-03 0.0056 0.0731 14.4128 -0.4529 0.980 0.015 0.269 0.448 20 0.14 7.098E-03 0.0054 0.0458 14.4128 -0.4533 0.980 0.017 0.269 0.448 30 0.21 7.098E-03 0.0057 0.0399 14.4128 -0.4537 0.980 0.016 0.269 0.448 40 0.28 7.098E-03 0.0059 0.0455 14.4127 -0.4541 0.980 0.013 0.269 0.448 Simulation finished after 41 time-steps Writing final snapshot at time t = 0.29100461362395436 Wall clock time [hours] = 3.604E-03 (+/- 2.7778E-07) Wall clock time/timestep/meshpoint [microsec] = 9.656208 Wall clock time/timestep/local meshpoint [microsec] = 9.656208 Rhs wall clock time/timestep/local meshpoint [microsec] = 56.47613 Maximum used memory per cpu [MBytes] = 10.871 Maximum used memory [MBytes] = 10.871 real 78.80 user 74.52 sys 0.04