Get constants (e.g., g, R_Earth, Ω_Earth) from Oceananigans.defaults (#658)

* get constants from Oceananigans.defaults

* Update ClimaSeaIce compatibility version to 0.3.9

this captures CliMA/ClimaSeaIce.jl#94

Author: navidcy

Project.toml

@@ -45,7 +45,7 @@ ClimaOceanReactantExt = "Reactant"
 Adapt = "4"
 CFTime = "0.1, 0.2"
 CUDA = "4, 5"
-ClimaSeaIce = "0.3.8"
+ClimaSeaIce = "0.3.9"
 CondaPkg = "0.2.28"
 CubicSplines = "0.2"
 DataDeps = "0.7"
@@ -58,7 +58,7 @@ KernelAbstractions = "0.9"
 MPI = "0.20"
 MeshArrays = "0.3"
 NCDatasets = "0.12, 0.13, 0.14"
-Oceananigans = "0.100"
+Oceananigans = "0.100.3"
 OffsetArrays = "1.14"
 PrecompileTools = "1"
 PythonCall = "0.9"

src/OceanSeaIceModels/InterfaceComputations/coefficient_based_turbulent_fluxes.jl

@@ -1,5 +1,6 @@
 using DocStringExtensions
-using Oceananigans.BuoyancyFormulations: g_Earth
+
+g_Earth = Oceananigans.defaults.gravitational_acceleration
 
 """
     struct CoefficientBasedFluxes{CD, CH, CQ, S}

src/OceanSeaIceModels/InterfaceComputations/component_interfaces.jl

@@ -16,20 +16,19 @@ using ..OceanSeaIceModels.PrescribedAtmospheres:
 
 using ClimaSeaIce: SeaIceModel
 
-using Oceananigans.BuoyancyFormulations: g_Earth
 using Oceananigans: HydrostaticFreeSurfaceModel, architecture
 using Oceananigans.Grids: inactive_node, node, topology
 using Oceananigans.BoundaryConditions: fill_halo_regions!
-
 using Oceananigans.Fields: ConstantField, interpolate, FractionalIndices
 using Oceananigans.Utils: launch!, Time, KernelParameters
-
 using Oceananigans.Operators: ℑxᶜᵃᵃ, ℑyᵃᶜᵃ, ℑxᶠᵃᵃ, ℑyᵃᶠᵃ
 
 using KernelAbstractions: @kernel, @index
 
 import Oceananigans.Simulations: initialize!
 
+g_Earth = Oceananigans.defaults.gravitational_acceleration
+
 #####
 ##### Container for organizing information related to fluxes
 #####
@@ -189,11 +188,11 @@ function atmosphere_ocean_interface(atmos,
     downwelling_longwave  = Field{Center, Center, Nothing}(ocean.model.grid)
     downwelling_shortwave = Field{Center, Center, Nothing}(ocean.model.grid)
 
-    ao_fluxes = (; latent_heat, 
-                   sensible_heat, 
-                   water_vapor, 
+    ao_fluxes = (; latent_heat,
+                   sensible_heat,
+                   water_vapor,
                    x_momentum,
-                   y_momentum, 
+                   y_momentum,
                    friction_velocity,
                    temperature_scale,
                    water_vapor_scale,

src/OceanSeaIceModels/InterfaceComputations/similarity_theory_turbulent_fluxes.jl

@@ -1,6 +1,5 @@
 using Oceananigans.Utils: prettysummary
 using Oceananigans.Grids: AbstractGrid
-using Oceananigans.BuoyancyFormulations: g_Earth
 
 using Adapt
 using Thermodynamics: Liquid
@@ -14,6 +13,8 @@ using Statistics: norm
 import Thermodynamics as AtmosphericThermodynamics
 import Thermodynamics.Parameters: Rv_over_Rd
 
+g_Earth = Oceananigans.defaults.gravitational_acceleration
+
 #####
 ##### Bulk turbulent fluxes based on similarity theory
 #####

src/OceanSimulations/OceanSimulations.jl

@@ -9,6 +9,7 @@ using Oceananigans.Advection: FluxFormAdvection
 using Oceananigans.BoundaryConditions: DefaultBoundaryCondition
 using Oceananigans.ImmersedBoundaries: immersed_peripheral_node, inactive_node, MutableGridOfSomeKind
 using Oceananigans.OrthogonalSphericalShellGrids
+using Oceananigans.Operators
 
 using Oceananigans.TurbulenceClosures.TKEBasedVerticalDiffusivities:
     CATKEVerticalDiffusivity,
@@ -17,9 +18,8 @@ using Oceananigans.TurbulenceClosures.TKEBasedVerticalDiffusivities:
 
 using SeawaterPolynomials.TEOS10: TEOS10EquationOfState
 
-using Oceananigans.BuoyancyFormulations: g_Earth
-using Oceananigans.Coriolis: Ω_Earth
-using Oceananigans.Operators
+g_Earth = Oceananigans.defaults.gravitational_acceleration
+Ω_Earth = Oceananigans.defaults.planet_rotation_rate
 
 struct Default{V}
     value :: V

src/SeaIceSimulations.jl

@@ -9,9 +9,6 @@ using Oceananigans.Advection: FluxFormAdvection
 using Oceananigans.BoundaryConditions: DefaultBoundaryCondition
 using Oceananigans.ImmersedBoundaries: immersed_peripheral_node, inactive_node
 using Oceananigans.OrthogonalSphericalShellGrids
-
-using Oceananigans.BuoyancyFormulations: g_Earth
-using Oceananigans.Coriolis: Ω_Earth
 using Oceananigans.Operators
 
 using ClimaSeaIce
@@ -22,6 +19,9 @@ using ClimaSeaIce.Rheologies: IceStrength, ElastoViscoPlasticRheology
 
 using ClimaOcean.OceanSimulations: Default
 
+g_Earth = Oceananigans.defaults.gravitational_acceleration
+Ω_Earth = Oceananigans.defaults.planet_rotation_rate
+
 function sea_ice_simulation(grid, ocean=nothing;
                             Δt = 5minutes,
                             ice_salinity = 4, # psu