SUBROUTINE  PLTKP (Z1,ZE,Z2,NZ,KX,NX,KY,NY,PL)

C     [CONTOUR PLOT]
C     FUNCTION VALUES TAKEN FROM AN ARRAY
C     ESSENTIALLY PLTKO, BUT WITHOUT FRAMING AND EJECT
C     ZE(NX,NY)  ARRAY TO BE CONTOURED
C     (KX,KY)    NUMBER OF X- AND Y-SUBDIVISIONS
C     (Z1,Z2)    CONTOURING INTERVAL
C     NZ         NUMBER OF CONTOUR LEVELS
C     PL         COORDINATE CONVERSION SUBROUTINE
C     [16-FEB-74]

      EXTERNAL    PL
      DIMENSION   ZE(1)

      I0=MAX0(NX/(2*KX),1)
      J0=MAX0(NY/(2*KY),1)
      IX=2*I0
      IY=2*J0
      DZ=(Z2-Z1)/FLOAT(NZ-1)

      II=I0+1
      JJ=J0+1
   10 JA=MAX0(JJ-J0,1)
      JB=MIN0(JJ+J0,NY)
   20 IA=MAX0(II-I0,1)
      IB=MIN0(II+I0,NX)
      Z=Z1
      DO 30 K=1,NZ
      CALL KONSC (Z,0.0,0.0,IA,IB,JA,JB,ZE,NX,NY,PL)
   30 Z=Z+DZ
      II=II+IX
      IF ((II-I0.LT.NX).AND.(II+I0.GT.1)) GO TO 20
      IX=-IX
      II=II+IX
      JJ=JJ+IY
      IF (JJ-J0.LT.NY) GO TO 10
      RETURN

      END