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 2-D run
nxgrid, nygrid, nzgrid= 64 1 64
Lx, Ly, Lz= 4.00000000 0.00000000 3.20000005
Vbox= 12.8000002
rsnap: read snapshot var.dat in 0.0000000000000000 seconds
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 2.95406263E-02 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)
Note: rho_const or lnrho_const are not used
initialize_density: no need to read initial stratification for lanti_shockdiffusion=F.
select_eos_variable: Using lnrho and ss
initialize_energy: hcond0 given by mixinglength_flux= 3.99999972E-03 9.99999978E-03
initialize_energy: cs2top,cs2cool,cs2cool2= 0.333333045 0.00000000 0.00000000
initialize_energy: now set cs2cool=cs2top
initialize_energy: nheatc_max,iheatcond= 4 K-profile nothing nothing nothing
heat conduction: K-profile
cooling_profile,z2,wcool,cs2cool=cubic_step 1.00000000 0.100000001 0.333333045
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.
| Put lenforce_maux_check=F in run_pars to overide this.!
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: 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$
calc_pencils_energy: max(advec_cs2) = 1321.81567
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>
duu_dt: max(advec_uu) = 1.28062749E-02
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= 1.12515295 2.05379200 1.00000000
calc_heatcond: hcond0= 3.99999972E-03
calc_heatcond: lgravz= T
calc_heatcond: Fbot,Ftop= 9.99999978E-03 0.00000000
calc_heatcond: added thdiff
calc_heat_cool: lgravz, lgravr, lgravx, lspherical_coords= T F F F
cooling_profile,z2,wcool,cs2cool=cubic_step 1.00000000 0.100000001 0.333333045
---it------t-------dt------urms----rhom----ssm----dtc---dtnu-dtchi-
0 0.00 6.841E-03 0.0005 1.8457 -0.1250 0.995 0.088 0.418
10 0.07 6.842E-03 0.0038 1.8457 -0.1278 0.995 0.088 0.419
20 0.14 6.843E-03 0.0069 1.8457 -0.1286 0.995 0.088 0.419
30 0.21 6.844E-03 0.0084 1.8457 -0.1289 0.995 0.088 0.419
40 0.27 6.844E-03 0.0094 1.8456 -0.1291 0.995 0.088 0.419
Simulation finished after 41 time-steps
Writing final snapshot at time t = 0.28055941185448319
Wall clock time [hours] = 8.681E-05 (+/- 2.7778E-07)
Wall clock time/timestep/meshpoint [microsec] = 1.860828
Wall clock time/timestep/local meshpoint [microsec] = 1.860828
Rhs wall clock time/timestep/local meshpoint [microsec] = 2.232993
Maximum used memory per cpu [MBytes] = 14.875
Maximum used memory [MBytes] = 14.875
real 0.40
user 0.39
sys 0.00