C $Header: /u/gcmpack/MITgcm/pkg/generic_advdiff/gad_osc_hat_y.F,v 1.1 2016/03/13 01:44:02 jmc Exp $
C $Name:  $

#     include "GAD_OPTIONS.h"

C--  File gad_osc_hat_y.F: Routines ???
C--   Contents
C--   o GAD_OSC_LOC_Y
C--   o GAD_OSC_HAT_Y

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

      SUBROUTINE GAD_OSC_LOC_Y(iy,mask,fbar,ohat)

          implicit none

C     =============================================== global variables
#         include "SIZE.h"

C     ====================================================== arguments
          integer iy
          _RL mask(1-OLy:sNy+OLy)
          _RL fbar(1-OLy:sNy+OLy)
          _RL ohat(1:2,
     &             1-OLy:sNy+OLy)

C     ====================================================== variables
          _RL floc(-2:+2)

          if (iy.gt. +1-OLy .and.
     &        iy.lt.sNy+OLy) then

C     ================ mask local stencil: expand from centre outwards

          floc(+0) = fbar(+0+iy)

          floc(-1) = floc(+0) +
     &      mask(iy-1)*(fbar(iy-1)-floc(+0))
          floc(+1) = floc(+0) +
     &      mask(iy+1)*(fbar(iy+1)-floc(+0))

C     ================ calc. 1st & 2nd derivatives over masked stencil

          ohat(+1,iy) = floc(+1)*0.25 _d 0
     &                - floc(-1)*0.25 _d 0

          ohat(+2,iy) = floc(+1)*0.25 _d 0
     &                - floc(+0)*0.50 _d 0
     &                + floc(-1)*0.25 _d 0

          else

          if (iy.eq. +1-OLy) then

C     ================ mask local stencil: expand from centre outwards

          floc(+0) = fbar(+0+iy)

          floc(+1) = floc(+0) +
     &      mask(iy+1)*(fbar(iy+1)-floc(+0))
          floc(+2) = floc(+1) +
     &      mask(iy+2)*(fbar(iy+2)-floc(+1))

C     ================ calc. 1st & 2nd derivatives over masked stencil

          ohat(+1,iy) = floc(+1)*0.50 _d 0
     &                - floc(+0)*0.50 _d 0

          ohat(+2,iy) = floc(+2)*0.25 _d 0
     &                - floc(+1)*0.50 _d 0
     &                + floc(+0)*0.25 _d 0

          end


if if (iy.eq.sNy+OLy) then C ================ mask local stencil: expand from centre outwards floc(+0) = fbar(+0+iy) floc(-1) = floc(+0) + & mask(iy-1)*(fbar(iy-1)-floc(+0)) floc(-2) = floc(-1) + & mask(iy-2)*(fbar(iy-2)-floc(-1)) C ================ calc. 1st & 2nd derivatives over masked stencil ohat(+1,iy) = floc(+0)*0.50 _d 0 & - floc(-1)*0.50 _d 0 ohat(+2,iy) = floc(+0)*0.25 _d 0 & - floc(-1)*0.50 _d 0 & + floc(-2)*0.25 _d 0 end


if end


if return c end subroutine GAD_OSC_LOC_Y end


C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| SUBROUTINE GAD_OSC_HAT_Y(bi,bj,kk,ix, & mask,fbar,ohat,myThid) C |================================================================| C | OSC_HAT_Y: compute WENO oscillation derivatives in Y. | C |================================================================| implicit none C =============================================== global variables # include "SIZE.h" C ====================================================== arguments integer bi,bj,kk,ix _RL mask(1-OLy:sNy+OLy) _RL fbar(1-OLy:sNy+OLy) _RL ohat(1:2, & 1-OLy:sNy+OLy) integer myThid C ====================================================== variables integer iy C ================================ derivatives for WENO indicators do iy = 1-OLy+0, sNy+OLy-0 CALL GAD_OSC_LOC_Y(iy,mask,fbar,ohat) end


do return c end subroutine GAD_OSC_HAT_Y end