New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 8085 for branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limmp.F90 – NEMO

Ignore:
Timestamp:
2017-05-29T16:58:42+02:00 (7 years ago)
Author:
vancop
Message:

More on melt ponds in LIM

File:
1 edited

Legend:

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

    r8063 r8085  
    101101         WRITE(numout,*)'    Activate melt ponds                                         ln_pnd        = ', ln_pnd 
    102102         WRITE(numout,*)'    Type of melt pond scheme =0 presc, =1 empirical = 2 topo    nn_pnd_scheme = ', nn_pnd_scheme 
    103          WRITE(numout,*)'    Type of melt pond coupling =0 pass., =1 full, =2 rad, 3=fw  nn_pnd_cpl    = ', nn_pnd_cp 
     103         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      ENDIF 
     106 
     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 ' 
     110         nn_pnd_scheme = 0 
     111         nn_pnd_cpl    = 0 
     112         rn_apnd       = 0._wp 
    105113      ENDIF 
    106114      ! 
     
    134142      SELECT CASE ( nn_pnd_scheme ) 
    135143 
    136           CASE (1) 
    137  
    138              CALL lim_mp_cesm ! empirical melt ponds 
    139  
    140           CASE (2) 
    141  
    142              CALL lim_mp_topo    (at_i, a_i,                                       & 
    143                        &          vt_i, v_i, v_s,            t_i, s_i, a_ip_frac,  & 
    144                        &          h_ip,     t_su) 
     144         CASE (1) 
     145 
     146            CALL lim_mp_cesm ! empirical melt ponds 
     147 
     148         CASE (2) 
     149 
     150            CALL lim_mp_topo    (at_i, a_i,                                       & 
     151                      &          vt_i, v_i, v_s,            t_i, s_i, a_ip_frac,  & 
     152                      &          h_ip,     t_su) 
    145153 
    146154      END SELECT 
     
    193201       REAL(wp)                            :: z1_rhofw          ! inverse freshwater density 
    194202       REAL(wp)                            :: z1_zpnd_aspect    ! inverse pond aspect ratio 
    195        REAL(wp)                            :: zvold             ! dummy pond volume 
     203       REAL(wp)                            :: zvpold            ! dummy pond volume 
    196204 
    197205       INTEGER                             :: ji, jj, jl, ij    ! loop indices 
     
    261269            ELSE                        !--- Update pond characteristics 
    262270 
    263                !--- Add retained melt water 
     271               !--- Add retained melt water to melt ponds 
    264272               v_ip(ji,jj,jl)      = v_ip(ji,jj,jl) + zrfrac(ji,jj,jl) * z1_rhofw * zwfx_mlw(ji,jj) * a_i(ji,jj,jl) * rdt_ice 
    265273 
Note: See TracChangeset for help on using the changeset viewer.