Changeset 10518


Ignore:
Timestamp:
2019-01-15T18:06:06+01:00 (18 months ago)
Author:
clem
Message:

solve ticket #1997 as proposed

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/BDY/bdydyn2d.F90

    r10425 r10518  
    270270            zcoef1 = bdytmask(ii-1,ij  ) +  bdytmask(ii+1,ij  ) 
    271271            zcoef2 = bdytmask(ii  ,ij-1) +  bdytmask(ii  ,ij+1) 
    272             IF ( zcoef1+zcoef2 == 0 ) THEN 
    273                ! corner 
    274 !               zcoef = tmask(ii-1,ij,1) + tmask(ii+1,ij,1) +  tmask(ii,ij-1,1) +  tmask(ii,ij+1,1) 
    275 !               zssh(ii,ij) = zssh(ii-1,ij  ) * tmask(ii-1,ij  ,1) + & 
    276 !                 &           zssh(ii+1,ij  ) * tmask(ii+1,ij  ,1) + & 
    277 !                 &           zssh(ii  ,ij-1) * tmask(ii  ,ij-1,1) + & 
    278 !                 &           zssh(ii  ,ij+1) * tmask(ii  ,ij+1,1) 
    279                zcoef = bdytmask(ii-1,ij) + bdytmask(ii+1,ij) +  bdytmask(ii,ij-1) +  bdytmask(ii,ij+1) 
    280                zssh(ii,ij) = zssh(ii-1,ij  ) * bdytmask(ii-1,ij  ) + & 
    281                  &           zssh(ii+1,ij  ) * bdytmask(ii+1,ij  ) + & 
    282                  &           zssh(ii  ,ij-1) * bdytmask(ii  ,ij-1) + & 
    283                  &           zssh(ii  ,ij+1) * bdytmask(ii  ,ij+1) 
     272            IF ( zcoef1+zcoef2 == 0 ) THEN   ! corner 
     273               zcoef = bdytmask(ii-1,ij-1) + bdytmask(ii+1,ij+1) + bdytmask(ii+1,ij-1) + bdytmask(ii-1,ij+1) 
     274               zssh(ii,ij) = zssh( ii-1, ij-1 ) * bdytmask( ii-1, ij-1) + & 
     275                 &           zssh( ii+1, ij+1 ) * bdytmask( ii+1, ij+1) + & 
     276                 &           zssh( ii+1, ij-1 ) * bdytmask( ii+1, ij-1) + & 
     277                 &           zssh( ii-1, ij+1 ) * bdytmask( ii-1, ij+1) 
    284278               zssh(ii,ij) = ( zssh(ii,ij) / MAX( 1, zcoef) ) * tmask(ii,ij,1) 
    285279            ELSE 
Note: See TracChangeset for help on using the changeset viewer.