Example
using PML, Ipaper
+Calibrate PMLV2 model using FLUXNET data · PML.jl Calibrate PMLV2 model using FLUXNET data
using PML, Ipaper
df_out, df, par = deserialize(file_FLUXNET_CRO)
df.GPP_obs = df.GPPobs
@@ -29,31 +29,13 @@
(; var="ET", GOF(df.ET_obs, df_out.ET)...),
(; var="GPP", GOF(df.GPP_obs, df_out.GPP)...)] |> DataFrame
DataFrame(gof)
2×10 DataFrameRow var NSE R2 KGE R RMSE MAE bias bias_perc n_valid String Float64 Float64 Float64 Float64 Float64 Float64 Float64 Float64 Int64 1 ET 0.663915 0.664915 0.759815 0.815423 1.25573 0.887616 -0.0139188 -0.463569 263 2 GPP 0.535 0.61344 0.46745 0.783224 3.33344 2.65305 1.02043 27.1677 260
绘图
using Plots
-gr(framestyle=:box)
-t = df.date
275-element Vector{Dates.Date}:
- 2009-01-01
- 2009-01-09
- 2009-01-17
- 2009-01-25
- 2009-02-02
- 2009-02-10
- 2009-02-18
- 2009-02-26
- 2009-03-06
- 2009-03-14
- ⋮
- 2014-10-24
- 2014-11-01
- 2014-11-09
- 2014-11-17
- 2014-11-25
- 2014-12-03
- 2014-12-11
- 2014-12-19
- 2014-12-27
inds = 1:length(t)
inds = 1:46
+gr(framestyle=:box, titlefontsize=12)
+t = df.date
+inds = 1:46*1
-function plot_var(var; label=string(var), title=string(var), kw...)
- plot(t[inds], df_out[inds, var]; label, title, kw...)
+function plot_var(var; label=string(var), title=string(var),
+ data=df_out, scale=1, kw...)
+ plot(t[inds], data[inds, var] * scale; label, title, kw...)
end
function plot_var!(p, var; label=string(var), kw...)
plot!(p, t[inds], df_out[inds, var]; label, kw...)
@@ -67,10 +49,12 @@
plot(
p1,
plot_var(:GPP; title="GPP (gC m-2 d-1)", label="GPP"),
+ plot_var(:Eeq; title="Eeq (mm/d)", label="Eeq"),
plot_var(:Gc_w; title="Conductance (m s-1)", label="Gc"),
plot_var(:fval_soil; title="β_soil", label="β_soil"),
+ plot_var(:VPD; data=df, scale=-1, title="-VPD (kPa)", label="-VPD"),
plot_var(:Pi; title="P - Ei (mm/d)"),
plot_var(:Es_eq; title="Es_eq (mm/d)"),
size=(800, 700),
- layout=(3, 2)
-)
This page was generated using Literate.jl.
Settings
This document was generated with Documenter.jl version 1.7.0 on Tuesday 24 September 2024. Using Julia version 1.10.5.
+ layout=(4, 2)
+)
This page was generated using Literate.jl.