add a simulation

Author: simone-silvestri

src/OceanSeaIceModels/InterfaceComputations/assemble_net_fluxes.jl

@@ -21,7 +21,7 @@ using ClimaOcean.OceanSeaIceModels: sea_ice_concentration
     return zero(Iˢʷ)
 end
 
-@inline get_radiative_forcing(ocean::Simulation{<:HydrostaticFreeSurfaceModel}) = get_radiative_forcing(ocean.model.forcing.T)
+@inline get_radiative_forcing(ocean::OceananigansSimulation) = get_radiative_forcing(ocean.model.forcing.T)
 @inline get_radiative_forcing(FT) = FT
 
 @inline function get_radiative_forcing(FT::MultipleForcings)

src/OceanSeaIceModels/InterfaceComputations/atmosphere_ocean_fluxes.jl

@@ -4,7 +4,7 @@ using ClimaOcean.OceanSeaIceModels.PrescribedAtmospheres: thermodynamics_paramet
                                                           surface_layer_height,
                                                           boundary_layer_height
 
-@inline get_ocean_state(ocean::Simulation{<:HydrostaticFreeSurfaceModel}, coupled_model) = 
+@inline get_ocean_state(ocean::OceananigansSimulation, coupled_model) = 
                   (u = ocean.model.velocities.u,
                    v = ocean.model.velocities.v,
                    T = ocean.model.tracers.T,

src/OceanSeaIceModels/InterfaceComputations/component_interfaces.jl

@@ -167,7 +167,7 @@ Base.summary(crf::ComponentInterfaces) = "ComponentInterfaces"
 Base.show(io::IO, crf::ComponentInterfaces) = print(io, summary(crf))
 
 atmosphere_ocean_interface(::Nothing, args...) = nothing
-atmosphere_ocean_interface(ocean::Simulation{<:HydrostaticFreeSurfaceModel}, args...) = 
+atmosphere_ocean_interface(ocean::OceananigansSimulation, args...) = 
     atmosphere_ocean_interface(ocean.model.grid, args...)
 
 function atmosphere_ocean_interface(grid::AbstractGrid, 
@@ -296,7 +296,7 @@ function default_ao_specific_humidity(ocean)
     return ImpureSaturationSpecificHumidity(phase, x_H₂O)
 end
 
-function ocean_surface_fluxes(ocean::Simulation{<:HydrostaticFreeSurfaceModel})
+function ocean_surface_fluxes(ocean::OceananigansSimulation)
     τx = surface_flux(ocean.model.velocities.u)
     τy = surface_flux(ocean.model.velocities.v)
     tracers = ocean.model.tracers

src/OceanSeaIceModels/OceanSeaIceModels.jl

@@ -41,6 +41,7 @@ const default_gravitational_acceleration = 9.80665
 const default_freshwater_density = 1000
 
 const SeaIceSimulation = Simulation{<:SeaIceModel}
+const OceananigansSimulation = Simulation{<:HydrostaticFreeSurfaceModel}
 
 sea_ice_thickness(::Nothing) = ZeroField()
 sea_ice_thickness(sea_ice::SeaIceSimulation) = sea_ice.model.ice_thickness