[ncl-talk] How to correctly read in the merged CMORPH data usig NCL

[wen]

Hi. The following script is able to read the data. I hope it work for you. 
 


 setfileoption("bin","ReadByteOrder","littleEndian")
 mfile = "1.grd"
  rec = 0
  crain = fbindirread(mfile,rec,(/ny,nx/),"float")
  rec = 1
  gsamp = fbindirread(mfile,rec,(/ny,nx/),"float")





--

Best regards，
Guanhuan Wen

在 2017-04-30 12:07:44，"易路" <dg1225033 at smail.nju.edu.cn> 写道：

Hi all,
I downloaed the merged CMORPH data from China Meteorological Administration to get the hourly precipitation imformation, the suffix of the data file is grd which is a standard format of the tool GrADS. According to the associated ctl file of the grd file,  I tried to write the NCL script to read it but can't get the right data, and I couldn't find out the prolbems...
 
So would you please help me to find out the problems ? Thanks very mcuh for your precious time and kindness!
 
Sincerely
Yi Lu
 
;;;;;;;;;;; the part of NCL script to read in the precipitation ;;;;;;;;;; 
tim = "1.grd"
setfileoption("bin","ReadByteOrder","LittleEndian")
rain = fbindirread(tim,0,(/nlat,nlon/),"double")

nlat = 440
nlon = 700
r=new((/nlat,nlon/),"double")
k=0
do i  = 0, nlat-1
  do j = 0, nlon-1
   r(i,j) = rain(nlat-1-i,j)
        k = k+1
  end do
end do

lon = fspan(70.05,139.95,nlon)
lon!0 = "lon"
lon at units = "degrees_east"
lon at long_name = "longitude"

lat = fspan(15.05,58.95,nlat)
lat!0 = "lat"
lat at long_name = "latitude"
lat at units = "degrees_north"

r!0 = "lat"
r!1 = "lon"
r&lat = lat
r&lon = lon
r at units = "mm/hour"
r at long_name = "merged CMORPH data "
r at _FillValue = -999.00

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





| |
易路

南大邮件系统/学生/博士生/12级博士生

南京市汉口路22号 

|
 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170503/7c193563/attachment.html