Changeset 1413


Ignore:
Timestamp:
2009-04-24T12:03:11+02:00 (12 years ago)
Author:
flavoni
Message:

correct bug for Meridionnal Stream function in diaptr.F90 when using mpi, see ticket #421 and #361

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/DIA/diaptr.F90

    r1397 r1413  
    125125         DO jj = 2, jpjm1 
    126126            DO ji = fs_2, fs_jpim1   ! Vector opt. 
    127                p_fval(jj) = p_fval(jj) + pva(ji,jj,jk) * tmask_i(ji,jj+1) * tmask_i(ji,jj)  
     127               p_fval(jj) = p_fval(jj) + pva(ji,jj,jk) * tmask_i(ji,jj)  
    128128            END DO 
    129129         END DO 
     
    160160      DO jj = 2, jpjm1 
    161161         DO ji = nldi, nlei   ! No vector optimisation here. Better use a mask ? 
    162             p_fval(jj) = p_fval(jj) + pva(ji,jj) * tmask_i(ji,jj+1) * tmask_i(ji,jj) 
     162            p_fval(jj) = p_fval(jj) + pva(ji,jj) * tmask_i(ji,jj) 
    163163         END DO 
    164164      END DO 
     
    199199               DO ji =  nldi, nlei   ! No vector optimisation here. Better use a mask ? 
    200200                  p_fval(jj,jk) = p_fval(jj,jk) + pva(ji,jj,jk) * e1v(ji,jj) * fse3v(ji,jj,jk)   & 
    201                      &                                          * tmask_i(ji,jj+1) * tmask_i(ji,jj) & 
    202                      &                                          * bmask(ji,jj) 
     201                     &                                          * tmask_i(ji,jj) * bmask(ji,jj) 
    203202               END DO 
    204203            END DO 
     
    209208               DO ji =  nldi, nlei   ! No vector optimisation here. Better use a mask ? 
    210209                  p_fval(jj,jk) = p_fval(jj,jk) + pva(ji,jj,jk) * e1v(ji,jj) * fse3v(ji,jj,jk)   & 
    211                      &                                          * tmask_i(ji,jj+1) * tmask_i(ji,jj) 
     210                     &                                          * tmask_i(ji,jj) 
    212211               END DO 
    213212            END DO 
Note: See TracChangeset for help on using the changeset viewer.