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 and more and more...

Another run ended

Saturn1000 was the 12eV GS1m run. It has been analysed.

Demonstrating also happened, 2 hours of waves optics and accoustics in the maths problem class. That is now over and done with eternally. No more problem classes!!!!!!!!!

Anasuya sent another problem paper through alan, proving some things never end. Not alan...

Runs

News of more marking from Anasuya...

Two runs finished overnight and have been analysed. All runs are either going or done, so no recoding and resetting off anymore.

SaturnX was the 12eV GS100m run. This was analysed, completing the GS100m set.

Saturn50000 was the 12eV BSR0.25 run. This was also analysed, changing 12eV's grade to a pass and completing the BSR0.25 set.

This leaves 12 runs filling fourteen slots (as the middle run of the three sets is the control run, which is the same for all three).

Have looked at thesis. Looks the same as always

Ok, back onto things for the paper (remember that). Grodent gave me examples of his stuff for maxwellian energy distributions of characteristic energy: 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, 30, 40, 50, 100, 200, 500 keV, equivalent to mean energies of: 2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 30, 40, 60, 80, 100, 200, 400, 1000 keV. For it to work best, I'll keep the mean at 40keV at the top, 2 at the bottom and 10 and 20 in the middle - ie use the 1, 5, 10 and 20 keV maxwellian results to run this.

These need to be normalised to 50 erg cm-1 s-1 (although I then need to pull slowly away from cgi units...). For each of them, I need to generate an electron distribution, with appropriate normalisation. I then need to extract the atmosphere from each of the Grodent datafiles and then run the generated electron spectrum into the atmosphere using RIDE.

Have sorted out some emails that needed deleting. What's that about avoidance?

Electron numbers now calculated.

Last runs set off

For the sensitivities thing. Excluding all the required reruns that'll follow in April, of course. Decided to wait no longer and have everything going continuously until they relent and finish.

Speaking of relentless. Marking time. Had a late script off Howarth and still have stuff from anasuya to do.

Howarth done. All astro/medical Anasuya scripts done. Just leaves the physicists behind. As ever.

Another quick note

One run ended today

Saturn1000 doing the final 50eV GS run (GS1m). It was recoded and set off as GS1m for 12eV. Analysis revealed no grade changes.

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...

Marking, marking, marking.

Raman's stuff remains...

Just got some spam from an oscilloscope maker...!

A run finished. SaturnX doing the 12eV GS10k run. That has been analysed, recoded and reset off as 12eV GS100m. No grade changes yet, but GS10km is now utterly complete. And it is another 12eV run down.

Marking done! Unless Raman has more hidden away... Anasuya's next...

More

...of the same.

Saturn1000 finished the last of the 100keV GS runs, was analysed with no surprises or grade changes (though GS is now at the 90% complete mark along with BSR, whilst SD languishes at just above 70%). Saturn1000 is now running the 50eV GS1m run.

Now, Martian cross-sections, where was I?

Ionisation cross-sections now laboriously transferred.

Dissociation cross-sections now laboriously transferred. I think my fingers are now half their previous length...

eCheck and eLevel now adjusted (the arrays which say which species I have excitation cross-sections for (echeck) and what energy of excitation level this corresponds to (elevel)). The excitation cross-section array has also been suitably annexed. Cross now belongs to the Martians...

Module state of play:
randgen - ready
atom - ready
cross - ready
electron - needs conversion
strike - needs conversion
proton - needs conversion
proton1 - needs conversion
selectron - needs conversion
recoil - needs conversion
back - needs conversion
hit - needs conversion

...and the inputs are ok too. Also have a paper with energy spectrum to use if and when this gets ready to roll. Should I get the mag field, then the pitch angle spectrum is also available.

So, in the dying embers of the day, lets see what can be dealt with:
electron - Adjusted common blocks, declarations, code, compiled, done!
strike - adjusted common blocks, declarations, coded, compiled, done!
proton - adjusted common blocks, declarations, coded, compiled, done!
proton1 - adjusted common blocks, declarations, coded, compiled, done!
selectron - adjusted common blocks, declarations, coded, compile, done!
recoil - adjusted common blocks, declarations, coded, compiled, done!
back - adjusted common blocks, declarations, coded, compiled, done!
hit - adjusted common blocks, declarations, coded, compiled, done!

So, in theory a complete Martian model. Of course, a theoretically complete fortran model means only one thing - Segmentation faults lie in wait... we shall see.

And so we have... bastard thing. The debugger's no good too.

Right, trying Keter out, debugger works for once...

Model running

Model still running

Model still... WHY!!!! Did I wrong Mars in another life?

Suppose I'd better leave it there if it's going to run and run. Well, it won't go too long as I have imposed a twenty four hour limit on operations. We'll see if this thing really does work...

Have been told Anasuya's marking is now due!

Mars is finished. Really, finished... end of day.

Another run dun

The penultimate one of these fast runs.

Saturn1000 finished the 100keV GS100m. It was analysed, recoded as GS10km and set off for one final time. No grade changes etc.

Have got all the dissociation and ionisation potentials from NIST. Their electron impact ionisation cross-sections match with mine too. Still no word on scattering though! This means atom.f is finally finished.

Hurray! Found one! By typing in "total cross-section for CO" into google, I of course got the... well... CO2 one, but hey, it's still one of the three I need...

Now got CO as well... just need nitric oxide

Have extrapolated and interpolated my way through CO and CO2

and now have mined, interpolated and extrapolated NO as well. That's the easy bit done, now to stick them all into cross whilst managing not to delete any of the existing cross-sectionts that I intend to keep...

Scattering cross-sections now transferred laboriously to cross.f

One run to finish off with Saturn100 completed the final 100eV run (SD-2). It has since been set off as 50eV SD-2. Analysis was undertaken confirming 100eV is 100% done. No other grade changes.

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.

More and more and more...

Another block of maths demonstrating out of the way.

One run has finished. Saturn1000 slot 100keV GS1m run. Analysis was as expected. Recoded as GS10 version of same energy and set off again. Have decided to kill off all the fast runs asap rather than do them concurrently with the slow runs, then do all the slow runs in one big block - asides from those currently being done and their heirs...

Not that any 'fast' runs are left asides fro the 100keV ones. All medium to slow.

Also had a quick look at Makenzie's newly produced profiles. She hasn't produced the further out ones, which are the important ones yet, so will wait until tomorrow for discussions. Or later. Who knows...

...one last thing. References for Martian cross-sections kept from Europlanet, this list I will be going through for obvious reasons:
Straub et al 1996 - already got
Lindsay et al 2000 - now got
Manyan et al 2000 - ah, ManGan et al... ahem, now got...
Deutsch et al 2002 - now got
Haider et al 1992 (modelling) - JGR, dammit!
Erdman and Zopt 1983 (cameron band?) - ZIPF, what sort of a name is that? Far prefer Zopt... Anyway, not included in UCL's subscription to Planetary and space science, so no-go despite importance...
Fox and Dalgarno 1979 - JGR - though did get hold of three 1979 papers by Dalgarno relating to proton processes I was in dire need of for DICE...
Lummerzhien and Lilensten 1994 - Lummerzheim, given I own half his back catalogue, it sometimes seems, I should know how to spell his name... plus annales has withdrawn the PDF! Did they know I wanted it?
Tran and Vidal 1998 - TIan... now got, oh!, there's another, got it too. Waaaaait. A third useful Tian and Vidal 1998 has appeared with cross-sections, my they had a busy year - nope, here's number four, maybe Tian and Vidal 1998 is a generic reference to their entire career...
Torr and Torr 1985 (Photon)
Avalgan et al 1998 (Photon)
Hitchcock et al 1980 (Photon)
Masuoka 1994 (Photon)

Will deal with the photon stuff another time. Mostly have them, but no real prospect of a photon thing to deal with yet as Tommi is otherwise engaged.

I make no garrentees that any of those names are spelt correctly, or in any way resemble the correct spellings.

Couple of runs

Much of the day has been taken up with Mill Hill and taking up Outreach leaflets to the observatory, but beforehand two runs did finish, the analysis of which was kept until afterwards.

The first was Saturn1000 finishing off the 1keV GS runs with GS1m. This was recorded and then recoded as 100keV GS1m run and set off. Later analysis was as expected. No grade changes.

The second was Saturn12000, doing a 1keV SD run (SD+1). This was recorded, recoded as SD+C and set off, making it the penultimate 1keV run. Analysis as expected. SD grade changes to a first. It is 70% done, GS is more than 80% done and BSR is 90% done. Overall, just over 80% done.

More and more and more...

Saturn40 has finished the final 40eV GS run. This has been analysed, recoded and set off as an SD run. Analysis raises the 40eV grade to a first.

Some research done on Martian aurora. Altitude is about 129km, cf Martim altitude range of about 60-250km. Have got hold of a quick 1D Martian atmosphere to do a few runs on, but more cross-sections are required, as ever. The species involved are CO2, N2, O2, O, CO, NO and Ar. From the Earth, cross-sections are provided for N2, O2, O and Ar. This leaves CO2, CO and NO in need of some data, slightly unfortunate as the main emitting species is CO2 and other bands of interest include CO bands. Nevertheless, a quick browse about has netted a few cross-sections and as I know the names of the bands I need to work with, I can grab them - or at least have a better idea of where to look. A current and energy distribution have also been uncovered, suggesting a number flux can be determined, meaning I can then simulate a few emissions, see if they lie in the right area and are of the correct magnitude, then a more complex look should provide interesting things. Interesting things to look at include band emission ratios and the effect of pitch angle on all this, as field aligned electrons seem not to be dominant, apparently. We shall see.

The Martian atmosphere is a lot thinner than the Earth's, this means the number densities cover a very small altitude range, which in turn means a fast running program, hopefully. The low energy and flux of electrons involved in the Martian aurora should also help in that...

Atmosphere extrapolated. Tomorrow, I will create the program base.

A bit done

Problem class has taken up a couple of hours. Demonstrating mathematics with a half written answer sheet. How I love UCL...

The Thesis has also had a bit of attention. History of observations of the aurora, just the timeline, individual bits on particular observations to go in later. I knew that long winded talk I did last year would come in handy at some point...

Saturn500 finished the 40eV GS100m run. This is the last of the 40eV GS runs not being dealt with, so moved it onto the SD runs. 40eV might actually get finished! Analysis of the run complete, no grade changes, but 80% of GS runs now complete, along with 90% BSR runs, 65.7... etc % SD runs and 78.5... etc % of total runs.

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.

Quick note

Saturn1000 finished the 1keV GS10m run. It was analysed, and the GS1m run set off.