pro read_eg_ecmwf,lat,lon,t,q
;Code to read example binary ECMWF fields
;Create arrays to read data into
lat=fltarr(161)
lon=fltarr(281)
q=fltarr(281,161,21,4);21 levels, 4 times
t=fltarr(281,161,21,4)
;Read ECMWF data (readu is read binary)
; Temperature, Kelvin
openr,12,'/nfs/see-fs-02_users/earnadr/public_html/Teaching/IDL_course/data/ECMWF/20070626_t.dat'
readu,12,t
close,12
;kg/kg water vapour mixing ratio
openr,12,'/nfs/see-fs-02_users/earnadr/public_html/Teaching/IDL_course/data/ECMWF/20070626_q.dat'
readu,12,q
close,12
openr,12,'/nfs/see-fs-02_users/earnadr/public_html/Teaching/IDL_course/data/ECMWF/20070626_lat.dat'
readu,12,lat
close,12
openr,12,'/nfs/see-fs-02_users/earnadr/public_html/Teaching/IDL_course/data/ECMWF/20070626_lon.dat'
readu,12,lon
close,12
END
