compilation test...

Randgen - compiles fine of course (I never wrote it...)
atom - one array too big, sorted
cross - one array to small, sorted
electron - statement number error, sorted
strike - fine
proton1 - fine
hit - fine
back - fine
selectron - same statement error as electron, sorted
recoil - fine
proton - fine
Jupiter - fine (aledgedly...)

Just a couple of quick runs to test it
10keV, pitch angle zero (field aligned), 1000 electrons

Loads of questions of course - how much computational time will be taken by the larger atmosphere? Is the atmosphere in the right way up for this program? Does it actually work? All the usual ones...

well, at least I know they can run for ages... lets see if they can finish...

Meanwhile, back at the office, the new keter submission system is alledgedly up and working. The maximum runtime for the fast server is 62 days, the max for the medium one is 124 and the slow servers are unlimited. Great, except for one minor quibble. The maximum runtime it is possible to specifiy on the runscript is 999hours 59mins 59 seconds, or just over 41 days... Dontcha love basic maths... of course, it may be that the numbers can be varied... that had better be the case, at least. Lets see if the sensitivities will go back up (all seven remaining ones...)

Saturn2 - BSR/GS/SD 0.5 12eV now off on new V890 (or queued at least)
Saturn100000 - BSR1.0 100KeV now off on new V890 (or queued at least)
Saturn1 - GS10m 12eV now off on new V880
Saturn3 - BSR0.75 now off on new V880
Saturn40 - SD+C 12eV now off on new V880
Saturn50 - SD-1 12eV now off on new V880
Saturn12000 - SD+2 12eV now off on new V880

All seem to be ambling along (those that are set off that is...). We shall see if they work...

Jupiter run finished. Atmosphere is upside down. Sigh.

Sorted, attempt 2

That worked fine. Now need to do two things, the first is run an energy of known height, then known energy influx.

What else needs to be done?

Validation of terrestrial electrons - need to get another model and go through this again
Validation of terrestrial protons - as above
Photon validation - got 'model', so ready
Ovation stuff - once validated, the electron ionisation results for this can be used
Hbeta runs - when validated and when data arrives, this can be done
Photon runs - can be done during validation
Makenzie stuff - when data arrives etc - connects on from this current stuff
MICdev stuff - improvements should I have time

So, jobs:

Get proton model
Get electron model
Validate jovian electrons (being done)
complete sensitivities (being done)
Validate protons
Validate electrons
Validate photons
Ovation ionisation rates
Ovation CSDA model
Ovation/TIROS comparisons
Extend Grodent profile
Extend Grodent atmosphere
Fit to Makenzie data
Experiment with proton Hbeta doppler profile parameters
Fit to HBeta data
Run photon scattering models
Write up!

There. Easy...

Just been battling latex again, identifying where the missing style files are. Between starlink and apl, I might be able to muster all the required files... we shall see...

More stuff

One run finished overnight:

Saturn50 50eV SD-1. This was analysed, recoded as 12eV SD-1 and set off again. No grade changes.

Mill Hill

Three runs - and marking

Three runs have finished overnight (well, they didn't finish in an overnight burst, they just happened to finish, after a long time, between me leaving last night and arriving today... but you know what I mean.. don't you?)

Saturn12 is the 12eV BSR1.0 run, finished with a day to spare... It was recoded, analysed and reset off as BSR0.9. Analysis gave no grade changes, but it did finally finish off the BSR1.0 set of runs, the first of the BSRs to finish.

Saturn50 is the 50eV GS10m run. This was analysed and then recoded as SD-1 (since the last 50eV GS run is already being done) and set off, meaning all six slots that are destined to attack the SD 12eV runs are in the SD block, ready to roll once they've done their stuff on 40, 50 and 1000 eV. No grade changes.

Saturn50000 was the 12eV BSR0.1 run. This was analysed and recoded as BSR0.25. Analysis revealed no grade chages, but BSR0.1 became the second BSR set to finish.

A good day for 12eV!

...and some marking...

Couple 'a runs

Since Steve isn't in after all...

Two runs have finished overnight

Saturn50 finished the 50eV GS100m run. It was analysed, recoded as GS10m and set off again. No surprises or grade changes.

Saturn1000 finished the 100keV GS10m run. It was analysed, recoded as GS100m and set off again. No surprises or grade changes again.

Dealing with those wonderful cross-sections:
Straub et al 1996 - no further information to extract
Lindsay et al 2000 - dissociation cross-sections for NO, plus further references - mined
Mangan et al 2000 - CO dissociation cross-sections - mined
Deutsch et al 2002 - relevant to deeper model (ionisation of an ion)
Butler and Dalgarno 1979 - protons
Butler et al 1979 - protons
Bieniek and Dalgarno 1979 - protons
Tian and Vidal 1998 - information obtained
Tian and Vidal 1998a - relevant to deeper model, up to quadruple ionisation/dissociation
Tian and Vidal 1998b - relevant to deeper model with C2H2 and Ch4 cross-secs as well as CO
Tian and Vidal 1998c - CO2 dissociation - mined

Good, that's dissociation done as well as ionisation. Only require elastic scattering and maybe a few of the excitations to finish with.

Its amazing how large the NO dissociation cross-section is compared to its ionisation cross-section, compare to the disparity between say N2 or O2. Must be something to do with the heterogenaity of the atoms redistributing electron density unevenly across the bond. As ever.

Yep, same thing with the CO one.

Even worse for CO2, guess its because its polyatomic too.

So now have all appropriate dissociation and ionisation cross-sections, just missing elastic and (eventually) excitations, though the latter aren't necessary for immediate evaluations of the profiles.

Another quick note

Two runs completed today:

Saturn50 finished the 50eV GS10km run. This was analysed, with no grade changes, and the next run set up and off.

Saturn12000 finished the 1keV SD-1 run. This was also anaylsed, set up and off again.