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 3283 for branches/2011/dev_NEMO_MERGE_2011/NEMOGCM – NEMO

Ignore:
Timestamp:
2012-01-27T14:15:50+01:00 (12 years ago)
Author:
rblod
Message:

Array out of bound in diafwb, see ticket #919

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/DIA/diafwb.F90

    r3168 r3283  
    176176         END SELECT 
    177177         !  
    178          DO ji = mi0(ii0), mi1(ii1) 
     178         DO ji = mi0(ii0), MIN(mi1(ii1),jpim1) 
    179179            DO jj = mj0(ij0), mj1(ij1) 
    180180               DO jk = 1, jpk  
    181181                  zt = 0.5 * ( tsn(ji,jj,jk,jp_tem) + tsn(ji+1,jj,jk,jp_tem) ) 
    182182                  zs = 0.5 * ( tsn(ji,jj,jk,jp_sal) + tsn(ji+1,jj,jk,jp_sal) ) 
    183                   zu = un(ji,jj,jk) * fse3t(ji,jj,jk) * e2u(ji,jj) 
     183                  zu = un(ji,jj,jk) * fse3t(ji,jj,jk) * e2u(ji,jj) * tmask_i(ji,jj) 
    184184 
    185185                  IF( un(ji,jj,jk) > 0.e0 ) THEN  
    … …  
    223223         END SELECT 
    224224         !  
    225          DO ji = mi0(ii0), mi1(ii1) 
     225         DO ji = mi0(ii0), MIN(mi1(ii1),jpim1) 
    226226            DO jj = mj0(ij0), mj1(ij1) 
    227227               DO jk = 1, jpk  
    228228                  zt = 0.5 * ( tsn(ji,jj,jk,jp_tem) + tsn(ji+1,jj,jk,jp_tem) ) 
    229229                  zs = 0.5 * ( tsn(ji,jj,jk,jp_sal) + tsn(ji+1,jj,jk,jp_sal) ) 
    230                   zu = un(ji,jj,jk) * fse3t(ji,jj,jk) * e2u(ji,jj) 
     230                  zu = un(ji,jj,jk) * fse3t(ji,jj,jk) * e2u(ji,jj) * tmask_i(ji,jj) 
    231231                   
    232232                  IF( un(ji,jj,jk) > 0.e0 ) THEN  
    … …  
    270270         END SELECT 
    271271         !  
    272          DO ji = mi0(ii0), mi1(ii1) 
     272         DO ji = mi0(ii0), MIN(mi1(ii1),jpim1) 
    273273            DO jj = mj0(ij0), mj1(ij1) 
    274274               DO jk = 1, jpk  
    275275                  zt = 0.5 * ( tsn(ji,jj,jk,jp_tem) + tsn(ji+1,jj,jk,jp_tem) ) 
    276276                  zs = 0.5 * ( tsn(ji,jj,jk,jp_sal) + tsn(ji+1,jj,jk,jp_sal) ) 
    277                   zu = un(ji,jj,jk) * fse3t(ji,jj,jk) * e2u(ji,jj) 
     277                  zu = un(ji,jj,jk) * fse3t(ji,jj,jk) * e2u(ji,jj) * tmask_i(ji,jj) 
    278278                   
    279279                  IF( un(ji,jj,jk) > 0.e0 ) THEN  
    … …  
    317317         END SELECT 
    318318         !  
    319          DO ji = mi0(ii0), mi1(ii1) 
     319         DO ji = mi0(ii0), MIN(mi1(ii1),jpim1) 
    320320            DO jj = mj0(ij0), mj1(ij1) 
    321321               DO jk = 1, jpk 
    322322                  zt = 0.5 * ( tsn(ji,jj,jk,jp_tem) + tsn(ji+1,jj,jk,jp_tem) ) 
    323323                  zs = 0.5 * ( tsn(ji,jj,jk,jp_sal) + tsn(ji+1,jj,jk,jp_sal) ) 
    324                   zu = un(ji,jj,jk) * fse3t(ji,jj,jk) * e2u(ji,jj) 
     324                  zu = un(ji,jj,jk) * fse3t(ji,jj,jk) * e2u(ji,jj) * tmask_i(ji,jj) 
    325325                   
    326326                  IF( un(ji,jj,jk) > 0.e0 ) THEN  
Note: See TracChangeset for help on using the changeset viewer.