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.
Ticket Diff – NEMO

Changes between Initial Version and Version 1 of Ticket #2237


Ignore:
Timestamp:
2019-02-07T16:36:35+01:00 (5 years ago)
Author:
nicolasmartin
Comment:

Legend:

Unmodified
Added
Removed
Modified

  • Ticket #2237 – Description

    initial v1  
    1212Calculate from tmask: 
    1313 
     14{{{#!f 
    1415               fmask(ji,jj,jk) = tmask(ji,jj  ,jk) * tmask(ji+1,jj  ,jk)   & 
    1516                  &            * tmask(ji,jj+1,jk) * tmask(ji+1,jj+1,jk) 
     17}}} 
    1618 
    1719Apply slip condition: 
    1820 
     21{{{#!f 
    1922      IF( rn_shlat /= 0 ) THEN      ! Not free-slip lateral boundary condition 
    2023         ! 
    … …  
    3134               END DO 
    3235            END DO 
     36}}} 
    3337 
    3438This selects fmask points at the topography-sea boundary by looking for points where fmask is zero but has a neighbouring fmask which is one. But this criterion doesn't work for single-point channels. A better criterion would be to check the neighbouring umask and vmask values: 
    3539 
    36 Instead of:   MAX( zwf(ji+1,jj), zwf(ji,jj+1),zwf(ji-1,jj), zwf(ji,jj-1) ) 
     40Instead of:    
     41{{{#!f 
     42MAX( zwf(ji+1,jj), zwf(ji,jj+1),zwf(ji-1,jj), zwf(ji,jj-1) ) 
     43}}} 
    3744 
    38 this would be better:   MAX( umask(ji,jj,jk), umask(ji,jj+1,jk), vmask(ji,jj,jk), vmask(ji+1,jj,jk) ) 
     45this would be better: 
     46{{{#!f 
     47MAX( umask(ji,jj,jk), umask(ji,jj+1,jk), vmask(ji,jj,jk), vmask(ji+1,jj,jk) ) 
     48}}} 
    3949 
    4050Maybe this was deliberate? But if you want no-slip with free-slip in channels it is probably better to do that explicitly via a spatially varying shlat parameter.