program cfget_somellor
c
c   extracts cf from cfl3d.prout file
c   This program is HARDWIRED for the specific So Mellor
c   testcases somellor_sa.inp and somellor_sarc.inp
c
      open(2,file='cfl3d.prout',form='formatted',status='old')
      open(3,file='cf.dat',form='formatted',status='unknown')
c  read front matter that we do not need
      do n=1,269
        read(2,*)
      enddo
      sswall=0.
c  read cf
      write(3,'(''variables="s","cf"'')')
      do n=1,255
        read(2,*) idum,idum,idum,x,y,z,dum,dum,dum,cf
        if (n .eq. 1) then
          dist=x
        else
          dist=sqrt((x-xold)**2+(z-zold)**2)
        end if
        sswall=sswall+dist
        if (cf .lt. 0.) cf=-cf
        write(3,'(2e15.5)') sswall,cf
        xold=x
        zold=z
      enddo
      write(6,'('' tecplot file output to cf.dat'')')
      stop
      end