SVN: ------- v. ( ) $Id$
WARNING: find_namelist: namelist "special_init_pars" is missing!!
SVN: start.in,v v. 1.1 (2009-09-12 16:05:39) wlyra
SVN: run.in,v v. 1.2 (2009-09-24 08:59:37) wlyra
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 16
Lx, Ly, Lz= 0.40000000000000002 1.6000000000000001 0.20000000000000001
Vbox= 0.12800000000000003
rsnap: read snapshot var.dat in 1.6261199999999999E-004 seconds
setup_slices: slice_position = p
setup_slices: ix_loc,iy_loc,iz_loc,iz2_loc, (video files) = 4 4 4 19
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= 3.4363488390850127E-009
units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 -3.6888794541139358 0.10000000000000001 7.1428571428571461E-003 0.28571428571428570
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
diffusion: mesh hyperdiffusion
diffusion: shock diffusion
WARNING: initialize_density: For diffusion energy/momentum correction should use lmassdiff_fix=T!
initialize_density: no need to read initial stratification for lanti_shockdiffusion=F.
select_eos_variable: Using rho and ss
initialize_energy: nheatc_max,iheatcond= 4 hyper3-mesh shock nothing nothing
heat conduction: hyperdiffusivity of ss
heat conduction: shock
viscous force: nu_hyper3_mesh/pi^5 *(Deltav)^6/Deltaq
viscous force: nu_shock*(XXXXXXXXXXX)
sh: 1: Syntax error: Bad fd number
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 rho ( 13) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil grho ( 16) 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 cv1 ( 73) 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) = 139.18980252797598
WARNING: calc_pencils_viscosity: viscous heating not implemented for lvisc_hyper3_mesh!
duu_dt: SOLVE
Bcs for ux, x: < she>, y: < p>, z: < p>
Bcs for uy, x: < she>, y: < p>, z: < p>
Bcs for uz, x: < she>, y: < p>, z: < p>
coriolis_cartesian: add Coriolis force; Omega= 1.0000000000000000
duu_dt: max(advec_uu) = 0.0000000000000000
dlnrho_dt: SOLVE
Bcs for lnrho, x: < she>, y: < p>, z: < p>
dlnrho_dt: diffrho_shock= 1.0000000000000000
dlnrho_dt: diffrho_hyper3_mesh= 5.0000000000000000
dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000
dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000
denergy_dt: SOLVE denergy_dt
Bcs for ss, x: < she>, y: < p>, z: < p>
denergy_dt: lnTT,cs2,cp1= -3.6946269386926391 9.9426900061287231E-003 1.1728390000000000
calc_heatcond_shock: use shock diffusion
calc_heatcond_shock: chi_shock= 1.0000000000000000
calc_heatcond_shock: added thdiff
calc_heatcond_hyper3_mesh: chi_hyper3= 0.0000000000000000
calc_heatcond_hyper3: added thdiff
calc_heat_cool: lgravz, lgravr, lgravx, lspherical_coords= F F F F
dspecial_dt: SOLVE dSPECIAL_dt
shearing: Sshear,Sshear1= -1.5000000000000000 -1.5000000000000000
shearing: qshear,qshear0= 1.5000000000000000 0.0000000000000000
---it-------t--------dt-------uxm-------uym-------rhom-------rho2m---rhomin---rhomax-------ppm-------ethm--------csm--------TTm-----TTmax--TTmin--
0 0.000 1.921E-02 0.00E+00 0.00E+00 9.983E-01 9.991E-01 0.73 1.189E+00 7.143E-03 1.786E-02 1.002E-01 2.511E-02 0.034 0.0210
10 0.192 1.920E-02 8.15E-06 -5.18E-07 9.983E-01 9.989E-01 0.73 1.186E+00 7.143E-03 1.786E-02 1.002E-01 2.510E-02 0.034 0.0212
20 0.384 1.919E-02 1.20E-05 -1.42E-06 9.983E-01 9.986E-01 0.75 1.181E+00 7.144E-03 1.786E-02 1.002E-01 2.510E-02 0.033 0.0213
30 0.576 1.920E-02 1.70E-05 -2.49E-06 9.983E-01 9.984E-01 0.78 1.174E+00 7.144E-03 1.786E-02 1.002E-01 2.509E-02 0.032 0.0215
40 0.768 1.921E-02 2.13E-05 -3.87E-06 9.983E-01 9.982E-01 0.81 1.164E+00 7.144E-03 1.786E-02 1.002E-01 2.509E-02 0.030 0.0216
50 0.960 1.920E-02 2.43E-05 -5.59E-06 9.983E-01 9.981E-01 0.82 1.151E+00 7.144E-03 1.786E-02 1.002E-01 2.509E-02 0.031 0.0217
60 1.152 1.920E-02 2.75E-05 -7.61E-06 9.983E-01 9.980E-01 0.82 1.142E+00 7.144E-03 1.786E-02 1.002E-01 2.508E-02 0.031 0.0219
70 1.344 1.924E-02 2.94E-05 -9.98E-06 9.983E-01 9.978E-01 0.82 1.130E+00 7.144E-03 1.786E-02 1.002E-01 2.508E-02 0.030 0.0221
80 1.537 1.926E-02 3.05E-05 -1.25E-05 9.983E-01 9.977E-01 0.83 1.123E+00 7.144E-03 1.786E-02 1.002E-01 2.508E-02 0.030 0.0222
90 1.729 1.925E-02 3.10E-05 -1.51E-05 9.983E-01 9.977E-01 0.84 1.126E+00 7.144E-03 1.786E-02 1.001E-01 2.507E-02 0.029 0.0222
Simulation finished after 91 time-steps
Writing final snapshot at time t = 1.7483417219707393
Wall clock time [hours] = 7.400E-04 (+/- 2.7778E-13)
Wall clock time/timestep/meshpoint [microsec] = 1.786845
Wall clock time/timestep/local meshpoint [microsec] = 7.147378
Rhs wall clock time/timestep/local meshpoint [microsec] = 14.77104
Maximum used memory per cpu [MBytes] = 31.055
Maximum used memory [MBytes] = 123.340
real 7.94
user 30.27
sys 0.43