pro pmaiac_daily_hdf

dir='/misc/mcc20/common/satellite/maiac/DataRelease/CMG/CMG_5km_Comb/AOT5km/2002/'
f0='MAIACCMG.2002183'
file_name=f0+'.hdf'

; Choose a variable: AOT, WV, FineModeFraction, CloudFraction
sds_name='AOT'

; Open the file
sd_id=hdf_sd_start(dir+file_name,/read)

; Find the index of the sds to read using its name
sds_index=hdf_sd_nametoindex(sd_id,sds_name)

; Select it
sds_id=hdf_sd_select(sd_id,sds_index)

; Get data set information including dimension information
hdf_sd_getinfo,sds_id,name=sds_name,natts=num_attributes,ndim=num_dims, $
dims=dim_sizes
nlat=dim_sizes[1]
nlon=dim_sizes[0]
print,'nlon=',nlon
print,'nlat=',nlat

; Read the data
hdf_sd_getdata,sds_id,data
print,'max min ',max(data),min(data)
data=data*0.0001
print,'max min ',max(data(where(data gt 0))),min(data(where(data gt 0)))
mk=where(data lt 0)
if total(mk) gt 0. then begin
 o=where(data lt 0.)
 data(o)=-999.9
endif
print,size(data)
print,'max min ',max(data),min(data)

; End, close
hdf_sd_endaccess,sds_id
hdf_sd_end,sd_id

; lon, lat
lon=findgen(nlon)*0.05-180.
lat=findgen(nlat)*0.05- 90.

aod=data
; Reverse latitude direction from N->S to S->N:
for j=0,nlat-1 do begin
 jx=nlat-1-j 
 aod(*,j)=data(*,jx)
endfor

; plot ----
lvl=findgen(10)*0.04
;lvl=[0,0.1,0.2,0.3,0.4,0.5,0.6,0.8,1.0,1.5]
plt='gif'
pmapgif,aod,lon,lat,levels=lvl,figname=f0
print,f0+'.'+plt

stop
end