;===================================================
;Dust2netcdf.pro
;transfer GASp binary data format to netcdf format
;written by Jian Zeng
;2006-06-21
;===================================================
@read_dust_binary_all.pro
@convert_julianday.pro


PRO Dust2netcdf_2,yr,jday,imonth,iday,time,dir_in,dir_out

;infile1='MYD.comb.20120320.hr19.v6.3.4.dat'
infile1='MYD.comb.'+yr+imonth+iday+'.hr'+time+'.v6.3.4.dat'
;dir_in='/data/test/Dust/output/20120320/Output/'
infile=dir_in+infile1
print,'infile=',infile

f1=file_search(infile,count=nfl)

if nfl gt 0 then begin

nx_g=601
ny_g=251

;baod=read_dust_binary_all(dir_in+infile,nx_g=nx_g,ny_g=ny_g,lon=lon,lat=lat,combaod=combaod,$
;               combconc=combconc)

;read_dust_binary_all(dir_in+infile,nx_g=nx_g,ny_g=ny_g,lon=lon,lat=lat,combaod=combaod,combconc=combconc)
;read_dust_binary_all
nx=601L
ny=251L
lon=fltarr(nx,ny)
lat=fltarr(nx,ny)
combaod=fltarr(nx,ny)
combconc=fltarr(nx,ny)
openr,iun,infile,/get_lun,/f77_unformatted
readu,iun,nx,ny,lon,lat,combaod,combconc
close,iun
free_lun,iun

outfile1='MYD.comb.'+yr+imonth+iday+'.hr'+time+'.v6.3.4.ncdf'
;outfile='MYD.comb.20120320.hr19.v6.3.4.dat.ncdf'
;yr='2012'
;jday='080'
;time='19'
;convert_julianday,yr,jday,mon,day
;imonth='03'
;iday='20'

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

;define global attributes
ncdf_attput,fid,'Title','DUST Modis ',/global
ncdf_attput,fid,'Platform','Modis',/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',601)
yid=ncdf_dimdef(fid,'y',251)

;define variables
;vid=ncdf_vardef(fid,'combaod',[xid,yid],/byte)
vid=ncdf_vardef(fid,'combaod',[xid,yid],/float)
offset=-0.5
scale_factor=1./100.
ncdf_attput,fid,vid,'long_name','DUST Modis combaods'
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,combaod

ncdf_control,fid,/redef
;

vid=ncdf_vardef(fid,'combconc',[xid,yid],/float)
offset=0.
scale_factor=1./100.
ncdf_attput,fid,vid,'long_name','DUST Modis combconc '
ncdf_attput,fid,vid,'offset',offset,/float
ncdf_attput,fid,vid,'scale_factor',scale_factor,/float
ncdf_attput,fid,vid,'unit','Kilogram/Cube Kilometre'
ncdf_control,fid,/endef
ncdf_varput,fid,vid,combconc

ncdf_control,fid,/redef



ncdf_close,fid

endif
end