Changeset 13008

Timestamp:

2020-06-03T08:30:25+02:00 (3 years ago)

Author:

rrenshaw

Message:

revert Sarah's boundary conditions change for v5 reanalysis

File:

• branches/UKMO/AMM15_v3_6_STABLE_package_reanalysis4/NEMOGCM/NEMO/TOP_SRC/trcbdy.F90 (modified) (2 diffs)

branches/UKMO/AMM15_v3_6_STABLE_package_reanalysis4/NEMOGCM/NEMO/TOP_SRC/trcbdy.F90

@@ -9 +9 @@
    !!            3.5  !  2012     (S. Mocavero, I. Epicoco) Optimization of BDY communications
    !!            3.6  !  2015     (T. Lovato) Adapt BDY for tracers in TOP component
+! RJR: Sarah's boundary changes reverted here to see if they fix the BGC blow-ups
    !!----------------------------------------------------------------------
 #if defined key_bdy && key_top
@@ -109 +110 @@
             ii = idx%nbi(ib,igrd)
             ij = idx%nbj(ib,igrd)
-            zwgt = idx%nbw(ib,igrd) * rdttrc(ik) / 86400.d0  ! damping with a timescale of day
+! RJR: Sarah's boundary changes reverted here to see if they fix the BGC blow-ups
+!           zwgt = idx%nbw(ib,igrd) * rdttrc(ik) / 86400.d0  ! damping with a timescale of day
+            zwgt = idx%nbw(ib,igrd)
             tra(ii,ij,ik,jn) = ( tra(ii,ij,ik,jn) + zwgt * ( ( dta%trc(ib,ik) * dta%rn_fac)  & 
                         &  - tra(ii,ij,ik,jn) ) ) * tmask(ii,ij,ik)
          END DO
       END DO 
-! tra is overwritten at the boundary so damping doesn't work here - need neumann in addition
-! to duplicate the internal value at the boundary
-      CALL bdy_trc_nmn( jn,  idx, dta, kt )
+! ! tra is overwritten at the boundary so damping doesn't work here - need neumann in addition
+! ! to duplicate the internal value at the boundary
+!       CALL bdy_trc_nmn( jn,  idx, dta, kt )