#include "dims.h"
      SUBROUTINE ADVECL(NX,S,K)
      use cons_module,only: imaxp2,dlamda_2div3,dlamda_1div12,
     |  dphi_2div3,dphi_1div12,re_inv
      implicit none
#include "params.h"
#include "fcom.h"
#include "buff.h"
#include "index.h"
#include "phys.h"
!
! Args:
!     DIMENSION S(1)
      real,intent(out) :: s(zimxp)
      integer,intent(in) :: nx,k
!
! Local:
      integer :: i,nuk,nxjp2k,nxjp1k,nxk,nxjm1k,nxjm2k,nvk,nvjp2k,
     |  nvjp1k,nvjm1k,nvjm2k
!
      NUK=NJ+NU+(K-1)
      NXJP2K=NJP2+NX
      NXJP1K=NJP1+NX
      NXK=NJ+NX
      NXJM1K=NJM1+NX
      NXJM2K=NJM2+NX
      NVK=NV+K-1
      NVJP2K=NJP2+NVK
      NVJP1K=NJP1+NVK
      NVJM1K=NJM1+NVK
      NVJM2K=NJM2+NVK
      DO I=3,IMAXP2
        S(I)=.5*RACS*(dlamda_2div3*(F(I+1,NXK)-F(I-1,NXK))*(F(I+1,NUK)+
     |    F(I-1,NUK))-dlamda_1div12*(F(I+2,NXK)-F(I-2,NXK))*(F(I+2,NUK)+
     |    F(I-2,NUK)))+.5*re_inv*(dphi_2div3*(F(I,NXJP1K)-F(I,NXJM1K))*
     |    (F(I,NVJP1K)+F(I,NVJM1K))-dphi_1div12*
     |    (F(I,NXJP2K)-F(I,NXJM2K))*(F(I,NVJP2K)+F(I,NVJM2K)))
      enddo
      RETURN
      END
C