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 13492 – NEMO

Changeset 13492


Ignore:
Timestamp:
2020-09-18T17:27:27+02:00 (4 years ago)
Author:
smasson
Message:

r4.0-HEAD: avoid to use undefined values in GLS. pass all sette etsts in debug mode

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/releases/r4.0/r4.0-HEAD/src/OCE/ZDF/zdfgls.F90

    r13284 r13492  
    445445         END DO 
    446446      END DO 
    447       DO jk = 2, jpk                               ! Second recurrence : Lk = RHSk - Lk / Dk-1 * Lk-1 
     447      DO jk = 2, jpkm1                             ! Second recurrence : Lk = RHSk - Lk / Dk-1 * Lk-1 
    448448         DO jj = 2, jpjm1 
    449449            DO ji = fs_2, fs_jpim1    ! vector opt. 
     
    452452         END DO 
    453453      END DO 
    454       DO jk = jpk-1, 2, -1                         ! thrid recurrence : Ek = ( Lk - Uk * Ek+1 ) / Dk 
     454      DO jk = jpkm1, 2, -1                         ! thrid recurrence : Ek = ( Lk - Uk * Ek+1 ) / Dk 
    455455         DO jj = 2, jpjm1 
    456456            DO ji = fs_2, fs_jpim1    ! vector opt. 
     
    673673         END DO 
    674674      END DO 
    675       DO jk = 2, jpk                               ! Second recurrence : Lk = RHSk - Lk / Dk-1 * Lk-1 
     675      DO jk = 2, jpkm1                             ! Second recurrence : Lk = RHSk - Lk / Dk-1 * Lk-1 
    676676         DO jj = 2, jpjm1 
    677677            DO ji = fs_2, fs_jpim1    ! vector opt. 
     
    680680         END DO 
    681681      END DO 
    682       DO jk = jpk-1, 2, -1                         ! Third recurrence : Ek = ( Lk - Uk * Ek+1 ) / Dk 
     682      DO jk = jpkm1, 2, -1                         ! Third recurrence : Ek = ( Lk - Uk * Ek+1 ) / Dk 
    683683         DO jj = 2, jpjm1 
    684684            DO ji = fs_2, fs_jpim1    ! vector opt. 
Note: See TracChangeset for help on using the changeset viewer.