计算nino指数 | 与冬季海温异常的相关/回归系数 nino=wgt_areaave(anmDJF,-5,5,-170,-120) ninoSD=nino/nino.std(dim='time') # simultaneous cor0 = xr.corr(ninoSD, anmDJF, dim="time") reg0 = xr.cov(ninoSD, anmDJF, dim="time")/nino
# == seasonal mean anmS = anm.rolling(time=3, center=True).mean('time') anmDJF=anmS.sel(time=slice(f'{ystr}-01-01',f'{yend}-12-01',12)) print(anmDJF) def wgt_areaave(indat, latS, latN, lonW, lonE): lat=indat.lat lon=indat.lon if ( ((lonW < 0) or (lonE < 0...