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 2374 – NEMO

Changeset 2374


Ignore:
Timestamp:
2010-11-10T22:45:27+01:00 (13 years ago)
Author:
rblod
Message:

Correct fix for obcfla on the trunk

File:
1 edited

Legend:

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

    r2270 r2374  
    8686               &            * ( ( sshn_e(ji,jj) + sshn_e(ji+1,jj) ) * 0.5        & 
    8787               &            - sshfoe(jj) )  ) * uemsk(jj,1) 
    88          END DO 
    89          DO jj = 1, jpj 
    9088            sshfoe_b(ji,jj) = sshfoe_b(ji,jj) + sqrt( grav*hur(ji,jj) )     & 
    9189               &             * ( ( sshn_e(ji,jj) + sshn_e(ji+1,jj) ) * 0.5  & 
    92                &                 - sshfoe(jj) ) * uemsk(jj,1) 
     90               &            - sshfoe(jj) )  ) * uemsk(jj,1) 
     91         END DO 
     92      END DO 
     93      DO ji = nie0p1, nie1p1 
     94         DO jj = 1, jpj 
    9395            ssha_e(ji,jj) = ssha_e(ji,jj) * ( 1. - temsk(jj,1) ) & 
    9496               &            + temsk(jj,1) * sshfoe(jj) 
     
    153155               &            * ( ( sshn_e(ji,jj) + sshn_e(ji,jj+1) ) * 0.5         & 
    154156               &                - sshfon(ji) ) ) * vnmsk(ji,1) 
    155          END DO 
    156          DO ji = 1, jpi 
    157157            sshfon_b(ji,jj) = sshfon_b(ji,jj) + sqrt( grav * hvr(ji,jj) )  & 
    158158               &              * ( ( sshn_e(ji,jj) + sshn_e(ji,jj+1) ) * 0.5    & 
    159159               &                  - sshfon(ji) ) * vnmsk(ji,1) 
     160         END DO 
     161      END DO 
     162      DO jj = njn0p1, njn1p1 
     163         DO ji = 1, jpi 
    160164            ssha_e(ji,jj) = ssha_e(ji,jj) * ( 1. - tnmsk(ji,1) ) & 
    161165               &            + sshfon(ji) * tnmsk(ji,1) 
Note: See TracChangeset for help on using the changeset viewer.