Changeset 5324


Ignore:
Timestamp:
2015-06-01T10:54:30+02:00 (5 years ago)
Author:
pabouttier
Message:

Fix a mistake in jpkm1 loop in dynkeg.F90

File:
1 edited

Legend:

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

    r5323 r5324  
    9898      zhke(:,:) = 0._wp 
    9999 
    100       SELECT CASE ( kscheme )             !== Horizontal kinetic energy at T-point  ==! 
    101       ! 
    102       CASE ( nkeg_C2 )                          !--  Standard scheme  --! 
    103          DO jk = 1, jpkm1 
     100      DO jk = 1, jpkm1 
     101 
     102         SELECT CASE ( kscheme )             !== Horizontal kinetic energy at T-point  ==! 
     103         ! 
     104         CASE ( nkeg_C2 )                          !--  Standard scheme  --! 
     105          
    104106            DO jj = 2, jpj 
    105107               DO ji = fs_2, jpi   ! vector opt. 
     
    111113               END DO 
    112114            END DO 
    113          END DO 
    114115         ! 
    115       CASE ( nkeg_HW )                          !--  Hollingsworth scheme  --! 
    116          DO jk = 1, jpkm1 
     116         CASE ( nkeg_HW )                          !--  Hollingsworth scheme  --! 
     117 
    117118            DO jj = 2, jpjm1 
    118119               DO ji = fs_2, jpim1   ! vector opt. 
     
    129130               END DO 
    130131            END DO 
    131          END DO 
    132          CALL lbc_lnk( zhke, 'T', 1. ) 
     132            CALL lbc_lnk( zhke, 'T', 1. ) 
     133            ! 
     134         END SELECT 
    133135         ! 
    134       END SELECT 
    135       ! 
    136       DO jk = 1, jpkm1                    !==  grad( KE ) added to the general momentum trends  ==! 
    137136         DO jj = 2, jpjm1 
    138137            DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    141140            END DO 
    142141         END DO 
     142 
    143143      END DO 
    144144      ! 
Note: See TracChangeset for help on using the changeset viewer.