Its ok to iterate on land a fixed number of times (#376)

Author: glwagner

src/OceanSeaIceModels/InterfaceComputations/atmosphere_ocean_fluxes.jl

@@ -116,7 +116,12 @@ end
     qₛ = saturation_specific_humidity(q_formulation, ℂₐ, 𝒬ₐ.ρ, Tᵢ, Sᵢ) 
     initial_interface_state = InterfaceState(u★, u★, u★, uᵢ, vᵢ, Tᵢ, Sᵢ, qₛ)
 
-    if inactive_node(i, j, kᴺ, grid, Center(), Center(), Center())
+    # Don't use convergence criteria in an inactive cell
+    stop_criteria = turbulent_flux_formulation.solver_stop_criteria
+    needs_to_converge = stop_criteria isa ConvergenceStopCriteria
+    not_water = inactive_node(i, j, kᴺ, grid, Center(), Center(), Center())
+
+    if needs_to_converge && not_water
         interface_state = zero_interface_state(FT)
     else
         interface_state = compute_interface_state(turbulent_flux_formulation,
@@ -129,6 +134,9 @@ end
                                                   ocean_properties)
     end
 
+    # In the case of FixedIterations, make sure interface state is zero'd
+    interface_state = ifelse(not_water, zero_interface_state(FT), interface_state)
+
     u★ = interface_state.u★
     θ★ = interface_state.θ★
     q★ = interface_state.q★

src/OceanSeaIceModels/InterfaceComputations/atmosphere_sea_ice_fluxes.jl

@@ -129,10 +129,13 @@ end
     # Guess
     Sₛ = zero(FT) # what should we use for interface salinity?
     initial_interface_state = InterfaceState(u★, u★, u★, uᵢ, vᵢ, Tₛ, Sₛ, convert(FT, qₛ))
-    land = inactive_node(i, j, kᴺ, grid, Center(), Center(), Center())
+    not_water = inactive_node(i, j, kᴺ, grid, Center(), Center(), Center())
     ice_free = ℵᵢ == 0
 
-    if (land | ice_free)
+    stop_criteria = turbulent_flux_formulation.solver_stop_criteria
+    needs_to_converge = stop_criteria isa ConvergenceStopCriteria
+
+    if needs_to_converge && not_water || ice_free
         interface_state = InterfaceState(zero(FT), zero(FT), zero(FT), uᵢ, vᵢ, Tᵢ, Sₛ, zero(FT))
     else
         interface_state = compute_interface_state(turbulent_flux_formulation,