Pencil1: iapp, nprocs, ncpus= 0 4 4 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 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! diffusion: nothing (i.e. no mass diffusion) initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. 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! 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 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 gYYk ( 13) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil lambda1 ( 18) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil S0_R ( 23) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cs2 ( 24) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil glnpp ( 25) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil del2pp ( 26) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil glnrho2 ( 38) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil del2lnrho ( 39) 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.651E-05 1.199E-05 8.772E-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.652E-05 1.199E-05 8.775E-04 1.637E-05 2.095E-02 4.242E-03 4.096E-02 3.361E-05 6 6.00E-08 1.00E-08 1.24E+02 952.2 8.653E-05 1.199E-05 8.779E-04 1.635E-05 2.095E-02 4.242E-03 4.096E-02 3.363E-05 8 8.00E-08 1.00E-08 1.24E+02 952.2 8.654E-05 1.199E-05 8.782E-04 1.632E-05 2.095E-02 4.242E-03 4.096E-02 3.366E-05 10 1.00E-07 1.00E-08 1.24E+02 952.2 8.654E-05 1.199E-05 8.785E-04 1.629E-05 2.095E-02 4.242E-03 4.096E-02 3.368E-05 Simulation finished after 11 time-steps Writing final snapshot at time t = 1.1000000000000003E-007 Wall clock time [hours] = 4.568E-04 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 18.25125 Maximum used memory per cpu [MBytes] = 42.812 Maximum used memory [MBytes] = 168.906 real 24.32 user 94.41 sys 0.38