fidz3gm = ncdf_open('data_dump_Z3GM_MillstoneHill.nc')
fideDens = ncdf_open('data_dump_eDens_MillstoneHill.nc')
ncdf_varget, fidz3gm, 'data', Z3GM
ncdf_varget, fideDens, 'data', eDens

uthours = findgen(288) / 12

levels = 42
clevels = findgen(levels) / 13.65 + 9
psfile = 'eclipse_Millstone.ps'

set_plot, 'ps'
loadct, 25
DEVICE, /COLOR, FILENAME = psfile
CONTOUR, alog10(eDens[*,0:100]*1.E06), uthours, Z3GM[*,0:100]/1000,/fill,levels=clevels,yrange=[100,500],xrange=[0,24],xstyle=1
device, /close