(注:本文仅为个人笔记)
多文件读写
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