Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified 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 WARNING: run: quad precision suppressed! 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.5708200000000000E-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! 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.637E-04 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 10.53317 Wall clock time/timestep/local meshpoint [microsec] = 42.13269 Rhs wall clock time/timestep/local meshpoint [microsec] = 266.1853 Maximum used memory per cpu [MBytes] = 52.879 Maximum used memory [MBytes] = 208.922 real 14.61 user 57.14 sys 0.19