SUBROUTINE COUTLI(Z,S,JX) 3

      IMPLICIT NONE
      INCLUDE 'champ.cdk'
      INCLUDE 'chinc.cdk'
      INCLUDE 'chobs.cdk'

      COMPLEX Z(IKTX,IKTY),S(IKTX,IKTY)

      REAL KX,KY,KW,JX

      DO IKX=1,IKTX
         KX = FLOAT(IKX-1) 
         DO IKY=1,IKTY
            KY = FLOAT(IKY - KTY - 1)
            NZK(IKX,IKY) = CMPLX(0.,0.)
            NSK(IKX,IKY) = CMPLX(0.,0.)
            IF ( L(IKX+KTX,IKY).eq.1 ) then
               KW = MAX(KX*KX+KY*KY ,0.001 ) 
               NZK(IKX,IKY) = (ZT(IKX,IKY)-ZOBS(IKX,IKY,NT+2)+Z(IKX,IKY))
     .              *  PZ(IKX,IKY,NT+2)
C     NSK(IKX,IKY) = (ST(IKX,IKY)-SOBS(IKX,IKY,NT+2)+S(IKX,IKY))
c     .              *  PS(IKX,IKY,NT+2)
               JX = JX + 0.5*REAL(NZK(IKX,IKY)*CONJG(NZK(IKX,IKY)))/KW
c     JX = JX + 0.5*REAL(NSK(IKX,IKY)*CONJG(NSK(IKX,IKY)))
            end if
         ENDDO
      ENDDO

      RETURN
      END