Changeset 11978


Ignore:
Timestamp:
2019-11-27T12:51:32+01:00 (2 months ago)
Author:
dford
Message:

Update BDY scheme to remove dead zone.

File:
1 edited

Legend:

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

    r10270 r11978  
    109109            ii = idx%nbi(ib,igrd) 
    110110            ij = idx%nbj(ib,igrd) 
    111             zwgt = idx%nbw(ib,igrd) 
     111            zwgt = idx%nbw(ib,igrd) * rdttrc(ik) / 86400.d0  ! damping with a timescale of day 
    112112            tra(ii,ij,ik,jn) = ( tra(ii,ij,ik,jn) + zwgt * ( ( dta%trc(ib,ik) * dta%rn_fac)  &  
    113113                        &  - tra(ii,ij,ik,jn) ) ) * tmask(ii,ij,ik) 
    114114         END DO 
    115115      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 
     120 
    116121      ! 
    117122      IF( kt .eq. nit000 ) CLOSE( unit = 102 ) 
Note: See TracChangeset for help on using the changeset viewer.