Description Usage Arguments Value Note Author(s) References Examples
View source: R/F_plot_meas_sim.R
Plots curves for comparison between measures and simulations (hourly interpolations)
1 2 3 | plot_meas_sim(meas, sim, series = NULL, chart.start,
chart.end, date.format = "ymd", missing_code = NA,
wait = 1, plot.leg = TRUE, leg.pos = "bottomright")
|
meas |
measured hourly values file (table), where the first column is the series' ID |
sim |
simulated hourly list (output of Th_int_series) |
series |
names of the seris to plot. If |
chart.start |
start date for the plotting. Format example: "1Jan2000" |
chart.end |
end date for the plotting. Format example: "1Jan2000" |
date.format |
input date format for measurements
(formats for function |
missing_code |
code (either real or character) for
missing values in measurements. Default is |
wait |
lag time (seconds) between plot appearance on the screen (default is 1 second) |
plot.leg |
logical: if |
leg.pos |
position of legends (only if
|
A plot with two curves: measured values and hourly interpolations
If daily minimum and maximum values are the absolute ones for each day, the interpolated curve will be generally either higher (maximum) or lower (minimum) than hourly measurements, which are mean hourly values (hence, slightly lower and higher than the absolute daily ones, respectively). This may explain an apparent mismatch between the two curves.
Emanuele Eccel, Emanuele Cordano emanuele.eccel@iasma.it
Eccel, E., 2010: What we can ask to hourly temperature recording. Part II: hourly interpolation of temperatures for climatology and modelling. Italian Journal of Agrometeorology XV(2):45-50 http://www.agrometeorologia.it/documenti/Rivista2010_2/AIAM%202-2010_pag45.pdf,www.agrometeorologia.it
See also: Eccel, E., 2010: What we can ask to hourly temperature recording. Part I: statistical vs. meteorological meaning of minimum temperature. Italian Journal of Agrometeorology XV(2):41-43. http://www.agrometeorologia.it/documenti/Rivista2010_2/AIAM%202-2010_pag41.pdf
1 2 3 4 5 | data(Trentino_hourly_T)
stations <- c("T0001","T0010","T0129")
plot_meas_sim(meas = h_d_t, sim = Th_int_list, series=stations,
missing_code = -999.9, chart.start = "1Feb2004",
chart.end = "29Feb2004", leg.pos = "top")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.