- Timestamp:
- 2015-11-11T19:11:01+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_bdy_sponge_temp/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_iso.F90
r5876 r5878 32 32 USE wrk_nemo ! Memory Allocation 33 33 USE timing ! Timing 34 #if defined key_bdy 35 USE bdy_oce 36 #endif 34 37 35 38 IMPLICIT NONE … … 110 113 REAL(wp), POINTER, DIMENSION(:,: ) :: z2d 111 114 REAL(wp), POINTER, DIMENSION(:,:,:) :: zdkt, zdk1t, zdit, zdjt, ztfw 115 REAL(wp), DIMENSION(jpi,jpj) :: zfactor ! multiplier for diffusion 112 116 !!---------------------------------------------------------------------- 113 117 ! … … 124 128 ENDIF 125 129 ! 130 #if defined key_bdy 131 zfactor(:,:) = sponge_factor(:,:) 132 #else 133 zfactor(:,:) = 1.0 134 #endif 126 135 ! ! =========== 127 136 DO jn = 1, kjpt ! tracer loop … … 200 209 DO jj = 1 , jpjm1 201 210 DO ji = 1, fs_jpim1 ! vector opt. 202 zabe1 = ( fsahtu(ji,jj,jk) + pahtb0 ) * re2u_e1u(ji,jj) * fse3u_n(ji,jj,jk)203 zabe2 = ( fsahtv(ji,jj,jk) + pahtb0 ) * re1v_e2v(ji,jj) * fse3v_n(ji,jj,jk)211 zabe1 = zfactor(ji,jj) * ( fsahtu(ji,jj,jk) + pahtb0 ) * re2u_e1u(ji,jj) * fse3u_n(ji,jj,jk) 212 zabe2 = zfactor(ji,jj) * ( fsahtv(ji,jj,jk) + pahtb0 ) * re1v_e2v(ji,jj) * fse3v_n(ji,jj,jk) 204 213 ! 205 214 zmsku = 1. / MAX( tmask(ji+1,jj,jk ) + tmask(ji,jj,jk+1) &
Note: See TracChangeset
for help on using the changeset viewer.