Changeset 8098 for branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limmp.F90

Timestamp:

2017-05-30T13:17:53+02:00 (7 years ago)

Author:

vancop

Message:

Further melt pond work

File:

• branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limmp.F90 (modified) (5 diffs)

branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limmp.F90

@@ -83 +83 @@
       !!-------------------------------------------------------------------
       INTEGER  ::   ios                 ! Local integer output status for namelist read
-      NAMELIST/namicemp/  ln_pnd, nn_pnd_scheme, nn_pnd_cpl, rn_apnd
+      NAMELIST/namicemp/  ln_pnd, nn_pnd_scheme, nn_pnd_cpl, rn_apnd, rn_hpnd
       !!-------------------------------------------------------------------
 
@@ -103 +103 @@
          WRITE(numout,*)'    Type of melt pond coupling =0 pass., =1 full, =2 rad, 3=fw  nn_pnd_cpl    = ', nn_pnd_cpl
          WRITE(numout,*)'    Prescribed pond fraction                                    rn_apnd       = ', rn_apnd
+         WRITE(numout,*)'    Prescribed pond depth                                       rn_hpnd       = ', rn_hpnd
       ENDIF
 
-      IF ( ln_pnd == .FALSE. ) THEN
-         WRITE(numout) ' Melt ponds are not activated '
-         WRITE(numout) ' nn_pnd_scheme, nn_pnd_cpl and rn_apnd are set to zero '
+      IF ( .NOT. ln_pnd ) THEN
+         WRITE(numout,*) ' Melt ponds are not activated '
+         WRITE(numout,*) ' nn_pnd_scheme, nn_pnd_cpl, rn_apnd and rn_hpnd are set to zero '
          nn_pnd_scheme = 0
          nn_pnd_cpl    = 0
          rn_apnd       = 0._wp
+         rn_hpnd       = 0._wp
       ENDIF
+
       !
    END SUBROUTINE lim_mp_init
-
+   
 
 
@@ -142 +145 @@
       SELECT CASE ( nn_pnd_scheme )
 
+         CASE (0)
+
+            CALL lim_mp_cstt       ! staircase melt ponds
+
          CASE (1)
 
-            CALL lim_mp_cesm ! empirical melt ponds
+            CALL lim_mp_cesm       ! empirical melt ponds
 
          CASE (2)
 
-            CALL lim_mp_topo    (at_i, a_i,                                       &
-                      &          vt_i, v_i, v_s,            t_i, s_i, a_ip_frac,  &
+            CALL lim_mp_topo   &   ! topographic melt ponds 
+                      &          (at_i, a_i,                                       &
+                      &          vt_i, v_i, v_s,            t_i, s_i, a_ip_frac,   &
                       &          h_ip,     t_su)
 
@@ -162 +170 @@
    END SUBROUTINE lim_mp 
 
-
+   SUBROUTINE lim_mp_cstt 
+       !!-------------------------------------------------------------------
+       !!                ***  ROUTINE lim_mp_cstt  ***
+       !!
+       !! ** Purpose    : Compute melt pond evolution
+       !!
+       !! ** Method     : Melt pond fraction and thickness are prescribed 
+       !!                 to non-zero values when t_su = 0C
+       !!
+       !! ** Tunable parameters : pond fraction (rn_apnd), pond depth (rn_hpnd)
+       !!                
+       !! ** Note       : Coupling with such melt ponds is only radiative
+       !!                 Advection, ridging, rafting... are bypassed
+       !!
+       !! ** References : Bush, G.W., and Trump, D.J. (2017)
+       !!    
+       !!-------------------------------------------------------------------
+       INTEGER                             :: ji, jj, jl
+       REAL(wp)                            :: z1_jpl            ! 1/jpl
+       !!-------------------------------------------------------------------
+
+       z1_jpl     = 1. / FLOAT(jpl)
+
+       WHERE ( ( a_i > epsi10 ) .AND. ( t_su >= rt0-epsi06 ) ) 
+          a_ip      = rn_apnd * z1_jpl 
+          h_ip      = rn_hpnd    
+          v_ip      = a_ip * h_ip 
+          a_ip_frac = a_ip / a_i
+       ELSE WHERE
+          a_ip      = 0._wp
+          h_ip      = 0._wp
+          v_ip      = 0._wp
+          a_ip_frac = 0._wp
+       END WHERE
+
+       wfx_pnd(:,:) = 0._wp
+
+   END SUBROUTINE lim_mp_cstt
 
    SUBROUTINE lim_mp_cesm
@@ -172 +217 @@
        !! ** Method     : Empirical method. A fraction of meltwater is accumulated 
        !!                 in pond volume. It is then released exponentially when
-       !!                 surface is freezingAccumulation of meltwater and exponential release
+       !!                 surface is freezing.
        !!
-       !! ** Tunable parameters :
-       !!                
+       !! ** Tunable parameters : (no expertise yet)
        !! 
        !! ** Note       : Stolen from CICE for quick test of the melt pond
        !!                 radiation and freshwater interfaces
+       !!                 Coupling can be radiative AND freshwater
+       !!                 Advection, ridging, rafting are called
        !!
        !! ** References : Holland, M. M. et al (J Clim 2012)