setprob.f.html CLAWPACK  
 Source file:   setprob.f
 Directory:   /Users/rjl/git/clawpack/classic/examples/acoustics_1d_heterogeneous
 Converted:   Wed Dec 28 2016 at 22:28:20   using clawcode2html
 This documentation file will not reflect any later changes in the source file. 

 
      subroutine setprob

      implicit none

      common /cqinit/ beta,ic
      integer ic
      double precision beta

      common /comaux/ Zl, cl, Zr, cr
      double precision rhol, Zl, cl, rhor, Zr, cr
c
c     # Set the material parameters for the acoustic equations
c
      character*25 fname

      fname = 'setprob.data'
      call opendatafile(7, fname)

c     # choice of initial data:
      read(7,*) ic
c     # beta for initial conditions:
      read(7,*) beta
c
c     # Piecewise constant medium with single interface at x=0
c     # Density and sound speed to left and right:
      read(7,*) rhol
      read(7,*) cl
      Zl = rhol*cl

      read(7,*) rhor
      read(7,*) cr
      Zr = rhor*cr

      return
      end