- Timestamp:
- 2019-11-27T12:51:32+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/AMM15_v3_6_STABLE_package_collate_utils305/NEMOGCM/NEMO/TOP_SRC/trcbdy.F90
r10270 r11978 109 109 ii = idx%nbi(ib,igrd) 110 110 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 112 112 tra(ii,ij,ik,jn) = ( tra(ii,ij,ik,jn) + zwgt * ( ( dta%trc(ib,ik) * dta%rn_fac) & 113 113 & - tra(ii,ij,ik,jn) ) ) * tmask(ii,ij,ik) 114 114 END DO 115 115 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 116 121 ! 117 122 IF( kt .eq. nit000 ) CLOSE( unit = 102 )
Note: See TracChangeset
for help on using the changeset viewer.