New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 2236 for branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/DYN/dynnxt.F90 – NEMO

Ignore:
Timestamp:
2010-10-12T20:49:32+02:00 (13 years ago)
Author:
cetlod
Message:

First guess of NEMO_v3.3

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/DYN/dynnxt.F90

    r2148 r2236  
    1515   !!            2.3  !  2007-07  (D. Storkey) Calls to BDY routines.  
    1616   !!            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 
    1718   !!------------------------------------------------------------------------- 
    1819   
     
    3435   USE bdydta          ! unstructured open boundary conditions 
    3536   USE bdydyn          ! unstructured open boundary conditions 
    36    USE agrif_opa_update 
    37    USE agrif_opa_interp 
    3837   USE in_out_manager  ! I/O manager 
    3938   USE lbclnk          ! lateral boundary condition (or mpp link) 
    4039   USE prtctl          ! Print control 
     40#if defined key_agrif 
     41   USE agrif_opa_update 
     42   USE agrif_opa_interp 
     43#endif 
    4144 
    4245   IMPLICIT NONE 
     
    5053   !! NEMO/OPA 3.2 , LOCEAN-IPSL (2009)  
    5154   !! $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)  
    5356   !!------------------------------------------------------------------------- 
    5457 
     
    171174# elif defined key_bdy  
    172175      !                                !* BDY open boundaries 
    173       !RB all this part should be in a specific routine 
    174176      IF( lk_dynspg_exp .OR. lk_dynspg_ts ) THEN       ! except for filtered option 
    175          ! 
    176177         CALL bdy_dyn_frs( kt ) 
    177          ! 
    178          IF( ln_bdy_dyn_fla ) THEN 
    179             ua_e(:,:) = 0.e0 
    180             va_e(:,:) = 0.e0 
    181             ! Set these variables for use in bdy_dyn_fla 
    182             hur_e(:,:) = hur(:,:) 
    183             hvr_e(:,:) = hvr(:,:) 
    184             DO jk = 1, jpkm1   !! Vertically integrated momentum trends 
    185                ua_e(:,:) = ua_e(:,:) + fse3u(:,:,jk) * umask(:,:,jk) * ua(:,:,jk) 
    186                va_e(:,:) = va_e(:,:) + fse3v(:,:,jk) * vmask(:,:,jk) * va(:,:,jk) 
    187             END DO 
    188             ua_e(:,:) = ua_e(:,:) * hur(:,:) 
    189             va_e(:,:) = va_e(:,:) * hvr(:,:) 
    190             DO jk = 1 , jpkm1 
    191                ua(:,:,jk) = ua(:,:,jk) - ua_e(:,:) 
    192                va(:,:,jk) = va(:,:,jk) - va_e(:,:) 
    193             END DO 
    194             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 updated 
    197             CALL lbc_lnk( va_e, 'V', -1. )   ! 
    198             DO jk = 1 , jpkm1 
    199                ua(:,:,jk) = ( ua(:,:,jk) + ua_e(:,:) ) * umask(:,:,jk) 
    200                va(:,:,jk) = ( va(:,:,jk) + va_e(:,:) ) * vmask(:,:,jk) 
    201             END DO 
    202          ENDIF 
    203          ! 
    204178      ENDIF 
    205179# endif 
Note: See TracChangeset for help on using the changeset viewer.