...continuing...

BSR:
PH2 - done!
FWHM - done!
Integrated - done!
E per ion-electron pair - done!
Peak/integrated - done!
Peak rate position - done!

GSD:
PR - done!
PH1 - done!
PH2 - done!
FWHM - done!
Integrated - done!
E per ion-electron pair - done!
Peak/integrated - done!
Peak rate position - done!

SDC:
PR - done!
PH1 - done!
PH2 - done!
FWHM - done!
Integrated - done!
E per ion-electron pair - done!
Peak/integrated - done!
Peak rate position - done!

...all results now in files. The stats in the spreadsheet can either be ignored or added into tables later. Certainly, there should at least be some comment as to which results show nonlinear trends and which seem to be more noise than actual deviation from a simple baseline.

Now, of course, I have to turn the present example plot into a general purpose multiplot... first thing's first - correct the axes... done!

Sensitivities

Finally mastered 3d graphs for idl (that was a harsh half hour) and will be plotting out the survived results of sensitivities (ie all except 100k BSR 0.0 results).

There are three categories of graphs - BSR, SD and GS representing backscatter ratios, errors in atmosphere/cross-sections and Grid Size (atmospheric contraction/expansion induced changes) respectively.

Within each category, there are eight graphs, representing Peak ionisation rate, peak height measured by which altitude saw the maximum rate (monte carlo models have inherent random effects that can make this uncertain) peak rate as measured through an average of five or so levels rather than single levels (biased towards higher altitudes due to asymmetrical shape of the Bragg curve), FWHM, integrated ionisation, average energy expended per ion-electron pair created (assuming all energy goes into ion electron creation), peak as a fraction of integrated rate, peak rate position as a fraction of the distance from the lower half width to the higher half width. This gives twenty four sets of graphs in total.

Within each graph, the points are ordered by initial electron energy. There are ten electron energies, providing a total of 240 points to litter the appendix with. P.Multi may be used...

I can fit around eight graphs to each page, which gives me three pages of graphs... not the most exciting tract of the thing, but then again could be for not very exciting people... actually, maybe six pages, four graphs to a page, these are pretty big things once the annotations go on (as tested now).

Also need to decide which example graphs go in the main text to complete the figures in chapter four (I say complete, despite some figures not yet being in there and others mislabelled due to the removal of Uranus).

Now transferring some of the datasets:

BSR:
pr - done
ph1 - done

Revising four

Due to the Uranus crash on keter, chapter four is getting pruned. I have converted all necessary text in that chapter and will also alter the planout. Chapter three has also been returned with scribbled corrections to be dealt with.

Original plan:
Text:
Results
-Ion production and direct heating - needs results - saturn already done, uranus scrapped...
-Chemistry (further)? - needs results
-Heating (further - Joule heating etc)? - not required
-Photons - needs further results (photon chemistry)
-Sensitivities - needs results
Discussion - needs above results
Conclusions - needs above results

Figures:
4.6 - Validation monoruns - runs ongoing (getting reset off again...)
4.13 - Uranus North ion production - ongoing - scrapped
4.14 - Uranus South ion production - ongoing - scrapped
4.16 - Uranus heating - ongoing - scrapped
4.19 - Peak rate altitudes, Uranus - needs runs - scrapped
4.24 - H3+ densities for three cases Saturn - needs runs - ready to be run
4.25 - H3+ densities for three cases Uranus - needs runs - scrapped
4.26-8 - Sensitivity graphs for Saturn, three whole 3D graphs - runs done (asides from one that died)

Tables:
4.2 - Ion production Uranus - scrapped
4.4 - Heating Uranus - runs ongoing - scrapped
4.5 - H3+ Saturn - needs runs - runs ready to go
4.6 - H3+ Uranus - needs runs - scrapped

Modelling required:
Validation runs (three cases - highly intensive due to powerful aurora) - ongoing
Uranus runs (two cases - very low intensity aurora) - ongoing - scrapped
Chemistry module - coding to update for certain Steve based things, then runs using model atmospheres and ion production
Sensitivities - you can't miss these runs, they've been going for a year now... - ongoing - done, asides from dead one

New plan:
Text:
Results
-Chemistry (further)? - needs results
-Photons - needs further results (photon chemistry)
-Sensitivities - needs results
Discussion - needs above results
Conclusions - needs above results

Figures:
4.6 - Validation monoruns - runs ongoing (getting reset off again...)
4.24 - H3+ densities for three cases Saturn - needs runs - ready to be run
4.26-8 - Sensitivity graphs for Saturn, three whole 3D graphs - runs done (asides from one that died)

Tables:
4.2 - H3+ Saturn - needs runs - runs ready to go

Modelling required:
Validation runs (three cases - highly intensive due to powerful aurora) - ongoing
Chemistry module - coding to update for certain Steve based things, then runs using model atmospheres and ion production

Please note, where the word 'scrapped' has been used, it is short for "still going, but will be used in a paper/project afterwards". I think there is as much mileage in the Uranus ideas as in the rest of the thesis, all things considered. This does help me rationalise chapter two a little too.

...also chatted to Alan about three...

I have been running and rerunning and debugging and, well, more debugging

MicroRide is just so small its difficult to know when its working or not...

More microRide

Have adjusted declarations in all modules to reflect the new input files...

Back.f has had its indices changed to reflect this as well as recording systems etc and beta
Hit.f has also been adjusted
Proton1.f has been adjusted
Electron.f has been adjusted
Recoil.f has been adjusted
Selectron.f has been adjusted
Proton.f has been added and adjusted

Slight battle with the forces of IS, more with the forces of email inertia giving (partially) wrong impressions of IS

All modules compiled first time, setting off 10keV field aligned stuff.

...which of course failed.

...now debugging...

Been investigating some of the doom laden keter runs. One finished yesterday, Grodent1. Dead. Will be set off again soon...

But back to microride:
Have adjusted the magnetic field strength and dip angle to the correct number of altitude levels.
Also transposed all the number densities of the species as well as neutral temperatures and replaced the original files with the ones used for the Rees/Tiros/Rees2 stuff. These files (Magnetic, number densities and temperatures) will be used for both microRide and microDice.

Have obtained pixel numbers of the MicDev and converted these to central wavelengths. The FWHM has been revealed to be roughly 2.5 pixels. I can therefore bin the wavelengths by their shifts into pixels and then describe the output as a series of Lorenzians as before in a 2D array similar to what I have at the moment. This 2D array of pixel intensities then needs to be convolved with the completed Lorentzian and then it is done.