pro tecplt,ps=ps

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)

xll = .15 & yll = .35
xur = .95 & yur = .85

 loadct,5  ;  5 -> Std Gamma II
;loadct,15 ; 15 -> Stearn Special
;blue_white_red,ps=ps

if (keyword_set(ps)) then begin
 set_plot,'ps'
 device,/color
 xll = .15 & yll = .25
 xur = .95 & yur = .85
 c_charsize=1.0 
 charsize=1.0
 xcharsize=1.0
 ycharsize=1.0
 print,'Set plot device for ps'
endif else begin
 set_plot,'X'
 c_charsize=1.2
 charsize=1.4
 xcharsize=1.3
 ycharsize=1.3
 print,'Set plot device X'
endelse
c_color = intarr(17) & c_color[*]=!d.table_size-1

 mappos = [xll,yll,xur,yur]
 map_set,/cyl,/noborder,position=mappos

 levels = findgen(17)*10
 c_labels = intarr(17) & c_labels[*]=1
 print,'levels=',levels

 contour,tec,x,y,/noerase,position=mappos,xstyle=5,ystyle=5,$
   /fill,levels=levels
 contour,tec,x,y,/noerase,xstyle=1,$
   ystyle=1,yrange=[-90.,90.],xtitle='Local Time',ytitle='Latitude',$
   yticks=6,yminor=3,xticks=4,xminor=6,position=mappos,title=$
  'Total Electron Content (December 21, 12 UT, Solar Maximum)',$
   /follow,levels=levels,c_labels=c_labels,c_charsize=c_charsize,$
   charsize=charsize,xcharsize=xcharsize,ycharsize=ycharsize,$
   c_color=c_color

;map_set,/cyl,/noborder,position=mappos,/noerase
;map_continents,/noerase,color=!d.table_size-1

;Make magnetic equator thick

 contour,alat,x,y,/overplot,levels=0.,thick=4.,c_color=[!d.table_size-1]

if (keyword_set(ps)) then device,/close

end