Figures:
Fig 1a: denslice,37,’x’,5.9,/io; density slice for low-density model
Fig1b: denslice,93,’x’,5.9,/io,vmax=5,vmin=1; density slice for high-density model 
Fig 1c: alfspeedslice,37,’x’,5.9/io; VA for low-density
Fig1d: alfspeedslice,93,’x’,5.9/io; VA for high-density
Fig2a: poyntslice,80,'x',5.9,600,600
Fig2b: poyntslice,88,'x',5.9,600,600
Fig3a: poyntslice,80,’x’,5.9,600,600,vmax=1
Fig3b: fieldslice,80,’ex’,’l’,5.9,600,600,vmax=5
Fig4a: fieldslice,80,’jz’,’k’,377,600,600
Fig4b: fieldslice,88,’jz’,’k’,377,600,600
Fig5a: poyntslice,81,’x’,6.1,600,600
Fig5b: poyntslice,89,’x’,6.1,600,600
Fig6a: fieldslice,81,’jz’,’k’,377,600,600
Fig6b: fieldslice,89,’jz’,’k’,377,600,600
Fig6c: fieldslice,81,’jz’,’k’,7,600,600
Fig6d: fieldslice,89,’jz’,’k’,7,600,600
Fig7a: poyntslice,90,'x',5.9,600,600
Fig7b: poyntslice,91,'x',5.9,600,600
Fig8a: fieldslice,90,’jz’,’k’,377,600,600
Fig8b: fieldslice,91,’jz’,’k’,377,600,600
Fig9: jphihist,94,5.88374,83,1,1,600
Fig10a: fieldslice,94,’jz’,’k’,377,600,600
Fig10b: phiparslice, 94,600,600,1,vmax=25

Density models:

In these models, the 5 columns are:
Radial distance, Latitude, Number density, Ionosphere mass density, Torus mass density

1:  Frank Crary hi-den case
2:  Based on Licia model, but dene symmetric
3:  Fixed dene, now organized by centrifugal lat, model 2, mode 3, n0=2e5, h=1200 
4:  same as 3 but increase scale height to 3000
5:  hi-den with larger scale height