program blasius

c     Writes out the Blasius solution for incompressible, 
c     laminar flow past a flat plate.

      parameter ( njmax = 45 )

      dimension eta (njmax)
      dimension f   (njmax)
      dimension fp  (njmax)

c.....Read in eta, f, f'.

      open ( unit=7, file='fplam.blasius' )

      do  n = 1, 4
        read ( 7, '(a72)' ) cline
      enddo

      do  n = 1, 41
        read (7,*) eta(n), f(n), fp(n)
      enddo

c.....Write out solution for u-velocity.

      open ( unit=8, file='u.blasius', form='formatted' )

      do  n = 1, 41
        write (8,'(2(f12.5))') fp(n), eta(n)
      enddo

c.....Write out solution for v-velocity.

      open ( unit=8, file='v.blasius', form='formatted' )

      do  n = 1, 41
        vnon = eta(n) * fp(n) - f(n)
        write (8,'(2(f12.5))') vnon, eta(n) 
      enddo

c.....Write out skin friction coefficient.

      open ( unit=9, file='cf.blasius', form='formatted' )

      ufs    = 129.6974
      dynvis = 6.5044372E-04

      do  i = 2, 30
        x   = ( i - 1.0 ) / 29.0
        rex = ufs * x / dynvis
        cf  = 0.664 / sqrt( rex )
        write (9,'(2(f12.5))') x, cf 
      enddo

      stop
      end