Changeset 8142 for branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/OPA_SRC/SBC/albedo.F90

Timestamp:

2017-06-06T09:39:43+02:00 (7 years ago)

Author:

vancop

Message:

Melt pond interfaces practically operational

File:

• branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/OPA_SRC/SBC/albedo.F90 (modified) (2 diffs)

branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/OPA_SRC/SBC/albedo.F90

@@ -110 +110 @@
       !!---------------------------------------------------------------------
 
-      WRITE(numout,*) ' albedo_ice '
-      WRITE(numout,*) ' ld_pnd     : ', ld_pnd
-      WRITE(numout,*) ' nn_ice_alb : ', nn_ice_alb
-
       ijpl = SIZE( pt_ice, 3 )                     ! number of ice categories
       
@@ -326 +322 @@
 
           ! Depth-dependent ponded ice albedo
-          zhref_pnd = 0.05        ! Characteristic length scale for thickness dependence of ponded ice albedo, Lecomte et al (2015)
-          zalb_pnd  = ralb_dp - ( ralb_dp - zalb_ice ) * EXP( - ph_pnd / zhref_pnd ) 
+          IF ( ld_pnd ) THEN
+             zhref_pnd = 0.05        ! Characteristic length scale for thickness dependence of ponded ice albedo, Lecomte et al (2015)
+             zalb_pnd  = ralb_dp - ( ralb_dp - zalb_ice ) * EXP( - ph_pnd / zhref_pnd ) 
+          ELSE
+             zalb_pnd  = ralb_dp
+          ENDIF
 
           ! Surface albedo is weighted mean of snow, ponds and bare ice contributions