Skip to content

Increase Charnock parameter to see effect on example global simulations #541

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
wants to merge 8 commits into
base: main
Choose a base branch
from
Open

Increase Charnock parameter to see effect on example global simulations #541

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
8 commits
Select commit Hold shift + click to select a range
f7f067f
Update roughness_lengths.jl
glwagner May 23, 2025
517c43d
Update roughness_lengths.jl
glwagner May 23, 2025
8104f35
Update ocean_simulation.jl
glwagner May 24, 2025
6e2f56a
Merge branch 'main' into glw/increase-charnock
navidcy May 30, 2025
36bcf30
Merge branch 'main' into glw/increase-charnock
navidcy May 30, 2025
c66c1db
Merge branch 'main' into glw/increase-charnock
glwagner Jun 5, 2025
8b982b3
Update one_degree_simulation.jl
glwagner Jun 5, 2025
d846b9f
Merge branch 'main' into glw/increase-charnock
simone-silvestri Jul 14, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion examples/one_degree_simulation.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down
4 changes: 2 additions & 2 deletions src/OceanSeaIceModels/InterfaceComputations/roughness_lengths.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down Expand Up @@ -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
Expand Down
2 changes: 1 addition & 1 deletion src/OceanSimulations/ocean_simulation.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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;
Expand Down
Loading