Archer

[josuemtzmo]

This PR works for the following versions of the packages:

pkg> status
Status `/data/hpcdata/users/josnez/Oceananigans.jl/ice-eddy/RUN/Project.toml`
  [6ba0ff68] ClimaSeaIce v0.5.0
  [85f8d34a] NCDatasets v0.14.15
  [904d977b] NumericalEarth v0.4.2 `/users/josnez/github/NumericalEarth.jl#main`
  [9e8cae18] Oceananigans v0.108.0

Additionally, some changes in the source code of the NumericalEarth package are required to run the model:

Change in the NumericalEarth is to add ocean_reference_density within the parameters of the function.

function sea_ice_dynamics(grid, ocean=nothing;
                          sea_ice_ocean_drag_coefficient = 3.24e-3,
                          rheology = ElastoViscoPlasticRheology(),
                          coriolis = default_coriolis(ocean),
                          free_drift = nothing,
                          ocean_reference_density = ocean_reference_density(ocean),
                          solver = default_solver(grid, ocean))

    SSU, SSV = ocean_surface_velocities(ocean)
    FT = eltype(grid)
    sea_ice_ocean_drag_coefficient = convert(FT, sea_ice_ocean_drag_coefficient)
    ρₑ = convert(FT, ocean_reference_density)

Additionally, changes described in NumericalEarth/NumericalEarth.jl#269 allow the model to work without atmosphere and radiation.