-
Notifications
You must be signed in to change notification settings - Fork 0
/
ch_read_init.f90
58 lines (49 loc) · 1.59 KB
/
ch_read_init.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
subroutine ch_read_init
use basin_module
use input_file_module
use maximum_data_module
use channel_data_module
use sd_channel_module
implicit none
character (len=80) :: titldum ! |title of file
character (len=80) :: header ! |header of file
integer :: eof ! |end of file
integer :: imax !units |description
logical :: i_exist ! |check to determine if file exists
integer :: ich !none |counter
eof = 0
imax = 0
inquire (file=in_cha%init, exist=i_exist)
if (.not. i_exist .or. in_cha%init == "null") then
allocate (ch_init(0:0))
allocate (sd_init(0:0))
else
do
open (105,file=in_cha%init)
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_init = imax
allocate (ch_init(0:imax))
allocate (sd_init(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_init
read (105,*,iostat=eof) ch_init(ich)
if (eof < 0) exit
end do
close (105)
exit
enddo
endif
return
end subroutine ch_read_init