Changeset 6777


Ignore:
Timestamp:
2016-07-04T13:16:12+02:00 (4 years ago)
Author:
timgraham
Message:

Bug fixes in updateTS

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/NEMO/NST_SRC/agrif_opa_update.F90

    r6454 r6777  
    263263      !!--------------------------------------------- 
    264264      INTEGER, INTENT(in) :: i1,i2,j1,j2,k1,k2,n1,n2 
    265       REAL(wp),DIMENSION(i1:i2,j1:j2,k1:k2,n1:n2+1), INTENT(inout) :: ptab 
     265      REAL(wp),DIMENSION(i1:i2,j1:j2,k1:k2,n1:n2), INTENT(inout) :: ptab 
    266266      LOGICAL, INTENT(in) :: before 
    267267      !! 
     
    280280      IF (before) THEN 
    281281         zrho_xy = Agrif_rhox() * Agrif_rhoy()  
    282          DO jn = n1,n2 
     282         DO jn = n1,n2-1 
    283283            DO jk=k1,k2 
    284284               DO jj=j1,j2 
     
    292292            DO jj=j1,j2 
    293293               DO ji=i1,i2 
    294                   ptab(ji,jj,jk,n2+1) = zrho_xy * e1e2t(ji,jj) * e3t_n(ji,jj,jk) * tmask(ji,jj,jk) 
     294                  ptab(ji,jj,jk,n2) = tmask(ji,jj,jk) * zrho_xy * e1e2t(ji,jj) * e3t_n(ji,jj,jk)  
    295295               END DO 
    296296            END DO 
     
    305305           N_in = 0 
    306306           DO jk=k1,k2 !k2 = jpk of child grid 
    307              IF (ptab(ji,jj,jk,n2+1) == 0) EXIT 
     307             IF (ptab(ji,jj,jk,n2) == 0) EXIT 
    308308             N_in = N_in + 1 
    309              tabin(jk,:) = ptab(ji,jj,jk,n1:n2)/ptab(ji,jj,jk,n2+1) 
    310              h_in(N_in) = ptab(ji,jj,jk,n2+1)/e1e2t(ji,jj) 
     309             tabin(jk,:) = ptab(ji,jj,jk,n1:n2-1)/ptab(ji,jj,jk,n2) 
     310             h_in(N_in) = ptab(ji,jj,jk,n2)/e1e2t(ji,jj) 
    311311           ENDDO 
    312312           N_out = 0 
     
    331331!               h_out(N_out) = - h_diff 
    332332             ENDIF 
    333              DO jn=n1,n2 
     333             DO jn=n1,n2-1 
    334334               CALL reconstructandremap(tabin(1:N_in,jn),h_in(1:N_in),ptab_child(ji,jj,1:N_out,jn),h_out(1:N_out),N_in,N_out) 
    335335             ENDDO 
     
    345345         IF (.NOT.(lk_agrif_fstep.AND.(neuler==0))) THEN 
    346346            ! Add asselin part 
    347             DO jn = n1,n2 
     347            DO jn = n1,n2-1 
    348348               DO jk=1,jpk 
    349349                  DO jj=j1,j2 
     
    359359            ENDDO 
    360360         ENDIF 
    361          DO jn = n1,n2 
     361         DO jn = n1,n2-1 
    362362            DO jk=1,jpk 
    363363               DO jj=j1,j2 
Note: See TracChangeset for help on using the changeset viewer.