;
; reader for GOCART aerosol aot output from online ccm
;Huisheng Bian, 12/14/2011
;

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

;dir_in = '/misc/mcc20/common/arctas/GEOS5/HaoHe/'
dir_in = 'where is your data'

;cycle for months
;cycle for days

infile = dir_in+'bR_control.tavg2d_aer_x.'+yr+mm(nm)+dd(nd)+'_0900z.nc4'
;
   ;Does the filename exist?
    if(file_search(infile) eq '') then begin
     print, infile+': not found, continuing'
     continue
    endif
    print,infile

   ;Open nc4 file
   cdfid = ncdf_open(infile)

constid = ncdf_varid(cdfid, 'lon')
ncdf_varget,cdfid,constid,lon
lon = reform(lon)
nlon = n_elements(lon)

constid = ncdf_varid(cdfid, 'lat')
ncdf_varget,cdfid,constid,lat
lat = reform(lat)
nlat = n_elements(lat)

constid = ncdf_varid(cdfid, 'DUEXTTAU') ; dust AOD at 550nm 
ncdf_varget,cdfid,constid,duexttau
duexttau= reform(duexttau)
print,'duexttau=',max(duexttau(*,*,*)), min(duexttau(*,*,*))

constid = ncdf_varid(cdfid, 'DUEXTTFM') ; dust AOD at 550nm for PM 1.0um
ncdf_varget,cdfid,constid,duexttfm
duexttfm= reform(duexttfm)
print,'duexttfm=',max(duexttfm(*,*,*)), min(duexttfm(*,*,*))

constid = ncdf_varid(cdfid, 'DUEXTT25') ; dust AOD at 550nm for PM 2.5um
ncdf_varget,cdfid,constid,duextt25
duextt25= reform(duextt25)
print,'duextt25=',max(duextt25(*,*,*)), min(duextt25(*,*,*))

constid = ncdf_varid(cdfid, 'DUSMASS') ; dust surface mass concentration kg/m3
ncdf_varget,cdfid,constid,dusmass
dusmass= reform(dusmass)
print,'dusmass=',max(dusmass(*,*,*)), min(dusmass(*,*,*))

constid = ncdf_varid(cdfid, 'DUCMASS') ; dust column mass density kg/m2
ncdf_varget,cdfid,constid,ducmass
ducmass= reform(ducmass)
print,'ducmass=',max(ducmass(*,*,*)), min(ducmass(*,*,*))

ncdf_close,cdfid

; end cycle days
; end cycle months

end