目录

函数addfile和addfiles

创建NetCDF文件

函数addfile和addfiles

f &＃61; addfile(filename, status)

当status为"c"&＃xff0c;即创建一个新文件时&＃xff0c;如果同路径下已经存在同名文件&＃xff0c;则会出现报错信息

对于此情况&＃xff0c;须在使用函数addfile前将同名文件删除&＃xff0c;可通过函数systemfunc调用Linux命令rm实现

filename &＃61; "dataT.nc"
systemfunc("rm -f"&＃43;filename)
g &＃61; addfile(filename,"c")

当用户需要一次读入多个文件时&＃xff0c;则使用函数addfiles&＃xff0c;例如&＃xff1a;

list_of_files &＃61; systemfunc("ls ./*.nc")
file_list &＃61; addfiles(list_of_files, status)

在此操作后&＃xff0c;可通过[:]读取列表中多个文件的同一个变量

t &＃61; f[:]->TEMP

注意&＃xff0c;NCL有两种方式读取同一个变量&＃xff0c;分别是“cat"和"join”

默认是cat方式&＃xff0c;若采用join方式&＃xff0c;则要求所有文件中的同名变量数组结构必须相同

ListSetType(f, "join")
x &＃61; f[:]->TEMP

同时&＃xff0c;函数addfiles打开多个文件后&＃xff0c;也可读取某些指定文件中的变量。例如&＃xff1a;用户可能需要从列表变量中每隔一个文件读取变量&＃xff1a;
 

T_sec &＃61; f[0::2]->T

---

创建NetCDF文件

fin &＃61; addfile("tem_C.nc","r")
fout &＃61; addfile("tem_K.nc","c")
filedimdef(fout, "time", -1, True) ; 定义文件中的"time”维大小不固定
tK &＃61; fin->tsurf
tC &＃61; tK
tC &＃61; tK-273.15
tC&＃64;units &＃61; "degC";--输出变量tC至文件fout中&＃xff0c;并将该新写入的变量命名为tc
fout->tc &＃61; tC