Changeset 2723


Ignore:
Timestamp:
2011-04-06T12:21:48+02:00 (10 years ago)
Author:
rblod
Message:

Correction in dynnxt for bdy, see ticket #811

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynnxt.F90

    r2715 r2723  
    1616   !!            3.2  !  2009-06  (G. Madec, R.Benshila)  re-introduce the vvl option 
    1717   !!            3.3  !  2010-09  (D. Storkey, E.O'Dea) Bug fix for BDY module 
     18   !!            3.3  !  2011-03  (P. Oddo) Bug fix for time-splitting+(BDY-OBC) and not VVL 
    1819   !!------------------------------------------------------------------------- 
    1920   
     
    175176# elif defined key_bdy  
    176177      !                                !* BDY open boundaries 
    177       IF( .NOT. lk_dynspg_flt )   CALL bdy_dyn_frs( kt ) 
     178      IF( .NOT. lk_dynspg_flt ) THEN 
     179         CALL bdy_dyn_frs( kt ) 
     180#  if ! defined key_vvl 
     181         ua_e(:,:) = 0.e0 
     182         va_e(:,:) = 0.e0 
     183         ! Set these variables for use in bdy_dyn_fla 
     184         hur_e(:,:) = hur(:,:) 
     185         hvr_e(:,:) = hvr(:,:) 
     186         DO jk = 1, jpkm1   !! Vertically integrated momentum trends 
     187            ua_e(:,:) = ua_e(:,:) + fse3u(:,:,jk) * umask(:,:,jk) * ua(:,:,jk) 
     188            va_e(:,:) = va_e(:,:) + fse3v(:,:,jk) * vmask(:,:,jk) * va(:,:,jk) 
     189         END DO 
     190         ua_e(:,:) = ua_e(:,:) * hur(:,:) 
     191         va_e(:,:) = va_e(:,:) * hvr(:,:) 
     192         DO jk = 1 , jpkm1 
     193            ua(:,:,jk) = ua(:,:,jk) - ua_e(:,:) 
     194            va(:,:,jk) = va(:,:,jk) - va_e(:,:) 
     195         END DO 
     196         CALL bdy_dta_fla( kt+1, 0,2*nn_baro) 
     197         CALL bdy_dyn_fla( sshn_b ) 
     198         CALL lbc_lnk( ua_e, 'U', -1. )   ! Boundary points should be updated 
     199         CALL lbc_lnk( va_e, 'V', -1. )   ! 
     200         DO jk = 1 , jpkm1 
     201            ua(:,:,jk) = ( ua(:,:,jk) + ua_e(:,:) ) * umask(:,:,jk) 
     202            va(:,:,jk) = ( va(:,:,jk) + va_e(:,:) ) * vmask(:,:,jk) 
     203         END DO 
     204#  endif 
     205      ENDIF 
    178206# endif 
    179207      ! 
Note: See TracChangeset for help on using the changeset viewer.