SVN: ------- v. ( ) $Id$ SVN: start.in,v v. 1.3 (2009-03-17 10:59:27) nbabkovs SVN: run.in,v v. 1.3 (2009-03-17 10:59:27) nbabkovs 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 1-D run nxgrid, nygrid, nzgrid= 128 1 1 Lx, Ly, Lz= 0.50000000000000000 0.10000000000000001 0.10000000000000001 Vbox= 0.50000000000000000 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: AR ind_glob= 14 ind_chem= 9 CO ind_glob= 17 ind_chem= 12 CO2 ind_glob= 18 ind_chem= 13 H ind_glob= 9 ind_chem= 4 H2 ind_glob= 6 ind_chem= 1 H2O ind_glob= 8 ind_chem= 3 H2O2 ind_glob= 13 ind_chem= 8 HE ind_glob= 16 ind_chem= 11 HO2 ind_glob= 12 ind_chem= 7 N2 ind_glob= 15 ind_chem= 10 O ind_glob= 10 ind_chem= 5 O2 ind_glob= 7 ind_chem= 2 OH ind_glob= 11 ind_chem= 6 Number of reactions= 25 Number of species= 13 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 cp ( 6) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cp1 ( 7) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil gradnu ( 12) 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 cs2 ( 24) 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 timestep_LSODE: Second order symmetric splitting procedure (Strang splitting) 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 = 3854762023.0274343 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! 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--------rhom---------TTm-------Y1m-------Y2m-------Y3m-------Y4m-------Y5m-------Y6m-------Y7m-------Y8m----------urms-----------divu2m---- 0 0.00E+00 1.00E-08 0.000223 1625.0000 1.147E-02 1.365E-01 1.025E-01 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 1.0175000000E+03 0.0000000E+00 5 5.00E-08 1.00E-08 0.000223 1624.8094 1.147E-02 1.366E-01 1.025E-01 1.504E-08 7.696E-09 5.332E-07 4.181E-07 7.966E-10 1.0174992587E+03 8.4906395E+04 10 1.00E-07 1.00E-08 0.000223 1624.6184 1.147E-02 1.366E-01 1.024E-01 2.697E-08 2.452E-08 1.144E-06 8.089E-07 1.634E-09 1.0175064873E+03 4.6354455E+05 15 1.50E-07 1.00E-08 0.000223 1624.4273 1.148E-02 1.366E-01 1.024E-01 3.658E-08 4.745E-08 1.819E-06 1.174E-06 2.425E-09 1.0175170873E+03 5.8707781E+05 20 2.00E-07 1.00E-08 0.000223 1624.2359 1.148E-02 1.366E-01 1.024E-01 4.446E-08 7.435E-08 2.549E-06 1.516E-06 3.168E-09 1.0175307551E+03 4.4123695E+05 25 2.50E-07 1.00E-08 0.000223 1624.0443 1.148E-02 1.366E-01 1.024E-01 5.104E-08 1.037E-07 3.325E-06 1.835E-06 3.866E-09 1.0175490797E+03 3.7377528E+05 30 3.00E-07 1.00E-08 0.000223 1623.8526 1.148E-02 1.367E-01 1.024E-01 5.665E-08 1.343E-07 4.141E-06 2.134E-06 4.521E-09 1.0175715582E+03 3.0056701E+05 Simulation finished after 31 time-steps Writing final snapshot at time t = 3.0999999999999984E-007 Wall clock time [hours] = 1.547E-04 (+/- 2.7778E-07) Wall clock time/timestep/meshpoint [microsec] = 140.3730 Maximum used memory per cpu [MBytes] = 17.773 Maximum used memory [MBytes] = 17.773 real 3.23 user 3.18 sys 0.01