|
setprob.f.html |
|
|
Source file: setprob.f
|
|
Directory: /Users/rjl/git/clawpack/amrclaw/examples/advection_2d_swirl
|
|
Converted: Wed Dec 28 2016 at 22:38:04
using clawcode2html
|
|
This documentation file will
not reflect any later changes in the source file.
|
subroutine setprob
implicit double precision (a-h,o-z)
character(len=25) fname
common /compsi/ pi
common /comvt/ tperiod,pi2
c
c # compute pi, used in psi.f
pi = 4.d0 * datan(1.d0)
c
c # save 2*pi and tperiod in common block for use in b4step2:
c
pi2 = 2.d0*pi
c
iunit = 7
fname = 'setprob.data'
c # open the unit with new routine from Clawpack 4.4 to skip over
c # comment lines starting with #:
call opendatafile(iunit, fname)
read(7,*) tperiod
return
end