From 0214a522a35937dbc13ae61598df8c5f515268d5 Mon Sep 17 00:00:00 2001 From: leneoMET Date: Thu, 15 Sep 2022 08:19:50 +0200 Subject: [PATCH] Verification scripts --- .../daily_obs_to_sqlite.R | 0 {Rscripts => verification}/forecast_max_min.R | 42 +++++++++---------- 2 files changed, 21 insertions(+), 21 deletions(-) rename {Rscripts => verification}/daily_obs_to_sqlite.R (100%) rename {Rscripts => verification}/forecast_max_min.R (62%) diff --git a/Rscripts/daily_obs_to_sqlite.R b/verification/daily_obs_to_sqlite.R similarity index 100% rename from Rscripts/daily_obs_to_sqlite.R rename to verification/daily_obs_to_sqlite.R diff --git a/Rscripts/forecast_max_min.R b/verification/forecast_max_min.R similarity index 62% rename from Rscripts/forecast_max_min.R rename to verification/forecast_max_min.R index 3faa580..b9326c5 100644 --- a/Rscripts/forecast_max_min.R +++ b/verification/forecast_max_min.R @@ -39,25 +39,25 @@ T2m_to_min_max_mean <- function(fcst, fcst_func){ #################### ## EXAMPLE OF USE ## #################### - -# READ FORECAST TEMPERATURE (3 HOURLY) -t2m <- read_point_forecast( - start_date = 20200401, - end_date = 20200531, - lead_time = seq(0,72,3), - fcst_model = "ec", - fcst_type = "det", - parameter = "T2m", - by = "1d", - file_path = "~/Documents/BMD2022/data/FCTABLE/deterministic" # must correspond to path.out in fcsttosqlite.R -) - -## CONVERT FORECAST TEMPERATURE FROM KELVIN TO DEGREES CELSIUS -t2m <- scale_point_forecast(t2m, -273.15, new_units = "degC") - -## REMOVE MISSING DATA FROM EC -t2m$ec <- t2m$ec %>% filter(ec_det<10^36) - -## FIND MAXIMUM TEMPERATURE -tmax.fcst <- T2m_to_min_max_mean(t2m, "max") +# +# # READ FORECAST TEMPERATURE (3 HOURLY) +# t2m <- read_point_forecast( +# start_date = 20200401, +# end_date = 20200531, +# lead_time = seq(0,72,3), +# fcst_model = "ec", +# fcst_type = "det", +# parameter = "T2m", +# by = "1d", +# file_path = "~/Documents/BMD2022/data/FCTABLE/deterministic" # must correspond to path.out in fcsttosqlite.R +# ) +# +# ## CONVERT FORECAST TEMPERATURE FROM KELVIN TO DEGREES CELSIUS +# t2m <- scale_point_forecast(t2m, -273.15, new_units = "degC") +# +# ## REMOVE MISSING DATA FROM EC +# t2m$ec <- t2m$ec %>% filter(ec_det<10^36) +# +# ## FIND MAXIMUM TEMPERATURE +# tmax.fcst <- T2m_to_min_max_mean(t2m, "max")