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 13540 for NEMO/branches/2020/r12377_ticket2386/src/OCE/DYN/wet_dry.F90 – NEMO

Ignore:
Timestamp:
2020-09-29T12:41:06+02:00 (4 years ago)
Author:
andmirek
Message:

Ticket #2386: update to latest trunk

Location:
NEMO/branches/2020/r12377_ticket2386
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/r12377_ticket2386

    • Property svn:externals
      •  

        old new  
        33^/utils/build/mk@HEAD         mk 
        44^/utils/tools@HEAD            tools 
        5 ^/vendors/AGRIF/dev@HEAD      ext/AGRIF 
         5^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS      ext/AGRIF 
        66^/vendors/FCM@HEAD            ext/FCM 
        77^/vendors/IOIPSL@HEAD         ext/IOIPSL 
        88 
        99# SETTE 
        10 ^/utils/CI/sette@HEAD         sette 
         10^/utils/CI/sette@13507        sette 
  • NEMO/branches/2020/r12377_ticket2386/src/OCE/DYN/wet_dry.F90

    r12511 r13540  
    3333   !! * Substitutions 
    3434#  include "do_loop_substitute.h90" 
     35#  include "domzgr_substitute.h90" 
    3536   !!---------------------------------------------------------------------- 
    3637   !! critical depths,filters, limiters,and masks for  Wetting and Drying 
     
    173174      ! 
    174175      wdmask(:,:) = 1._wp 
    175       DO_2D_01_01 
     176      DO_2D( 0, 1, 0, 1 ) 
    176177         ! 
    177178         IF( tmask(ji,jj,1)        < 0.5_wp )   CYCLE    ! we don't care about land cells 
     
    197198      wdramp(:,:) = min((ht_0(:,:) + psshb1(:,:) - rn_wdmin1)/(rn_wdmin0 - rn_wdmin1),1.0_wp) 
    198199      !jth assume don't need a lbc_lnk here 
    199       DO_2D_10_10 
     200      DO_2D( 1, 0, 1, 0 ) 
    200201         wdrampu(ji,jj) = MIN( wdramp(ji,jj) , wdramp(ji+1,jj) ) 
    201202         wdrampv(ji,jj) = MIN( wdramp(ji,jj) , wdramp(ji,jj+1) ) 
     
    210211         jflag = 0     ! flag indicating if any further iterations are needed 
    211212         ! 
    212          DO_2D_01_01 
     213         DO_2D( 0, 1, 0, 1 ) 
    213214            IF( tmask(ji, jj, 1) < 0.5_wp )   CYCLE  
    214215            IF( ht_0(ji,jj)      > zdepwd )   CYCLE 
     
    241242            ENDIF 
    242243         END_2D 
    243          CALL lbc_lnk_multi( 'wet_dry', zwdlmtu, 'U', 1., zwdlmtv, 'V', 1. ) 
     244         CALL lbc_lnk_multi( 'wet_dry', zwdlmtu, 'U', 1.0_wp, zwdlmtv, 'V', 1.0_wp ) 
    244245         ! 
    245246         CALL mpp_max('wet_dry', jflag)   !max over the global domain 
     
    257258      ! 
    258259!!gm TO BE SUPPRESSED ?  these lbc_lnk are useless since zwdlmtu and zwdlmtv are defined everywhere ! 
    259       CALL lbc_lnk_multi( 'wet_dry', puu(:,:,:,Kmm)  , 'U', -1., pvv(:,:,:,Kmm)  , 'V', -1. ) 
    260       CALL lbc_lnk_multi( 'wet_dry', uu_b(:,:,Kmm), 'U', -1., vv_b(:,:,Kmm), 'V', -1. ) 
     260      CALL lbc_lnk_multi( 'wet_dry', puu(:,:,:,Kmm)  , 'U', -1.0_wp, pvv(:,:,:,Kmm)  , 'V', -1.0_wp ) 
     261      CALL lbc_lnk_multi( 'wet_dry', uu_b(:,:,Kmm), 'U', -1.0_wp, vv_b(:,:,Kmm), 'V', -1.0_wp ) 
    261262!!gm 
    262263      ! 
     
    306307      zwdlmtv(:,:) = 1._wp 
    307308      ! 
    308       DO_2D_01_01 
     309      DO_2D( 0, 1, 0, 1 )      ! Horizontal Flux in u and v direction 
    309310         ! 
    310311         IF( tmask(ji, jj, 1 ) < 0.5_wp) CYCLE   ! we don't care about land cells 
     
    332333         jflag = 0     ! flag indicating if any further iterations are needed 
    333334         ! 
    334          DO_2D_01_01 
     335         DO_2D( 0, 1, 0, 1 ) 
    335336            ! 
    336337            IF( tmask(ji, jj, 1 ) < 0.5_wp )   CYCLE  
     
    366367         END_2D 
    367368         ! 
    368          CALL lbc_lnk_multi( 'wet_dry', zwdlmtu, 'U', 1., zwdlmtv, 'V', 1. ) 
     369         CALL lbc_lnk_multi( 'wet_dry', zwdlmtu, 'U', 1.0_wp, zwdlmtv, 'V', 1.0_wp ) 
    369370         ! 
    370371         CALL mpp_max('wet_dry', jflag)   !max over the global domain 
     
    378379      ! 
    379380!!gm THIS lbc_lnk is useless since it is already done at the end of the jk1-loop 
    380       CALL lbc_lnk_multi( 'wet_dry', zflxu, 'U', -1., zflxv, 'V', -1. ) 
     381      CALL lbc_lnk_multi( 'wet_dry', zflxu, 'U', -1.0_wp, zflxv, 'V', -1.0_wp ) 
    381382!!gm end 
    382383      ! 
Note: See TracChangeset for help on using the changeset viewer.