Changeset 15574 for NEMO/branches/2021/dev_r14318_RK3_stage1/src/ICE/icealb.F90

Timestamp:

2021-12-03T20:32:50+01:00 (3 years ago)

Author:

techene

Message:

#2605 #2715 trunk merged into dev_r14318_RK3_stage1

Location:

NEMO/branches/2021/dev_r14318_RK3_stage1

Files:

• . (modified) (1 prop)
• src/ICE/icealb.F90 (modified) (9 diffs)

NEMO/branches/2021/dev_r14318_RK3_stage1

@@ -9 +9 @@
 
 # SETTE
-^/utils/CI/sette@14244        sette
+^/utils/CI/sette@HEAD        sette
+

@@ -9 +9 @@
 
 # SETTE
-^/utils/CI/sette@14244        sette
+^/utils/CI/sette@HEAD        sette
+

NEMO/branches/2021/dev_r14318_RK3_stage1/src/ICE/icealb.F90

@@ -30 +30 @@
    PUBLIC   ice_alb        ! called in icesbc.F90 and iceupdate.F90
 
-   REAL(wp), PUBLIC, PARAMETER ::   rn_alb_oce = 0.066   !: ocean or lead albedo (Pegau and Paulson, Ann. Glac. 2001)
+   REAL(wp), PUBLIC, PARAMETER ::   rn_alb_oce = 0.066_wp   !: ocean or lead albedo (Pegau and Paulson, Ann. Glac. 2001)
    !
    !                             !!* albedo namelist (namalb)
@@ -38 +38 @@
    REAL(wp) ::   rn_alb_imlt      ! bare puddled ice albedo
    REAL(wp) ::   rn_alb_dpnd      ! ponded ice albedo
+   REAL(wp) ::   rn_alb_hpiv      ! pivotal ice thickness in meters (above which albedo is constant)
 
    !! * Substitutions
@@ -59 +60 @@
       !!                     which are an update of Allison et al. (JGR 1993) ; Brandt et al. 1999
       !!                     0-5cm  : linear function of ice thickness
-      !!                     5-150cm: log    function of ice thickness
-      !!                     > 150cm: constant
+      !!                     5-100cm: log    function of ice thickness
+      !!                     > 100cm: constant
       !!                  2) Albedo dependency on snow thickness follows the findings from Grenfell & Perovich (2004)
       !!                     i.e. it increases as -EXP(-snw_thick/0.02) during freezing and -EXP(-snw_thick/0.03) during melting
@@ -115 +116 @@
       IF( ln_timing )   CALL timing_start('icealb')
       !
-      z1_href_pnd = 1. / 0.05
-      z1_c1 = 1. / ( LOG(1.5) - LOG(0.05) ) 
-      z1_c2 = 1. / 0.05
-      z1_c3 = 1. / 0.02
-      z1_c4 = 1. / 0.03
+      z1_href_pnd = 1._wp / 0.05_wp
+      z1_c1 = 1._wp / ( LOG(rn_alb_hpiv) - LOG(0.05_wp) ) 
+      z1_c2 = 1._wp / 0.05_wp
+      z1_c3 = 1._wp / 0.02_wp
+      z1_c4 = 1._wp / 0.03_wp
       !
       CALL ice_var_snwfra( ph_snw, za_s_fra )   ! calculate ice fraction covered by snow
       !
       DO jl = 1, jpl
-         DO_2D( 1, 1, 1, 1 )
+         DO_2D( nn_hls, nn_hls, nn_hls, nn_hls )   ! palb_ice used over the full domain in icesbc
             !
             !---------------------------------------------!
@@ -140 +141 @@
             !--- Albedos ---!
             !---------------!               
-            !                       !--- Bare ice albedo (for hi > 150cm)
+            !                       !--- Bare ice albedo (for hi > 100cm)
             IF( ld_pnd_alb ) THEN
                zalb_ice = rn_alb_idry
@@ -147 +148 @@
                ELSE                                                                ;   zalb_ice = rn_alb_idry   ;   ENDIF
             ENDIF
-            !                       !--- Bare ice albedo (for hi < 150cm)
-            IF( 0.05 < ph_ice(ji,jj,jl) .AND. ph_ice(ji,jj,jl) <= 1.5 ) THEN      ! 5cm < hi < 150cm
-               zalb_ice = zalb_ice    + ( 0.18 - zalb_ice   ) * z1_c1 * ( LOG(1.5) - LOG(ph_ice(ji,jj,jl)) )
-            ELSEIF( ph_ice(ji,jj,jl) <= 0.05 ) THEN                               ! 0cm < hi < 5cm
-               zalb_ice = rn_alb_oce  + ( 0.18 - rn_alb_oce ) * z1_c2 * ph_ice(ji,jj,jl)
+            !                       !--- Bare ice albedo (for hi < 100cm)
+            IF( 0.05 < ph_ice(ji,jj,jl) .AND. ph_ice(ji,jj,jl) <= rn_alb_hpiv ) THEN      ! 5cm < hi < 100cm
+               zalb_ice = zalb_ice    + ( 0.18_wp - zalb_ice   ) * z1_c1 * ( LOG(rn_alb_hpiv) - LOG(ph_ice(ji,jj,jl)) )
+            ELSEIF( ph_ice(ji,jj,jl) <= 0.05_wp ) THEN                                    ! 0cm < hi < 5cm
+               zalb_ice = rn_alb_oce  + ( 0.18_wp - rn_alb_oce ) * z1_c2 * ph_ice(ji,jj,jl)
             ENDIF
             !
@@ -166 +167 @@
             zalb_os = ( zafrac_snw * zalb_snw + zafrac_pnd * zalb_pnd + zafrac_ice * zalb_ice ) * tmask(ji,jj,1)
             !
-            zalb_cs = zalb_os - ( - 0.1010 * zalb_os * zalb_os  &
-               &                  + 0.1933 * zalb_os - 0.0148 ) * tmask(ji,jj,1)
+            zalb_cs = zalb_os - ( - 0.1010_wp * zalb_os * zalb_os  &
+               &                  + 0.1933_wp * zalb_os - 0.0148_wp ) * tmask(ji,jj,1)
             !
             ! albedo depends on cloud fraction because of non-linear spectral effects
@@ -191 +192 @@
       INTEGER ::   ios   ! Local integer output status for namelist read
       !!
-      NAMELIST/namalb/ rn_alb_sdry, rn_alb_smlt, rn_alb_idry, rn_alb_imlt, rn_alb_dpnd
+      NAMELIST/namalb/ rn_alb_sdry, rn_alb_smlt, rn_alb_idry, rn_alb_imlt, rn_alb_dpnd, rn_alb_hpiv
       !!----------------------------------------------------------------------
       !
@@ -210 +211 @@
          WRITE(numout,*) '      albedo of bare puddled ice           rn_alb_imlt = ', rn_alb_imlt
          WRITE(numout,*) '      albedo of ponded ice                 rn_alb_dpnd = ', rn_alb_dpnd
+         WRITE(numout,*) '      pivotal ice thickness (m)            rn_alb_hpiv = ', rn_alb_hpiv
       ENDIF
       !

@@ -9 +9 @@
 
 # SETTE
-^/utils/CI/sette@14244        sette
+^/utils/CI/sette@HEAD        sette
+