Commit ea536ffa authored by DORAY's avatar DORAY
Browse files

sat pelgas res update

parent 7260e55e
......@@ -202,8 +202,6 @@ cat("Raster stack saved as:", filei2, "\n")
#sschla.rasterStack=stack(filei2)
#rm('sschla.rasterStack')
# # Resample by week
# plot(sschla.rasterStack)
# #rasteri2=aggregate(sschla.rasterStack)
......@@ -238,6 +236,9 @@ if (exists('sst.rasterStack')){
lidik2=NULL
}
rm('sst.rasterStack')
lidik2=NULL
jj=seq(from=1,to=366,by=1)
jj=seq(from=131,to=366,by=1)
year=2000:2021
......@@ -312,14 +313,14 @@ filei2 = paste(path.export.grid.sat,'SST/raster/',
writeRaster(sst.rasterStack,filei2,overwrite=TRUE)
# import raster
echosonde='/run/user/1000/gvfs/smb-share:domain=ifr,server=echosonde,share=data,user=mdoray/'
donnees2='/run/user/1000/gvfs/smb-share:domain=ifr,server=nantes,share=donnees2,user=mdoray/'
path.export.grid.sat=paste(echosonde,'Satellite/Atlantic/',sep='')
year=2000:2017
idi=paste('SST',paste(min(year),max(year),sep='-'),sep = "_")
filei2 = paste(path.export.grid.sat,'SST/raster/',paste('rasterStack',idi,sep='_'),'.grd',sep='')
# read in a raster stack
sst.rasterStack=stack(filei2)
# echosonde='/run/user/1000/gvfs/smb-share:domain=ifr,server=echosonde,share=data,user=mdoray/'
# donnees2='/run/user/1000/gvfs/smb-share:domain=ifr,server=nantes,share=donnees2,user=mdoray/'
# path.export.grid.sat=paste(echosonde,'Satellite/Atlantic/',sep='')
# year=2000:2017
# idi=paste('SST',paste(min(year),max(year),sep='-'),sep = "_")
# filei2 = paste(path.export.grid.sat,'SST/raster/',paste('rasterStack',idi,sep='_'),'.grd',sep='')
# # read in a raster stack
# sst.rasterStack=stack(filei2)
# 4. Average rasters over time -------------
# ************************************************
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment