;===================================================
;GASP2netcdf.pro
;transfer GASp binary data format to netcdf format
;written by Jian Zeng
;2006-06-21
;===================================================
@/net/orbit047L/home/jzeng/Netcdf/read_GASP_unscl.pro
@/net/orbit047L/home/jzeng/Tools/tools/convert_julianday.pro

infile='2006113191515_i16_US.all.aod.gz'
dir_in='/net/orbit111L/disk2/pub/zeng/GOES/AOD/'

lat_infile='/net/orbit111l/disk2/pub/aprados/ref/lat.dat'
lon_infile='/net/orbit111l/disk2/pub/aprados/ref/lon.dat'
nx_g=2000
ny_g=850

baod=read_GASP_unscl(dir_in+infile,nx_g=nx_g,ny_g=ny_g,baodstd=baodstd,$
               bsfc=bsfc,bch1=bch1,bcld=bcld,bmos=bmos,bsig=bsig,bmsk=bmsk,bcls=bcls)

outfile='2006113191515_i16_US.all.aod.ncdf'
yr='2006'
jday='113'
time='191515'
convert_julianday,yr,jday,mon,day
imonth='04'
iday='23'

fid=ncdf_create(outfile,/clobber)
ncdf_control,fid,/fill

;define global attributes
ncdf_attput,fid,'Title','GASP Aerosol Optical Thickness',/global
ncdf_attput,fid,'Platform','GOES-12',/global
ncdf_attput,fid,'Year',yr,/global
ncdf_attput,fid,'Month',imonth,/global
ncdf_attput,fid,'Day',iday,/global
ncdf_attput,fid,'Universal_Time',time,/global

;define dimensions
xid=ncdf_dimdef(fid,'x',2000)
yid=ncdf_dimdef(fid,'y',850)

;define variables
vid=ncdf_vardef(fid,'baod',[xid,yid],/byte)
offset=-0.5
scale_factor=1./100.
ncdf_attput,fid,vid,'long_name','Aerosol optical depth of all pixels scaled bytes'
ncdf_attput,fid,vid,'offset',offset,/float
ncdf_attput,fid,vid,'scale_factor',scale_factor,/float
ncdf_attput,fid,vid,'unit','none'
ncdf_control,fid,/endef
ncdf_varput,fid,vid,baod

ncdf_control,fid,/redef

vid=ncdf_vardef(fid,'baodstd',[xid,yid],/byte)
offset=0.
scale_factor=1./100.
ncdf_attput,fid,vid,'long_name','Aerosol optical depth standard deviation of 2 by 2'
ncdf_attput,fid,vid,'offset',offset,/float
ncdf_attput,fid,vid,'scale_factor',scale_factor,/float
ncdf_attput,fid,vid,'unit','none'
ncdf_control,fid,/endef
ncdf_varput,fid,vid,baodstd

ncdf_control,fid,/redef



ncdf_close,fid

end