+ Now report total energies (summed over phases) + And also total summed over type (thermal, kinetic, magnetic)
… up yet The total energies, calculated from the `en` cubes are of order `1e+58` whereas the sum of the kinetic, thermal, and magnetic energies are only of order `1e+47`, so something is wrong.
Not sure why these were done - it was in December and I don't understand why they never got committed.
…this will fix problems with some PDF viewers
+ mhdcuts-Bkey.py - make key bigger 128x128 and do some variants. Most are not used, except for the ones called "sat" and "val". In the "sat", we have the saturation falling off towards the middle of the disk, corresponding to increasingly prominent out-of-plane component of the field. But unlike in the original version, we also reduce the value slightly in the middle (to 0.8), so it doesn't look so washed out there. In "val", the value is reduced towrads the center, representing reduced B field. + mhdcuts-B30krum-stitchup-nolabels.py and krumx too - use the "sat" and "val" variants of the color key. Put them at the top of the top panels, instead of between top and bottom panels.
…tes for myfitsutils.py + projected-bfield-zoom.py: - Lots of new command line options for tweaking the appearance: --scale, --gamma, --bscale, --molecular-bscale - Disable warnings, so that PyX and PyXgraph don't write so much rubbish + myfitsutils.py: FitsImage clas has new public attributes min, max
…-zoom.py + makerotbmaps.f90: now does integration of B-field projected along q,u axes before converting back to x,y axes at the end. This deals with the theta, theta+180 degeneracy. Final method does it with angles, not with sqrts. + projected-bfield-zoom.py: by default uses the "unsigned" vector field. Now uses data arrays directly, rather than going via the parametric function. Still need to port this back to projected-bfield.py
+ Use new runs B30krum, B30krumx + Include the gas velocity at the fronts
+ Now in glorious black and white + Keys no longer overlap lines + Corrections to labels + Larger font + Densities of all 3 phases + Mass and volume fracs of all 3 phases
…f90) + cubedenstats.f90: Now distinguish between molecular and neutral gas + molfrac.f90: Add double precision version of the function
…ages (mhd-pressures-rgb.py) + mhd-pressures-rgb.py: Increase gamma of the n-B and pram-pmag plots + plotstats_turb_cmp_m.py: Add graphs for O star runs
… gas dynamics better
Extracted generic bivariate stuff into bivar.py, which is can now be used by the programs like mhd-pressres-rgb.py, although it isn't yet. New file vel-distro-rgb.py that does bivariate distribuion of radial velocity and radius for molecular, neutral, and ionized gas. This has problems on my laptop, and I haven't tested it on the linux servers yet. Line graph program plotstats_turb_cmp_m.py now does molecular, neutral, and ionized separately. This is what gave a negative velocity for the molecular gas, which started all these problems.
* makerotbmaps.f90 - save the column densities (ionized, neutral, molecular) as well as the integrated B fields * projected-bfield.py - now do 4 panels without any text on the axes. Labels are moved inside the box. Three panels of the B field for each phase, plus one panel of RGB column density map.
… graphs of them * makerotbmaps.f90 - this makes the maps * projected-bfield.py - and this makes the graphs - grayscale of LOS component plus vectors of POS components
- bivar.py is the bivariate PDF image plotting stuff that was previously used in programs such as mhd-pressures-rgb.py - molfrac.py contains a single function molfrac(AV) that calculates the fake molecular fraction These modules will get used in subsequent programs
…tats, cubedenstats). This eliminates the annoying ~10% discrepancy that was creeping into the mass fractions.
- Make sure the fractions are positive - Correct error with mass-weighted radial velocities - Volume fractions for neutral/molecular - Mean radius of dissociation front