Ignore:
Timestamp:
2020-07-09T12:39:43+02:00 (3 months ago)
Author:
clem
Message:

merge with r4.0-HEAD at r13278

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icedyn_rhg_evp.F90

    r13064 r13279  
    160160      REAL(wp), DIMENSION(jpi,jpj) ::   zmsk01x, zmsk01y                ! dummy arrays 
    161161      REAL(wp), DIMENSION(jpi,jpj) ::   zmsk00x, zmsk00y                ! mask for ice presence 
    162       REAL(wp), DIMENSION(jpi,jpj) ::   zfmask, zwf                     ! mask at F points for the ice 
     162      REAL(wp), DIMENSION(jpi,jpj) ::   zfmask                          ! mask at F points for the ice 
    163163 
    164164      REAL(wp), PARAMETER          ::   zepsi  = 1.0e-20_wp             ! tolerance parameter 
     
    202202 
    203203      ! Lateral boundary conditions on velocity (modify zfmask) 
    204       zwf(:,:) = zfmask(:,:) 
    205204      DO jj = 2, jpjm1 
    206205         DO ji = fs_2, fs_jpim1   ! vector opt. 
    207206            IF( zfmask(ji,jj) == 0._wp ) THEN 
    208                zfmask(ji,jj) = rn_ishlat * MIN( 1._wp , MAX( zwf(ji+1,jj), zwf(ji,jj+1), zwf(ji-1,jj), zwf(ji,jj-1) ) ) 
     207               zfmask(ji,jj) = rn_ishlat * MIN( 1._wp , MAX( umask(ji,jj,1), umask(ji,jj+1,1), & 
     208                  &                                          vmask(ji,jj,1), vmask(ji+1,jj,1) ) ) 
    209209            ENDIF 
    210210         END DO 
     
    212212      DO jj = 2, jpjm1 
    213213         IF( zfmask(1,jj) == 0._wp ) THEN 
    214             zfmask(1  ,jj) = rn_ishlat * MIN( 1._wp , MAX( zwf(2,jj), zwf(1,jj+1), zwf(1,jj-1) ) ) 
     214            zfmask(1  ,jj) = rn_ishlat * MIN( 1._wp , MAX( vmask(2,jj,1), umask(1,jj+1,1), umask(1,jj,1) ) ) 
    215215         ENDIF 
    216216         IF( zfmask(jpi,jj) == 0._wp ) THEN 
    217             zfmask(jpi,jj) = rn_ishlat * MIN( 1._wp , MAX( zwf(jpi,jj+1), zwf(jpim1,jj), zwf(jpi,jj-1) ) ) 
    218          ENDIF 
     217            zfmask(jpi,jj) = rn_ishlat * MIN( 1._wp , MAX( umask(jpi,jj+1,1), vmask(jpim1,jj,1), umask(jpi,jj-1,1) ) ) 
     218        ENDIF 
    219219      END DO 
    220220      DO ji = 2, jpim1 
    221221         IF( zfmask(ji,1) == 0._wp ) THEN 
    222             zfmask(ji,1  ) = rn_ishlat * MIN( 1._wp , MAX( zwf(ji+1,1), zwf(ji,2), zwf(ji-1,1) ) ) 
     222            zfmask(ji, 1 ) = rn_ishlat * MIN( 1._wp , MAX( vmask(ji+1,1,1), umask(ji,2,1), vmask(ji,1,1) ) ) 
    223223         ENDIF 
    224224         IF( zfmask(ji,jpj) == 0._wp ) THEN 
    225             zfmask(ji,jpj) = rn_ishlat * MIN( 1._wp , MAX( zwf(ji+1,jpj), zwf(ji-1,jpj), zwf(ji,jpjm1) ) ) 
     225            zfmask(ji,jpj) = rn_ishlat * MIN( 1._wp , MAX( vmask(ji+1,jpj,1), vmask(ji-1,jpj,1), umask(ji,jpjm1,1) ) ) 
    226226         ENDIF 
    227227      END DO 
Note: See TracChangeset for help on using the changeset viewer.