Changeset 12725 for NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icealb.F90

Timestamp:

2020-04-09T15:48:28+02:00 (4 years ago)

Author:

clem

Message:

fix previous revision. Make sure a_ip and a_ip_eff are correct

File:

• NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icealb.F90 (modified) (3 diffs)

NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icealb.F90

@@ -96 +96 @@
       LOGICAL , INTENT(in   )                   ::   ld_pnd_alb   !  effect of melt ponds on albedo
       REAL(wp), INTENT(in   ), DIMENSION(:,:,:) ::   pafrac_pnd   !  melt pond relative fraction (per unit ice area)
-                                                                  !  This is the effective fraction not covered up by a pond lid
       REAL(wp), INTENT(in   ), DIMENSION(:,:,:) ::   ph_pnd       !  melt pond depth
       REAL(wp), INTENT(  out), DIMENSION(:,:,:) ::   palb_cs      !  albedo of ice under clear    sky
@@ -120 +119 @@
          DO jj = 1, jpj
             DO ji = 1, jpi
-               !                       !--- Specific snow, ice and pond fractions (for now, we prevent melt ponds and snow at the same time)
-               IF( ph_snw(ji,jj,jl) == 0._wp ) THEN
+               !---------------------------------------------!
+               !--- Specific snow, ice and pond fractions ---!
+               !---------------------------------------------!               
+               IF( ph_snw(ji,jj,jl) == 0._wp ) THEN   !--- no snow : we prevent melt ponds and snow at the same time (for now)
                   zafrac_snw = 0._wp
                   IF( ld_pnd_alb ) THEN
@@ -130 +131 @@
                   zafrac_ice = 1._wp - zafrac_pnd
                ELSE
-                  zafrac_snw = 1._wp      ! Snow fully "shades" melt ponds and ice
+                  zafrac_snw = 1._wp                  !--- snow : fully "shades" melt ponds and ice
                   zafrac_pnd = 0._wp
                   zafrac_ice = 0._wp
                ENDIF
                !
+               !---------------!
+               !--- Albedos ---!
+               !---------------!               
                !                       !--- Bare ice albedo (for hi > 150cm)
                IF( ld_pnd_alb ) THEN