Use AdvectiveFormulation for the skew fluxes (#628)

* Update eddy closure to use new import syntax

Refactor eddy closure initialization to use imported types.

* Update one_degree_simulation.jl

* update packages

* update packages

---------

Co-authored-by: Navid C. Constantinou <[email protected]>

Author: simone-silvestri

Project.toml

@@ -46,7 +46,7 @@ ClimaOceanReactantExt = "Reactant"
 Adapt = "4"
 CFTime = "0.1, 0.2"
 CUDA = "4, 5"
-ClimaSeaIce = "0.3.5"
+ClimaSeaIce = "0.3.7"
 CondaPkg = "0.2.28"
 CubicSplines = "0.2"
 DataDeps = "0.7"

examples/one_degree_simulation.jl

@@ -51,7 +51,9 @@ grid = ImmersedBoundaryGrid(underlying_grid, GridFittedBottom(bottom_height);
 # eddy fluxes. For vertical mixing at the upper-ocean boundary layer we include the CATKE
 # parameterization.
 
-eddy_closure = Oceananigans.TurbulenceClosures.IsopycnalSkewSymmetricDiffusivity(κ_skew=1e3, κ_symmetric=1e3)
+using Oceananigans.TurbulenceClosures: IsopycnalSkewSymmetricDiffusivity, AdvectiveFormulation
+
+eddy_closure = IsopycnalSkewSymmetricDiffusivity(κ_skew=1e3, κ_symmetric=1e3, skew_flux_formulation=AdvectiveFormulation()) 
 vertical_mixing = ClimaOcean.OceanSimulations.default_ocean_closure()
 
 # ### Ocean simulation