correction (#492)

Author: simone-silvestri

src/OceanSeaIceModels/InterfaceComputations/assemble_net_fluxes.jl

@@ -183,7 +183,8 @@ function compute_net_sea_ice_fluxes!(coupled_model)
     kernel_parameters = interface_kernel_parameters(grid)
 
     sea_ice_surface_temperature = coupled_model.interfaces.atmosphere_sea_ice_interface.temperature
-
+    ice_concentration = sea_ice_concentration(sea_ice)
+    
     launch!(arch, grid, kernel_parameters,
             _assemble_net_sea_ice_fluxes!,
             top_fluxes,
@@ -193,6 +194,7 @@ function compute_net_sea_ice_fluxes!(coupled_model)
             atmosphere_sea_ice_fluxes,
             sea_ice_ocean_fluxes,
             freshwater_flux,
+            ice_concentration,
             sea_ice_surface_temperature,
             downwelling_radiation,
             sea_ice_properties,
@@ -208,6 +210,7 @@ end
                                                atmosphere_sea_ice_fluxes,
                                                sea_ice_ocean_fluxes,
                                                freshwater_flux, # Where do we add this one?
+                                               ice_concentration,
                                                surface_temperature,
                                                downwelling_radiation,
                                                sea_ice_properties,
@@ -220,7 +223,8 @@ end
     @inbounds begin
         Ts = surface_temperature[i, j, kᴺ]
         Ts = convert_to_kelvin(sea_ice_properties.temperature_units, Ts)
-
+        ℵi = ice_concentration[i, j, 1]
+        
         Qs = downwelling_radiation.Qs[i, j, 1]
         Qℓ = downwelling_radiation.Qℓ[i, j, 1]
         Qc = atmosphere_sea_ice_fluxes.sensible_heat[i, j, 1] # sensible or "conductive" heat flux
@@ -239,7 +243,7 @@ end
     Qu = upwelling_radiation(Ts, σ, ϵ)
     Qd = net_downwelling_radiation(i, j, grid, time, α, ϵ, Qs, Qℓ)
 
-    ΣQt = Qd + Qu + Qc + Qv
+    ΣQt = (Qd + Qu + Qc + Qv) * ℵi # If ℵi == 0 there is no heat flux from the top!
     ΣQb = Qf + Qi
 
     # Mask fluxes over land for convenience