diff --git a/mon_mean_areaavg.r b/mon_mean_areaavg.r
new file mode 100644
index 0000000000000000000000000000000000000000..6f246d57909caf132cf16149fafe4b3d180f5486
--- /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 70cacb07df53ea818d41a9e6f790dfcd2e03602e..48e9bfd06d6b904f91fe43f04f5e581288810601 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/'