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 13008 for branches/UKMO – NEMO

Changeset 13008 for branches/UKMO


Ignore:
Timestamp:
2020-06-03T08:30:25+02:00 (4 years ago)
Author:
rrenshaw
Message:

revert Sarah's boundary conditions change for v5 reanalysis

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/AMM15_v3_6_STABLE_package_reanalysis4/NEMOGCM/NEMO/TOP_SRC/trcbdy.F90

    r11639 r13008  
    99   !!            3.5  !  2012     (S. Mocavero, I. Epicoco) Optimization of BDY communications 
    1010   !!            3.6  !  2015     (T. Lovato) Adapt BDY for tracers in TOP component 
     11! RJR: Sarah's boundary changes reverted here to see if they fix the BGC blow-ups 
    1112   !!---------------------------------------------------------------------- 
    1213#if defined key_bdy && key_top 
     
    109110            ii = idx%nbi(ib,igrd) 
    110111            ij = idx%nbj(ib,igrd) 
    111             zwgt = idx%nbw(ib,igrd) * rdttrc(ik) / 86400.d0  ! damping with a timescale of day 
     112! RJR: Sarah's boundary changes reverted here to see if they fix the BGC blow-ups 
     113!           zwgt = idx%nbw(ib,igrd) * rdttrc(ik) / 86400.d0  ! damping with a timescale of day 
     114            zwgt = idx%nbw(ib,igrd) 
    112115            tra(ii,ij,ik,jn) = ( tra(ii,ij,ik,jn) + zwgt * ( ( dta%trc(ib,ik) * dta%rn_fac)  &  
    113116                        &  - tra(ii,ij,ik,jn) ) ) * tmask(ii,ij,ik) 
    114117         END DO 
    115118      END DO  
    116 ! tra is overwritten at the boundary so damping doesn't work here - need neumann in addition 
    117 ! to duplicate the internal value at the boundary 
    118       CALL bdy_trc_nmn( jn,  idx, dta, kt ) 
     119! ! tra is overwritten at the boundary so damping doesn't work here - need neumann in addition 
     120! ! to duplicate the internal value at the boundary 
     121!       CALL bdy_trc_nmn( jn,  idx, dta, kt ) 
    119122 
    120123 
Note: See TracChangeset for help on using the changeset viewer.