<div>Hi all,</div><div>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...</div><div> </div><div>So would you please help me to find out the problems ? Thanks very mcuh for your precious time and kindness!</div><div> </div><div><div>Sincerely </div><div>Yi Lu</div><div> </div><div>;;;;;;;;;;; the part of NCL script to read in the precipitation ;;;;;;;;;; </div></div><div>tim = "1.grd"<br>setfileoption("bin","ReadByteOrder","LittleEndian")<br>rain =
fbindirread(tim,0,(/nlat,nlon/),"double")<br><br>nlat = 440<br>nlon =
700<br>r=new((/nlat,nlon/),"double")<br>k=0<br>do i = 0, nlat-1<br> do j = 0,
nlon-1<br> r(i,j) = rain(nlat-1-i,j)<br> k = k+1<br> end do<br>end
do<br><br>lon = fspan(70.05,139.95,nlon)<br>lon!0 = "lon"<br><a href="mailto:lon@units" target="_blank"><font color="#0066cc">lon@units</font></a> =
"degrees_east"<br><a href="mailto:lon@long_name" target="_blank"><font color="#0066cc">lon@long_name</font></a> = "longitude"<br><br>lat =
fspan(15.05,58.95,nlat)<br>lat!0 = "lat"<br><a href="mailto:lat@long_name" target="_blank">lat@long_name</a> = "latitude"<br><a href="mailto:lat@units" target="_blank">lat@units</a> = "degrees_north"<br><br>r!0 = "lat"<br>r!1 =
"lon"<br>r&lat = lat<br>r&lon = lon<br><a href="mailto:r@units" target="_blank">r@units</a> = "mm/hour"<br><a href="mailto:r@long_name" target="_blank">r@long_name</a> = "merged CMORPH data "<br><a href="mailto:r@_FillValue" target="_blank">r@_FillValue</a> = -999.00<br></div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div> </div><div> </div><div><sign signid="98"><div style="color: rgb(144, 144, 144); font-family: Arial Narrow; font-size: 12px;"><br><br><br><br></div><div style="color: rgb(0, 0, 0); font-family: Verdana; font-size: 14px;"><div style="overflow: hidden;"><table class="preview_table"><tbody><tr><td valign="middle"><img id="preview_qrcode" src="https://exmail.qq.com/cgi-bin/setting_qrcode_card?action=qrcode_card&alias=dg1225033@smail.nju.edu.cn&key=33b1e5708e6a55b5"></td><td valign="middle" style="color: rgb(160, 160, 160); font-size: 12px;"><h4 class="name" style="margin: 0px; color: rgb(0, 0, 0); line-height: 28px; font-size: 14px; font-weight: bold; -ms-zoom: 1;">易路</h4><p class="department" style="margin: 0px;">南大邮件系统/学生/博士生/12级博士生</p><p class="addr" style="margin: 0px; color: rgb(160, 160, 160); line-height: 22px;">南京市汉口路22号 </p></td></tr></tbody></table></div></div></sign></div><div> </div><div><includetail><!--<![endif]--></includetail></div>