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 11205 for NEMO/branches/2019/dev_r10973_AGRIF-01_jchanut_small_jpi_jpj/src/NST/agrif_oce_interp.F90 – NEMO

Ignore:
Timestamp:
2019-07-02T14:25:46+02:00 (5 years ago)
Author:
jchanut
Message:

#2199
1) Make sponge independent of sub-domain size. Partially masked open boundary segments are not taken into account anymore. To do so, sponge coefficients should be read in a file for realistic applications (then nesting tools need to be modified accordingly).
2) Replace East-West-North-South barotropic data arrays by a global 2d array. Then apply barotropic open boundary conditions thanks to mi0/mi1, mj0/mj1 indexes.
3) Call AGRIF bdy update one more time in dynspg_ts during extrapolation phase. This removes a dozen lines of code in dynspg_ts routine.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2019/dev_r10973_AGRIF-01_jchanut_small_jpi_jpj/src/NST/agrif_oce_interp.F90

    r10068 r11205  
    401401      !! 
    402402      INTEGER :: ji, jj 
     403      INTEGER :: istart, iend, jstart, jend 
    403404      !!----------------------------------------------------------------------   
    404405      ! 
    405406      IF( Agrif_Root() )   RETURN 
    406407      ! 
    407       IF((nbondi == -1).OR.(nbondi == 2)) THEN 
     408      !--- West ---! 
     409      istart = 2 
     410      iend   = nbghostcells+1 
     411      DO ji = mi0(istart), mi1(iend) 
    408412         DO jj=1,jpj 
    409             va_e(2:nbghostcells+1,jj) = vbdy_w(1:nbghostcells,jj) * hvr_e(2:nbghostcells+1,jj) 
     413            va_e(ji,jj) = vbdy(ji,jj) * hvr_e(ji,jj) 
    410414            ! Specified fluxes: 
    411             ua_e(2:nbghostcells+1,jj) = ubdy_w(1:nbghostcells,jj) * hur_e(2:nbghostcells+1,jj) 
    412             ! Characteristics method (only if ghostcells=1): 
    413             !alt            ua_e(2,jj) = 0.5_wp * ( ubdy_w(jj) * hur_e(2,jj) + ua_e(3,jj) & 
    414             !alt                       &           - sqrt(grav * hur_e(2,jj)) * (sshn_e(3,jj) - hbdy_w(jj)) ) 
    415          END DO 
    416       ENDIF 
    417       ! 
    418       IF((nbondi == 1).OR.(nbondi == 2)) THEN 
     415            ua_e(ji,jj) = ubdy(ji,jj) * hur_e(ji,jj) 
     416         END DO 
     417      END DO 
     418      ! Characteristics method (only at boundary point): 
     419!      istart = 2 
     420!      iend   = 2 
     421!      DO ji = mi0(istart), mi1(iend) 
     422!         DO jj=1,jpj 
     423!            ua_e(ji,jj) = 0.5_wp * ( ubdy(ji,jj) * hur_e(ji,jj) + ua_e(ji+1,jj) & 
     424!                        &            - sqrt(grav * hur_e(ji,jj)) * (sshn_e(ji+1,jj) - hbdy(ji,jj)) ) 
     425!         END DO 
     426!      END DO 
     427      ! 
     428      !--- East ---! 
     429      istart = jpiglo-nbghostcells 
     430      iend   = jpiglo-1 
     431      DO ji = mi0(istart), mi1(iend) 
    419432         DO jj=1,jpj 
    420             va_e(nlci-nbghostcells:nlci-1,jj)   = vbdy_e(1:nbghostcells,jj) * hvr_e(nlci-nbghostcells:nlci-1,jj) 
     433            va_e(ji,jj) = vbdy(ji,jj) * hvr_e(ji,jj) 
     434         END DO 
     435      END DO 
     436      istart = jpiglo-nbghostcells-1 
     437      iend   = jpiglo-2 
     438      DO ji = mi0(istart), mi1(iend) 
     439         DO jj=1,jpj 
     440            ua_e(ji,jj) = ubdy(ji,jj) * hur_e(ji,jj) 
     441         END DO 
     442      END DO 
     443      ! Characteristics method (only at boundary point): 
     444!      istart = jpiglo-2 
     445!      iend   = jpiglo-2 
     446!      DO ji = mi0(istart), mi1(iend) 
     447!         DO jj=1,jpj 
     448!            ua_e(ji,jj) = 0.5_wp * ( ubdy(ji,jj) * hur_e(ji,jj) + ua_e(ji-1,jj) & 
     449!                        &            + sqrt(grav * hur_e(ji,jj)) * (sshn_e(ji,jj) - hbdy(ji+1,jj)) ) 
     450!         END DO 
     451!      END DO 
     452      ! 
     453      !--- South ---! 
     454      jstart = 2 
     455      jend   = nbghostcells+1 
     456      DO jj = mj0(jstart), mj1(jend) 
     457         DO ji=1,jpi 
     458            ua_e(ji,jj) = ubdy(ji,jj) * hur_e(ji,jj) 
    421459            ! Specified fluxes: 
    422             ua_e(nlci-nbghostcells-1:nlci-2,jj) = ubdy_e(1:nbghostcells,jj) * hur_e(nlci-nbghostcells-1:nlci-2,jj) 
    423             ! Characteristics method (only if ghostcells=1): 
    424             !alt            ua_e(nlci-2,jj) = 0.5_wp * ( ubdy_e(jj) * hur_e(nlci-2,jj) + ua_e(nlci-3,jj) & 
    425             !alt                            &           + sqrt(grav * hur_e(nlci-2,jj)) * (sshn_e(nlci-2,jj) - hbdy_e(jj)) ) 
    426          END DO 
    427       ENDIF 
    428       ! 
    429       IF((nbondj == -1).OR.(nbondj == 2)) THEN 
     460            va_e(ji,jj) = vbdy(ji,jj) * hvr_e(ji,jj) 
     461         END DO 
     462      END DO 
     463      ! Characteristics method (only at boundary point): 
     464!      jstart = 2 
     465!      jend   = 2 
     466!      DO jj = mj0(jstart), mj1(jend) 
     467!         DO ji=1,jpi 
     468!            va_e(ji,jj) = 0.5_wp * ( vbdy(ji,jj) * hvr_e(ji,jj) + va_e(ji,jj+1) & 
     469!                        &            - sqrt(grav * hvr_e(ji,jj)) * (sshn_e(ji,jj+1) - hbdy(ji,jj)) ) 
     470!         END DO 
     471!      END DO 
     472      ! 
     473      !--- North ---! 
     474      jstart = jpjglo-nbghostcells 
     475      jend   = jpjglo-1 
     476      DO jj = mj0(jstart), mj1(jend) 
    430477         DO ji=1,jpi 
    431             ua_e(ji,2:nbghostcells+1) = ubdy_s(ji,1:nbghostcells) * hur_e(ji,2:nbghostcells+1) 
    432             ! Specified fluxes: 
    433             va_e(ji,2:nbghostcells+1) = vbdy_s(ji,1:nbghostcells) * hvr_e(ji,2:nbghostcells+1) 
    434             ! Characteristics method (only if ghostcells=1): 
    435             !alt            va_e(ji,2) = 0.5_wp * ( vbdy_s(ji) * hvr_e(ji,2) + va_e(ji,3) & 
    436             !alt                       &           - sqrt(grav * hvr_e(ji,2)) * (sshn_e(ji,3) - hbdy_s(ji)) ) 
    437          END DO 
    438       ENDIF 
    439       ! 
    440       IF((nbondj == 1).OR.(nbondj == 2)) THEN 
     478            ua_e(ji,jj) = ubdy(ji,jj) * hur_e(ji,jj) 
     479         END DO 
     480      END DO 
     481      jstart = jpjglo-nbghostcells-1 
     482      jend   = jpjglo-2 
     483      DO jj = mj0(jstart), mj1(jend) 
    441484         DO ji=1,jpi 
    442             ua_e(ji,nlcj-nbghostcells:nlcj-1)   = ubdy_n(ji,1:nbghostcells) * hur_e(ji,nlcj-nbghostcells:nlcj-1) 
    443             ! Specified fluxes: 
    444             va_e(ji,nlcj-nbghostcells-1:nlcj-2) = vbdy_n(ji,1:nbghostcells) * hvr_e(ji,nlcj-nbghostcells-1:nlcj-2) 
    445             ! Characteristics method (only if ghostcells=1): 
    446             !alt            va_e(ji,nlcj-2) = 0.5_wp * ( vbdy_n(ji) * hvr_e(ji,nlcj-2)  + va_e(ji,nlcj-3) & 
    447             !alt                            &           + sqrt(grav * hvr_e(ji,nlcj-2)) * (sshn_e(ji,nlcj-2) - hbdy_n(ji)) ) 
    448          END DO 
    449       ENDIF 
     485            va_e(ji,jj) = vbdy(ji,jj) * hvr_e(ji,jj) 
     486         END DO 
     487      END DO 
     488      ! Characteristics method (only at boundary point): 
     489!      jstart = jpjglo-2 
     490!      jend   = jpjglo-2 
     491!      DO jj = mj0(jstart), mj1(jend) 
     492!         DO ji=1,jpi 
     493!            va_e(ji,jj) = 0.5_wp * ( vbdy(ji,jj) * hvr_e(ji,jj) + va_e(ji,jj-1) & 
     494!                        &            + sqrt(grav * hvr_e(ji,jj)) * (sshn_e(ji,jj) - hbdy(ji,jj+1)) ) 
     495!         END DO 
     496!      END DO 
    450497      ! 
    451498   END SUBROUTINE Agrif_dyn_ts 
    … …  
    485532      ELSE ! Linear interpolation 
    486533         bdy_tinterp = 0 
    487          ubdy_w(:,:) = 0._wp   ;   vbdy_w(:,:) = 0._wp  
    488          ubdy_e(:,:) = 0._wp   ;   vbdy_e(:,:) = 0._wp  
    489          ubdy_n(:,:) = 0._wp   ;   vbdy_n(:,:) = 0._wp  
    490          ubdy_s(:,:) = 0._wp   ;   vbdy_s(:,:) = 0._wp 
     534         ubdy(:,:) = 0._wp   ;   vbdy(:,:) = 0._wp  
    491535         CALL Agrif_Bc_variable( unb_id, procname=interpunb ) 
    492536         CALL Agrif_Bc_variable( vnb_id, procname=interpvnb ) 
    … …  
    503547      INTEGER, INTENT(in) ::   kt 
    504548      ! 
    505       INTEGER  :: ji, jj, indx, indy 
     549      INTEGER  :: ji, jj 
     550      INTEGER  :: istart, iend, jstart, jend 
    506551      !!----------------------------------------------------------------------   
    507552      ! 
    … …  
    516561      ! 
    517562      ! --- West --- ! 
    518       IF((nbondi == -1).OR.(nbondi == 2)) THEN 
    519          indx = 1+nbghostcells 
     563      istart = 2 
     564      iend   = 1 + nbghostcells 
     565      DO ji = mi0(istart), mi1(iend) 
    520566         DO jj = 1, jpj 
    521             DO ji = 2, indx 
    522                ssha(ji,jj) = hbdy_w(ji-1,jj) 
    523             ENDDO 
     567            ssha(ji,jj) = hbdy(ji,jj) 
    524568         ENDDO 
    525       ENDIF 
     569      ENDDO 
    526570      ! 
    527571      ! --- East --- ! 
    528       IF((nbondi == 1).OR.(nbondi == 2)) THEN 
    529          indx = nlci-nbghostcells 
     572      istart = jpiglo - nbghostcells 
     573      iend   = jpiglo - 1 
     574      DO ji = mi0(istart), mi1(iend) 
    530575         DO jj = 1, jpj 
    531             DO ji = indx, nlci-1 
    532                ssha(ji,jj) = hbdy_e(ji-indx+1,jj) 
    533             ENDDO 
     576            ssha(ji,jj) = hbdy(ji,jj) 
    534577         ENDDO 
    535       ENDIF 
     578      ENDDO 
    536579      ! 
    537580      ! --- South --- ! 
    538       IF((nbondj == -1).OR.(nbondj == 2)) THEN 
    539          indy = 1+nbghostcells 
    540          DO jj = 2, indy 
    541             DO ji = 1, jpi 
    542                ssha(ji,jj) = hbdy_s(ji,jj-1) 
    543             ENDDO 
     581      jstart = 2 
     582      jend   = 1 + nbghostcells 
     583      DO jj = mj0(jstart), mj1(jend) 
     584         DO ji = 1, jpi 
     585            ssha(ji,jj) = hbdy(ji,jj) 
    544586         ENDDO 
    545       ENDIF 
     587      ENDDO 
    546588      ! 
    547589      ! --- North --- ! 
    548       IF((nbondj == 1).OR.(nbondj == 2)) THEN 
    549          indy = nlcj-nbghostcells 
    550          DO jj = indy, nlcj-1 
    551             DO ji = 1, jpi 
    552                ssha(ji,jj) = hbdy_n(ji,jj-indy+1) 
    553             ENDDO 
     590      jstart = jpjglo - nbghostcells 
     591      jend   = jpjglo - 1 
     592      DO jj = mj0(jstart), mj1(jend) 
     593         DO ji = 1, jpi 
     594            ssha(ji,jj) = hbdy(ji,jj) 
    554595         ENDDO 
    555       ENDIF 
     596      ENDDO 
    556597      ! 
    557598   END SUBROUTINE Agrif_ssh 
    … …  
    564605      INTEGER, INTENT(in) ::   jn 
    565606      !! 
    566       INTEGER :: ji, jj, indx, indy 
    567       !!----------------------------------------------------------------------   
    568       !! clem ghost (starting at i,j=1 is important I think otherwise you introduce a grad(ssh)/=0 at point 2) 
     607      INTEGER :: ji, jj 
     608      INTEGER  :: istart, iend, jstart, jend 
     609      !!----------------------------------------------------------------------   
    569610      ! 
    570611      IF( Agrif_Root() )   RETURN 
    571612      ! 
    572613      ! --- West --- ! 
    573       IF((nbondi == -1).OR.(nbondi == 2)) THEN 
    574          indx = 1+nbghostcells 
     614      istart = 2 
     615      iend   = 1+nbghostcells 
     616      DO ji = mi0(istart), mi1(iend) 
    575617         DO jj = 1, jpj 
    576             DO ji = 2, indx 
    577                ssha_e(ji,jj) = hbdy_w(ji-1,jj) 
    578             ENDDO 
     618            ssha_e(ji,jj) = hbdy(ji,jj) 
    579619         ENDDO 
    580       ENDIF 
     620      ENDDO 
    581621      ! 
    582622      ! --- East --- ! 
    583       IF((nbondi == 1).OR.(nbondi == 2)) THEN 
    584          indx = nlci-nbghostcells 
     623      istart = jpiglo - nbghostcells 
     624      iend   = jpiglo - 1 
     625      DO ji = mi0(istart), mi1(iend) 
    585626         DO jj = 1, jpj 
    586             DO ji = indx, nlci-1 
    587                ssha_e(ji,jj) = hbdy_e(ji-indx+1,jj) 
    588             ENDDO 
     627            ssha_e(ji,jj) = hbdy(ji,jj) 
    589628         ENDDO 
    590       ENDIF 
     629      ENDDO 
    591630      ! 
    592631      ! --- South --- ! 
    593       IF((nbondj == -1).OR.(nbondj == 2)) THEN 
    594          indy = 1+nbghostcells 
    595          DO jj = 2, indy 
    596             DO ji = 1, jpi 
    597                ssha_e(ji,jj) = hbdy_s(ji,jj-1) 
    598             ENDDO 
     632      jstart = 2 
     633      jend   = 1+nbghostcells 
     634      DO jj = mj0(jstart), mj1(jend) 
     635         DO ji = 1, jpi 
     636            ssha_e(ji,jj) = hbdy(ji,jj) 
    599637         ENDDO 
    600       ENDIF 
     638      ENDDO 
    601639      ! 
    602640      ! --- North --- ! 
    603       IF((nbondj == 1).OR.(nbondj == 2)) THEN 
    604          indy = nlcj-nbghostcells 
    605          DO jj = indy, nlcj-1 
    606             DO ji = 1, jpi 
    607                ssha_e(ji,jj) = hbdy_n(ji,jj-indy+1) 
    608             ENDDO 
     641      jstart = jpjglo - nbghostcells 
     642      jend   = jpjglo - 1 
     643      DO jj = mj0(jstart), mj1(jend) 
     644         DO ji = 1, jpi 
     645            ssha_e(ji,jj) = hbdy(ji,jj) 
    609646         ENDDO 
    610       ENDIF 
     647      ENDDO 
    611648      ! 
    612649   END SUBROUTINE Agrif_ssh_ts 
    … …  
    848885   END SUBROUTINE interptsn 
    849886 
    850    SUBROUTINE interpsshn( ptab, i1, i2, j1, j2, before, nb, ndir ) 
     887   SUBROUTINE interpsshn( ptab, i1, i2, j1, j2, before ) 
    851888      !!---------------------------------------------------------------------- 
    852889      !!                  ***  ROUTINE interpsshn  *** 
    … …  
    855892      REAL(wp), DIMENSION(i1:i2,j1:j2), INTENT(inout) ::   ptab 
    856893      LOGICAL                         , INTENT(in   ) ::   before 
    857       INTEGER                         , INTENT(in   ) ::   nb , ndir 
    858       ! 
    859       LOGICAL :: western_side, eastern_side,northern_side,southern_side 
     894      ! 
    860895      !!----------------------------------------------------------------------   
    861896      ! 
    … …  
    863898         ptab(i1:i2,j1:j2) = sshn(i1:i2,j1:j2) 
    864899      ELSE 
    865          western_side  = (nb == 1).AND.(ndir == 1) 
    866          eastern_side  = (nb == 1).AND.(ndir == 2) 
    867          southern_side = (nb == 2).AND.(ndir == 1) 
    868          northern_side = (nb == 2).AND.(ndir == 2) 
    869          !! clem ghost 
    870          IF(western_side)  hbdy_w(1:nbghostcells,j1:j2) = ptab(i1:i2,j1:j2) * tmask(i1:i2,j1:j2,1) 
    871          IF(eastern_side)  hbdy_e(1:nbghostcells,j1:j2) = ptab(i1:i2,j1:j2) * tmask(i1:i2,j1:j2,1) 
    872          IF(southern_side) hbdy_s(i1:i2,1:nbghostcells) = ptab(i1:i2,j1:j2) * tmask(i1:i2,j1:j2,1)  
    873          IF(northern_side) hbdy_n(i1:i2,1:nbghostcells) = ptab(i1:i2,j1:j2) * tmask(i1:i2,j1:j2,1) 
     900         hbdy(i1:i2,j1:j2) = ptab(i1:i2,j1:j2) * tmask(i1:i2,j1:j2,1) 
    874901      ENDIF 
    875902      ! 
    … …  
    10451072   END SUBROUTINE interpvn 
    10461073 
    1047    SUBROUTINE interpunb( ptab, i1, i2, j1, j2, before, nb, ndir ) 
     1074   SUBROUTINE interpunb( ptab, i1, i2, j1, j2, before) 
    10481075      !!---------------------------------------------------------------------- 
    10491076      !!                  ***  ROUTINE interpunb  *** 
    … …  
    10521079      REAL(wp), DIMENSION(i1:i2,j1:j2), INTENT(inout) ::   ptab 
    10531080      LOGICAL                         , INTENT(in   ) ::   before 
    1054       INTEGER                         , INTENT(in   ) ::   nb , ndir 
    10551081      ! 
    10561082      INTEGER  ::   ji, jj 
    10571083      REAL(wp) ::   zrhoy, zrhot, zt0, zt1, ztcoeff 
    1058       LOGICAL  ::   western_side, eastern_side,northern_side,southern_side 
    10591084      !!----------------------------------------------------------------------   
    10601085      ! 
    … …  
    10621087         ptab(i1:i2,j1:j2) = e2u(i1:i2,j1:j2) * hu_n(i1:i2,j1:j2) * un_b(i1:i2,j1:j2) 
    10631088      ELSE 
    1064          western_side  = (nb == 1).AND.(ndir == 1) 
    1065          eastern_side  = (nb == 1).AND.(ndir == 2) 
    1066          southern_side = (nb == 2).AND.(ndir == 1) 
    1067          northern_side = (nb == 2).AND.(ndir == 2) 
    10681089         zrhoy = Agrif_Rhoy() 
    10691090         zrhot = Agrif_rhot() 
    … …  
    10821103         ENDIF 
    10831104         !    
    1084          IF(western_side)   ubdy_w(1:nbghostcells,j1:j2) = ubdy_w(1:nbghostcells,j1:j2) + ztcoeff * ptab(i1:i2,j1:j2)   
    1085          IF(eastern_side)   ubdy_e(1:nbghostcells,j1:j2) = ubdy_e(1:nbghostcells,j1:j2) + ztcoeff * ptab(i1:i2,j1:j2)   
    1086          IF(southern_side)  ubdy_s(i1:i2,1:nbghostcells) = ubdy_s(i1:i2,1:nbghostcells) + ztcoeff * ptab(i1:i2,j1:j2) 
    1087          IF(northern_side)  ubdy_n(i1:i2,1:nbghostcells) = ubdy_n(i1:i2,1:nbghostcells) + ztcoeff * ptab(i1:i2,j1:j2)  
     1105         ubdy(i1:i2,j1:j2) = ubdy(i1:i2,j1:j2) + ztcoeff * ptab(i1:i2,j1:j2) 
    10881106         !             
    10891107         IF( bdy_tinterp == 0 .OR. bdy_tinterp == 2) THEN 
    1090             IF(western_side)   ubdy_w(1:nbghostcells,j1:j2) = ubdy_w(1:nbghostcells,j1:j2) / (zrhoy*e2u(i1:i2,j1:j2)) * umask(i1:i2,j1:j2,1) 
    1091             IF(eastern_side)   ubdy_e(1:nbghostcells,j1:j2) = ubdy_e(1:nbghostcells,j1:j2) / (zrhoy*e2u(i1:i2,j1:j2)) * umask(i1:i2,j1:j2,1) 
    1092             IF(southern_side)  ubdy_s(i1:i2,1:nbghostcells) = ubdy_s(i1:i2,1:nbghostcells) / (zrhoy*e2u(i1:i2,j1:j2)) * umask(i1:i2,j1:j2,1) 
    1093             IF(northern_side)  ubdy_n(i1:i2,1:nbghostcells) = ubdy_n(i1:i2,1:nbghostcells) / (zrhoy*e2u(i1:i2,j1:j2)) * umask(i1:i2,j1:j2,1) 
     1108            ubdy(i1:i2,j1:j2) = ubdy(i1:i2,j1:j2) / (zrhoy*e2u(i1:i2,j1:j2)) * umask(i1:i2,j1:j2,1) 
    10941109         ENDIF 
    10951110      ENDIF 
    … …  
    10981113 
    10991114 
    1100    SUBROUTINE interpvnb( ptab, i1, i2, j1, j2, before, nb, ndir ) 
     1115   SUBROUTINE interpvnb( ptab, i1, i2, j1, j2, before ) 
    11011116      !!---------------------------------------------------------------------- 
    11021117      !!                  ***  ROUTINE interpvnb  *** 
    … …  
    11051120      REAL(wp), DIMENSION(i1:i2,j1:j2), INTENT(inout) ::   ptab 
    11061121      LOGICAL                         , INTENT(in   ) ::   before 
    1107       INTEGER                         , INTENT(in   ) ::   nb , ndir 
    11081122      ! 
    11091123      INTEGER  ::   ji,jj 
    11101124      REAL(wp) ::   zrhox, zrhot, zt0, zt1, ztcoeff    
    1111       LOGICAL  ::   western_side, eastern_side,northern_side,southern_side 
    11121125      !!----------------------------------------------------------------------   
    11131126      !  
    … …  
    11151128         ptab(i1:i2,j1:j2) = e1v(i1:i2,j1:j2) * hv_n(i1:i2,j1:j2) * vn_b(i1:i2,j1:j2) 
    11161129      ELSE 
    1117          western_side  = (nb == 1).AND.(ndir == 1) 
    1118          eastern_side  = (nb == 1).AND.(ndir == 2) 
    1119          southern_side = (nb == 2).AND.(ndir == 1) 
    1120          northern_side = (nb == 2).AND.(ndir == 2) 
    11211130         zrhox = Agrif_Rhox() 
    11221131         zrhot = Agrif_rhot() 
    … …  
    11331142            ztcoeff = 1 
    11341143         ENDIF 
    1135          !! clem ghost 
    1136          IF(western_side)   vbdy_w(1:nbghostcells,j1:j2) = vbdy_w(1:nbghostcells,j1:j2) + ztcoeff * ptab(i1:i2,j1:j2)   
    1137          IF(eastern_side)   vbdy_e(1:nbghostcells,j1:j2) = vbdy_e(1:nbghostcells,j1:j2) + ztcoeff * ptab(i1:i2,j1:j2)    
    1138          IF(southern_side)  vbdy_s(i1:i2,1:nbghostcells) = vbdy_s(i1:i2,1:nbghostcells) + ztcoeff * ptab(i1:i2,j1:j2) 
    1139          IF(northern_side)  vbdy_n(i1:i2,1:nbghostcells) = vbdy_n(i1:i2,1:nbghostcells) + ztcoeff * ptab(i1:i2,j1:j2)  
     1144         vbdy(i1:i2,j1:j2) = vbdy(i1:i2,j1:j2) + ztcoeff * ptab(i1:i2,j1:j2) 
    11401145         !             
    11411146         IF( bdy_tinterp == 0 .OR. bdy_tinterp == 2) THEN 
    1142             IF(western_side)   vbdy_w(1:nbghostcells,j1:j2) = vbdy_w(1:nbghostcells,j1:j2) / (zrhox*e1v(i1:i2,j1:j2)) * vmask(i1:i2,j1:j2,1) 
    1143             IF(eastern_side)   vbdy_e(1:nbghostcells,j1:j2) = vbdy_e(1:nbghostcells,j1:j2) / (zrhox*e1v(i1:i2,j1:j2)) * vmask(i1:i2,j1:j2,1) 
    1144             IF(southern_side)  vbdy_s(i1:i2,1:nbghostcells) = vbdy_s(i1:i2,1:nbghostcells) / (zrhox*e1v(i1:i2,j1:j2)) * vmask(i1:i2,j1:j2,1) 
    1145             IF(northern_side)  vbdy_n(i1:i2,1:nbghostcells) = vbdy_n(i1:i2,1:nbghostcells) / (zrhox*e1v(i1:i2,j1:j2)) * vmask(i1:i2,j1:j2,1) 
    1146          ENDIF 
     1147            vbdy(i1:i2,j1:j2) = vbdy(i1:i2,j1:j2) / (zrhox*e1v(i1:i2,j1:j2)) * vmask(i1:i2,j1:j2,1) 
     1148         ENDIF             
    11471149      ENDIF 
    11481150      ! 
    … …  
    11501152 
    11511153 
    1152    SUBROUTINE interpub2b( ptab, i1, i2, j1, j2, before, nb, ndir ) 
     1154   SUBROUTINE interpub2b( ptab, i1, i2, j1, j2, before ) 
    11531155      !!---------------------------------------------------------------------- 
    11541156      !!                  ***  ROUTINE interpub2b  *** 
    … …  
    11571159      REAL(wp), DIMENSION(i1:i2,j1:j2), INTENT(inout) ::   ptab 
    11581160      LOGICAL                         , INTENT(in   ) ::   before 
    1159       INTEGER                         , INTENT(in   ) ::   nb , ndir 
    11601161      ! 
    11611162      INTEGER  ::   ji,jj 
    1162       REAL(wp) ::   zrhot, zt0, zt1,zat 
    1163       LOGICAL  ::   western_side, eastern_side,northern_side,southern_side 
     1163      REAL(wp) ::   zrhot, zt0, zt1, zat 
    11641164      !!----------------------------------------------------------------------   
    11651165      IF( before ) THEN 
    … …  
    11701170         ENDIF 
    11711171      ELSE 
    1172          western_side  = (nb == 1).AND.(ndir == 1) 
    1173          eastern_side  = (nb == 1).AND.(ndir == 2) 
    1174          southern_side = (nb == 2).AND.(ndir == 1) 
    1175          northern_side = (nb == 2).AND.(ndir == 2) 
    11761172         zrhot = Agrif_rhot() 
    11771173         ! Time indexes bounds for integration 
    … …  
    11811177         zat = zrhot * (  zt1**2._wp * (-2._wp*zt1 + 3._wp)    & 
    11821178            &           - zt0**2._wp * (-2._wp*zt0 + 3._wp)    )  
    1183          !! clem ghost 
    1184          IF(western_side ) ubdy_w(1:nbghostcells,j1:j2) = zat * ptab(i1:i2,j1:j2)   
    1185          IF(eastern_side ) ubdy_e(1:nbghostcells,j1:j2) = zat * ptab(i1:i2,j1:j2)   
    1186          IF(southern_side) ubdy_s(i1:i2,1:nbghostcells) = zat * ptab(i1:i2,j1:j2) 
    1187          IF(northern_side) ubdy_n(i1:i2,1:nbghostcells) = zat * ptab(i1:i2,j1:j2)  
     1179         ! 
     1180         ubdy(i1:i2,j1:j2) = zat * ptab(i1:i2,j1:j2)  
    11881181      ENDIF 
    11891182      !  
    … …  
    11911184    
    11921185 
    1193    SUBROUTINE interpvb2b( ptab, i1, i2, j1, j2, before, nb, ndir ) 
     1186   SUBROUTINE interpvb2b( ptab, i1, i2, j1, j2, before ) 
    11941187      !!---------------------------------------------------------------------- 
    11951188      !!                  ***  ROUTINE interpvb2b  *** 
    … …  
    11981191      REAL(wp), DIMENSION(i1:i2,j1:j2), INTENT(inout) ::   ptab 
    11991192      LOGICAL                         , INTENT(in   ) ::   before 
    1200       INTEGER                         , INTENT(in   ) ::   nb , ndir 
    12011193      ! 
    12021194      INTEGER ::   ji,jj 
    12031195      REAL(wp) ::   zrhot, zt0, zt1,zat 
    1204       LOGICAL ::   western_side, eastern_side,northern_side,southern_side 
    12051196      !!----------------------------------------------------------------------   
    12061197      ! 
    … …  
    12121203         ENDIF 
    12131204      ELSE       
    1214          western_side  = (nb == 1).AND.(ndir == 1) 
    1215          eastern_side  = (nb == 1).AND.(ndir == 2) 
    1216          southern_side = (nb == 2).AND.(ndir == 1) 
    1217          northern_side = (nb == 2).AND.(ndir == 2) 
    12181205         zrhot = Agrif_rhot() 
    12191206         ! Time indexes bounds for integration 
    … …  
    12241211            &           - zt0**2._wp * (-2._wp*zt0 + 3._wp)    )  
    12251212         ! 
    1226          IF(western_side )   vbdy_w(1:nbghostcells,j1:j2) = zat * ptab(i1:i2,j1:j2)   
    1227          IF(eastern_side )   vbdy_e(1:nbghostcells,j1:j2) = zat * ptab(i1:i2,j1:j2)   
    1228          IF(southern_side)   vbdy_s(i1:i2,1:nbghostcells) = zat * ptab(i1:i2,j1:j2) 
    1229          IF(northern_side)   vbdy_n(i1:i2,1:nbghostcells) = zat * ptab(i1:i2,j1:j2)  
     1213         vbdy(i1:i2,j1:j2) = zat * ptab(i1:i2,j1:j2) 
    12301214      ENDIF 
    12311215      !       
Note: See TracChangeset for help on using the changeset viewer.