Changeset 9901


Ignore:
Timestamp:
2018-07-09T13:32:57+02:00 (2 years ago)
Author:
acc
Message:

Bug fix in traadv_fct.F90. See ticket #2110

File:
1 edited

Legend:

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

    r9598 r9901  
    548548!!gm   
    549549      ! 
     550      IF ( ln_isfcav ) THEN            ! set level two values which may not be set in ISF case 
     551         zwd(:,:,2) = 1._wp  ;  zwi(:,:,2) = 0._wp  ;  zws(:,:,2) = 0._wp  ;  zwrm(:,:,2) = 0._wp 
     552      END IF 
     553      ! 
    550554      DO jj = 2, jpjm1                 ! 2nd order centered at top & bottom 
    551555         DO ji = fs_2, fs_jpim1 
     
    556560            zwi (ji,jj,ikt) = 0._wp 
    557561            zws (ji,jj,ikt) = 0._wp 
    558             zwrm(ji,jj,ikt) = 0.5_wp * ( pt_in(ji,jj,jk-1) + pt_in(ji,jj,jk) ) 
     562            zwrm(ji,jj,ikt) = 0.5_wp * ( pt_in(ji,jj,ikt-1) + pt_in(ji,jj,ikt) ) 
    559563            ! 
    560564            zwd (ji,jj,ikb) = 1._wp          ! bottom 
    561565            zwi (ji,jj,ikb) = 0._wp 
    562566            zws (ji,jj,ikb) = 0._wp 
    563             zwrm(ji,jj,ikb) = 0.5_wp * ( pt_in(ji,jj,jk-1) + pt_in(ji,jj,jk) )             
     567            zwrm(ji,jj,ikb) = 0.5_wp * ( pt_in(ji,jj,ikb-1) + pt_in(ji,jj,ikb) )             
    564568         END DO 
    565569      END DO    
Note: See TracChangeset for help on using the changeset viewer.