model_calib passed

jl-pkgs · Sep 24, 2024 · 292186f · 292186f
1 parent ddeb4a8
commit 292186f
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 4 deletions.
diff --git a/scripts/calib_IGBP.jl b/scripts/calib_IGBP.jl
@@ -5,10 +5,6 @@ df = data[data.IGBP.==IGBP, :]
 printstyled("[IGBP=$IGBP] \n", bold=true, color=:green, underline=true)
 
 # @time theta, goal, flag = m_calib(df; IGBPcode, maxn=2500);
-@time _theta, goal, flag = model_calib(df, par0; IGBPcode=df.IGBPcode[1], maxn=2500);
-
-model_goal(df, _theta; verbose=true)
-
 r = PMLV2_sites(df; par=theta2par(_theta))
 r = cbind(df[:, [:site, :date, :ET_obs, :GPP_obs]], r)
 

diff --git a/test/test-PMLV2.jl b/test/test-PMLV2.jl
@@ -47,3 +47,12 @@ end
   @test GOF(df_out.ET_sim, r.ET).MAE <= 0.015
   @test GOF(df_out.GPP_sim, r.GPP).MAE <= 1E-8
 end
+
+
+@testset "model_calib" begin
+  df.GPP_obs = df.GPPobs
+  df.ET_obs = df.ETobs
+  @time _theta, goal, flag = model_calib(df, par0; IGBPcode=df.IGBPcode[1], maxn=2500)
+  goal = model_goal(df, _theta; verbose=true)
+  @test goal > 0.55 # mean(KGE_GPP, KGE_ET)
+end