循环遍历 R 中的 nc 文件

Question

Good morning everyone, I am currently using the code written by Antonio Olinto Avila-da-Silva on this link: https://oceancolor.gsfc.nasa.gov/forum/oceancolor/topic_show.pl?tid=5954 It allows me to extract data of type sst/chlor_a from nc file. It uses a loop to create an excel file with all the data. Unfortunately, I noticed that the function only takes the first data file in the loop. Thus, I find myself with 20 times the same data in a row in my excel file. Does anyone have a solution to make this loop work properly?

Answer 1

I would first check out that these two lines contain all the files you are expecting:

(f <- list.files(".", pattern="*.L3m_MO_SST_sst_9km.nc",full.names=F))
(lf<-length(f))

And then there's a bug in the for-loop. This line:

data<-nc_open(f)

Needs to reference the iterator i , so change it to something like this:

data<-nc_open(f[[i]])

It appears both scripts have this same bug.