Back to J

Correcting the small increments thing in the Jupiter model. It is in there, but not recorded, bizzarely enough.

Subroutines to alter:
electron - checked, declared, changed - done!
strike - checked, declared, changed - done!
proton - declarations - done!
proton1 - declarations - done!
hit - checked against Earth model, changed declarations, changed recording - done!
back - checked, declared, changed - done!
selectron - checked, declared, changed - done!
recoil - checked, declared, changed - done!

All done.

Compiled.

Running.

Why did I do that? All the hard work had been done, just hadn't redeclared the results arrays as real variables, instead leaving them as integers. Makes it easier to analyse, I suppose, but artificially inflates results in denser or energy degraded areas. I know this has little effect in Gas Giant atmospheres and I may have removed it to reduce computational time and keep things looking pretty, but the calculations were left in, so no time was reduced and things won't look pretty if they're recording incorrectly...

Ok, results are in. That's better. Enough of this madness, at last. Shape is looking good, peak position great, magnitude very good.

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.