ncl批量读取多个nc文件

diri = "/mnt/g/data/3/"
files = systemfunc("ls " + diri + "*nc4") ;批量读取文件路径
f1 = addfiles(files,"r")
ListSetType (f1, "cat")
lat = f1[0]->lat({20:30})
lon = f1[0]->lon({95:110})
time = f1[0]->time
pre = f1[:]->precipitation(:,{95:110},{20:30})