Btukfyl

I am trying to subset 8 netcdf files (one of them is here) in a foor loop according to a shorter period of time and then to save them as a new netcdf file.

I saw that other people already asked on how to subset netcdf files according to different time period (here or here) but once I do it on my netcdf files the "for loop" keeps on running without finishing (not even the first netcdf file) and I can't figure out why.

Here the code I use:

library(raster)

library(netcdf4)

library(lubridate)



#setting wd=indir containing netcdf files



setwd(indir)



files=list.files(pattern="nc")



for (j in seq_along(files)){



#setting wd containing netcdf files in the loop

setwd(indir)



b<-brick(files[j])



nc<-nc_open(files[j])



#variable

varname<-names(nc[['var']][3])

varunits <- ncatt_get(nc,varname,"units")[[2]]



lon<-ncvar_get(nc,"lon")



lat<-ncvar_get(nc,"lat", verbose = F)



time<-ncvar_get(nc, "time")



tunits <- ncatt_get(nc, "time", "units")[[2]]



dlname <- ncatt_get(nc, varname,"long_name")[[2]]



nc_close(nc)



#assigning a crs

proj4string(b)<-"+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"



#setting time as.Date 

tm<-ymd(getZ(b))



#setting time to rasterBrick

b<-setZ(b, tm)



# subsetting

b2<-subset(b, which(tm < as.Date('2006-01-01')))      



#setting wd where I want to save the "new" netcdf files

setwd(outdir)



writeRaster(b2, filename = paste0(varname, "_1971_2006_Noce.nc"),

format="CDF", varname=varname, varunit=varunits, longname=dlname,

xname="lon", yname="lat", zname="time", zunit=tunits, overwrite=TRUE)

}

Any help on how to get the loop working would be very much appreciated!