Skip to content

Commit

Permalink
Update arctic_night.jl for latest Oceananigans (#325)
Browse files Browse the repository at this point in the history
* Update arctic_night.jl

* Update idealized_single_column_simulation.jl

* Update idealized_single_column_simulation.jl

---------

Co-authored-by: Navid C. Constantinou <[email protected]>
  • Loading branch information
glwagner and navidcy authored Jan 27, 2025
1 parent c5108fc commit aa475d8
Show file tree
Hide file tree
Showing 2 changed files with 4 additions and 5 deletions.
4 changes: 2 additions & 2 deletions examples/arctic_night.jl
Original file line number Diff line number Diff line change
Expand Up @@ -29,8 +29,8 @@ parent(atmosphere.downwelling_radiation.shortwave) .= Qℓ # W
grid = RectilinearGrid(size=400, z=(-400, 0), topology=(Flat, Flat, Bounded))
ocean = ocean_simulation(grid, coriolis=FPlane(; f))

eos = ocean.model.buoyancy.model.equation_of_state
g = ocean.model.buoyancy.model.gravitational_acceleration
eos = ocean.model.buoyancy.formulation.equation_of_state
g = ocean.model.buoyancy.formulation.gravitational_acceleration
α = SeawaterPolynomials.thermal_expansion(T₀, S₀, 0, eos)
dTdz =/* g)
Tᵢ(z) = max(-1.8, T₀ + dTdz * z)
Expand Down
5 changes: 2 additions & 3 deletions examples/idealized_single_column_simulation.jl
Original file line number Diff line number Diff line change
Expand Up @@ -29,8 +29,8 @@ parent(atmosphere.downwelling_radiation.shortwave) .= Qℓ # W
grid = RectilinearGrid(size=100, z=(-400, 0), topology=(Flat, Flat, Bounded))
ocean = ocean_simulation(grid, coriolis=FPlane(; f))

eos = ocean.model.buoyancy.model.equation_of_state
g = ocean.model.buoyancy.model.gravitational_acceleration
eos = ocean.model.buoyancy.formulation.equation_of_state
g = ocean.model.buoyancy.formulation.gravitational_acceleration
α = SeawaterPolynomials.thermal_expansion(T₀, S₀, 0, eos)
dTdz =/* g)
Tᵢ(z) = T₀ + dTdz * z
Expand Down Expand Up @@ -115,4 +115,3 @@ record(fig, "idealized_atmosphere.mp4", 1:Nt, framerate=24) do nn
@info "Drawing frame $nn of $Nt..."
n[] = nn
end

0 comments on commit aa475d8

Please sign in to comment.