===================================================================== WARNING: using SINGLE PRECISION, you'd better know what you're doing! ===================================================================== SVN: ------- v. ( ) $Id$ SVN: start.in,v v. 1.1 (2010/10/09 19:44:40) brandenb SVN: run.in,v v. 1.3 (2010/10/20 19:46:23) brandenb initialize_mpicomm: enabled MPI 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= 14 (ldebug= F ) This is a 3-D run nxgrid, nygrid, nzgrid= 32 32 32 Lx, Ly, Lz= 6.28318548 6.28318548 6.28318548 Vbox= 248.050232 input_persist_forcing: 0.00000000 0.00000000 0.00000000 input_persist_forcing: -10.0000000 rsnap: read snapshot var.dat in 1.1126639000000001E-002 seconds setup_slices: slice_position = p setup_slices: ix_loc,iy_loc,iz_loc, (video files) = 4 4 4 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= 1.2027220775801925E-008 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.00000000 0.00000000 1.00000000 1.00000000 1.00000000 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 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 cs2 resistivity: constant eta Number of passive scalars = 1 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 sij2 ( 114) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil uxb ( 189) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil diva ( 221) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil gLam ( 249) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cc ( 260) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil gLamRA ( 285) 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_hydro: upwinding advection term calc_pencils_hydro: max(advec_uu) = 0.00000000 calc_pencils_energy: max(advec_cs2) = 77.8146667 duu_dt: SOLVE Bcs for ux, x: < p>, y: < p>, z: < p> Bcs for uy, x: < p>, y: < p>, z: < p> Bcs for uz, x: < p>, y: < p>, z: < p> dlnrho_dt: SOLVE Bcs for lnrho, x: < p>, y: < p>, z: < p> dlnrho_dt: max(diffus_diffrho ) = 0.00000000 dlnrho_dt: max(diffus_diffrho3) = 0.00000000 daa_dt: SOLVE Bcs for Ax, x: < p>, y: < p>, z: < p> Bcs for Ay, x: < p>, y: < p>, z: < p> Bcs for Az, x: < p>, y: < p>, z: < p> daa_dt: iresistivity=eta-const daa_dt: use upwinding in advection term daa_dt: max(diffus_eta) = 0.778146625 daa_dt: max(diffus_eta2) = 0.00000000 daa_dt: max(diffus_eta3) = 0.00000000 SOLVE dlncc_dt Bcs for cc1, x: < p>, y: < p>, z: < p> dlncc_dt: pscalar_diff= 9.99999978E-03 dspecial_dt: SOLVE dSPECIAL_dt Bcs for LamRA, x: < p>, y: < p>, z: < p> ---it--------t---------dt-------urms------brms-----abrms------abm-------ubm-------oum-------epsK------epsM-----apbrms-----ucm---- 0 0.000 8.73E-02 0.000E+00 1.352E-03 1.35E-07 -1.57E-09 0.00E+00 0.00E+00 0.000E+00 8.087E-06 1.35E-07 0.00E+00 10 0.853 8.34E-02 2.399E-02 3.252E-04 9.44E-09 -1.15E-10 -4.81E-09 1.90E-03 6.263E-05 1.788E-07 9.43E-09 -8.55E-05 20 1.679 8.21E-02 3.104E-02 1.521E-04 2.58E-09 5.41E-12 1.75E-08 3.08E-03 9.852E-05 1.897E-08 2.57E-09 -2.63E-04 30 2.495 8.09E-02 3.597E-02 9.426E-05 1.18E-09 1.17E-11 3.13E-08 4.04E-03 1.268E-04 4.568E-09 1.18E-09 -5.06E-04 40 3.296 7.99E-02 4.050E-02 6.736E-05 6.90E-10 7.60E-12 3.99E-08 5.16E-03 1.627E-04 1.708E-09 6.85E-10 -8.04E-04 Simulation finished after 41 time-steps Writing final snapshot at time t = 3.3762768153101206 Wall clock time [hours] = 7.901E-04 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 2.117043 Wall clock time/timestep/local meshpoint [microsec] = 8.468174 Rhs wall clock time/timestep/local meshpoint [microsec] = 50.42670 Maximum used memory per cpu [MBytes] = 33.625 Maximum used memory [MBytes] = 132.137 real 25.23 user 97.86 sys 0.89