Changeset 13775 for NEMO/branches/UKMO

Timestamp:

2020-11-11T13:54:38+01:00 (3 years ago)

Author:

dancopsey

Message:

Correct calculation of a_ip_eff so a_ip_eff=0 when h_il=zhl_max and a_ip_eff=a_ip_frac when h_il=zhl_min.

File:

• NEMO/branches/UKMO/NEMO_4.0.3_fix_a_ip_eff/src/ICE/icevar.F90 (modified) (1 diff)

NEMO/branches/UKMO/NEMO_4.0.3_fix_a_ip_eff/src/ICE/icevar.F90

@@ -244 +244 @@
       ELSEWHERE( h_il(:,:,:) >= zhl_max )  ;   a_ip_eff(:,:,:) = 0._wp                  ! lid is very thick. Cover all the pond up with ice and snow
       ELSEWHERE                            ;   a_ip_eff(:,:,:) = a_ip_frac(:,:,:) * &   ! lid is in between. Expose part of the pond
-         &                                                       ( h_il(:,:,:) - zhl_min ) / ( zhl_max - zhl_min )
+         &                                                       ( zhl_max - h_il(:,:,:) ) / ( zhl_max - zhl_min )
       END WHERE
       !