Changeset 4389


Ignore:
Timestamp:
2014-01-31T15:55:04+01:00 (7 years ago)
Author:
hliu
Message:

add a missed loop in wadlmt.F90

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_r4050_NOC_WaD/NEMOGCM/NEMO/OPA_SRC/DYN/wadlmt.F90

    r4385 r4389  
    5252      REAL(wp) ::   zzflxp, zzflxn      ! local scalars 
    5353      REAL(wp) ::   ztmp                ! local scalars 
    54       INTEGER,  POINTER,  DIMENSION(:,:) ::   iwdflg                   !: W/D cell flag 
     54      INTEGER,  POINTER,  DIMENSION(:,:) ::   iwdflg                   !: W/D cell flag (should become global) 
    5555      REAL(wp), POINTER,  DIMENSION(:,:) ::   zwdlmtu, zwdlmtv         !: W/D flux limiters 
    5656      REAL(wp), POINTER,  DIMENSION(:,:) ::   zflxp,  zflxn            ! specific 2D workspace 
     
    158158           
    159159           IF(jk1 == nn_waditr) THEN 
    160                
    161               IF(zflxu1(ji,  jj  ) >= 0._wp) zwdlmtu(ji,  jj  ) = 0._wp 
    162               IF(zflxu1(ji-1,jj  ) <  0._wp) zwdlmtu(ji-1,jj  ) = 0._wp 
    163               IF(zflxv1(ji,  jj  ) >= 0._wp) zwdlmtv(ji,  jj  ) = 0._wp 
    164               IF(zflxv1(ji,  jj-1) <  0._wp) zwdlmtv(ji,  jj-1) = 0._wp 
    165  
    166               CALL lbc_lnk( zwdlmtu, 'U', 1. ) 
    167               CALL lbc_lnk( zwdlmtv, 'V', 1. ) 
     160              DO jj = 2, jpjm1 
     161                 DO ji = fs_2, fs_jpim1   ! vector opt. 
     162                  
     163                 IF(zflxu1(ji,  jj  ) >= 0._wp) zwdlmtu(ji,  jj  ) = 0._wp 
     164                 IF(zflxu1(ji-1,jj  ) <  0._wp) zwdlmtu(ji-1,jj  ) = 0._wp 
     165                 IF(zflxv1(ji,  jj  ) >= 0._wp) zwdlmtv(ji,  jj  ) = 0._wp 
     166                 IF(zflxv1(ji,  jj-1) <  0._wp) zwdlmtv(ji,  jj-1) = 0._wp 
     167             
     168                 CALL lbc_lnk( zwdlmtu, 'U', 1. ) 
     169                 CALL lbc_lnk( zwdlmtv, 'V', 1. ) 
     170                 END DO ! ji loop 
     171              END DO  ! jj loop 
    168172           END IF 
    169173 
Note: See TracChangeset for help on using the changeset viewer.