-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathbasin_reservoir_output.f90
104 lines (92 loc) · 3.9 KB
/
basin_reservoir_output.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
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
subroutine basin_reservoir_output
use time_module
use basin_module
use reservoir_module
use hydrograph_module
use water_body_module
implicit none
integer :: ires !none |counter
real :: const
!! zero daily variables
bres = resmz
bres_in_d = resmz
bres_out_d = resmz
bres_wat_d = wbodz
!! sum all reservoirs to get basin totals
do ires = 1, sp_ob%res
bres = bres + res(ires)
bres_in_d = bres_in_d + res_in_d(ires)
res_in_d(ires) = resmz
bres_out_d = bres_out_d + res_out_d(ires)
res_out_d(ires) = resmz
bres_wat_d = bres_wat_d + res_wat_d(ires)
!res_wat_d(ires) = wbodz
end do
!! sum monthly variables
bres_in_m = bres_in_m + bres_in_d
bres_out_m = bres_out_m + bres_out_d
bres_wat_m = bres_wat_m + bres_wat_d
!! daily print - RESERVOIR
if (pco%day_print == "y" .and. pco%int_day_cur == pco%int_day) then
if (pco%res_bsn%d == "y") then
write (2100,100) time%day, time%mo, time%day_mo, time%yrc, ires, " 1", bsn%name, bres_wat_d, &
bres, bres_in_d, bres_out_d
if (pco%csvout == "y") then
write (2104,'(*(G0.3,:","))') time%day, time%mo, time%day_mo, time%yrc, ires, " 1", bsn%name, &
bres_wat_d, bres, bres_in_d, bres_out_d
end if
end if
end if
!! monthly print - RESERVOIR
if (time%end_mo == 1) then
const = float (ndays(time%mo + 1) - ndays(time%mo))
bres_in_y = bres_in_y + bres_in_m
bres_out_y = bres_out_y + bres_out_m
bres_wat_m = bres_wat_m // const
bres_wat_y = bres_wat_y + bres_wat_m
if (pco%res_bsn%m == "y") then
write (2101,100) time%day, time%mo, time%day_mo, time%yrc, ires, " 1", bsn%name, bres_wat_m, bres, &
bres_in_m, bres_out_m
if (pco%csvout == "y") then
write (2105,'(*(G0.3,:","))') time%day, time%mo, time%day_mo, time%yrc, ires, " 1", bsn%name, &
bres_wat_m, bres, bres_in_m, bres_out_m
endif
end if
bres_in_m = resmz
bres_out_m = resmz
bres_wat_m = wbodz
end if
!! yearly print - RESERVOIR
if (time%end_yr == 1) then
bres_in_a = bres_in_a + bres_in_y
bres_out_a = bres_out_a + bres_out_y
bres_wat_y = bres_wat_y // 12.
bres_wat_a = bres_wat_a + bres_wat_y
if (pco%res_bsn%y == "y") then
write (2102,100) time%day, time%mo, time%day_mo, time%yrc, ires, " 1", bsn%name, bres_wat_y, &
bres, bres_in_y, bres_out_y
if (pco%csvout == "y") then
write (2106,'(*(G0.3,:","))') time%day, time%mo, time%day_mo, time%yrc, ires, " 1", bsn%name, &
bres_wat_y, bres, bres_in_y, bres_out_y
end if
end if
!! zero yearly variables
bres_in_y = resmz
bres_out_y = resmz
bres_wat_y = wbodz
end if
!! average annual print - RESERVOIR
if (time%end_sim == 1 .and. pco%res_bsn%a == "y") then
bres_in_a = bres_in_a / time%yrs_prt
bres_out_a = bres_out_a / time%yrs_prt
bres_wat_a = bres_wat_a / time%yrs_prt
write (2103,100) time%day, time%mo, time%day_mo, time%yrc, ires, " 1", bsn%name, bres_wat_a, bres, &
bres_in_a, bres_out_a
if (pco%csvout == "y") then
write (2107,'(*(G0.3,:","))') time%day, time%mo, time%day_mo, time%yrc, ires, " 1", bsn%name, &
bres_wat_a, bres, bres_in_a, bres_out_a
end if
end if
return
100 format (4i6,i8,2x,a,2x,a17,f14.4,59(1x,e14.4))
end subroutine basin_reservoir_output