- Timestamp:
- 2010-10-12T20:49:32+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/DYN/dynnxt.F90
r2148 r2236 15 15 !! 2.3 ! 2007-07 (D. Storkey) Calls to BDY routines. 16 16 !! 3.2 ! 2009-06 (G. Madec, R.Benshila) re-introduce the vvl option 17 !! 3.3 ! 2010-09 (D. Storkey, E.O'Dea) Bug fix for BDY module 17 18 !!------------------------------------------------------------------------- 18 19 … … 34 35 USE bdydta ! unstructured open boundary conditions 35 36 USE bdydyn ! unstructured open boundary conditions 36 USE agrif_opa_update37 USE agrif_opa_interp38 37 USE in_out_manager ! I/O manager 39 38 USE lbclnk ! lateral boundary condition (or mpp link) 40 39 USE prtctl ! Print control 40 #if defined key_agrif 41 USE agrif_opa_update 42 USE agrif_opa_interp 43 #endif 41 44 42 45 IMPLICIT NONE … … 50 53 !! NEMO/OPA 3.2 , LOCEAN-IPSL (2009) 51 54 !! $Id$ 52 !! Software is governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)55 !! Software is governed by the CeCILL licence (NEMOGCM/License_CeCILL.txt) 53 56 !!------------------------------------------------------------------------- 54 57 … … 171 174 # elif defined key_bdy 172 175 ! !* BDY open boundaries 173 !RB all this part should be in a specific routine174 176 IF( lk_dynspg_exp .OR. lk_dynspg_ts ) THEN ! except for filtered option 175 !176 177 CALL bdy_dyn_frs( kt ) 177 !178 IF( ln_bdy_dyn_fla ) THEN179 ua_e(:,:) = 0.e0180 va_e(:,:) = 0.e0181 ! Set these variables for use in bdy_dyn_fla182 hur_e(:,:) = hur(:,:)183 hvr_e(:,:) = hvr(:,:)184 DO jk = 1, jpkm1 !! Vertically integrated momentum trends185 ua_e(:,:) = ua_e(:,:) + fse3u(:,:,jk) * umask(:,:,jk) * ua(:,:,jk)186 va_e(:,:) = va_e(:,:) + fse3v(:,:,jk) * vmask(:,:,jk) * va(:,:,jk)187 END DO188 ua_e(:,:) = ua_e(:,:) * hur(:,:)189 va_e(:,:) = va_e(:,:) * hvr(:,:)190 DO jk = 1 , jpkm1191 ua(:,:,jk) = ua(:,:,jk) - ua_e(:,:)192 va(:,:,jk) = va(:,:,jk) - va_e(:,:)193 END DO194 CALL bdy_dta_bt( kt+1, 0)195 CALL bdy_dyn_fla( sshn_b )196 CALL lbc_lnk( ua_e, 'U', -1. ) ! Boundary points should be updated197 CALL lbc_lnk( va_e, 'V', -1. ) !198 DO jk = 1 , jpkm1199 ua(:,:,jk) = ( ua(:,:,jk) + ua_e(:,:) ) * umask(:,:,jk)200 va(:,:,jk) = ( va(:,:,jk) + va_e(:,:) ) * vmask(:,:,jk)201 END DO202 ENDIF203 !204 178 ENDIF 205 179 # endif
Note: See TracChangeset
for help on using the changeset viewer.