setdambreak_geo.f.html | ![]() |
Source file: setdambreak_geo.f | |
Directory: /Users/rjl/git/rjleveque/clawpack-4.6.3/geoclaw/2d/lib | |
Converted: Mon Jan 21 2013 at 20:15:57 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