|
我有很多nc文件,每个文件中我需要提取的是lswt这个变量的数值,对于单个nc文件,我用来提取lswt的代码是这样的:
> library(ncdf)
> arclake=open.ncdf("ALID0151_PLREC9D_TS366LM.nc")
> lswt=get.var.ncdf(arclake,"lswt")
> library(xlsx)
> write.xlsx(lswt-273,"result.xlsx") #nc给出的lswt的单位是开氏温度,-273是为了划算成摄氏度,最后将结果输出成xlsx.
结果是正确的。
不过我现在有很多个不同的nc文件,都要进行相同的提取lswt的操作,并将每个nc文件提取出来的结果放在同一个xlsx中,应该怎样进行操作呢?
希望大家指点>>>.......
谢谢.
|
|