All qualified stuff today
ie, looking at programs, program cycles etc and deciding what to do.
RIDE produces heating rates and ion production rates. The effect on this on the atmosphere needs to be quantified. The initial problem is the H3+ column densities supported by auroral events. This means taking the production rates and sticking them into chemical rate programs before iterating through to produce H3+ densities, then integrating along the line of sight. RIDE produces production rates, chemical reaction rates module has come from KIM, and can be used with improvements from JIM and other sources. But this is just stage one...
The heating through chemistry and auroral events needs to be transported, as does the new set of chemicals. The heating also needs to iterate back into the rates. One part of the greater production rate at the peak of the ion production profile is that there is also a close by peak of the heating rates, allowing a greater number of endothermic reactions to take place.
When did I get into atmospheric chemistry? Or physics for that matter...?
To do this, need to understand the chemical reaction rates routine... make it my own... hehe
Inputs:
dtimeplas,iwrite,in,is, - related to position in flux tube, need to be rewritten for altitudes
He,H,H2,CH4, - neutral densities from RIDE, need to be expanded, eventually
Hp,H3p,Hep,elec, - ionic densities from RIDE - as above
prodhp,prodh3p,prodhep, - production rates of certain things from ride
betaChp,betaCh3p,betaChep,betaRhp,betaRh3p,betaRhep, - reaction rates, initial
CRC,vib_h2_sw,vib_h2_param, - related parameters
tarr,rch3p, - recombination rates
temp,recomb,chex - temps, plus charge exchange rate?
Usual declarations and inclusions. Most the things to alter from flux tubes to altitudes here.
Check for presence of specific things
All rates are zeroed to begin with
Rates are then calculated as a power law. Not all rates included, but pointers to remaining JIM rates are given where appropriate, and obviously they set up a pattern that can then be used for other rates beyond the original scope of this thing.
Rates then stuck into tubes and corrections for underflows etc made.
Have also been digging around for other things such as photochemistry et al, which provide either useful tools or contain bits of code that make the chemistry more complete - eg changing production and loss rates into actual chemical changes.
Opened thesis, didn't edit it though...
RIDE produces heating rates and ion production rates. The effect on this on the atmosphere needs to be quantified. The initial problem is the H3+ column densities supported by auroral events. This means taking the production rates and sticking them into chemical rate programs before iterating through to produce H3+ densities, then integrating along the line of sight. RIDE produces production rates, chemical reaction rates module has come from KIM, and can be used with improvements from JIM and other sources. But this is just stage one...
The heating through chemistry and auroral events needs to be transported, as does the new set of chemicals. The heating also needs to iterate back into the rates. One part of the greater production rate at the peak of the ion production profile is that there is also a close by peak of the heating rates, allowing a greater number of endothermic reactions to take place.
When did I get into atmospheric chemistry? Or physics for that matter...?
To do this, need to understand the chemical reaction rates routine... make it my own... hehe
Inputs:
dtimeplas,iwrite,in,is, - related to position in flux tube, need to be rewritten for altitudes
He,H,H2,CH4, - neutral densities from RIDE, need to be expanded, eventually
Hp,H3p,Hep,elec, - ionic densities from RIDE - as above
prodhp,prodh3p,prodhep, - production rates of certain things from ride
betaChp,betaCh3p,betaChep,betaRhp,betaRh3p,betaRhep, - reaction rates, initial
CRC,vib_h2_sw,vib_h2_param, - related parameters
tarr,rch3p, - recombination rates
temp,recomb,chex - temps, plus charge exchange rate?
Usual declarations and inclusions. Most the things to alter from flux tubes to altitudes here.
Check for presence of specific things
All rates are zeroed to begin with
Rates are then calculated as a power law. Not all rates included, but pointers to remaining JIM rates are given where appropriate, and obviously they set up a pattern that can then be used for other rates beyond the original scope of this thing.
Rates then stuck into tubes and corrections for underflows etc made.
Have also been digging around for other things such as photochemistry et al, which provide either useful tools or contain bits of code that make the chemistry more complete - eg changing production and loss rates into actual chemical changes.
Opened thesis, didn't edit it though...
Labels: 1D Hydrogen Atmosphere, Chemical rates
0 Comments:
Post a Comment
<< Home