-
Notifications
You must be signed in to change notification settings - Fork 3
/
HMC_Module_Args.f90
141 lines (116 loc) · 7.22 KB
/
HMC_Module_Args.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
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
!--------------------------------------------------------------------------------
! File: HMC_Module_Args.f90
! Author(s): Fabio Delogu, Francesco Silvestro, Simone Gabellani
! Created on February 15, 2017, 4:40 PM
!
! Module to read argument(s) defined by command-line
!--------------------------------------------------------------------------------
!--------------------------------------------------------------------------------
! Module Args
module HMC_Module_Args
!--------------------------------------------------------------------------------
! External module(s) and implicit none
use HMC_Module_Tools_Debug ! to import global variable(s) declaration
implicit none
!--------------------------------------------------------------------------------
contains
!--------------------------------------------------------------------------------
! Subroutine to read argument(s)
subroutine HMC_Args_Read(dUc, dUh, dCt, dCf, dCPI, dWTableHbr, dKSatRatio, dSlopeMax, &
dCN, dFrac, dWS, dWDL, &
sDomainName, &
sFileInfo, iArgsType)
!------------------------------------------------------------------------------------------
! Variable(s)
integer(kind = 4) :: iArgsN, iArgsType
real(kind = 4) :: dUc, dUh, dCt, dCf, dCPI, dWTableHbr, dKSatRatio, dSlopeMax
real(kind = 4) :: dCN, dFrac, dWS, dWDL
character(len = 700) :: sLineBuffer
character(len = 256) :: sDomainName
character(len = 700) :: sFileInfo
!------------------------------------------------------------------------------------------
!------------------------------------------------------------------------------------------
! Variable(s) initialization
iArgsN = -9999; iArgsType = -9999
dUc = 0.0; dUh = 0.0; dCt = 0.0; dCf = 0.0; dCPI = 0.0;
dWTableHbr = 0.0; dKSatRatio = 0.0; dSlopeMax = 0.0; sDomainName = "";
dCN = 0.0; dFrac = 0.0; dWS = 0.0; dWDL = 0.0;
!------------------------------------------------------------------------------------------
!------------------------------------------------------------------------------------------
! Get number of argument(s)
iArgsN = iargc()
!------------------------------------------------------------------------------------------
!------------------------------------------------------------------------------------------
! Get argument(s) from command line [VERSION < 3.x.x]
if (iArgsN == 9) then
!------------------------------------------------------------------------------------------
! Get friction coefficient in channels [m^0.5 s^-1] (Uc)
call getarg(1, sLineBuffer); read(sLineBuffer,*) dUc
! Get flow motion coefficient in hillslopes [s^-1] (Uh)
call getarg(2, sLineBuffer); read(sLineBuffer,*) dUh
! Get mean field capacity [-] (Ct)
call getarg(3, sLineBuffer); read(sLineBuffer,*) dCt
! Get infiltration capacity at saturation [-] (Cf)
call getarg(4, sLineBuffer); read(sLineBuffer,*) dCf
! Get domain name
call getarg(5, sLineBuffer); sDomainName = sLineBuffer
! Get soil humidity initial condition (CPI)
call getarg(6, sLineBuffer); read(sLineBuffer,*) dCPI
! Get maximum water capacity of the aquifer [mm] (WTableHbr)
call getarg(7, sLineBuffer); read(sLineBuffer,*) dWTableHbr
! Get anisotropy between the vertical and horizontal saturated conductivity, and to soil porosity [-]
call getarg(8, sLineBuffer); read(sLineBuffer,*) dKSatRatio
! Get subsoil maximum slope [%]
call getarg(9, sLineBuffer); read(sLineBuffer,*) dSlopeMax
! Define information file name
sFileInfo = trim(sDomainName)//'.info.txt'
! Define argument(s) type
iArgsType = 1
!------------------------------------------------------------------------------------------
! Get argument(s) from command line [VERSION >= 3.x.x]
elseif (iArgsN == 13) then
!------------------------------------------------------------------------------------------
! Get friction coefficient in channels [m^0.5 s^-1] (Uc)
call getarg(1, sLineBuffer); read(sLineBuffer,*) dUc
! Get flow motion coefficient in hillslopes [s^-1] (Uh)
call getarg(2, sLineBuffer); read(sLineBuffer,*) dUh
! Get mean field capacity [-] (Ct)
call getarg(3, sLineBuffer); read(sLineBuffer,*) dCt
! Get infiltration capacity at saturation [-] (Cf)
call getarg(4, sLineBuffer); read(sLineBuffer,*) dCf
! Get domain name
call getarg(5, sLineBuffer); sDomainName = sLineBuffer
! Get soil humidity initial condition (CPI)
call getarg(6, sLineBuffer); read(sLineBuffer,*) dCPI
! Get maximum water capacity of the aquifer [mm] (WTableHbr)
call getarg(7, sLineBuffer); read(sLineBuffer,*) dWTableHbr
! Get anisotropy between the vertical and horizontal saturated conductivity, and to soil porosity [-]
call getarg(8, sLineBuffer); read(sLineBuffer,*) dKSatRatio
! Get subsoil maximum slope [%]
call getarg(9, sLineBuffer); read(sLineBuffer,*) dSlopeMax
! Get curve number average initialization [-]
call getarg(10, sLineBuffer); read(sLineBuffer,*) dCN
! Get fracturation average initialization [-]
call getarg(11, sLineBuffer); read(sLineBuffer,*) dFrac
! Get watertable sources average initialization [1/t]
call getarg(12, sLineBuffer); read(sLineBuffer,*) dWS
! Get watertable deep losses average initialization [1/t]
call getarg(13, sLineBuffer); read(sLineBuffer,*) dWDL
! Define information file name
sFileInfo = trim(sDomainName)//'.info.txt'
! Define argument(s) type
iArgsType = 1
!------------------------------------------------------------------------------------------
elseif (iArgsN == 1) then
!------------------------------------------------------------------------------------------
! Get information file name
call getarg(1, sLineBuffer); sFileInfo = sLineBuffer
! Define argument(s) type
iArgsType = 2
!------------------------------------------------------------------------------------------
endif
!------------------------------------------------------------------------------------------
end subroutine HMC_Args_Read
!--------------------------------------------------------------------------------
end module HMC_Module_Args
!--------------------------------------------------------------------------------