Wednesday, June 20, 2007

Even four

Text:
Results
-Ion production and direct heating - needs results
-Chemistry (further)? - needs coding
-Heating (further - Joule heating etc)? - needs coding
-Photons - needs results
-Sensitivities - needs results
Discussion - needs results
Conclusions - needs results

Figures:
4.6 - Validation monoruns - runs ongoing
4.13 - Uranus North ion production - ongoing
4.14 - Uranus South ion production - ongoing
4.16 - Uranus heating - ongoing
4.19 - Peak rate altitudes, Uranus - needs runs
4.24 - H3+ densities for three cases Saturn(?) - needs coding, runs
4.25 - H3+ densities for three cases Uranus(?) - needs coding, runs
4.26-8 - Sensitivity graphs for Saturn, three whole 3D graphs - runs ongoing

Tables:
4.2 - Ion production Uranus - runs ongoing
4.4 - Heating Uranus - runs ongoing
4.5 - H3+ Saturn(?) - needs coding, runs
4.6 - H3+ Uranus(?) - needs coding, runs

Modelling required:
Validation runs (three cases - highly intensive due to powerful aurora) - ongoing
Uranus runs (two cases - very low intensity aurora) - ongoing
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

Two bits of coding remain as well as loads of calculations etc. Once these two are done and set off, I can start reformatting graphs on the paper version. I have two to compress and one to create, which depends on one of the set of runs listed here...

More questions from BAS... sigh. If they'd only wait a couple of days, I'll be on Chapter five anyway...

Testrun numbers calculated for Grodent cases (ok, updated...). All cases running (all four that is, as one extra has crept in). This means keter now has Uranus, Validation and Sensitivities runs going. C3 has alpha runs, except one that rimmer has.

Soooo. Chemistry. Chemical rates are heirarchical, which means you must get the major species right before thinking about the minor species, else you'll be iterating massively forever. H2 is the most major species, followed by He and then H (neutrals). The largest charged species is e-, followed initially by H2+, later to become H3+, then He+ and H+. I could go further, inlude hydrocarbons, He++ etc, but for the initial check not necessary - though Steve wants rotational level chemistry, which means sticking in all 27 levels of H3+ seperately and letting them decay into appropriate populations...

I have devised simple balance equations, using the production, loss and initial density terms for each thing. Next will be to iterate over and again until this lot reach new densities. Highest danger will be to the photoionisation rates at the top, which will mean considerable ionisation of some low density species, low densities fluctuate, hence so will the photoionisation rates. Similar considerations are required for auroral ionisation later... The reason for using this order isn't just because the big things changing changes the little things so much, but also because the little things change the big things so little...

But first, bits and pieces of graphs. Ok, after battles with p!multi, extra y-axes, positioning things etc, I finally have two graphs in one eps with both pressure and altitude scales attached. The other graphs are less essential. Or perhaps I'm lazy, either way they won't get 'em. Have adjusted captions and official review response. There now really is only one thing to do - the validations, which are cooking on Keter as we speak.

Labels: , , , , ,

Thursday, April 26, 2007

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

Labels: , , , , , , , , , , , , , ,

Friday, March 09, 2007

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.

Labels: , , , , , , , ,