Changeset 1775


Ignore:
Timestamp:
2009-11-26T12:02:18+01:00 (12 years ago)
Author:
rblod
Message:

Correct arrays out of bounds in diaptr, see ticket #619

File:
1 edited

Legend:

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

    r1715 r1775  
    2727   USE in_out_manager 
    2828   USE lib_mpp 
     29   USE lbclnk 
    2930 
    3031   IMPLICIT NONE 
     
    398399            ! T times V on T points (include bolus velocities) 
    399400#if defined key_diaeiv  
    400             DO jj = 1, jpj 
     401            DO jj = 2, jpj 
    401402               DO ji = 1, jpi 
    402403                  vt(ji,jj,:) = tn(ji,jj,:) * ( vn(ji,jj,:) + vn(ji,jj-1,:) + u_eiv(ji,jj,:) + u_eiv(ji,jj-1,:) )*0.5 
     
    405406            END DO 
    406407#else 
    407             DO jj = 1, jpj 
     408            DO jj = 2, jpj 
    408409               DO ji = 1, jpi 
    409410                  vt(ji,jj,:) = tn(ji,jj,:) * ( vn(ji,jj,:) + vn(ji,jj-1,:) )*0.5 
    410                   vs(ji,jj,:) = sn(ji,jj,:) * ( vn(ji,jj,:) + vn(ji,jj-1,:)  )*0.5 
     411                  vs(ji,jj,:) = sn(ji,jj,:) * ( vn(ji,jj,:) + vn(ji,jj-1,:) )*0.5 
    411412               END DO 
    412413            END DO 
    413414#endif  
     415            CALL lbc_lnk( vs, 'V', -1. )   ;   CALL lbc_lnk( vt, 'V', -1. ) 
    414416 
    415417            ht_glo(:) = SUM( ptr_vjk( vt(:,:,:)), 2 ) 
Note: See TracChangeset for help on using the changeset viewer.