;
c_month = ' '
;
print, ' Enter the TWO-DIGIT month to plot'
read, c_month
;
set_plot,'ps'
;;;device,filename='c:\sep_1987.ps',/landscape,/palatino,/inches, $
;;;    xsize=7.0, ysize=10.0, xoffset=0.5, yoffset=0.5
device,filename='c:\isccp\month_'+c_month+'.ps',/landscape,/palatino
;;;    /inches, ysize=10.0, yoffset=0.5
;;;window, 0, xsize=900, ysize=500
openr,3,'c:\isccp\month_'+c_month+'.isc'
isccp = fltarr(864,432)
lat_idx = 0
lon_idx = 0
cloud = 0.0
rlat=0.0
rlon=0.0
lon_ii = 0
lat_jj = 0
while not eof(3) do begin
   readf, 3, rlat, rlon, cloud, format='(1x, f6.2, 1x, f7.2, 1x, f9.3)'
   lat_idx = 6 * (round((91.25 + rlat) / 2.5) - 1)
   lon_idx = 6 * (round((181.25 + rlon) / 2.5) - 1)
   for ii = 0, 5 do begin
      lon_ii = lon_idx + ii
      for jj = 0, 5 do begin
         lat_jj = lat_idx + jj
         isccp(lon_ii,lat_jj) = 0.0
         if cloud gt 0.0 then isccp(lon_ii,lat_jj) = cloud
      endfor
   endfor
endwhile
;
close,3
isccp(0,0) = 0.0
isccp(863,431) = 1023.0
isccp_im = bytscl(isccp)
;;;tv, isccp_im, 18, 34
;;;TV, ISCCP_IM
;
title_month = ' '
print,' Enter the month to use for the title'
read, title_month
;
map_set,0.0,0.0,0.0,/cylindrical,limit=[-88.75,-180.0,88.75,180.0], $
  /noerase, /noborder, xmargin=[-68.0,-68.0], ymargin=[0.0,0.0], $
  title=title_month+' 1984-1990 Average Cloud Pressure Heights'
NEWIM = MAP_IMAGE(ISCCP_IM, 18, 34, /WHOLE_MAP, $
                  LATMIN=-88.75,LATMAX=88.75,LONMIN=-180.0,LONMAX=180.0)
TV, NEWIM, 18, 34
map_grid, color=255, latdel=10.0, londel=10.0
map_continents, color=255
device, /close_file
stop
end