SVN: ------- v. ( ) $Id$
SVN: start.in,v v. 1.1 (2010-05-19 14:36:03) nilshau
Read inlet data from file!
Read inlet data from file!
Read inlet data from file!
Read inlet data from file!
SVN: run.in,v v. 1.1 (2010-05-19 14:36:03) nilshau
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$
Write pc_constants.pro in chemistry.f90
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= 512 16 1
Lx, Ly, Lz= 0.50000000000000000 0.50000000000000000 0.10000000000000001
Vbox= 0.25000000000000000
rsnap: read snapshot var.dat in 8.8693099999999998E-004 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.5449077018110318
units_eos: chem.imp is found! Now cp, cv, gamma, mu are pencils ONLY!
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 lnTT
heat conduction: nothing
heat conduction: nothing
heat conduction: nothing
tran.in/dat file with transport data is found.
the following species are found in tran.in/dat:
CH2O ind_glob= 20 ind_chem= 15
CH3 ind_glob= 18 ind_chem= 13
CH3O ind_glob= 21 ind_chem= 16
CH4 ind_glob= 14 ind_chem= 9
CO ind_glob= 15 ind_chem= 10
CO2 ind_glob= 16 ind_chem= 11
H ind_glob= 8 ind_chem= 3
H2 ind_glob= 7 ind_chem= 2
H2O ind_glob= 12 ind_chem= 7
H2O2 ind_glob= 17 ind_chem= 12
HCO ind_glob= 19 ind_chem= 14
HO2 ind_glob= 13 ind_chem= 8
N2 ind_glob= 6 ind_chem= 1
O ind_glob= 9 ind_chem= 4
O2 ind_glob= 10 ind_chem= 5
OH ind_glob= 11 ind_chem= 6
lewis.dat file with transport data is found.
Species diffusion coefficients calculated using constant Lewis numbers.
N2 ind_glob= 6 Lewis= 1.0000000000000000
H2 ind_glob= 7 Lewis= 0.29999999999999999
H ind_glob= 8 Lewis= 0.17999999999999999
O ind_glob= 9 Lewis= 0.69999999999999996
O2 ind_glob= 10 Lewis= 1.1100000000000001
OH ind_glob= 11 Lewis= 0.72999999999999998
H2O ind_glob= 12 Lewis= 0.82999999999999996
HO2 ind_glob= 13 Lewis= 1.1000000000000001
CH4 ind_glob= 14 Lewis= 0.96999999999999997
CO ind_glob= 15 Lewis= 1.1000000000000001
CO2 ind_glob= 16 Lewis= 1.3899999999999999
H2O2 ind_glob= 17 Lewis= 1.1200000000000001
CH3 ind_glob= 18 Lewis= 1.0000000000000000
HCO ind_glob= 19 Lewis= 1.2700000000000000
CH2O ind_glob= 20 Lewis= 1.2800000000000000
CH3O ind_glob= 21 Lewis= 1.3000000000000000
Number of reactions= 35
Number of species= 16
Reinitializing chemistry.
Reading initial conditions in file CH4_p01_0phi0_7000tu0300
initchem: nothing
viscous force: nu is calculated for a mixture
WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_chemistry.
The values provided are for one-atomic ideal gas. Use at own risk!
sh: 1: Syntax error: Bad fd number
pencil_consistency_check: checking pencil case
calc_for_chem_mixture: writing mix_quant.out file
get_reaction_rate: writing react.out file
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 cp1 ( 7) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil gXXk ( 9) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil lambda1 ( 19) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil g2XXk ( 21) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil S0_R ( 26) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil cs2 ( 27) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil glnpp ( 28) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil del2pp ( 29) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil glnrho2 ( 47) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil del2lnrho ( 48) 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$
duu_dt: SOLVE
Bcs for ux, x: < set:a2>, y: < p>, z: < p>
Bcs for uy, x: < set:a2>, y: < p>, z: < p>
Bcs for uz, x: < set:a2>, y: < p>, z: < p>
dlnrho_dt: SOLVE
Bcs for lnrho, x: < a2>, y: < p>, z: < p>
denergy_dt: SOLVE denergy_dt
Bcs for lnTT, x: < set:a2>, y: < p>, z: < p>
denergy_dt: cs2 = 1240076546.2489269
dchemistry_dt: SOLVE dchemistry_dt
WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_chemistry.
The values provided are for one-atomic ideal gas. Use at own risk!
inlet_profile: uniform
inlet_YY_profile: uniform,
WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_chemistry.
The values provided are for one-atomic ideal gas. Use at own risk!
WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_chemistry.
The values provided are for one-atomic ideal gas. Use at own risk!
---it-------t-------dt-------uxm----TTm-----Y2m-------Y3m-------Y6m-------Y8m-------Y9m-------Y10m------Y11m------Y13m---
0 0.00E+00 1.00E-08 1.23E+02 952.2 8.650E-05 1.199E-05 8.768E-04 1.643E-05 2.095E-02 4.242E-03 4.096E-02 3.355E-05
2 2.00E-08 1.00E-08 1.23E+02 952.2 8.650E-05 1.200E-05 8.765E-04 1.640E-05 2.095E-02 4.242E-03 4.096E-02 3.358E-05
4 4.00E-08 1.00E-08 1.24E+02 952.2 8.651E-05 1.201E-05 8.763E-04 1.638E-05 2.095E-02 4.242E-03 4.096E-02 3.360E-05
6 6.00E-08 1.00E-08 1.24E+02 952.2 8.651E-05 1.202E-05 8.760E-04 1.635E-05 2.095E-02 4.242E-03 4.096E-02 3.362E-05
8 8.00E-08 1.00E-08 1.24E+02 952.2 8.651E-05 1.203E-05 8.758E-04 1.632E-05 2.095E-02 4.242E-03 4.096E-02 3.365E-05
10 1.00E-07 1.00E-08 1.24E+02 952.2 8.651E-05 1.204E-05 8.755E-04 1.630E-05 2.095E-02 4.242E-03 4.096E-02 3.367E-05
Simulation finished after 11 time-steps
Writing final snapshot at time t = 1.1000000000000003E-007
Wall clock time [hours] = 2.786E-04 (+/- 2.7778E-13)
Wall clock time/timestep/meshpoint [microsec] = 11.13094
Wall clock time/timestep/local meshpoint [microsec] = 44.52376
Rhs wall clock time/timestep/local meshpoint [microsec] = 261.2351
Maximum used memory per cpu [MBytes] = 43.668
Maximum used memory [MBytes] = 171.762
real 15.34
user 59.73
sys 0.19