Bits and bobs

Group meeting occured on the topics of the Move and the Rolling Grant, plus a little bit about computing.
The Move: The meeting took place in the new PhD office. People have been assigned to the office and a further meeting was to happen afterwards to decide desk allocation. This has happened, room has been rearranged, and asides from grabbing a couple of sofas and a coffee table, all is ready.
The Rolling Grant: All I heard was "moan, moan, moan" which translates as "they've changed the way the grants are assessed". Should be fun...
Computing: After yesterday's events, APL has threatened to withdraw funding from UCL's supercomputers by using the Miracle grant. We'll see how negotiations with getting serial jobs back on the fast processors go. Steve is all for trying to renegotiate everything back. Alan wants to pull us out...

Other things.

Howarth handed me more marking and left a message with my officemates about him mailing me the model answers later. This was not conveyed until Howarth returned and made them convey it...

Anasuya got her marking back in the group meeting. I wasn't missing a script.

One run finished, or rather hit the time limit... The only one that was ever likely to. It has been reset off on Rimmer, where it joins a long CMAT2 run already going, which means both of Rimmer's processors are now in use. Where's Alun...

Antionio Hales is our newest Dr

Grodent stuff moving along slowly.

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.