NCL新手 插值出错

葡萄不是哭包

我在ncl交互模式下根据官网的例子，尝试插值，未果。。
报错内容是：
fatal:No coordinate variable exists for dimension (lon) in variable (bcavgday)
fatal:(lon) is not coordinate variable in variable(bcavgday).
fatal:["Execute.c":8638]:Execute: Error occurred at or near line 38

具体过程如下：

Ncl 0 > f=addfile("/home/zhangyj/MERRA2_300.tavg1_2d_aer_Nx.20071201.SUB.nc","r")

ncl 1> bc=f->BCEXTTAU

ncl 2>  N = num(.not.ismissing(bc))

ncl 8> bcT=transpose(bc)  %转置以后的结果，方便使用dim_avg

                               
登录/注册后可看大图

ncl 10> bcavgday=dim_avg(bcT)%对数组最右边的一维求平均值

ncl 19> bcavgday!0="lon" %这里明明对坐标变量命名了lon 和 lat为什么还会报错说是lon维度的坐标变量不存在呢？？
ncl 20> bcavgday!1="lat"

ncl 22> bcavgday@_FillValue=1e+15                                          
ncl 23> dimbcavgday = dimsizes(bcavgday)                                   
ncl 24> nlon=dimbcavgday(0)                                                
ncl 25> nlat=dimbcavgday(1)                                                
ncl 26> diro="./home/zhangyj/"                                             
ncl 27> filo="daily.dat"                                                   
ncl 28> ptho=diro+filo                                                     
ncl 29> mlon=576                                                           
ncl 30> mlat=361                                                           
ncl 31> klon=ispan(0,mlon-1,1)*0.625                                       
ncl 32> klat=ispan(0,mlat-1,1)*0.5                                         
ncl 33> model=linint2(bcavgday&lon,bcavgday&lat,bcavgday,True,klon,klat,0)
fatal:No coordinate variable exists for dimension (lon) in variable (bcavgday)
fatal:["Execute.c":8638]:Execute: Error occurred at or near line 33

ncl 34> printVarSummary(bcavgday)                   %查看该变量的信息，显示坐标变量名称正确了呀。。。                      

Variable: bcavgday
Type: float
Total Size: 831744 bytes
            207936 values
Number of Dimensions: 2
Dimensions and sizes:    [lon | 576] x [lat | 361]   
Coordinates:
Number Of Attributes: 1
  _FillValue :    1e+15

求好心人解答~~~~~~~

stepdance

print(lon)出来试试，看看有没有经度数据

我猜是第10行dim_avg出问题了，这个函数算平均的时候，不会保留原始数据的元信息，最好用dim_avg_Wrap

irides

没有对lon和lat赋值吧，里面没有经纬度信息

Adam-Lee

还是写成批处理脚本好看些，还能随便改

Yang.hahha

我新装了一下6.4版的ncl，在使用函数的时候也出现了这个8638这个错误，请问楼主解决了吗。然后你加个赋值bcavgday!0    。bcavgday&0=什么什么，就不会说不存在了

Yang.hahha

我新装了一下6.4版的ncl，在使用函数的时候也出现了这个8638这个错误，请问楼主解决了吗。然后你加个赋值bcavgday!0    。bcavgday&0=什么什么，就不会说不存在了