Changeset 14392

Timestamp:

2021-02-04T10:08:21+01:00 (3 years ago)

Author:

dancopsey

Message:

Include heat involved in sublimation in the qns_ice flux. Weight qsr_tot by the leads fraction (ziceld).

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 and melting
+         ! Get the sea ice non solar heat flux from conductive, melting and sublimation fluxes
          IF( TRIM(sn_rcv_iceflx%cldes) == 'coupled' ) THEN
-            zqns_ice(:,:,:) = qml_ice(:,:,:) + qcn_ice(:,:,:)
+            zqns_ice(:,:,:) = qml_ice(:,:,:) + qcn_ice(:,:,:) + evap_ice(:,:,:) * rLsub
          ELSE
             zqns_ice(:,:,:) = 0._wp
@@ -2032 +2032 @@
       !                                                      ! ========================= !
       CASE( 'oce only' )
-         zqsr_tot(:,:  ) = MAX( 0._wp , frcv(jpr_qsroce)%z3(:,:,1) )
+         zqsr_tot(:,:  ) = MAX( 0._wp , frcv(jpr_qsroce)%z3(:,:,1) ) * ziceld(:,:)
          ! For Met Office sea ice solar fluxes are already delt with by JULES so setting to zero
          ! here so the only flux is the ocean only one.