Changeset 12554


Ignore:
Timestamp:
2020-03-13T17:22:06+01:00 (8 months ago)
Author:
jcastill
Message:

Fixes to the velocity limiter and vertical interpolation

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/r12083_India_uncoupled/src/OCE/DYN/dynnxt.F90

    r12453 r12554  
    5858   PUBLIC    dyn_nxt   ! routine called by step.F90 
    5959 
     60   !! Substitution 
     61#  include "vectopt_loop_substitute.h90" 
    6062   !!---------------------------------------------------------------------- 
    6163   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
     
    406408  
    407409      DO jk = 1, jpkm1  
    408          DO jj = 1, jpjm1  
    409             DO ji = 1, jpim1  
     410         DO jj = 2, jpjm1  
     411            DO ji = fs_2, fs_jpim1   ! vect. opt. 
    410412               ! U direction  
    411413               zzu = un(ji,jj,jk)  
     
    425427               ! limit currents  
    426428               un(ji,jj,jk) = min ( zzu,zplim) * isp + max(zzu,zmlim) *ism  
     429 
    427430               ! V  direction  
    428431               zzu = vn(ji,jj,jk)  
     
    439442               zmlim = -zzcn * (e3t_n(ji,jj+1,jk) * e1t(ji,jj+1) * e2t(ji,jj+1)) / &  
    440443                                               (2.0*rdt * ze3e1)*vmask(ji,jj,jk)  
     444               ! limit currents  
    441445               vn(ji,jj,jk) = min ( zzu,zplim) * isp + max(zzu,zmlim) *ism  
    442446            ENDDO  
    443447         ENDDO  
    444448      ENDDO  
     449      CALL lbc_lnk_multi( 'dynnxt', un(:,:,:), 'U',  -1., vn(:,:,:), 'V',  -1. ) 
    445450  
    446451    END SUBROUTINE dyn_limit_velocity 
Note: See TracChangeset for help on using the changeset viewer.