wip

CliMA · Oct 23, 2023 · 74281df · 74281df
1 parent 506a1e3
commit 74281df
Show file tree

Hide file tree

Showing 3 changed files with 11 additions and 1 deletion.
diff --git a/examples/Project.toml b/examples/Project.toml
@@ -5,6 +5,7 @@ ClimaAtmos = "b2c96348-7fb7-4fe0-8da9-78d88439e717"
 Distributions = "31c24e10-a181-5473-b8eb-7969acd0382f"
 EnsembleKalmanProcesses = "aa8a2aa5-91d8-4396-bcef-d4f2ec43552d"
 JLD2 = "033835bb-8acc-5ee8-8aae-3f567f8a3819"
+NetCDF = "30363a11-5582-574a-97bb-aa9a979735b9"
 TOML = "fa267f1f-6049-4f14-aa54-33bafae1ed76"
 
 [compat]

diff --git a/examples/driver.sbatch b/examples/driver.sbatch
@@ -8,7 +8,7 @@ export JULIA_LOAD_PATH=":/examples/load_path"
 
 job_id=${1?Error: no job ID given}
 
-julia --project=examples -e "using Pkg; Pkg.instantiate(;verbose=true)"
+julia --project=examples -e 'using Pkg; Pkg.develop('.") Pkg.instantiate(;verbose=true)'
 
 # This is a bit messy - could be file
 julia --project=examples -e "import CalibrateAtmos; using Distributions; using EnsembleKalmanProcesses.ParameterDistributions \

diff --git a/examples/hs_obs_map.jl b/examples/hs_obs_map.jl
@@ -0,0 +1,9 @@
+using NetCDF
+
+ta = ncread("~/Downloads/default/ta_60.0d_average.nc", "ta");
+ta_third_height = ta[3:size(ta)[1],:,:,3]
+dims = (2, 3)
+
+avg_ta_third_height = dropdims(mean(ta_third_height; dims = dims), dims = dims)
+extrema_ta_third_height = [extrema(ta_third_height[i,:,:]) for i in 1:size(ta_third_height)[1]]
+