diff --git a/examples/one_degree_simulation.jl b/examples/one_degree_simulation.jl index 65e3acf6f..4a886b3d2 100644 --- a/examples/one_degree_simulation.jl +++ b/examples/one_degree_simulation.jl @@ -174,7 +174,7 @@ seaice.output_writers[:surface] = JLD2Writer(ocean.model, seaice_outputs; run!(simulation) simulation.Δt = 20minutes -simulation.stop_time = 365days +simulation.stop_time = 5 * 365days run!(simulation) # ### A pretty movie diff --git a/src/OceanSeaIceModels/InterfaceComputations/roughness_lengths.jl b/src/OceanSeaIceModels/InterfaceComputations/roughness_lengths.jl index 62845744a..4bb72aba0 100644 --- a/src/OceanSeaIceModels/InterfaceComputations/roughness_lengths.jl +++ b/src/OceanSeaIceModels/InterfaceComputations/roughness_lengths.jl @@ -88,7 +88,7 @@ function MomentumRoughnessLength(FT=Oceananigans.defaults.FloatType; gravitational_acceleration = default_gravitational_acceleration, maximum_roughness_length = 1, air_kinematic_viscosity = 1.5e-5, - wave_formulation = 0.02, + wave_formulation = 0.04, smooth_wall_parameter = 0.11) if wave_formulation isa Number @@ -207,7 +207,7 @@ ReynoldsScalingFunction(FT = Oceananigans.defaults.FloatType; A = 5.85e-5, b = 0 # Clip ℓ_max = ℓ.maximum_roughness_length - return min(ℓs, ℓ_max) + ℓu #return min(ℓs, ℓ_max) end # Convenience for users diff --git a/src/OceanSimulations/ocean_simulation.jl b/src/OceanSimulations/ocean_simulation.jl index 3786a689a..ecf632ee7 100644 --- a/src/OceanSimulations/ocean_simulation.jl +++ b/src/OceanSimulations/ocean_simulation.jl @@ -86,7 +86,7 @@ function default_ocean_closure(FT=Oceananigans.defaults.FloatType) end function default_radiative_forcing(grid) - ϵʳ = 0.6 # red fraction + ϵʳ = 0.0 # red fraction λʳ = 1 # red decay scale λᵇ = 16 # blue decay scale forcing = TwoColorRadiation(grid;