NCL多文件读写+单位改变+得出nino3.4

（注：本文仅为个人笔记）

多文件读写

fns = systemfunc(“ls /data/ACCESS_rlilp1_1/rcp45_???.nc”)
infile = addfiles(fns,“r”)
ListSetType(infile,“cat”) ;将所有年份rcp45资料整合到一个文件，维度不变
st = infile[:]->tos ;将整合好的文件中tos输出为st

单位改变

sst = st - 273.15 ;将K转变为C

sst!0 = “time” ;when doing simple calculation, build an array
sst&time = st&time
sst!1 = “lat”
sst&lat = st&latitude
sst!2 = “lon”
sst&lon = st&longitude

sst@units = “C”
sst@long_name = “sea surface temperature[C]”

Nino3.4

latS = -70.0
latN = 70.0
lonL = 30.0
lonR = 290.0
yrStrt = 2006
yrLast = 2100
yrClimStrt = 2010
yrClimLast = 2030
time = infile[:]->time
lat = in