Changeset 8597 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icealb.F90

Timestamp:

2017-10-05T16:44:46+02:00 (7 years ago)

Author:

clem

Message:

first step to make melt ponds compliant with the new code

File:

• branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icealb.F90 (modified) (5 diffs)

branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icealb.F90

@@ -49 +49 @@
 CONTAINS
 
-   SUBROUTINE ice_alb( pt_su, ph_ice, ph_snw, pafrac_pnd, ph_pnd, palb_cs, palb_os )
+   SUBROUTINE ice_alb( pt_su, ph_ice, ph_snw, ld_pnd_alb, pafrac_pnd, ph_pnd, palb_cs, palb_os )
       !!----------------------------------------------------------------------
       !!               ***  ROUTINE ice_alb  ***
@@ -98 +98 @@
       REAL(wp), INTENT(in   ), DIMENSION(:,:,:) ::   ph_ice       !  sea-ice thickness
       REAL(wp), INTENT(in   ), DIMENSION(:,:,:) ::   ph_snw       !  snow depth
+      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)
       REAL(wp), INTENT(in   ), DIMENSION(:,:,:) ::   ph_pnd       !  melt pond depth
@@ -125 +126 @@
                IF( ph_snw(ji,jj,jl) == 0._wp ) THEN
                   zafrac_snw = 0._wp
-                  zafrac_pnd = pafrac_pnd(ji,jj,jl)
+                  IF( ld_pnd_alb ) THEN
+                     zafrac_pnd = pafrac_pnd(ji,jj,jl)
+                  ELSE
+                     zafrac_pnd = 0._wp
+                  ENDIF
                   zafrac_ice = 1._wp - zafrac_pnd
                ELSE
@@ -134 +139 @@
                !
                !                       !--- Bare ice albedo (for hi > 150cm)
-               IF( zafrac_pnd > 0._wp ) THEN
+               IF( ld_pnd_alb ) THEN
                   zalb_ice = rn_alb_idry
                ELSE
@@ -154 +159 @@
                ENDIF
                !                       !--- Ponded ice albedo
-               IF( zafrac_pnd > 0._wp ) THEN
+               IF( ld_pnd_alb ) THEN
                   zalb_pnd = rn_alb_dpnd - ( rn_alb_dpnd - zalb_ice ) * EXP( - ph_pnd(ji,jj,jl) * z1_href_pnd ) 
                ELSE