pro tecplt
alat=fltarr(73,61)
alon=fltarr(73,61)
tec=fltarr(73,61)
xla=fltarr(61)
xlo=fltarr(61)
y=fltarr(61)
x=fltarr(73)
infile = 'mlatlon.1995'
tecfile = 'tecgrid'
openr,1,infile
for lo=0,72 do begin
	readf,1,format='(61f7.2)',xla,xlo
	alat(lo,0:60) = xla
	alon(lo,0:60) = xlo
endfor
close,1
openr,1,tecfile
for lo=0,72 do begin
	readf,1,format='(61e10.3)',xla
	tec(lo,0:60) = 1.e-16*xla
endfor
close,1
for i=0,60 do y(i)=3.*(i-30)
for i=0,72 do x(i)=5.*(i-36)
set_plot,'ps'
map_set,/continent,/cyl,title='Total Electron Content (December 21, 12 UT, Solar Maximum)',/noborder
!p.thick=2.
contour,tec,x,y,/overplot,levels=findgen(17)*10,c_labels=[1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]
;Make magnetic equator thick
!p.thick=4.
contour,alat,x,y,/overplot,levels=0.
!p.thick=1.
!x.margin=[1.,1.]
!y.margin=[1.,2.]
plot,[-180,180],[-90,90],/nodata,/noerase,xstyle=1,xrange=[0,24],ystyle=1,yrange=[-90.,90.]$
,xtitle='Local Time',ytitle='Latitude',yticks=6,yminor=3,xticks=4,xminor=6
device,/close
end