Merge pull request #1004 from CliMA/zs/tinyamip_test

change ocean roughness and ice albedo
CliMA · Oct 9, 2024 · e9aedc9 · e9aedc9
2 parents d72c420 + 746ad8f
commit e9aedc9
Show file tree

Hide file tree

Showing 3 changed files with 6 additions and 3 deletions.
diff --git a/config/nightly_configs/amip_coarse.yml b/config/nightly_configs/amip_coarse.yml
@@ -15,6 +15,7 @@ hourly_checkpoint_dt: 720
 land_albedo_type: "map_temporal"
 mode_name: "amip"
 mono_surface: false
+netcdf_interpolation_num_points: [90, 45, 31]
 netcdf_output_at_levels: true
 output_default_diagnostics: true
 rayleigh_sponge: true

diff --git a/experiments/ClimaEarth/components/ocean/prescr_seaice.jl b/experiments/ClimaEarth/components/ocean/prescr_seaice.jl
@@ -46,7 +46,7 @@ Base.@kwdef struct IceSlabParameters{FT <: AbstractFloat}
     z0b::FT = 1e-4          # roughness length for tracers [m]
     T_freeze::FT = 271.2    # freezing temperature of sea water [K]
     k_ice::FT = 2           # thermal conductivity of sea ice [W / m / K] (less in HM71)
-    α::FT = 0.8             # albedo of sea ice [0, 1]
+    α::FT = 0.65            # albedo of sea ice, roughly tuned to match observations
 end
 
 Interfacer.name(::IceSlabParameters) = "IceSlabParameters"

diff --git a/experiments/ClimaEarth/run_amip.jl b/experiments/ClimaEarth/run_amip.jl
@@ -297,8 +297,10 @@ if mode_name == "amip"
     ocean_sim = Interfacer.SurfaceStub((;
         T_sfc = SST_init,
         ρ_sfc = CC.Fields.zeros(boundary_space),
-        z0m = FT(1e-3),
-        z0b = FT(1e-3),
+        # ocean roughness follows GFDL model 
+        # (https://github.com/NOAA-GFDL/ice_param/blob/main/ocean_rough.F90#L47)
+        z0m = FT(5.8e-5),
+        z0b = FT(5.8e-5),
         beta = FT(1),
         α_direct = CC.Fields.ones(boundary_space) .* FT(0.06),
         α_diffuse = CC.Fields.ones(boundary_space) .* FT(0.06),