Ignore:
Timestamp:
2017-05-30T13:17:53+02:00 (4 years ago)
Author:
vancop
Message:

Further melt pond work

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limmp.F90

    r8085 r8098  
    8383      !!------------------------------------------------------------------- 
    8484      INTEGER  ::   ios                 ! Local integer output status for namelist read 
    85       NAMELIST/namicemp/  ln_pnd, nn_pnd_scheme, nn_pnd_cpl, rn_apnd 
     85      NAMELIST/namicemp/  ln_pnd, nn_pnd_scheme, nn_pnd_cpl, rn_apnd, rn_hpnd 
    8686      !!------------------------------------------------------------------- 
    8787 
     
    103103         WRITE(numout,*)'    Type of melt pond coupling =0 pass., =1 full, =2 rad, 3=fw  nn_pnd_cpl    = ', nn_pnd_cpl 
    104104         WRITE(numout,*)'    Prescribed pond fraction                                    rn_apnd       = ', rn_apnd 
     105         WRITE(numout,*)'    Prescribed pond depth                                       rn_hpnd       = ', rn_hpnd 
    105106      ENDIF 
    106107 
    107       IF ( ln_pnd == .FALSE. ) THEN 
    108          WRITE(numout) ' Melt ponds are not activated ' 
    109          WRITE(numout) ' nn_pnd_scheme, nn_pnd_cpl and rn_apnd are set to zero ' 
     108      IF ( .NOT. ln_pnd ) THEN 
     109         WRITE(numout,*) ' Melt ponds are not activated ' 
     110         WRITE(numout,*) ' nn_pnd_scheme, nn_pnd_cpl, rn_apnd and rn_hpnd are set to zero ' 
    110111         nn_pnd_scheme = 0 
    111112         nn_pnd_cpl    = 0 
    112113         rn_apnd       = 0._wp 
     114         rn_hpnd       = 0._wp 
    113115      ENDIF 
     116 
    114117      ! 
    115118   END SUBROUTINE lim_mp_init 
    116  
     119    
    117120 
    118121 
     
    142145      SELECT CASE ( nn_pnd_scheme ) 
    143146 
     147         CASE (0) 
     148 
     149            CALL lim_mp_cstt       ! staircase melt ponds 
     150 
    144151         CASE (1) 
    145152 
    146             CALL lim_mp_cesm ! empirical melt ponds 
     153            CALL lim_mp_cesm       ! empirical melt ponds 
    147154 
    148155         CASE (2) 
    149156 
    150             CALL lim_mp_topo    (at_i, a_i,                                       & 
    151                       &          vt_i, v_i, v_s,            t_i, s_i, a_ip_frac,  & 
     157            CALL lim_mp_topo   &   ! topographic melt ponds  
     158                      &          (at_i, a_i,                                       & 
     159                      &          vt_i, v_i, v_s,            t_i, s_i, a_ip_frac,   & 
    152160                      &          h_ip,     t_su) 
    153161 
     
    162170   END SUBROUTINE lim_mp  
    163171 
    164  
     172   SUBROUTINE lim_mp_cstt  
     173       !!------------------------------------------------------------------- 
     174       !!                ***  ROUTINE lim_mp_cstt  *** 
     175       !! 
     176       !! ** Purpose    : Compute melt pond evolution 
     177       !! 
     178       !! ** Method     : Melt pond fraction and thickness are prescribed  
     179       !!                 to non-zero values when t_su = 0C 
     180       !! 
     181       !! ** Tunable parameters : pond fraction (rn_apnd), pond depth (rn_hpnd) 
     182       !!                 
     183       !! ** Note       : Coupling with such melt ponds is only radiative 
     184       !!                 Advection, ridging, rafting... are bypassed 
     185       !! 
     186       !! ** References : Bush, G.W., and Trump, D.J. (2017) 
     187       !!     
     188       !!------------------------------------------------------------------- 
     189       INTEGER                             :: ji, jj, jl 
     190       REAL(wp)                            :: z1_jpl            ! 1/jpl 
     191       !!------------------------------------------------------------------- 
     192 
     193       z1_jpl     = 1. / FLOAT(jpl) 
     194 
     195       WHERE ( ( a_i > epsi10 ) .AND. ( t_su >= rt0-epsi06 ) )  
     196          a_ip      = rn_apnd * z1_jpl  
     197          h_ip      = rn_hpnd     
     198          v_ip      = a_ip * h_ip  
     199          a_ip_frac = a_ip / a_i 
     200       ELSE WHERE 
     201          a_ip      = 0._wp 
     202          h_ip      = 0._wp 
     203          v_ip      = 0._wp 
     204          a_ip_frac = 0._wp 
     205       END WHERE 
     206 
     207       wfx_pnd(:,:) = 0._wp 
     208 
     209   END SUBROUTINE lim_mp_cstt 
    165210 
    166211   SUBROUTINE lim_mp_cesm 
     
    172217       !! ** Method     : Empirical method. A fraction of meltwater is accumulated  
    173218       !!                 in pond volume. It is then released exponentially when 
    174        !!                 surface is freezingAccumulation of meltwater and exponential release 
     219       !!                 surface is freezing. 
    175220       !! 
    176        !! ** Tunable parameters : 
    177        !!                 
     221       !! ** Tunable parameters : (no expertise yet) 
    178222       !!  
    179223       !! ** Note       : Stolen from CICE for quick test of the melt pond 
    180224       !!                 radiation and freshwater interfaces 
     225       !!                 Coupling can be radiative AND freshwater 
     226       !!                 Advection, ridging, rafting are called 
    181227       !! 
    182228       !! ** References : Holland, M. M. et al (J Clim 2012) 
Note: See TracChangeset for help on using the changeset viewer.