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 work

Brief overview with the 'new era of Jovian auroral understanding' ie Steve's ideas on Makenzie's profiles. Yep, there are new ones, lots of them, with interesting things to look at.

One run finished (the last of the set that were running, meaning all I have to do is wait until the new regime is in place and set off the seven old ones again). Saturn5000 was the 12eV SD+1 run. It has been analysed.

Of course, first there is CTIP...

I have just had a look at putting the small increments thing in and it seems ok. Like everything else...

Sigh, "bus errror", should'a known it...

hmmm, will try a few runs on other computers using proton1...

Fails on rimmer

can't compile on C3

Going on keter... on the slow processors...

Trying out 31eV runs whilst waiting. It doesn't dip into the small increments level too deeply at this point, so may or may not avoid bus errors. On the other hand, ramping up the number of electrons fired does invite more errors to show their hand...

...bus error at 100k dammit!

Hmmm, difficult one. The lower energy stuff could probably quite happilly live without the small increments stuff, but the high energy stuff really needs it, otherwise there's the risk of crossing into an unrealistic pressure level. I'll run through the other energies and see if any of them can complete...

...one minor comfort, doing the 46eV lot, is that I do get ionisation in all pressure levels, as hoped.

...bus error at 100k

66eV
...100k again

96eV
...1M this time. Dammit... so close...

140eV
100k

206eV
100k

300eV
1M

440eV
1M

650eV
1M+ (hurray!)

950eV
1M

Don't yet have time to try out:
1380, 2030, 2980, 4400, 6400, 9400, 14000, 20800, 30500

I have actually remembered that since I have to divide all answers by a thousand to account for path length, 100k is more than enough anyway... and the answers (suggesting final peak rates from just shy of 1000 to just shy of 10,000 pcm3) is what is expected.

Must/may also do runs with pitch angle distributions, plus magnetic effects... after all...

...more and more and more...

Marking, 'bout halfway through

One run ended and was analysed. Saturn500 doing the 12eV SD-2 run. This leaves one run going on Keter, and six ready to make their reappearance once the queueing system has gone through...

Quick note

One run finished

Saturn100 finished the 12eV SD-C run. It was analysed etc.

Bit of stuff

One run completed today, the Saturn5000 40eV SD+1 run. This was the final 40eV run, and was analysed, recoded and set off as 12eV SD+1, meaning all 12eV SD runs are off, all 40eV runs are done, and, indeed, asides from the 100keV run, and 12eV runs, all is done. Sometime after the first week in March, the 12eV results should start rolling in, though this final run is taking a long while to set off due to TAMPA choosing today to launch mulitple copies of every program they have with no regard for proper submission protocols...

Also did all Howarth's marking.

Back to it

Meetings to do with seminar feedback today. Very positive.

Finished off Raman marking, then got some papers off Howarth, including one of the last lot, which has been dealt with, also saw anasuya about her marking, which needs swiftly shifting.

Problem class.

One run finished. Saturn500 finished the 40eV SD-1 run, was analysed recoded as 12eV SD-2 and reset off. This was the penultimate 40eV run remaining and means only two non 12eV runs remain to be completed. Also, only one of the SD 12eV slots is yet to be set off.

Run, marking, seminar

That's the plan.

The run that finished overnight was Saturn100, a 50eV SD-C run. It was analysed, recoded and set off as SD-C for 12eV. The analysis showed that it was the final 50eV run, meaning all energies except 12, 40 and 100,000 are done, with the remaining 100k and 40eV runs underway. Half of all the remaining to be completed and analysed 18 runs are SD runs.

Ok seminar. Three bits. Modelling, Photon, Auroral morphology - not in that order.

Begin with the Earth, starting at Bolearis, then shift outwards, until the Earth is just another planet. Show that visible light is just one of the emissions and the list of emissions on all eight planets (including the null of Mercury). Talk about ENAs at Mercury and why it hasn't an aurora. Talk about the dynamo at Venus, the swept aurora and how this is being applied to Mars. Move onto the Giant Planets, with Jupiter's internal dynamo, Saturn's state of affairs, and mention the other two. Then talk about how particle interactions change the current only analogy sometimes used. Have the photon as an anology. Then use it to introduce RIDE and shine up a result or two. Then finish and run off to Mill Hill to avoid questions...

Ok, most the way through it now. Finished off tomorrow.

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

Another quick note

One run:

Saturn12000 finished the 1keV SD+2 run, finishing off all 1keV runs. It was analysed and recoded as 12eV SD+2 and set off. Analysis revealed no grade changes beyond the end of 1keV (all energies asides from 100k, 12eV, 40eV and 50eV are now done). However, with only 11 runs left to be recorded as done, SD is finally down to one set of slots.

Quick note

Two completed runs:

Saturn5000 completed the 50eV SD+1 run, was analysed, recoded and set off as the 40eV SD+1 run (with a new 31 day runscript). Analysis saw no grade changes, but it did mean 22 runs were remaining to go - ie twice the number of available slots. Also, all 50eV runs and all 40eV runs were either done or going, indeed that is true for all energies except 12eV.

Saturn40 completed the 40eV SD+C run. This was analysed, recoded and reset off as the 12eV SD+C run. Analysis revealed no grade changes, but overall, we hit the 90% mark, so far as runs done goes.

Also spent half the (satur)day at Mill Hill, helping with the SPA tours.

Marking, marking, marking.

Marking

Ok, one run finished too. Saturn500 finished the 40eV SD-C run, was analysed, recoded and set off as SD-1. Analysis brought no grade changes, but SD finally hit the 80% complete mark. overall, just over 89% complete... 23 runs left to finish, eleven of which are currently going.

Marking

Have created a title for my seminar (under great duress) and grabbed a couple of images. Also tried to book the day off from Mill Hill...

Marking for Howarth done.

Quick note

One run and a bit of marking done so far today.

Saturn100 finished the 50eV SD-2 run, and was recoded and set off as the 50eV SD-C run, the only unfilled slot of that energy. Analysis revealed no grade changes.

Marking, marking, marking.

It hasn't stopped... Raman's stuff may be done and ready to be shipped off, but Anasuya's stuff remains... it becomes today's task. And tomorrow's etc etc etc...

One run has come to an end, Saturn12000 doing the 1keV SD+C thing. This was analysed, recoded and set off as the SD+2 1keV run, which will be the last of the 1keV runs.

Another quick note

Two runs finished today.

Saturn500 was the 40eV SD-2 run. This was analysed, recoded and set off as SD-C. This marked the 33 runs left point, meaning if all the runs not yet analysed were run in the eleven slots, they'd be used three times each. It also means that since eleven of them are running, that only twenty-two runs remain to be set off. No grade changes.

Saturn5000 was the 50eV SD+C run. This was analysed, recoded and set off as SD+1. Analysis brought 50eV up to a first, meaning all energies except 12eV are now on a first of one sort or another. 12eV remains a fail.

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

Quick note

No entry yesterday as most the day was spent going down a blind alley - a paper promising dissociation cross-sections for the three species that turned out not to quite be what I wanted - it was dissociation plus attachement at once. I could do a bit more work to extract something useful from there, and of course the process will be useful for other things eventually, however, not for the time being... sigh...

Today's delight was a completed run. The penultimate 100eV run from the Saturn100 slot. This was analysed, recoded and set off again. Analysis yielded no grade changes, though it did mark the point at which 80% of all the runs was completed. This was the SD-C run. The new run is SD-2 and will mark the completion (in about six days) of the 100eV runs.

Bits and bobs

Upgrade to titan's IDL install affected to allow running of Solar2000. Now works

Have been thinking again about the suicidal x-ray problem. Suicidal x-rays refer to the effect that an x-ray flux impacting on a neutral atmosphere will have. The atmosphere gets ionised by the impact, creating electrons. Electrons have a similar size to the wavelength of x-rays, which means Compton scattering cross-sections see their peak in the x-ray region. This, in turn means that x-rays get scattered more by the electrons they create. The effect of scattering is that x-rays lose energy and become EUV or less, meaning the x-ray flux kills itself. The chapman function method of calculating this is to apply the cf and have a compton scattering addition to the path length, ignoring the creation of EUV photons and concurrent destruction of x-rays, not to mention the directional effects of the scattering.

Another thought experiment has included space charges created in the aurorae and the currents created by the churning of background electrons to remove it. Though this effect is well known as it causes the induced currents that knock out powerlines and stations in powerful substorms, the heating effect in the upper atmosphere doesn't seem to be as well considered...

Of course, the majority of the work done today was on the Mars model, getting hold of, interpolating and extrapolating NO ionisation cross-sections. That is nitrogen monoxide, not "no" cross-sections... though the feeling did occur...

Saturn5000 slot run of 50eV in SD+2 mode completed. Analysis raised 50eV grade to 2i. Recoded as SD+C and set off again.

Another run dun

100eV finished the SD-1 run. This was analysed and the SD-C run set off.

Thesis work - the morphologies of aurorae on all (EIGHT) planets. Mercury's easy. Nothing seen, nothing expected... unless you count particle emission, ENAs, rather than photon emission... damnit. Time to alter the definition...

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.