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

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.

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.

Quick note

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