working version prior to rebase [skip ci]

Author: Katherine Deck

Artifacts.toml

@@ -20,10 +20,10 @@ git-tree-sha1 = "cbbc6b3752d9cb9b667ec33cfbeb46819f8db418"
     url = "https://caltech.box.com/shared/static/b3u4dv0dsoswvqgp8y63bzj7awbhwztd.gz"
 
 [soil_ic_2008_50m]
-git-tree-sha1 = "fd9cda235e203b4235136b6dcdcd07af788c00f6"
+git-tree-sha1 = "7b0f0572f4620a8ba2a7ae816a547d2dcfb9eb9b"
 
     [[soil_ic_2008_50m.download]]
-    sha256 = "22e3cee446ee3b34edd6956b305bace5f04d12b8c1f124c1891a93048bc91b34"
+    sha256 = "b0080d605585072b88871bc8987f615c9dfc35e64ad16e00a52087d879fc5c3a"
     url = "https://caltech.box.com/shared/static/atiztw8uu2i78g2vd02n8onb68cer9tn.gz"
 
 [soilgrids]

experiments/calibration/calibrate_land.jl

@@ -4,7 +4,7 @@ model_type = "land" # "land" or "bucket"
 const variable_list = ["swu"] # variable(s) you want to capture by adjusting your priors
 const n_iterations = 10 # 1 iterations takes ~ 1.5 hour with current settings ((50, 15) resolution, 2 year simulation)
 const spinup_period = Year(1)
-const start_date = DateTime(2008, 12, 01) # this is the start of the forward model spinup
+const start_date = DateTime(2000, 12, 01) # this is the start of the forward model spinup
 @assert month(start_date + spinup_period) == 12 "The start of your calibration period should be December."
 const nelements = (101, 15) # resolution - (horizontal elements (lon, lat), vertical elements (soil discretization))
 const dirname = "land_snow_zenith_hires" # ideally, describe your calibration in a few words

experiments/calibration/priors_land.jl

@@ -1,6 +1,6 @@
-prior_α_0 = EKP.constrained_gaussian("α_0", 0.65, 0.2, 0.2, 0.95);
+prior_α_0 = EKP.constrained_gaussian("α_0", 0.6, 0.2, 0.2, 0.95);
 prior_Δα = EKP.constrained_gaussian("Δα", 0.2, 0.1, 0.0, 1.0);
-prior_k = EKP.constrained_gaussian("k", 10, 5, 0.1, 25);
+prior_k = EKP.constrained_gaussian("k", 4, 2, 0.01, 25);
 #prior_α_0 = EKP.constrained_gaussian("α_0", 0.7, 0.1, 0.5, 0.95);
 #prior_beta_snow_cover = EKP.constrained_gaussian("beta_snow_cover", 1.77, 0.2, 0.1, 2.5);
 #prior_z0_snow_cover = EKP.constrained_gaussian("z0_snow_cover", 0.106, 0.03, 0.0, 0.3);