diff --git a/mon_mean.r b/mon_mean_onept.r similarity index 57% rename from mon_mean.r rename to mon_mean_onept.r index ab89b6494f1583707813a29fa9e1ae4d80e01685..5d6711a28e2fc8d0734bb9df221c8b2f2116d5cf 100644 --- a/mon_mean.r +++ b/mon_mean_onept.r @@ -1,12 +1,13 @@ -mon_mean <- function(var,lat, lon) { +mon_mean_onept <- function(var, lat, lon) { #read in daily mean data, calculate monthly mean -var_onept <- var[lon,lat,] +#'var' is the variable read from netCDF +var_onept <- var[lon, 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 08f3b831d8fea2515636f5f8fedf80cfb39e5621..70cacb07df53ea818d41a9e6f790dfcd2e03602e 100644 --- a/run.r +++ b/run.r @@ -10,6 +10,6 @@ tmp <- read_netcdf(var_name, path, file) lat <- 120 lon <- 255 -monmean <- mon_mean(tmp, lat, lon) +monmean <- mon_mean_onept(tmp, lat, lon) message('Monthly mean at [',lon,',',lat,']: ',monmean)