From 080e54a73455f5ad806118cb83a6cbaed9f5ace9 Mon Sep 17 00:00:00 2001
From: aho <an.ho@bsc.es>
Date: Mon, 25 Mar 2019 10:55:27 +0100
Subject: [PATCH] small change to run.r

---
 mon_mean_areaavg.r | 16 ++++++++++++++++
 run.r              |  2 +-
 2 files changed, 17 insertions(+), 1 deletion(-)
 create mode 100644 mon_mean_areaavg.r

diff --git a/mon_mean_areaavg.r b/mon_mean_areaavg.r
new file mode 100644
index 0000000..6f246d5
--- /dev/null
+++ b/mon_mean_areaavg.r
@@ -0,0 +1,16 @@
+mon_mean_areaavg <- function(var, latmin, latmax, lonmin, lonmax) {
+
+#read in daily mean data, calculate monthly mean
+#read in lat and lon range, calculate area mean
+#'var' is the variable read from netCDF
+var_sub <- var[lonmin:lonmax, latmin:latmax, ]
+lat_sub <- 
+cos_lat <- 
+monmean <- mean(var_onept)
+#monmean <- apply(as.matrix(var_onept), 3, mean)
+
+#Outputs
+return(monmean)
+
+}
+
diff --git a/run.r b/run.r
index 70cacb0..48e9bfd 100644
--- a/run.r
+++ b/run.r
@@ -1,6 +1,6 @@
 #Load function files
 source("read_netcdf.r")
-source("mon_mean.r")
+source("mon_mean_onept.r")
 
 var_name <- 'tas'
 path <- '/esarchive/recon/ecmwf/erainterim/daily_mean/tas/'
-- 
GitLab