pips

160d

140d

6 helems, 30 z, 120s

FriersonDiffusion

albedo

FriersonDiffusion

debug

debug

format

var albedo for slab ocean

Author: LenkaNovak

.buildkite/hierarchies/pipeline.yml

@@ -62,6 +62,16 @@ steps:
         command:
           - "julia --color=yes --project=experiments/ClimaEarth/ experiments/ClimaEarth/run_moist_held_suarez.jl"
         artifact_paths:  "moist_held_suarez/moist_held_suarez/clima_atmos/*"
+        agents:
+          queue: clima
+          slurm_mem: 20GB
+          slurm_gpus: 1
+          modules: common
+
+      - label: "Clima: GPU ClimaCoupler Cloudless Aquaplanet"
+        command:
+          - "julia --color=yes --project=experiments/ClimaEarth/ experiments/ClimaEarth/run_cloudless_aquaplanet.jl"
+        artifact_paths:  "cloudless_aquaplanet/cloudless_aquaplanet/clima_atmos/*"
         agents:
           queue: clima
           slurm_mem: 20GB

.buildkite/pipeline.yml

@@ -396,6 +396,15 @@ steps:
         agents:
           slurm_mem: 20GB
 
+      - label: ":construction: Cloudless Aquaplanet"
+        key: "cloudless_aquaplanet"
+        command:
+          - sed 's/t_end = "1000days"/t_end = "1days"/' experiments/ClimaEarth/run_cloudless_aquaplanet.jl > experiments/ClimaEarth/run_cloudless_aquaplanet_short.jl
+          - "julia --color=yes --project=experiments/ClimaEarth/ experiments/ClimaEarth/run_cloudless_aquaplanet_short.jl"
+        artifact_paths:  "cloudless_aquaplanet/cloudless_aquaplanet/clima_atmos/*"
+        agents:
+          slurm_mem: 20GB
+
   - group: "GPU integration tests"
     steps:
       # GPU RUNS: slabplanet

experiments/ClimaEarth/components/ocean/slab_ocean.jl

@@ -89,6 +89,8 @@ function ocean_init(
         thermo_params = thermo_params,
         # add dss_buffer to cache to avoid runtime dss allocation
         dss_buffer = CC.Spaces.create_dss_buffer(CC.Fields.zeros(space)),
+        α_direct = CC.Fields.ones(space) .* params.α,
+        α_diffuse = CC.Fields.ones(space) .* params.α,
     )
 
     ode_algo = CTS.ExplicitAlgorithm(stepper)
@@ -110,8 +112,8 @@ Interfacer.get_field(sim::SlabOceanSimulation, ::Val{:area_fraction}) = sim.inte
 Interfacer.get_field(sim::SlabOceanSimulation, ::Val{:beta}) = convert(eltype(sim.integrator.u), 1.0)
 Interfacer.get_field(sim::SlabOceanSimulation, ::Val{:roughness_buoyancy}) = sim.integrator.p.params.z0b
 Interfacer.get_field(sim::SlabOceanSimulation, ::Val{:roughness_momentum}) = sim.integrator.p.params.z0m
-Interfacer.get_field(sim::SlabOceanSimulation, ::Union{Val{:surface_direct_albedo}, Val{:surface_diffuse_albedo}}) =
-    sim.integrator.p.params.α
+Interfacer.get_field(sim::SlabOceanSimulation, ::Val{:surface_direct_albedo}) = sim.integrator.p.α_direct
+Interfacer.get_field(sim::SlabOceanSimulation, ::Val{:surface_diffuse_albedo}) = sim.integrator.p.α_diffuse
 Interfacer.get_field(sim::SlabOceanSimulation, ::Val{:surface_humidity}) = sim.integrator.p.q_sfc
 Interfacer.get_field(sim::SlabOceanSimulation, ::Val{:surface_temperature}) = sim.integrator.u.T_sfc
 Interfacer.get_field(sim::SlabOceanSimulation, ::Val{:water}) = nothing
@@ -137,6 +139,12 @@ end
 function Interfacer.update_field!(sim::SlabOceanSimulation, ::Val{:turbulent_energy_flux}, field)
     parent(sim.integrator.p.F_turb_energy) .= parent(field)
 end
+function Interfacer.update_field!(sim::SlabOceanSimulation, ::Val{:surface_direct_albedo}, field::CC.Fields.Field)
+    sim.integrator.p.α_direct .= field
+end
+function Interfacer.update_field!(sim::SlabOceanSimulation, ::Val{:surface_diffuse_albedo}, field::CC.Fields.Field)
+    sim.integrator.p.α_diffuse .= field
+end
 
 # extensions required by FieldExchanger
 Interfacer.step!(sim::SlabOceanSimulation, t) = Interfacer.step!(sim.integrator, t - sim.integrator.t, true)