关于求SSTA的问题

Tobyright

我是ncl初学者，这个问题能很简单，但还是困扰了自己不短的时间，希望各位可以帮助解答一下。

我先求了1979-2019所有经纬的平均值，得出一个二维（lat，lon）的avesst，然后用sst与avesst相减，但由于sst是三维的，因此两者相减时会报错。用过conform函数后，把avesst变成三维的,但这时候avesst只有格点数，没有坐标范围，因为自己后面还需要求对区域的平均，所以没有区域范围的话后面也还是会报错，希望大家解答一下我的问题。

以下是我的ncl

f = addfile("E:\Desktop\sst_mon_mean.nc","r")
yrS    = 1979
yrE    = 2019
t2     = f->time
yr2    = cd_calendar(t2,-1)/100                  
iyr2   = ind(yr2.ge.yrS .and. yr2.le.yrE)
sst    = f->sst(iyr2,:,:)
avesst  = dim_avg_n_Wrap(sst(:,:,:),0)
avesst1 = conform(sst,avesst,(/1,2/))

Tobyright

                               
登录/注册后可看大图

Adiligirl

可以给变量赋坐标维
avesst!0 = "time"
avesst!1 = "lat"
avesst!2 = "lon"
avesst&lat = f->lat
avesst&lon = f->lon

Tobyright

Adiligirl 发表于 2020-3-31 13:57
可以给变量赋坐标维
avesst!0 = "time"
avesst!1 = "lat"

根据你的方法，问题解决了，非常感谢

Adiligirl

Tobyright 发表于 2020-3-31 14:39
根据你的方法，问题解决了，非常感谢

不客气，共同学习！