Changeset 5061 for trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv_ubs.F90
- Timestamp:
- 2015-02-05T15:14:40+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv_ubs.F90
r4990 r5061 16 16 USE oce ! ocean dynamics and tracers 17 17 USE dom_oce ! ocean space and time domain 18 USE trd_oce ! trends: ocean variables 19 USE trddyn ! trend manager: dynamics 20 ! 18 USE trdmod ! ocean dynamics trends 19 USE trdmod_oce ! ocean variables trends 21 20 USE in_out_manager ! I/O manager 22 21 USE prtctl ! Print control 23 22 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 24 23 USE lib_mpp ! MPP library 25 USE wrk_nemo ! Memory Allocation26 USE timing ! Timing24 USE wrk_nemo ! Memory Allocation 25 USE timing ! Timing 27 26 28 27 IMPLICIT NONE … … 116 115 DO jj = 2, jpjm1 ! laplacian 117 116 DO ji = fs_2, fs_jpim1 ! vector opt. 117 ! 118 118 zlu_uu(ji,jj,jk,1) = ( ub (ji+1,jj,jk)-2.*ub (ji,jj,jk)+ub (ji-1,jj,jk) ) * umask(ji,jj,jk) 119 zlv_vv(ji,jj,jk,1) = ( vb (ji,jj+1,jk)-2.*vb (ji,jj,jk)+vb (ji,jj-1,jk) ) * vmask(ji,jj,jk) 120 zlu_uv(ji,jj,jk,1) = ( ub (ji,jj+1,jk)- 2.*ub (ji,jj,jk)+ub (ji,jj-1,jk) ) * umask(ji,jj,jk)121 zlv_vu(ji,jj,jk,1) = ( vb (ji+1,jj,jk)- 2.*vb (ji,jj,jk)+vb (ji-1,jj,jk) ) * vmask(ji,jj,jk)119 zlv_vv(ji,jj,jk,1) = ( vb (ji,jj+1,jk)-2.*vb (ji,jj,jk)+vb (ji,jj-1,jk) ) * vmask(ji,jj,jk) 120 zlu_uv(ji,jj,jk,1) = ( ub (ji,jj+1,jk)-ub (ji,jj,jk) ) * fmask(ji,jj,jk)-( ub (ji,jj,jk)-ub (ji,jj-1,jk) ) * fmask(ji,jj-1,jk) 121 zlv_vu(ji,jj,jk,1) = ( vb (ji+1,jj,jk)-vb (ji,jj,jk) ) * fmask(ji,jj,jk)-( vb (ji,jj,jk)-vb (ji-1,jj,jk) ) * fmask(ji-1,jj,jk) 122 122 ! 123 123 zlu_uu(ji,jj,jk,2) = ( zfu(ji+1,jj,jk)-2.*zfu(ji,jj,jk)+zfu(ji-1,jj,jk) ) * umask(ji,jj,jk) 124 124 zlv_vv(ji,jj,jk,2) = ( zfv(ji,jj+1,jk)-2.*zfv(ji,jj,jk)+zfv(ji,jj-1,jk) ) * vmask(ji,jj,jk) 125 zlu_uv(ji,jj,jk,2) = ( zfu (ji,jj+1,jk)-2.*zfu(ji,jj,jk)+zfu(ji,jj-1,jk) ) * umask(ji,jj,jk)126 zlv_vu(ji,jj,jk,2) = ( zfv (ji+1,jj,jk)-2.*zfv(ji,jj,jk)+zfv(ji-1,jj,jk) ) * vmask(ji,jj,jk)125 zlu_uv(ji,jj,jk,2) = ( zfu (ji,jj+1,jk)-zfu (ji,jj,jk) ) * fmask(ji,jj,jk)-( zfu (ji,jj,jk)-zfu (ji,jj-1,jk) ) * fmask(ji,jj-1,jk) 126 zlv_vu(ji,jj,jk,2) = ( zfv (ji+1,jj,jk)-zfv (ji,jj,jk) ) * fmask(ji,jj,jk)-( zfv (ji,jj,jk)-zfv (ji-1,jj,jk) ) * fmask(ji-1,jj,jk) 127 127 END DO 128 128 END DO … … 197 197 zfu_uw(:,:,:) = ua(:,:,:) - zfu_uw(:,:,:) 198 198 zfv_vw(:,:,:) = va(:,:,:) - zfv_vw(:,:,:) 199 CALL trd_ dyn( zfu_uw, zfv_vw, jpdyn_keg, kt )199 CALL trd_mod( zfu_uw, zfv_vw, jpdyn_trd_had, 'DYN', kt ) 200 200 zfu_t(:,:,:) = ua(:,:,:) 201 201 zfv_t(:,:,:) = va(:,:,:) … … 246 246 zfu_t(:,:,:) = ua(:,:,:) - zfu_t(:,:,:) 247 247 zfv_t(:,:,:) = va(:,:,:) - zfv_t(:,:,:) 248 CALL trd_ dyn( zfu_t, zfv_t, jpdyn_zad, kt )248 CALL trd_mod( zfu_t, zfv_t, jpdyn_trd_zad, 'DYN', kt ) 249 249 ENDIF 250 250 ! ! Control print
Note: See TracChangeset
for help on using the changeset viewer.