Changeset 14616

Timestamp:

2021-03-18T17:35:09+01:00 (3 years ago)

Author:

dancopsey

Message:

Add qevap_ice as the enthalpy in sublimating ice/snow.

File:

• NEMO/branches/UKMO/NEMO_4.0.4_penetrating_solar/src/OCE/SBC/sbccpl.F90 (modified) (2 diffs)

NEMO/branches/UKMO/NEMO_4.0.4_penetrating_solar/src/OCE/SBC/sbccpl.F90

@@ -1883 +1883 @@
       CASE( 'oce only' )         ! the required field is directly provided
 
-         ! Get the sea ice non solar heat flux from conductive, melting and sublimation fluxes
+         ! Get the sea ice non solar heat flux from conductive and melting fluxes
          IF( TRIM(sn_rcv_iceflx%cldes) == 'coupled' ) THEN
             zqns_ice(:,:,:) = qml_ice(:,:,:) + qcn_ice(:,:,:)
@@ -1958 +1958 @@
       ! --- heat content of evap over ice in W/m2 (to be used in 1D-thermo) --- !
       DO jl = 1, jpl
-         zqevap_ice(:,:,jl) = 0._wp ! should be -evap * ( ( Tice - rt0 ) * rcpi ) but atm. does not take it into account
+         ! zqevap_ice(:,:,jl) = 0._wp ! should be -evap * ( ( Tice - rt0 ) * rcpi ) but atm. does not take it into account
+
+         ! How much enthalpy is stored in sublimating snow and ice. At this stage we don't know if it is snow or ice that is
+         ! sublimating so we will use the combined snow and ice layer temperature t1_ice.
+         zqevap_ice(:,:,jl) = -zevap_ice(:,:,jl) * ( ( rt0 - t1_ice(:,:,jl) ) * rcpi + rLfus )
+
+         ! This enthalpy change should always be negative as only positive sublimation is used (energy loss from the ice).
+         WHERE( zqevap_ice(:,:,jl) > 0._wp ) zqevap_ice(:,:,jl) = 0._wp
       END DO