Ignore:
Timestamp:
2020-12-24T09:07:31+01:00 (9 months ago)
Author:
vancop
Message:

bug fixes in topo ponds for trunk

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/ICE/icethd_pnd.F90

    r14072 r14252  
    511511      !!------------------------------------------------------------------- 
    512512      REAL(wp), PARAMETER :: &   ! shared parameters for topographic melt ponds 
    513          zTd     = 0.15_wp       , & ! temperature difference for freeze-up (C) 
     513         zTd     = 273._wp       , & ! temperature difference for freeze-up (K) 
    514514         zvp_min = 1.e-4_wp          ! minimum pond volume (m) 
    515515 
     
    530530         zv_mlt          ! total amount of meltwater produced 
    531531 
    532       REAL(wp), DIMENSION(jpi,jpj) ::   zvolp, &     !! total melt pond water available before redistribution and drainage 
    533                                         zvolp_res    !! remaining melt pond water available after drainage 
     532      REAL(wp), DIMENSION(jpi,jpj) ::   zvolp_ini , &   !! total melt pond water available before redistribution and drainage 
     533                                        zvolp     , &   !! total melt pond water volume 
     534                                        zvolp_res       !! remaining melt pond water available after drainage 
    534535 
    535536      REAL(wp), DIMENSION(jpi,jpj,jpl) ::   z1_a_i 
     
    624625      END DO ! ji 
    625626 
     627      zvolp_ini(:,:) = zvolp(:,:) 
     628 
    626629      !-------------------------------------------------------------- 
    627630      ! Redistribute and drain water from ponds 
     
    637640         DO_2D( 1, 1, 1, 1 ) 
    638641 
    639             IF ( at_i(ji,jj) > 0.01 .AND. hm_i(ji,jj) > rn_himin .AND. vt_ip(ji,jj) > zvp_min * at_i(ji,jj) ) THEN 
     642            IF ( at_i(ji,jj) > 0.01 .AND. hm_i(ji,jj) > rn_himin .AND. zvolp_ini(ji,jj) > zvp_min * at_i(ji,jj) ) THEN 
    640643 
    641644               !-------------------------- 
     
    688691 
    689692                        ! differential growth of base of surface floating ice layer 
    690                         zdTice = MAX( - t_su(ji,jj,jl) - zTd , 0._wp ) ! > 0 
     693                        zdTice = MAX( - ( t_su(ji,jj,jl) - zTd ) , 0._wp ) ! > 0 
    691694                        zomega = rcnd_i * zdTice / zrhoi_L 
    692695                        zdHui  = SQRT( 2._wp * zomega * rDt_ice + ( v_il(ji,jj,jl) / a_i(ji,jj,jl) )**2 ) & 
Note: See TracChangeset for help on using the changeset viewer.