Changeset 14978


Ignore:
Timestamp:
2021-06-11T15:21:08+02:00 (6 months ago)
Author:
hadcv
Message:

#2640: correct tra_adv_qck results on the north fold

Location:
NEMO/trunk/src/OCE/TRA
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/TRA/traadv_qck.F90

    r14834 r14978  
    246246         IF (nn_hls==1) CALL lbc_lnk( 'traadv_qck', zfc(:,:,:), 'T', 1.0_wp , zfd(:,:,:), 'T', 1.0_wp, ld4only= .TRUE. )   ! Lateral boundary conditions 
    247247 
     248         ! Correct zfd on northfold after lbc_lnk; see #2640 
     249         IF( nn_hls == 1 .AND. l_IdoNFold .AND. ntej == Nje0 ) THEN 
     250            DO jk = 1, jpkm1 
     251               WHERE( tmask_i(ntsi:ntei,ntej:jpj) == 0._wp ) zfd(ntsi:ntei,ntej:jpj,jk) = zfc(ntsi:ntei,ntej:jpj,jk) 
     252            END DO 
     253         ENDIF 
    248254         ! 
    249255         ! Horizontal advective fluxes 
  • NEMO/trunk/src/OCE/TRA/traadv_qck_lf.F90

    r14834 r14978  
    111111      ! 
    112112      !        ! horizontal fluxes are computed with the QUICKEST + ULTIMATE scheme 
    113       CALL tra_adv_qck_i( kt, cdtype, p2dt, pU, Kbb, Kmm, pt, kjpt, Krhs ) 
    114       CALL tra_adv_qck_j( kt, cdtype, p2dt, pV, Kbb, Kmm, pt, kjpt, Krhs ) 
     113      CALL tra_adv_qck_i_lf( kt, cdtype, p2dt, pU, Kbb, Kmm, pt, kjpt, Krhs ) 
     114      CALL tra_adv_qck_j_lf( kt, cdtype, p2dt, pV, Kbb, Kmm, pt, kjpt, Krhs ) 
    115115 
    116116      !        ! vertical fluxes are computed with the 2nd order centered scheme 
    117       CALL tra_adv_cen2_k( kt, cdtype, pW, Kmm, pt, kjpt, Krhs ) 
     117      CALL tra_adv_cen2_k_lf( kt, cdtype, pW, Kmm, pt, kjpt, Krhs ) 
    118118      ! 
    119119   END SUBROUTINE tra_adv_qck_lf 
    120120 
    121121 
    122    SUBROUTINE tra_adv_qck_i( kt, cdtype, p2dt, pU, Kbb, Kmm, pt, kjpt, Krhs ) 
     122   SUBROUTINE tra_adv_qck_i_lf( kt, cdtype, p2dt, pU, Kbb, Kmm, pt, kjpt, Krhs ) 
    123123      !!---------------------------------------------------------------------- 
    124124      !! 
     
    193193      END DO 
    194194      ! 
    195    END SUBROUTINE tra_adv_qck_i 
    196  
    197  
    198    SUBROUTINE tra_adv_qck_j( kt, cdtype, p2dt, pV, Kbb, Kmm, pt, kjpt, Krhs ) 
     195   END SUBROUTINE tra_adv_qck_i_lf 
     196 
     197 
     198   SUBROUTINE tra_adv_qck_j_lf( kt, cdtype, p2dt, pV, Kbb, Kmm, pt, kjpt, Krhs ) 
    199199      !!---------------------------------------------------------------------- 
    200200      !! 
     
    273273      END DO 
    274274      ! 
    275    END SUBROUTINE tra_adv_qck_j 
    276  
    277  
    278    SUBROUTINE tra_adv_cen2_k( kt, cdtype, pW, Kmm, pt, kjpt, Krhs ) 
     275   END SUBROUTINE tra_adv_qck_j_lf 
     276 
     277 
     278   SUBROUTINE tra_adv_cen2_k_lf( kt, cdtype, pW, Kmm, pt, kjpt, Krhs ) 
    279279      !!---------------------------------------------------------------------- 
    280280      !! 
     
    323323      END DO 
    324324      ! 
    325    END SUBROUTINE tra_adv_cen2_k 
     325   END SUBROUTINE tra_adv_cen2_k_lf 
    326326 
    327327 
Note: See TracChangeset for help on using the changeset viewer.