-
Notifications
You must be signed in to change notification settings - Fork 0
/
calsoft_read_codes.f90
47 lines (40 loc) · 1.61 KB
/
calsoft_read_codes.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
subroutine calsoft_read_codes
use calibration_data_module
use plant_data_module
use input_file_module
use soil_module
use plant_module
use hydrograph_module
use hru_lte_module
use sd_channel_module
use organic_mineral_mass_module
use mgt_operations_module
use conditional_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 !none |check to determine if file exists
eof = 0
inquire (file=in_chg%codes_sft, exist=i_exist)
if (.not. i_exist .or. in_chg%codes_sft == "null") then
! allocate (cal_codes(0:0))
else
do
open (107,file=in_chg%codes_sft)
read (107,*,iostat=eof) titldum
if (eof < 0) exit
read (107,*,iostat=eof) header
if (eof < 0) exit
read (107,*,iostat=eof) cal_codes
if (eof < 0) exit
exit
enddo
if (cal_codes%hyd_hru /= "n" .or. cal_codes%hyd_hrul == "y".or. &
cal_codes%plt == "y" .or. cal_codes%sed == "y" .or. &
cal_codes%nut == "y" .or. cal_codes%chsed == "y" .or. &
cal_codes%chnut == "y" .or. cal_codes%res == "y") cal_soft = "y"
end if
close(107)
return
end subroutine calsoft_read_codes