setdambreak_geo.f.html
Source file:   setdambreak_geo.f
Directory:   /home/rjl/git/rjleveque/clawpack-4.x/geoclaw/2d/lib
Converted:   Sun May 15 2011 at 19:15:41   using clawcode2html
This documentation file will not reflect any later changes in the source file.

 
c=========================================================================
      subroutine setdambreak()
c=========================================================================

c     This was copied from settsunami_geo.f and needs to be fixed up
c     for dam break problems.
      use geoclaw_module
      
      implicit double precision (a-h,o-z)
      character*25 fname
      logical foundFile

      include "call.i"

      write(parmunit,*) ' '
      write(parmunit,*) '--------------------------------------------'
      write(parmunit,*) 'SETDAMBREAK:'
      write(parmunit,*) '------------'

c       # read user parameters from setdambreak.data

      fname  = 'setdambreak.data'
      inquire(file=fname,exist=foundFile)
      if (.not. foundFile) then
        write(*,*) 'You must provide a file ', fname
        stop
      endif

      iunit = 7
      call opendatafile(iunit, fname)

      read(7,*) sealevel
      read(7,*) drytolerance
      read(7,*) wavetolerance
      read(7,*) depthdeep
      read(7,*) maxleveldeep
      read(7,*) coeffmanning
      read(7,*) frictiondepth
      close(7)

      write(parmunit,*) '   drytolerance:',drytolerance
      write(parmunit,*) '   wavetolerance:',wavetolerance
      write(parmunit,*) '   maxleveldeep:', maxleveldeep
      write(parmunit,*) '   depthdeep:', depthdeep
      write(parmunit,*) '   Manning coefficient:',coeffmanning
      write(parmunit,*) '   frictiondepth:',frictiondepth

      return
      end