Changeset 1129 for trunk/NEMO/OPA_SRC/DYN/dynadv_ubs.F90
- Timestamp:
- 2008-06-24T08:45:37+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DYN/dynadv_ubs.F90
r719 r1129 17 17 USE dynspg_oce ! surface pressure gradient 18 18 USE in_out_manager ! I/O manager 19 USE dynspg_rl ! I/O manager19 USE dynspg_rl ! surface pressure gradient 20 20 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 21 USE trdmod ! ocean dynamics trends 22 USE trdmod_oce ! ocean variables trends 23 USE prtctl ! Print control 21 24 22 25 IMPLICIT NONE … … 97 100 zlv_vu(:,:,:,:) = 0.e0 98 101 99 102 IF( l_trddyn ) THEN ! Save ua and va trends 103 zfu_uw(:,:,:) = ua(:,:,:) 104 zfv_vw(:,:,:) = va(:,:,:) 105 ENDIF 106 100 107 ! ! =============== 101 108 DO jk = 1, jpkm1 ! Horizontal slab … … 191 198 ua(ji,jj,jk) = ua(ji,jj,jk) + zua 192 199 va(ji,jj,jk) = va(ji,jj,jk) + zva 193 #if defined key_trddyn194 utrd(ji,jj,jk,1) = zua ! save the horizontal advective trend of momentum195 vtrd(ji,jj,jk,1) = zva196 #endif197 200 END DO 198 201 END DO … … 201 204 ! ! =============== 202 205 206 IF( l_trddyn ) THEN ! save the horizontal advection trend for diagnostic 207 zfu_uw(:,:,:) = ua(:,:,:) - zfu_uw(:,:,:) 208 zfv_vw(:,:,:) = va(:,:,:) - zfv_vw(:,:,:) 209 CALL trd_mod( zfu_uw, zfv_vw, jpdyn_trd_had, 'DYN', kt ) 210 ENDIF 203 211 204 212 ! II. Vertical advection 205 213 ! ---------------------- 214 215 IF( l_trddyn ) THEN ! Save ua and va trends 216 zfu_t(:,:,:) = ua(:,:,:) 217 zfv_t(:,:,:) = va(:,:,:) 218 ENDIF 206 219 207 220 ! Second order centered tracer flux at w-point … … 251 264 END DO 252 265 266 IF( l_trddyn ) THEN ! save the vertical advection trend for diagnostic 267 zfu_t(:,:,:) = ua(:,:,:) - zfu_t(:,:,:) 268 zfv_t(:,:,:) = va(:,:,:) - zfv_t(:,:,:) 269 CALL trd_mod( zfu_t, zfv_t, jpdyn_trd_zad, 'DYN', kt ) 270 ENDIF 271 272 ! ! Control print 273 IF(ln_ctl) CALL prt_ctl( tab3d_1=ua, clinfo1=' ubs2 adv - Ua: ', mask1=umask, & 274 & tab3d_2=va, clinfo2= ' Va: ', mask2=vmask, clinfo3='dyn' ) 275 253 276 END SUBROUTINE dyn_adv_ubs 254 277
Note: See TracChangeset
for help on using the changeset viewer.