小白刚开始接触NCL，希望有好心人能解惑

songyy

代码：
begin
        f=addfile("d:/1960-2015/var1continous-summer-rain652stn_1960_2015_daily.nc","r")
        var = f->var3
       
        wks = gsn_open_wks("pdf","rain")
        res=True
        res@gsnLeftString=""
        res@gsnRightString=""

        res@pmTickMarkDisplayMode="Always"
       
        res@mpMinLatF=10
        res@mpMaxLatF=60
        res@mpMinLonF=90
        res@mpMaxLonF=180
       
        res@mpFillOn=False
        res@mpGridAndLimbOn=True
        res@mpGridLineThicknessF=0.5
        res@mpGridLineDashPattern=2
        res@mpGridSpacingF=10
       
        res@cnLinesOn=False
       
        res@cnLevelSelectionMode="ExplicitLevels"
        res@cnLevels=(/2,4,6,8,10/)
        res@cnFillColors=(/2,4,6,8,11,13/)
       
        plot =gsn_csm_contour_map_ce(wks,var,res)
end

杨松

没有坐标系？lon lat 你看看数据有没有问题

来自安卓APP客户端

songyy

杨松 发表于 2021-2-9 17:58
没有坐标系？lon lat 你看看数据有没有问题

数据没有问题。我在想是不是这么读变量有问题

穿越的小瓶子

你的数据读取没有错，但读出来设置一下var&var=var , var@units =" 输入单位 "，使读出来的变量有单位和坐标对应

songyy

穿越的小瓶子 发表于 2021-2-9 22:31
你的数据读取没有错，但读出来设置一下var&var=var , var@units =" 输入单位 "，使读出来的变量有单位和坐 ...

还是不行啊，按照您说的更改后会出现   fatal:(var) is not a named dimension in variable (var).

t_imba

这个是以站点数据存储的NC资料。不是格点形式的资料。你这写法读取不到经纬度信息，要想想怎么插值

songyy

t_imba 发表于 2021-2-10 14:54
这个是以站点数据存储的NC资料。不是格点形式的资料。你这写法读取不到经纬度信息，要想想怎么插值

啊，这么复杂吗，还要插值

songyy

t_imba 发表于 2021-2-10 14:54
这个是以站点数据存储的NC资料。不是格点形式的资料。你这写法读取不到经纬度信息，要想想怎么插值

我之前试过用mask_data插值，但是运行的时候运行不了，说是识别不了这个函数