## ISSUES/TODO
## if eg prev_pr is used than the map is moved one gridcell upwards
## vapor pressure is too high???
## swdown can only be done globally: fix this!
## check if start yday is correct
## check base_offset it is now right for 6hourly calculations but maybe not for others like 3 hourly
## unitconversion
## Cleanup
rm(list=ls(all=TRUE))
## Load the library
library(metGeneratoR)
## Override standard settings
mgsetLonlatbox(c(-179.75, 179.75, -89.75, 89.75))
# mgsetLonlatbox(c(92.25,110.25,7.25,36.25))
mgsetPeriod(startdate = "1998-1-1", enddate = "1998-1-1")
mgsetInVars(list(
pr = list(ncname = "pr", filename = "../example_data4mtclim/Global/pr_day_HadGEM2-ES_historical_r1i1p1_EWEMBI_landonly_1998.nc"),
tasmin = list(ncname = "tasmin", filename = "../example_data4mtclim/Global/tasmin_day_HadGEM2-ES_historical_r1i1p1_EWEMBI_landonly_1998.nc"),
tasmax = list(ncname = "tasmax", filename = "../example_data4mtclim/Global/tasmax_day_HadGEM2-ES_historical_r1i1p1_EWEMBI_landonly_1998.nc"),
psurf = list(ncname = "ps", filename = "../example_data4mtclim/Global/ps_day_HadGEM2-ES_historical_r1i1p1_EWEMBI_landonly_1998.nc"),
# relhum = list(ncname = "hurs", filename = "../example_data4mtclim/Global/hurs_day_HadGEM2-ES_historical_r1i1p1_EWEMBI_landonly_1998.nc"),
qair = list(ncname = "Qair", filename = "../example_data4mtclim/Qair_1998.nc"),
swdown = list(ncname = "rsds", filename = "../example_data4mtclim/Global/rsds_day_HadGEM2-ES_historical_r1i1p1_EWEMBI_landonly_1998.nc"),
lwdown = list(ncname = "rlds", filename = "../example_data4mtclim/Global/rlds_day_HadGEM2-ES_historical_r1i1p1_EWEMBI_landonly_1998.nc"),
wind = list(ncname = "sfcWind", filename = "../example_data4mtclim/Global/wind_day_HadGEM2-ES_historical_r1i1p1_EWEMBI_landonly_1998.nc")
))
mgsetInDt(inDt = 24)
mgsetOutDt(outDt = 6)
mgsetOutVars(c("swdown", "lwdown", "tas", "pr", "psurf", "wind", "vp"))
# mgsetOutVars(c( "swdown"))
mgsetOutVars(c( "tas"))
# mgsetOutVars(c( "pr"))
mgsetOutName("./output/<VAR>/<VAR>_6hourly2_HadGEM2-ES_historical_r1i1p1_EWEMBI_landonly_<SYEAR>_local.nc")
## Run metGen
metGenRun()
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.