program fppost2

c     Reads the skin friction data Genplot file output by CFPOST 
c     and creates skin friction along plate.

      parameter ( nimax = 200 ) 

      character cline * 72

      dimension x (nimax)
      dimension y (nimax)
      dimension u (nimax)
      dimension v (nimax)

c.....Read x and y

      do n = 1, 6
        read(*,'(a72)') cline
      enddo 

      read (*,*) ni

      do  i = 1, ni
        read(*,*) x(i), y(i)
      enddo

c.....Read u

      do n = 1, 6
        read(*,'(a72)') cline
      enddo 

      do  i = 1, ni
        read(*,*) xx, u(i)
      enddo

c.....Read v

      do n = 1, 6
        read(*,'(a72)') cline
      enddo 

      do  i = 1, ni
        read(*,*) xx, v(i)
      enddo

c.....Compute tau and cf.

      ufs    = 129.6974
      dynvis = 6.5044372E-04

      open ( unit=7, file='cf.d' )

      do  i = 1, ni
        if ( x(i) .gt. 0.0 ) then
          tau = u(i) / y(i)
          cf  = 2.0 * tau * dynvis / ( ufs * ufs )
          write (7,'(2(f12.5))') x(i), cf
        endif
      enddo

      stop
      end