Changeset 1347


Ignore:
Timestamp:
2009-03-30T17:29:37+02:00 (12 years ago)
Author:
rblod
Message:

Bug correction in lim2 with vector optimization, see ticket #380

Location:
trunk/NEMO
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/LIM_SRC_2/limdyn_2.F90

    r1228 r1347  
    152152          
    153153         DO jj = 1, jpjm1 
    154             DO ji = 1, fs_jpim1   ! vector opt. 
     154            DO ji = 1, jpim1   ! NO vector opt. 
    155155               zu_io(ji,jj) = 0.5 * ( ui_ice(ji+1,jj+1) + ui_ice(ji+1,jj  ) ) - ssu_m(ji,jj) 
    156156               zv_io(ji,jj) = 0.5 * ( vi_ice(ji+1,jj+1) + vi_ice(ji  ,jj+1) ) - ssv_m(ji,jj) 
     
    159159         ! frictional velocity at T-point 
    160160         DO jj = 2, jpjm1 
    161             DO ji = fs_2, fs_jpim1   ! vector opt. 
     161            DO ji = 2, jpim1   ! NO vector opt. 
    162162               ust2s(ji,jj) = 0.5 * cw                                                          & 
    163163                  &         * (  zu_io(ji,jj) * zu_io(ji,jj) + zu_io(ji-1,jj) * zu_io(ji-1,jj)   & 
  • trunk/NEMO/LIM_SRC_2/limrhg_2.F90

    r1156 r1347  
    243243         !-  Determination of zc1u, zc2u, zc1v and zc2v. 
    244244         DO jj = k_j1+1, k_jpj-1 
    245             DO ji = fs_2, fs_jpim1 
     245            DO ji = 2, fs_jpim1   ! NO vector opt. 
    246246               !* zc1u , zc2v 
    247247               zvis11 = 2.0 * zviseta (ji-1,jj-1) + dm 
     
    371371            DO jj = k_j1+1, k_jpj-1 
    372372!CDIR NOVERRCHK 
    373                DO ji = fs_2, fs_jpim1 
     373               DO ji = 2, fs_jpim1   ! NO vector opt. 
    374374                  ! 
    375375                  ze11 =   akappa(ji,jj-1,1,1) * zu_a(ji+1,jj) + akappa(ji,jj-1,1,2) * zv_a(ji+1,jj) 
  • trunk/NEMO/LIM_SRC_2/limwri_2.F90

    r1339 r1347  
    129129      zcmo(:,:, 1:jpnoumax ) = 0.e0  
    130130      DO jj = 2 , jpjm1 
    131          DO ji = fs_2 , fs_jpim1 
     131         DO ji = 1 , jpim1   ! NO vector opt. 
    132132            zindh  = MAX( zzero , SIGN( zone , hicif(ji,jj) * (1.0 - frld(ji,jj) ) - 0.10 ) ) 
    133133            zinda  = MAX( zzero , SIGN( zone , ( 1.0 - frld(ji,jj) ) - 0.10 ) ) 
  • trunk/NEMO/OPA_SRC/SBC/sbcblk_clio.F90

    r1275 r1347  
    491491!CDIR COLLAPSE 
    492492         DO jj = 2, jpj 
    493             DO ji = fs_2, jpi   ! vector opt. 
     493            DO ji = 2, jpi   ! B grid : no vector opt. 
    494494               p_taui(ji,jj) = 0.5 * ( utau(ji-1,jj  ) + utau(ji-1,jj-1) ) 
    495495               p_tauj(ji,jj) = 0.5 * ( vtau(ji  ,jj-1) + vtau(ji-1,jj-1) ) 
  • trunk/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90

    r1270 r1347  
    118118         !     B-grid dynamics :  I-point (F-point with sea-ice indexation) 
    119119         DO jj = 2, jpj 
    120             DO ji = fs_2, jpi   ! vector opt. 
     120            DO ji = 2, jpi   ! B grid : no vector opt. 
    121121               ui_oce(ji,jj) = 0.5 * ( ssu_m(ji-1,jj  ) + ssu_m(ji-1,jj-1) ) * tmu(ji,jj) 
    122122               vi_oce(ji,jj) = 0.5 * ( ssv_m(ji  ,jj-1) + ssv_m(ji-1,jj-1) ) * tmu(ji,jj) 
Note: See TracChangeset for help on using the changeset viewer.