More stuff

Text:
Results
-Ion production and direct heating - needs results
-Chemistry (further)? - needs results
-Heating (further - Joule heating etc)? - not required
-Photons - needs further results (photon chemistry)
-Sensitivities - needs results
Discussion - needs above results
Conclusions - needs above results

Figures:
4.6 - Validation monoruns - runs ongoing
4.13 - Uranus North ion production - ongoing
4.14 - Uranus South ion production - ongoing
4.16 - Uranus heating - ongoing
4.19 - Peak rate altitudes, Uranus - needs runs
4.24 - H3+ densities for three cases Saturn - needs runs
4.25 - H3+ densities for three cases Uranus - needs runs
4.26-8 - Sensitivity graphs for Saturn, three whole 3D graphs - runs ongoing

Tables:
4.2 - Ion production Uranus - runs ongoing
4.4 - Heating Uranus - runs ongoing
4.5 - H3+ Saturn - needs runs
4.6 - H3+ Uranus - needs runs

Modelling required:
Validation runs (three cases - highly intensive due to powerful aurora) - ongoing
Uranus runs (two cases - very low intensity aurora) - ongoing
Chemistry module - coding to update for certain Steve based things, then runs using model atmospheres and ion production
Sensitivities - you can't miss these runs, they've been going for a year now... - ongoing

The problem with vanishing hydrogen in the chemistry routine has itself vanished as the coding mistake was discovered. The problem with H3+ being depleted compared to Grodent is likely to be due to the philosophy of the code (me not including transport terms) than anything else as the rates are very consistent and believable. Of course, the possibility that not having the H+ or He+ production rates reduces the amount of H3+ due to the lesser production through eqn 9 (lesser as some H+ appears chemically, so equation 9 does produce a little H3+) or through production of H2+ through charge exchange with He+ (also a production route for H+, which then adds to 9). Then there's the problem of high density regions requiring tiny timesteps to stop the densities overpowering fortran, and low density regions requiring low timesteps so their changes don't become so tiny they vanish in a rounding error... This at least is being worked on and the chemistry routine is starting to look a little like it could work in the future... though the runtime will be massive, especially when Steve's stuff goes in...

Final proof of the above arrived. The H3+ densities conform almost entirely to the H2+ ion production rate, which agrees with Chris' assessment that all H2+ goes to H3+ (within a certain tolerance). So using eqn 9 and charge exchange will top up the H3+. Plus another source of ionisation has been divined... photoionisation, which I forgot... I have the rates, or will have in a sec, but haven't stuck them in. Tut.

OK, have been fitting in the extra bit - variable resolution to cope with the different densities, and the test worked really well. Have had to adjust it as H3+ was being ignored (it needs a few iterations to get going) but got it done in the end. The only real question is how many iterations to allow it to have... as usual a case of results versus time - though the initial results were so good that may not be of concern. We shall see...

Further heating now cancelled due to lack of spectral generators to act as coolant in any atmosphere (asides from the H3+ coolant code from JIM etc). Having said that, lessons from the chemical rates module have informed corrections to the thermal profile one.

Also answered questions for BAS