program fppost1

c     Reads the boundary layer Genplot file output by CFPOST 
c     and creates profiles in similarity coordinates.

      parameter ( njmax = 100 ) 

      character cline * 72

      dimension y (njmax)
      dimension u (njmax)
      dimension v (njmax)
      dimension t (njmax)

c.....Read y and x

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

      read (*,*) nj

      do  j = 1, nj
        read(*,*) y(j), xx
      enddo

c.....Read u

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

      do  j = 1, nj
        read(*,*) yy, u(j)
      enddo

c.....Read v

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

      do  j = 1, nj
        read(*,*) yy, v(j)
      enddo

c.....Read t

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

      do  j = 1, nj
        read(*,*) yy, t(j)
      enddo

c.....Convert y to eta and nondimensionalize u.

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

      ufs    = 129.6974
      dynvis = 6.5044372E-04

      do  j = 1, nj
        eta = y(j) * sqrt( 0.5 * ufs / ( dynvis * xx ) )
        if ( eta .le. 7.0 ) write (7,'(2(f12.5))') u(j)/ufs, eta
      enddo

c.....Convert y to eta and nondimensionalize v.

      open ( unit=8, file='v.d' ) 

      rex = ufs * xx / dynvis

      do  j = 1, nj
        eta  = y(j) * sqrt( 0.5 * ufs / ( dynvis * xx ) )
        vnon = v(j) * sqrt( 2.0 * rex ) / ufs
        if ( eta .le. 7.0 ) write (8,'(2(f12.5))') vnon, eta
      enddo

      stop
      end