From d36b3b9f126096a34b4045d79c113bc082044ce1 Mon Sep 17 00:00:00 2001 From: Colin Mahony Date: Sun, 5 Jan 2025 22:59:12 -0800 Subject: [PATCH] expand domain to -179E --- data_processing/Mosaic_NorAm_Implement.R | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/data_processing/Mosaic_NorAm_Implement.R b/data_processing/Mosaic_NorAm_Implement.R index d1f7b24..6e838cb 100644 --- a/data_processing/Mosaic_NorAm_Implement.R +++ b/data_processing/Mosaic_NorAm_Implement.R @@ -26,7 +26,7 @@ weightSurface <- function(lat, lon, ext){ # ext is c(lon.start, lon.end, lat.sta return(weights) } -studyarea <- ext(c(-170, -52, 14, 83.5)) +studyarea <- ext(c(-179.0625, -52, 14, 83.5)) # study area DEM (create, write, and read) # dir <- paste("//objectstore2.nrs.bcgov/ffec/Climatologies/PRISM_BC/", sep="") @@ -405,6 +405,8 @@ for(e in 1:3){ # Map showing the data sources #----------------------------- +plotarea <- ext(c(-170, -52, 14, 83.5)) + plot(gan.project, col="gray") comp <- dem values(comp) <- NA @@ -442,6 +444,7 @@ ocean[lat>47 & lat<49 & lon > (-123.5) & lon < (-122)] <- NA comp.final <- cover(ocean, comp) comp.final <- mask(comp.final, comp) +comp.final <- crop(comp.final, plotarea) plot(comp.final) # dir <- paste("//objectstore2.nrs.bcgov/ffec/Climatologies/PRISM_BC/", sep="")