【求助】ncl画grd数据问题

电波系主任

在网上下载了OMI全球NO2柱浓度观测数据，分辨率为0.125*0.125。没有ctl文件，但是grd文件开头有如下说明：
ncols 2880
nrows 1440
xllcorner -180.0
yllcorner  -90.0
cellsize 0.125
nodata_value -999
version 2.0
之后是个2880*1440的数组。我想画中国地区的浓度图
参考这个脚本 http://bbs.06climate.com/forum.php?mod=viewthread&tid=49716
我的ncl脚本中坐标定义如下：
path = "no2_201806.grd"
dims = (/2880,1440/)
f=fbindirread(path,0,dims,"float")
a=new((/2880,1440/),"float")
do i=0,2879
do j=0,1439
a(i,j)=where(f(i,j).eq.-999,0,f(i,j))
end do
end do   
a!0="lon"
a!1="lat"
a&lon=fspan(-180,180,2880)
a&lat=fspan(-90,90,1440)
a&lon@unit="degrees_east"
a&lat@unit="degrees_north"

printVarSummary(a)
wks = gsn_open_wks("x11","no2")


选取画图范围
res@gsnAddCyclic        = False
res@mpLimitMode       = "LatLon"
res@mpMinLatF         = 10      
res@mpMaxLatF         = 60
res@mpMinLonF         = 70
res@mpMaxLonF         = 136

在运行后没有报错，做图范围也是对的，但是提取出来的数据都在1以下，明显不对。
不知道是哪里出了问题，求大家帮忙看一下。

电波系主任

这个附件我不知道怎么设置下载不扣金钱，不好意思

覃建明

楼主搞定了么

有鱼吃不吃

楼主，你好，请问你的问题解决了吗？