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 13811 – NEMO

Changeset 13811


Ignore:
Timestamp:
2020-11-18T13:05:56+01:00 (3 years ago)
Author:
vancop
Message:

pre-test version

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/SI3-05_MP/src/ICE/icethd_pnd.F90

    r13810 r13811  
    5353      !!                
    5454      !! ** Purpose :   change melt pond fraction and thickness 
     55      !! 
     56      !! Note: Melt ponds affect currently affect radiative transfer  
     57      !!       They carry no heat, and the melt water they carry is not 
     58      !!       exchanged with the ocean 
     59      !! 
     60      !!                   This means freshwater is directly released after surface & bottom melt in ice_thd_dh 
     61      !! 
     62      !!       A wfx_pnd has been coded for diagnostic purposes 
     63      !!       Each time wfx_pnd is updated, wfx_sum / wfx_snw_sum must be updated 
     64      !! 
     65      !!!      The current diagnostic lacks a contribution from drainage 
    5566      !!                 
    5667      !!------------------------------------------------------------------- 
     
    226237            !--- mass flux ---! 
    227238            IF( zdv_mlt > 0._wp ) THEN 
    228                ! MV add comment on what that mass flux means 
    229                ! water removed from fw flux due to melt pond growth 
     239             ! MV add comment on what that mass flux means 
    230240               zfac = zdv_mlt * rhow * r1_rdtice                        ! melt pond mass flux < 0 [kg.m-2.s-1] 
    231241               wfx_pnd_1d(ji) = wfx_pnd_1d(ji) - zfac 
     
    277287             
    278288            ! Calculate the permeability of the ice (Assur 1958, see Flocco 2010) 
    279             ! zsbr = - ( t_i_1d(ji,jk) - rt0 ) / rTmlt 
    280289            DO jk = 1, nlay_i 
    281290               zsbr = - 1.2_wp                                  & 
     
    293302            zdv_flush   = MAX( zdv_flush, -v_ip_1d(ji) ) 
    294303            v_ip_1d(ji) = v_ip_1d(ji) + zdv_flush 
     304             
     305            ! MV --- why pond drainage does not give back water into freshwater flux ? 
    295306             
    296307            !--- Set new pond area and depth ---! assuming linear relation between h_ip and a_ip_frac 
     
    407418 
    408419       REAL (wp), PARAMETER :: & 
     420          ! MV ouate de phoque!!! constants hard coded ???? 
     421          ! 917 = rhoice 
     422          ! 0.334 =Lfus!!!! 
    409423          zr1_rlfus = 1._wp / 0.334e+6 / 917._wp , & ! (J/m^3) 
    410424          zTp       = -0.15_wp, & ! pond freezing temperature (C) 
     
    474488 
    475489      !js 23/04/19: rewired back to a fraction with a_i 
    476       zrfrac(:,:,:) = rn_pnd_fracmin + ( rn_pnd_fracmax - rn_pnd_fracmin ) * aicen(:,:,:) 
     490      !!! zrfrac(:,:,:) = rn_pnd_fracmin + ( rn_pnd_fracmax - rn_pnd_fracmin ) * aicen(:,:,:) 
     491      ! MV2020 
     492      zrfrac(:,:,:) = rn_apnd_min + ( rn_apnd_max - rn_apnd_min ) * at_i_1d(:) !  = ( 1 - r ) = fraction of melt water that is not flushed 
    477493      zwfx_tmp(:,:) = 0._wp 
    478  
    479 ! MV ---> use expression from level-ice ponds, clarify what is needed 
    480494 
    481495      !--- Add retained melt water to melt ponds 
     
    496510                                    * ( dh_i_pnd(:,:,jl)*rhoic + dh_s_pnd(:,:,jl)*rhosn )          &  
    497511                                    * z1_rhow * a_i(:,:,jl) 
     512! MV ---> use expression from level ice melt ponds (dv_mlt) 
     513                                             
    498514      END DO 
    499515 
    500 ! MV ---> use expression from level ice melt ponds (dv_mlt) 
    501  
     516 
     517!!! MV ---> rewrite this accounting for level ice melt ponds  
     518!!! 1) wfx_pnd_in is obsolete, ln_pnd_fw as well 
     519!!! 2)  
     520!!! wfx_sum 
     521      !!! 
     522      !!! 
    502523      !js 03/05/19: we truncate negative values after calculating zvolp, in a 
    503524      ! similar manner to the subroutine ice_thd_pnd_cesm. Variation dh_i_pnd and 
     
    531552         wfx_snw_sum(:,:) = ( 1._wp - zwfx_tmp(:,:) ) * wfx_snw_sum(:,:) 
    532553      ENDIF 
     554       
     555       
    533556 
    534557       !----------------------------------------------------------------- 
     
    568591      END DO                     ! jj 
    569592 
     593 
     594 
     595!!! MV sounds like this should be replaced by lid 
    570596 
    571597       DO ij = 1, icells 
Note: See TracChangeset for help on using the changeset viewer.