vtf-logo

fsi/sfc-amroc/ThinPlate/src/combl.f


cc     =====================================================
      subroutine combl()
c     =====================================================
c     
c     Create and initialize application specific common-blocks.
c
c     Copyright (C) 2003-2007 California Institute of Technology
c     Ralf Deiterding, ralf@amroc.net
c
      implicit double precision (a-h,o-z)
c
      include  "cuser.i"
      character *16 cwork
      parameter( lmechout = 13 )

      Wref    = 1.d6
      a2ref   = 1.d3
c     Air
      gamma  = 1.4d0  
      gamma1 = 0.4d0
      pi = 4.d0*atan(1.d0)
      Wk   = 29.08d-3/Wref
      cwork= 'AIR'
      RU = 0.83140000E+01/(Wref*a2ref)
      PA = 0.10132500E+06/a2ref
c
      open(unit=lmechout, status='unknown', form='formatted', 
     &     file='chem.dat')
      write (lmechout,400) RU
      write (lmechout,401) PA
      write (lmechout,402) gamma
      write (lmechout,403) cwork     
      write (lmechout,404) Wk        
      close (lmechout)
c
 400  format('RU       ',e16.8)
 401  format('PA       ',e16.8)
 402  format('Gamma    ',e16.8)
 403  format('Sp       ',a16)
 404  format('W        ',e16.8)
c
      idisc = 1
      alf   = 1.d0
      beta  = 0.d0
      x0    = 0.d0
      y0    = 0.d0
c       
      rhoamb = 1.d0
      uamb   = 0.d0
      vamb   = 0.d0
      wamb   = 0.d0
      pamb   = 1.d5
      eamb   = pamb/gamma1 + .5d0*rhoamb*(uamb**2+vamb**2+wamb**2)
c
      rhoshk = 1.d0
      ushk   = 0.d0
      vshk   = 0.d0
      wshk   = 0.d0
      pshk   = 11.d5
      eshk   = pshk/gamma1 + .5d0*rhoshk*(ushk**2+vshk**2+wshk**2)
c
      return
      end