vtf-logo

weno/applications/euler/3d/Homogeneous/src/combl.f

c     =====================================================
      subroutine combl()
c     =====================================================
c     
c     Create and initialize application specific common-blocks.
c
      implicit none
c
      include 'air.i'
      include 'cuser.i'
      include 'spiralsgs.i'
      
      integer lin
      parameter( lin = 5 )
      
c
      call stdcdat()
c
      open(unit=lin,status='old',form='formatted',file='init.dat')
      read (lin, *) ihomo_sgs, ihomo_str, ihomo_flag
      close (lin)
      
      Re = 375.0d0
      Mach = 0.3d0
      Pr = 0.7d0 
      alpha = Pr*Re*gamma1*Mach*Mach
      beta = gamma*Mach*Mach

      if ( ihomo_sgs .eq. 1 ) then
         call LES_Spiral_AlignmentModel(0)
      endif

      if ( ihomo_str .eq. 1 ) then
         dir = 1
      else
         dir = 0
      endif

      return
      end