-
Notifications
You must be signed in to change notification settings - Fork 0
/
ch_read_hyd.f90
67 lines (56 loc) · 1.98 KB
/
ch_read_hyd.f90
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
subroutine ch_read_hyd
use basin_module
use input_file_module
use maximum_data_module
use channel_data_module
implicit none
character (len=80) :: titldum ! |title of file
character (len=80) :: header ! |header of file
integer :: eof ! |end of file
logical :: i_exist ! |check to determine if file exists
integer :: imax !units |description
integer :: ich !none |counter
eof = 0
imax = 0
inquire (file=in_cha%hyd, exist=i_exist)
if (.not. i_exist .or. in_cha%hyd == "null") then
allocate (ch_hyd(0:0))
else
do
open (105,file=in_cha%hyd)
read (105,*,iostat=eof) titldum
if (eof < 0) exit
read (105,*,iostat=eof) header
if (eof < 0) exit
do while (eof == 0)
read (105,*,iostat=eof) titldum
if (eof < 0) exit
imax = imax + 1
end do
db_mx%ch_hyd = imax
allocate (ch_hyd(0:imax))
rewind (105)
read (105,*,iostat=eof) titldum
if (eof < 0) exit
read (105,*,iostat=eof) header
if (eof < 0) exit
do ich = 1, db_mx%ch_hyd
read (105,*,iostat=eof) titldum
if (eof < 0) exit
backspace (105)
read (105,*,iostat=eof) ch_hyd(ich)
if (eof < 0) exit
ch_hyd(ich)%alpha_bnk = Exp(-ch_hyd(ich)%alpha_bnk)
if (ch_hyd(ich)%s <= 0.) ch_hyd(ich)%s = .0001
if (ch_hyd(ich)%n <= 0.01) ch_hyd(ich)%n = .01
if (ch_hyd(ich)%n >= 0.70) ch_hyd(ich)%n = 0.70
if (ch_hyd(ich)%l <= 0.) ch_hyd(ich)%l = .0010
if (ch_hyd(ich)%wdr <= 0.) ch_hyd(ich)%wdr = 3.5
if (ch_hyd(ich)%side <= 1.e-6) ch_hyd(ich)%side = 2.0
end do
close (105)
exit
enddo
endif
return
end subroutine ch_read_hyd