NCL center_finite_diff_n 出错：Assignment type mismatch 跪求问大佬！！

猪猪2020

a = addfile("/home/yangyun/data_sets/pi/ne30/ocn/VVEL_regrid_.00101-50201.nc","r")
b = addfile("/home/yangyun/data_sets/pi/ne30/ocn/TEMP_regrid_.00101-50201.nc","r")
;47°-63°N 20°-45°W
VVEL0 = a->VVEL(75*12:375*12-1,0:19,46:64,314:341)
TEMP0 = b->TEMP(75*12:375*12-1,0:19,46:64,314:341)
lon = VVEL0&lon ; VVEL0和TEMP0两个nc文件dimsizes是相同的
lat = VVEL0&lat

;DJF
T=month_to_season (TEMP0, "DJF")
V=month_to_season (VVEL0, "DJF")

;dTdX
nlat=19
dlon = (lon(2)-lon(1))*0.0174533 ; convert to radians
                                    ; pre-allocate space
dTdX = new ( dimsizes(T), typeof(T), T@_FillValue)

do nl=0,nlat-1                    ; loop over each latitude
      dX = 6378388.*cos(0.0174533*lat(nl))*dlon  ; constant at this latitude
      dTdX(:,:,nl:nl,:) = center_finite_diff_n (T(:,:,nl:nl,:), dX , True,0,3)
end do

出错说：数据类型不匹配（T V 是海洋的栅格数据包含缺测值，即陆地上的值是缺测值）

流浪的心

把出错地方作用两边的变量用printVarSummary()函数输出，看看数据类型是否一致，若不一致可用数据类型转换的函数，转换一下就可以了！

猪猪2020

流浪的心 发表于 2022-7-5 08:27
把出错地方作用两边的变量用printVarSummary()函数输出，看看数据类型是否一致，若不一致可用数据类型转换 ...

大佬，我查看了一下数据类型，都是float

X1a0wu

dX也是float吗

流浪的心

主要是检查center_finite_diff_n()生产的是什么类型，你可以测试一下，即aaa = center_finite_diff_n (T(:,:,nl:nl,:), dX , True,0,3)，然后看一下aaa的类型是否是float！

guoguohh

数据类型一样的话可能你左右两边数据的维数不一样，输出dTdX(:,:,nl:nl,:)看看这个是不是3维的数据，你写了循环但是脚本没看到你的nlat用在哪里，最后你的dim又选择的是第四维，所以报错

井中月

好多数据的坐标变量会用double，你应该查看一下lon lat的类型

猪猪2020

具体原因不太知道，但是我写了下面的代码程序就不报错了dTdx(:,nl:nl,:) = tofloat(center_finite_diff_n (T(:,nl:nl,:), dx , True,0,2))
将输出的结果强制转换成float。目前没有发现结果有什么问题