New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 1694 for trunk/NEMO/LIM_SRC_2 – NEMO

Changeset 1694 for trunk/NEMO/LIM_SRC_2


Ignore:
Timestamp:
2009-10-30T17:06:23+01:00 (15 years ago)
Author:
smasson
Message:

avoid out of bounds access, see ticket:576

Location:
trunk/NEMO/LIM_SRC_2
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMO/LIM_SRC_2/limdyn_2.F90

    r1470 r1694  
    158158         ! frictional velocity at T-point 
    159159         DO jj = 2, jpjm1 
    160             DO ji = 2, jpim1   ! NO vector opt. 
     160            DO ji = 2, jpim1   ! NO vector opt. because of zu_io 
    161161               ust2s(ji,jj) = 0.5 * cw                                                          & 
    162162                  &         * (  zu_io(ji,jj) * zu_io(ji,jj) + zu_io(ji-1,jj) * zu_io(ji-1,jj)   & 

  • trunk/NEMO/LIM_SRC_2/limmsh_2.F90

    r1156 r1694  
    157157 
    158158      DO jj = 2, jpj 
    159          DO ji = 2, jpi 
     159         DO ji = 2, jpi   ! NO vector opt. 
    160160            zh1p  =  e1t(ji  ,jj  ) * wght(ji,jj,2,2)   & 
    161161               &   + e1t(ji-1,jj  ) * wght(ji,jj,1,2)   & 
    … …  
    234234      tmu(1,:) = 0.e0 
    235235      DO jj = 2, jpj               ! ice U.V-point: computed from ice T-point mask 
    236          DO ji = 2, jpim1 
     236         DO ji = 2, jpim1   ! NO vector opt. 
    237237            tmu(ji,jj) =  tms(ji,jj) * tms(ji-1,jj) * tms(ji,jj-1) * tms(ji-1,jj-1)             
    238238         END DO 

  • trunk/NEMO/LIM_SRC_2/limsbc_2.F90

    r1550 r1694  
    258258 
    259259         DO jj = 2, jpjm1 
    260             DO ji = 2, jpim1 
     260            DO ji = 2, jpim1   ! NO vector opt. 
    261261               ! ... ice-cover wheighted ice-ocean stress at U and V-points  (from I-point values) 
    262262               zutau  = 0.5 * ( ztio_u(ji+1,jj) + ztio_u(ji+1,jj+1) ) 

  • trunk/NEMO/LIM_SRC_2/limtrp_2.F90

    r1470 r1694  
    109109         zvbord = 1.0 + ( 1.0 - bound ) 
    110110         DO jj = 1, jpjm1 
    111             DO ji = 1, jpim1 
     111            DO ji = 1, jpim1   ! NO vector opt. 
    112112               zui_u(ji,jj) = ( u_ice(ji+1,jj  ) + u_ice(ji+1,jj+1) ) / ( MAX( tmu(ji+1,jj  ) + tmu(ji+1,jj+1), zvbord ) ) 
    113113               zvi_v(ji,jj) = ( v_ice(ji  ,jj+1) + v_ice(ji+1,jj+1) ) / ( MAX( tmu(ji  ,jj+1) + tmu(ji+1,jj+1), zvbord ) ) 

  • trunk/NEMO/LIM_SRC_2/limwri_dimg_2.h90

    r1470 r1694  
    9595    zcmo(:,:, 1:jpnoumax ) = 0.e0  
    9696    DO jj = 2 , jpjm1 
    97        DO ji = 2 , jpim1 
     97       DO ji = 2 , jpim1   ! NO vector opt. 
    9898          zindh  = MAX( zzero , SIGN( zone , hicif(ji,jj) * (1.0 - frld(ji,jj) ) - 0.10 ) ) 
    9999          zinda  = MAX( zzero , SIGN( zone , ( 1.0 - frld(ji,jj) ) - 0.10 ) ) 
    … …  
    138138          !  case of instantaneaous output rcmoy(:,:, 1:jpnoumax ) = 0.e0 
    139139          DO jj = 2 , jpjm1 
    140              DO ji = 2 , jpim1 
     140             DO ji = 2 , jpim1   ! NO vector opt. 
    141141                zindh  = MAX( zzero , SIGN( zone , hicif(ji,jj) * (1.0 - frld(ji,jj) ) - 0.10 ) ) 
    142142                zinda  = MAX( zzero , SIGN( zone , ( 1.0 - frld(ji,jj) ) - 0.10 ) ) 
Note: See TracChangeset for help on using the changeset viewer.