Changeset 7641


Ignore:
Timestamp:
2017-02-03T16:19:34+01:00 (4 years ago)
Author:
clem
Message:

bug fix for AGRIF+LIM2-VP, see ticket #1607

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/NST_SRC/agrif_lim2_update.F90

    r6204 r7641  
    177177         tabres = zrhox * tabres 
    178178      ELSE 
    179          DO jj=j1,j2 
    180             DO ji=i1,i2 
     179         DO jj=MAX(j1,2),j2 
     180            DO ji=MAX(i1,2),i2 
    181181               v_ice(ji,jj) = tabres(ji,jj) / (e1f(ji-1,jj-1)) 
    182182               v_ice(ji,jj) = v_ice(ji,jj) * tmu(ji,jj) 
     
    202202      IF( before ) THEN 
    203203         zrhoy = Agrif_Rhoy() 
    204          DO jj=MAX(j1,2),j2 
    205             DO ji=MAX(i1,2),i2 
     204         DO jj=j1,j2 
     205            DO ji=i1,i2 
    206206               tabres(ji,jj) = e2u(ji,jj) * u_ice(ji,jj) 
    207207            END DO 
     
    209209         tabres = zrhoy * tabres 
    210210      ELSE 
    211          DO jj=MAX(j1,2),j2 
    212             DO ji=MAX(i1,2),i2 
     211         DO jj=j1,j2 
     212            DO ji=i1,i2 
    213213               u_ice(ji,jj) = tabres(ji,jj) / (e2u(ji,jj)) 
    214214               u_ice(ji,jj) = u_ice(ji,jj) * tmu(ji,jj) 
     
    235235      IF( before ) THEN 
    236236         zrhox = Agrif_Rhox() 
    237          DO jj=MAX(j1,2),j2 
    238             DO ji=MAX(i1,2),i2 
     237         DO jj=j1,j2 
     238            DO ji=i1,i2 
    239239               tabres(ji,jj) = e1v(ji,jj) * v_ice(ji,jj) 
    240240            END DO 
Note: See TracChangeset for help on using the changeset viewer.