SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
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= 32 32 32
Lx, Ly, Lz= 1.0000000000000000 1.0000000000000000 2.0000000000000000
Vbox= 2.0000000000000000
rsnap: read snapshot var.dat in 0.12916684150695801 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.5449077018110318
units_eos: unit_temperature= 9.6217761721314082E-009
units_eos: cp, lnTT0, cs0, pp0, Rgas= 2.0000000000000000 -1.3862951936203205 0.57735000000000003 0.19999982149999285 0.79999995199999807
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.0000000000000000
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.33333302250000002 0.33333299999999999 0.0000000000000000
initialize_energy: now set cs2top=(cool*cs2cool + cool2*cs2cool2)/(cool+cool2)
initialize_energy: Calculated Fbot = 1.0000000600000059E-002
initialize_energy: nheatc_max,iheatcond= 4 K-profile nothing nothing nothing
heat conduction: K-profile
cooling_profile,z2,wcool,cs2cool=gaussian 1.0000000000000000 0.20000000000000001 0.33333299999999999
viscous force: nu*(del2u+graddivu/3+2S.glnrho)
sh: 1: Syntax error: Bad fd number
pencil_consistency_check: checking pencil case
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'
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 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) = 3071.9133269640456
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.10000000000000001
duu_dt: max(advec_uu) = 5.9250845432077044E-005
dlnrho_dt: SOLVE
Bcs for lnrho, x: < p>, y: < p>, z: < a2>
dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000
dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000
denergy_dt: SOLVE denergy_dt
Bcs for ss, x: < p>, y: < p>, z: < c1:cT>
denergy_dt: lnTT,cs2,cp1= 6.8312952739877485E-003 1.3424727748122127 0.50000000000000000
calc_heatcond: hcond0= 1.6000000000000000E-002
calc_heatcond: lgravz= T
calc_heatcond: Fbot,Ftop= 1.0000000600000059E-002 0.0000000000000000
calc_heatcond: added thdiff
calc_heat_cool: lgravz, lgravr, lgravx, lspherical_coords= T F F F
cooling_profile,z2,wcool,cs2cool=gaussian 1.0000000000000000 0.20000000000000001 0.33333299999999999
--it-----------t-------------------dt----------------urms--------------umax-------------rhom------------ssbycpm------------dtc-------------dtu-------------dtnu-----------dtchi------
0 0.0000000000000E+00 6.7934342668984E-03 0.00633494973307 0.09559838043254 14.4708006727069 -0.44600582855207 0.9782842680130 0.0251658508978 0.2072676794831 0.3454470813524
10 6.7934015851104E-02 6.7933327514499E-03 0.00555072815237 0.07230049197946 14.4708248641281 -0.44635726637377 0.9782849189789 0.0192098244316 0.2072645822467 0.3454423045229
20 1.3586674889854E-01 6.7932219627419E-03 0.00529708235339 0.04707023056431 14.4708460967685 -0.44673948348581 0.9782853113495 0.0191772568469 0.2072612020833 0.3454530326315
30 2.0379898731185E-01 6.7932487784274E-03 0.00561242542619 0.04134314153769 14.4708291937098 -0.44712008079588 0.9782834660600 0.0171112925353 0.2072620202298 0.3455191423222
40 2.7173167069761E-01 6.7932232805153E-03 0.00584389886042 0.04486324704298 14.4708024985952 -0.44747109877967 0.9782834560116 0.0134086078292 0.2072612422885 0.3456218264845
Simulation finished after 41 time-steps
Writing final snapshot at time t = 0.27852489397812757
Wall clock time [hours] = 2.495E-03 (+/- 2.7778E-10)
Wall clock time/timestep/meshpoint [microsec] = 6.685716
Wall clock time/timestep/local meshpoint [microsec] = 13.37143
Rhs wall clock time/timestep/local meshpoint [microsec] = 66.37233
Maximum used memory per cpu [MBytes] = 17.617
Maximum used memory [MBytes] = 32.848
real 51.36
user 71.12
sys 5.86