Pencil1: iapp, nprocs, ncpus= 0 4 4 SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) no cvsid is given i 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= 64 64 32 Lx, Ly, Lz= 0.69999999999999996 2.6179910000000000 1.5708000000000000 Vbox= 2.4286526417797338 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= 4.8108880860657041E-009 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 0.40546520810816927 1.0000000000000000 0.60000002400000085 0.39999997599999904 initialize_gravity: kepler x-grav, gravx= 3.0000000000000000 initialize_gravity: no y-gravity initialize_gravity: no z-gravity diffusion: nothing (i.e. no mass diffusion) initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. At line 3060 of file eos_idealgas.f90 Fortran runtime warning: An array temporary was created for argument 'f' of procedure 'getlnrho_2dyz' select_eos_variable: Using lnrho and ss initialize_energy: nheatc_max,iheatcond= 4 K-profile nothing nothing nothing heat conduction: K-profile resistivity: constant eta viscous force: nu*(del2u+graddivu/3+2S.glnrho) farray_check_maux: naux, maux = 0 4 WARNING: farray_check_maux: | The count of MAUXs is inconsistent with the actual number of | auxiliaries required. Please check your src/cparam.local. It is | also possible that someone unintentionally added *optional* | auxiliaries to MAUX CONTRIBUTION to one of the physics modules. | Please urge them to move it to src/cparam.local under their | specific application.! At line 3060 of file eos_idealgas.f90 Fortran runtime warning: An array temporary was created for argument 'f' of procedure 'getlnrho_2dyz' bc_ss_flux_x: FbotKbot= 6.1224589231721334 bc_ss_flux_x: FbotKbot= 6.1224589231721334 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$ bc_ss_flux_x: FbotKbot= 6.1224589231721334 calc_pencils_hydro: upwinding advection term calc_pencils_hydro: call gij_etc calc_pencils_energy: max(advec_cs2) = 14299.291360066360 duu_dt: SOLVE Bcs for ux, x: < a>, y: < s>, z: < p> Bcs for uy, x: < sfr>, y: < a>, z: < p> Bcs for uz, x: < sfr>, y: < sfr>, z: < p> coriolis_spherical: Omega= 0.44735680562864799 duu_dt: max(advec_uu) = 7.3049775177719866E-003 dlnrho_dt: SOLVE Bcs for lnrho, x: < a2>, y: < s>, z: < p> 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: < c1:cT>, y: < s>, z: < p> denergy_dt: lnTT,cs2,cp1= 0.41258263587703281 1.0071428168571432 1.0000000000000000 calc_heatcond: hcond0= 0.0000000000000000E+000 calc_heatcond: lgravz= F calc_headcond: "turbulent" entropy diffusion: chi_t= 2.5000000000000001E-004 calc_heatcond: added thdiff calc_heat_cool: lgravz, lgravr, lgravx, lspherical_coords= F F T T daa_dt: SOLVE Bcs for Ax, x: < s:a>, y: < a>, z: < p> Bcs for Ay, x: < a:nfr>, y: < s>, z: < p> Bcs for Az, x: < a:nfr>, y: < a>, z: < p> daa_dt: iresistivity=eta-const daa_dt: use upwinding in advection term daa_dt: max(diffus_eta) = 7.0989392570401577 daa_dt: max(diffus_eta2) = 0.0000000000000000E+000 daa_dt: max(diffus_eta3) = 0.0000000000000000E+000 ---it-------t---------dt-------urms------ux2m------uy2m------uz2m------orms-----umax------rhom------ssm-------eem-------ethm------mass------brms-----bx2m-----by2m-----bz2m-----urmsx-----brmsx-- 0 0.000 2.458E-03 4.98E-05 8.15E-10 8.29E-10 8.32E-10 4.23E-03 8.5E-05 0.100658 2.108304 1.250452 0.058424 0.244463 4.86E-04 2.3E-08 1.1E-07 1.0E-07 5.05E-05 9.3E-04 1 0.002 2.456E-03 2.41E-03 5.78E-06 7.76E-10 7.80E-10 4.08E-03 1.0E-02 0.100658 2.108149 1.250335 0.058424 0.244463 4.68E-04 2.2E-08 1.0E-07 9.7E-08 3.07E-04 8.9E-04 2 0.005 2.455E-03 4.76E-03 2.26E-05 7.22E-10 7.94E-10 3.95E-03 2.0E-02 0.100658 2.107981 1.250218 0.058423 0.244463 4.51E-04 2.0E-08 9.3E-08 9.0E-08 6.08E-04 8.6E-04 3 0.007 2.453E-03 7.03E-03 4.94E-05 6.72E-10 1.04E-09 3.88E-03 3.0E-02 0.100658 2.107801 1.250100 0.058422 0.244463 4.35E-04 1.9E-08 8.7E-08 8.3E-08 9.10E-04 8.3E-04 4 0.010 2.451E-03 9.21E-03 8.49E-05 6.26E-10 1.78E-09 3.90E-03 3.8E-02 0.100658 2.107609 1.249981 0.058422 0.244463 4.20E-04 1.8E-08 8.1E-08 7.8E-08 1.21E-03 8.0E-04 5 0.012 2.446E-03 1.13E-02 1.28E-04 5.86E-10 3.36E-09 4.08E-03 4.6E-02 0.100658 2.107409 1.249861 0.058421 0.244463 4.05E-04 1.7E-08 7.5E-08 7.2E-08 1.51E-03 7.7E-04 6 0.015 2.441E-03 1.33E-02 1.76E-04 5.51E-10 6.15E-09 4.44E-03 5.3E-02 0.100658 2.107202 1.249740 0.058421 0.244463 3.92E-04 1.5E-08 7.0E-08 6.8E-08 1.81E-03 7.5E-04 7 0.017 2.435E-03 1.51E-02 2.28E-04 5.19E-10 1.06E-08 4.98E-03 5.9E-02 0.100658 2.106990 1.249620 0.058421 0.244463 3.79E-04 1.5E-08 6.6E-08 6.3E-08 2.11E-03 7.2E-04 8 0.020 2.428E-03 1.68E-02 2.84E-04 4.90E-10 1.71E-08 5.68E-03 6.5E-02 0.100658 2.106776 1.249500 0.058420 0.244463 3.66E-04 1.4E-08 6.1E-08 5.9E-08 2.41E-03 7.0E-04 9 0.022 2.421E-03 1.85E-02 3.41E-04 4.64E-10 2.62E-08 6.51E-03 7.0E-02 0.100658 2.106560 1.249381 0.058420 0.244463 3.55E-04 1.3E-08 5.7E-08 5.5E-08 2.70E-03 6.8E-04 10 0.024 2.412E-03 2.00E-02 4.00E-04 4.41E-10 3.84E-08 7.44E-03 7.4E-02 0.100658 2.106346 1.249265 0.058420 0.244463 3.44E-04 1.2E-08 5.4E-08 5.2E-08 3.00E-03 6.6E-04 11 0.027 2.404E-03 2.14E-02 4.59E-04 4.21E-10 5.39E-08 8.45E-03 7.8E-02 0.100658 2.106135 1.249150 0.058420 0.244463 3.33E-04 1.2E-08 5.1E-08 4.9E-08 3.29E-03 6.4E-04 12 0.029 2.396E-03 2.27E-02 5.17E-04 4.03E-10 7.33E-08 9.51E-03 8.1E-02 0.100658 2.105928 1.249038 0.058420 0.244463 3.23E-04 1.1E-08 4.7E-08 4.6E-08 3.58E-03 6.2E-04 13 0.032 2.388E-03 2.40E-02 5.76E-04 3.87E-10 9.71E-08 1.06E-02 8.4E-02 0.100658 2.105728 1.248930 0.058421 0.244463 3.14E-04 1.0E-08 4.5E-08 4.3E-08 3.87E-03 6.0E-04 Simulation finished after 14 time-steps Writing final snapshot at time t = 3.4044819150705792E-002 Wall clock time [hours] = 6.911E-03 (+/- 2.7778E-10) Wall clock time/timestep/meshpoint [microsec] = 13.55908 Maximum used memory per cpu [MBytes] = 30.078 Maximum used memory [MBytes] = 118.766 real 39.96 user 58.89 sys 2.45