add exp reference

CliMA · Sep 19, 2024 · 7dbc323 · 7dbc323
1 parent 7753eb9
commit 7dbc323
Show file tree

Hide file tree

Showing 2 changed files with 17 additions and 4 deletions.
diff --git a/experiments/standalone/Vegetation/refs/exp_T_dt0.5_24.0hours_updateat3hours.txt b/experiments/standalone/Vegetation/refs/exp_T_dt0.5_24.0hours_updateat3hours.txt
@@ -0,0 +1,9 @@
+291.78999999999996
+291.5839907697842
+291.1769844057755
+293.47485047021405
+296.8251086150474
+297.1236273619955
+295.58039061588073
+291.027492059103
+291.6326663549803
diff --git a/experiments/standalone/Vegetation/timestep_test.jl b/experiments/standalone/Vegetation/timestep_test.jl
@@ -181,10 +181,14 @@ ref_dt = 0.5
 #     SciMLBase.solve(prob, ode_algo; dt = ref_dt, callback = cb, saveat = saveat);
 # ref_T = [parent(ref_sol.u[k].canopy.energy.T)[1] for k in 1:length(ref_sol.t)]
 
-# Read in solution from saved delimited file (experiment run explicitly with dt = 0.1s)
-# savedir = joinpath(pkgdir(ClimaLand), "experiments/standalone/Vegetation");
-# ref_file = joinpath(savedir, "refs", "exp_T_dt$(ref_dt)_$(sim_time)hours_updateat3hours.txt")
-# ref_T = vec(readdlm(ref_file, ','))
+# Read in solution from saved delimited file (experiment run explicitly)
+savedir = joinpath(pkgdir(ClimaLand), "experiments/standalone/Vegetation");
+ref_file = joinpath(
+    savedir,
+    "refs",
+    "exp_T_dt$(ref_dt)_$(sim_time)hours_updateat3hours.txt",
+)
+ref_T = vec(readdlm(ref_file, ','))
 
 mean_err = []
 p95_err = []