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

Changeset 10646


Ignore:
Timestamp:
2019-02-07T12:58:00+01:00 (5 years ago)
Author:
smueller
Message:

Enabling of the activation of the inverse barometer effect at open boundaries in configurations where both tidal forcing and the time-split surface pressure gradient computation have been selected; fix for the bug described in ticket #2223.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/releases/release-4.0/src/OCE/BDY/bdydta.F90

    r10529 r10646  
    371371      END DO  ! jbdy 
    372372 
     373      IF ( ln_apr_obc ) THEN 
     374         DO jbdy = 1, nb_bdy 
     375            IF (cn_tra(jbdy) /= 'runoff')THEN 
     376               igrd = 1                      ! meridional velocity 
     377               DO ib = 1, idx_bdy(jbdy)%nblenrim(igrd) 
     378                  ii   = idx_bdy(jbdy)%nbi(ib,igrd) 
     379                  ij   = idx_bdy(jbdy)%nbj(ib,igrd) 
     380                  dta_bdy(jbdy)%ssh(ib) = dta_bdy(jbdy)%ssh(ib) + ssh_ib(ii,ij) 
     381               END DO 
     382            ENDIF 
     383         END DO 
     384      ENDIF 
     385 
    373386      IF ( ln_tide ) THEN 
    374387         IF (ln_dynspg_ts) THEN      ! Fill temporary arrays with slow-varying bdy data                            
     
    389402      ENDIF 
    390403 
    391       IF ( ln_apr_obc ) THEN 
    392          DO jbdy = 1, nb_bdy 
    393             IF (cn_tra(jbdy) /= 'runoff')THEN 
    394                igrd = 1                      ! meridional velocity 
    395                DO ib = 1, idx_bdy(jbdy)%nblenrim(igrd) 
    396                   ii   = idx_bdy(jbdy)%nbi(ib,igrd) 
    397                   ij   = idx_bdy(jbdy)%nbj(ib,igrd) 
    398                   dta_bdy(jbdy)%ssh(ib) = dta_bdy(jbdy)%ssh(ib) + ssh_ib(ii,ij) 
    399                END DO 
    400             ENDIF 
    401          END DO 
    402       ENDIF 
    403404      ! 
    404405      IF( ln_timing )   CALL timing_stop('bdy_dta') 
Note: See TracChangeset for help on using the changeset viewer.