Update init topography profile + config file

	modified:   src/solver/type_getters.jl
	modified:   src/topography/topography.jl
+ Hughes2023 testcase + ci-plots

.buildkite/pipeline.yml

@@ -325,6 +325,13 @@ steps:
           julia --color=yes --project=examples examples/hybrid/driver.jl
           --config_file $CONFIG_PATH/sphere_baroclinic_wave_rhoe_topography_dcmip_rs.yml
         artifact_paths: "sphere_baroclinic_wave_rhoe_topography_dcmip_rs/output_active/*"
+
+      - label: ":computer: baroclinic wave (ρe) topography (hughes2023)"
+        command: >
+          julia --color=yes --project=examples examples/hybrid/driver.jl
+          --config_file $CONFIG_PATH/sphere_baroclinic_wave_rhoe_hughes2023.yml
+        artifact_paths: "sphere_baroclinic_wave_rhoe_hughes2023/output_active/*"
+
 
       - label: ":computer: held suarez (ρe) topography (dcmip)"
         command: >

config/model_configs/sphere_baroclinic_wave_rhoe_hughes2023.yml

@@ -0,0 +1,13 @@
+dt_save_state_to_disk: "2days"
+initial_condition: "DryBaroclinicWave"
+topography: "Hughes2023"
+topo_smoothing: false
+dt: "400secs"
+t_end: "10days"
+job_id: "sphere_baroclinic_wave_rhoe_hughes2023"
+diagnostics:
+  - short_name: [pfull, wa, va, rv, ke]
+    period: 1days
+  - short_name: [pfull, wa, va, rv, ke]
+    period: 1days
+    writer: h5

post_processing/ci_plots.jl

@@ -635,6 +635,7 @@ end
 
 function make_plots(
     ::Val{:sphere_baroclinic_wave_rhoe_topography_dcmip_rs},
+    ::Val{:sphere_baroclinic_wave_rhoe_hughes2023},
     output_paths::Vector{<:AbstractString},
 )
     simdirs = SimDir.(output_paths)

src/solver/type_getters.jl

@@ -141,9 +141,12 @@ function get_spaces(parsed_args, params, comms_ctx)
     bubble = parsed_args["bubble"]
     deep = parsed_args["deep_atmosphere"]
 
-    @assert topography in ("NoWarp", "DCMIP200", "Earth", "Agnesi", "Schar")
+    @assert topography in
+            ("NoWarp", "DCMIP200", "Earth", "Agnesi", "Schar", "Hughes2023")
     if topography == "DCMIP200"
         warp_function = topography_dcmip200
+    elseif topography == "Hughes2023"
+        warp_function = topography_hughes2023
     elseif topography == "Agnesi"
         warp_function = topography_agnesi
     elseif topography == "Schar"

src/topography/topography.jl

@@ -23,6 +23,39 @@ function topography_dcmip200(coords)
     return zₛ
 end
 
+"""
+    topography_hughes2023(λ,ϕ)
+λ = longitude (degrees)
+ϕ = latitude (degrees)
+Returns the surface elevation profile used in the baroclinic wave
+test problem defined by Hughes and Jablonowski (2023).
+"""
+function topography_hughes2023(coords)
+    λ, ϕ = coords.long, coords.lat
+    FT = eltype(λ)
+    h₀ = FT(2e3)
+    # Angles in degrees
+    ϕ₁ = FT(45)
+    ϕ₂ = FT(45)
+    λ_min = minimum(λ)
+    λ₁ = FT(72)
+    λ₂ = FT(140)
+    λₘ = FT(7)
+    ϕₘ = FT(40)
+    d = ϕₘ / 2 * (-log(0.1))^(-1 / 6)
+    c = λₘ / 2 * (-log(0.1))^(-1 / 2)
+    d₁ = @. (λ - λ_min) - λ₁
+    d₂ = @. (λ - λ_min) - λ₂
+    l₁ = @. λ - λ₁
+    l₂ = @. λ - λ₂
+    zₛ = @. FT(
+        h₀ * (
+            exp(-(((ϕ - ϕ₁) / d)^6 + (l₁ / c)^2)) +
+            exp(-(((ϕ - ϕ₂) / d)^6 + (l₂ / c)^2))
+        ),
+    )
+end
+
 function generate_topography_warp(earth_spline)
     function topography_earth(coords)
         λ, Φ = coords.long, coords.lat