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 2064 for trunk – NEMO

Changeset 2064 for trunk


Ignore:
Timestamp:
2010-08-27T13:22:40+02:00 (14 years ago)
Author:
rblod
Message:

Set correct value for uemsk and vemsk index, see ticket #708

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/OBC/obcfla.F90

    r1941 r2064  
    7878      !!------------------------------------------------------------------------------ 
    7979      !! * Local declaration 
    80       INTEGER ::   ji, jj, jk ! dummy loop indices 
     80      INTEGER ::   ji, jj ! dummy loop indices 
    8181      !!------------------------------------------------------------------------------ 
    8282 
     
    8585            ua_e(ji,jj) = (  ubtfoe(jj) * hur(ji,jj) + sqrt( grav*hur(ji,jj) )   & 
    8686               &            * ( ( sshn_e(ji,jj) + sshn_e(ji+1,jj) ) * 0.5        & 
    87                &            - sshfoe(jj) )  ) * uemsk(jj,jk) 
     87               &            - sshfoe(jj) )  ) * uemsk(jj,1) 
    8888         END DO 
    8989      END DO 
     
    9595            ssha_e(ji,jj) = ssha_e(ji,jj) * ( 1. - temsk(jj,1) ) & 
    9696               &            + temsk(jj,1) * sshfoe(jj) 
    97             va_e(ji,jj) = vbtfoe(jj) * hvr(ji,jj) * uemsk(jj,jk) 
     97            va_e(ji,jj) = vbtfoe(jj) * hvr(ji,jj) * uemsk(jj,1) 
    9898         END DO 
    9999      END DO 
     
    114114      !!------------------------------------------------------------------------------ 
    115115      !! * Local declaration 
    116       INTEGER ::   ji, jj, jk ! dummy loop indices 
     116      INTEGER ::   ji, jj ! dummy loop indices 
    117117      !!------------------------------------------------------------------------------ 
    118118 
     
    121121            ua_e(ji,jj) = ( ubtfow(jj) * hur(ji,jj) - sqrt( grav * hur(ji,jj) )   & 
    122122               &            * ( ( sshn_e(ji,jj) + sshn_e(ji+1,jj) ) * 0.5         & 
    123                &                - sshfow(jj) ) ) * uwmsk(jj,jk) 
    124             va_e(ji,jj) = vbtfow(jj) * hvr(ji,jj) * uwmsk(jj,jk) 
     123               &                - sshfow(jj) ) ) * uwmsk(jj,1) 
     124            va_e(ji,jj) = vbtfow(jj) * hvr(ji,jj) * uwmsk(jj,1) 
    125125         END DO 
    126126         DO jj = 1, jpj 
     
    147147      !!------------------------------------------------------------------------------ 
    148148      !! * Local declaration 
    149       INTEGER ::   ji, jj, jk ! dummy loop indices 
     149      INTEGER ::   ji, jj ! dummy loop indices 
    150150      !!------------------------------------------------------------------------------ 
    151151 
     
    154154            va_e(ji,jj) = ( vbtfon(ji) * hvr(ji,jj) + sqrt( grav * hvr(ji,jj) )   & 
    155155               &            * ( ( sshn_e(ji,jj) + sshn_e(ji,jj+1) ) * 0.5         & 
    156                &                - sshfon(ji) ) ) * vnmsk(ji,jk) 
     156               &                - sshfon(ji) ) ) * vnmsk(ji,1) 
    157157         END DO 
    158158      END DO 
     
    164164            ssha_e(ji,jj) = ssha_e(ji,jj) * ( 1. - tnmsk(ji,1) ) & 
    165165               &            + sshfon(ji) * tnmsk(ji,1) 
    166             ua_e(ji,jj) = ubtfon(ji) * hur(ji,jj) * vnmsk(ji,jk) 
     166            ua_e(ji,jj) = ubtfon(ji) * hur(ji,jj) * vnmsk(ji,1) 
    167167         END DO 
    168168      END DO 
     
    182182      !!------------------------------------------------------------------------------ 
    183183      !! * Local declaration 
    184       INTEGER ::   ji, jj, jk ! dummy loop indices 
     184      INTEGER ::   ji, jj ! dummy loop indices 
    185185 
    186186      !!------------------------------------------------------------------------------ 
     
    190190            va_e(ji,jj) = ( vbtfos(ji) * hvr(ji,jj) - sqrt( grav * hvr(ji,jj) )   & 
    191191               &            * ( ( sshn_e(ji,jj) + sshn_e(ji,jj+1) ) * 0.5         & 
    192                &                - sshfos(ji) ) ) * vsmsk(ji,jk) 
    193             ua_e(ji,jj) = ubtfos(ji) * hur(ji,jj) * vsmsk(ji,jk) 
     192               &                - sshfos(ji) ) ) * vsmsk(ji,1) 
     193            ua_e(ji,jj) = ubtfos(ji) * hur(ji,jj) * vsmsk(ji,1) 
    194194         END DO 
    195195         DO ji = 1, jpi 
Note: See TracChangeset for help on using the changeset viewer.