; STELAB density

 restore,'STELAB/stelab_v_d_1au.sav'

 dmax=10.
 dmin=0.
 mappart=d
 map=fltarr(360,180)
 for i = 29,148 do begin
  map(*,i+1)=mappart(*,i-29)
 endfor
 usetitle='STELAB Density (cm^-3) at 1.00 AU'
 plottitle='temp'
 loadct,3
 TVLCT, r, g, b, /Get
 TVLCT, Reverse(r), Reverse(g), Reverse(b)
 blackuse=255

; the following lines come from a subroutine used for all maps

 set_plot,'ps'

 !p.multi=[0,1,1]
 !p.font=0
 !y.minor=-1
 !x.minor=-1
 !x.margin=[10,3]
 !p.charsize=1.1
 device, bits=16, /color, /portrait, /encapsulated, $
 file='PLOTS/'+plottitle+'_whicarr.eps', xsize=24, ysize=12

 longitude=[ '0 ', ' ', ' ','60', ' ', ' ', '120', ' ', ' ', '180', $
            ' ', ' ', '240', ' ', ' ', '300', ' ' , ' ',  '   360' ]
 latitude=[-90,-45,0,45,90]

 nlev=20.
 deltlev=(dmax-dmin)/(nlev-1.)
 levuse=fltarr(nlev)
 for i = 0,nlev-1 do levuse(i)=dmin+i*deltlev

 contour,map,$
  color=blackuse,/fill,$
  title=usetitle,$
  yticks=4 , xticks=18, xtickname=longitude, $
  ticklen=0.02, $
  ytickname=latitude,$
  xs=1,ys=1,$
  levels=levuse

  colorbar,/vertical,position=[0.05,0.18,0.07,0.88],minrange=dmin,maxrange=dmax,color=blackuse,divisions=divuse

device,/close

end