3
This commit is contained in:
parent
6c7785fc6c
commit
d043f8e979
103
Fortran/1.json
103
Fortran/1.json
@ -1,8 +1,9 @@
|
|||||||
{
|
{
|
||||||
"fortranSourceFile": "D:\\\u5DE5\u4F5C2025\\\u7B97\u6CD5\u5E73\u53F0\\dll\\\u65B0\u5B89\u6C5F\u6A21\u578B\\XAJMX\\XAJMX\\XAJMX.f90",
|
"fortranSourceFile": "D:\\\u5DE5\u4F5C2025\\\u7B97\u6CD5\u5E73\u53F0\\dll\\\u65B0\u5B89\u6C5F\u6A21\u578B\\XAJMX\\XAJMX\\XAJMX - \u526F\u672C.f90",
|
||||||
"fortranFunctionName": "XAJMX",
|
"fortranFunctionName": "XAJMX",
|
||||||
"projectName": "FortranWebApi",
|
"projectName": "FortranWebApi",
|
||||||
"outputDirectory": "D:\\\u5DE5\u4F5C2025\\\u7B97\u6CD5\u5E73\u53F0\\\u8F93\u51FA\u5E93\\\u65B0\u5B89\u6C5F\u6A21\u578B",
|
"outputDirectory": "D:\\\u5DE5\u4F5C2025\\\u7B97\u6CD5\u5E73\u53F0\\\u8F93\u51FA\u5E93\\\u65B0\u5B89\u6C5F\u6A21\u578B",
|
||||||
|
"cmbLanguage": "Fortran",
|
||||||
"parameters": [
|
"parameters": [
|
||||||
{
|
{
|
||||||
"name": "N",
|
"name": "N",
|
||||||
@ -21,9 +22,105 @@
|
|||||||
"isSelected": true
|
"isSelected": true
|
||||||
},
|
},
|
||||||
{
|
{
|
||||||
"name": "PAR",
|
"name": "WUM",
|
||||||
"dataType": "Float",
|
"dataType": "Float",
|
||||||
"arrayType": "OneDimensional",
|
"arrayType": "Scalar",
|
||||||
|
"direction": "Input",
|
||||||
|
"description": "",
|
||||||
|
"isSelected": true
|
||||||
|
},
|
||||||
|
{
|
||||||
|
"name": "W1",
|
||||||
|
"dataType": "Float",
|
||||||
|
"arrayType": "Scalar",
|
||||||
|
"direction": "Input",
|
||||||
|
"description": "",
|
||||||
|
"isSelected": true
|
||||||
|
},
|
||||||
|
{
|
||||||
|
"name": "WDM",
|
||||||
|
"dataType": "Float",
|
||||||
|
"arrayType": "Scalar",
|
||||||
|
"direction": "Input",
|
||||||
|
"description": "",
|
||||||
|
"isSelected": true
|
||||||
|
},
|
||||||
|
{
|
||||||
|
"name": "KC",
|
||||||
|
"dataType": "Float",
|
||||||
|
"arrayType": "Scalar",
|
||||||
|
"direction": "Input",
|
||||||
|
"description": "",
|
||||||
|
"isSelected": true
|
||||||
|
},
|
||||||
|
{
|
||||||
|
"name": "C",
|
||||||
|
"dataType": "Float",
|
||||||
|
"arrayType": "Scalar",
|
||||||
|
"direction": "Input",
|
||||||
|
"description": "",
|
||||||
|
"isSelected": true
|
||||||
|
},
|
||||||
|
{
|
||||||
|
"name": "B",
|
||||||
|
"dataType": "Float",
|
||||||
|
"arrayType": "Scalar",
|
||||||
|
"direction": "Input",
|
||||||
|
"description": "",
|
||||||
|
"isSelected": true
|
||||||
|
},
|
||||||
|
{
|
||||||
|
"name": "IMP1",
|
||||||
|
"dataType": "Float",
|
||||||
|
"arrayType": "Scalar",
|
||||||
|
"direction": "Input",
|
||||||
|
"description": "",
|
||||||
|
"isSelected": true
|
||||||
|
},
|
||||||
|
{
|
||||||
|
"name": "SM",
|
||||||
|
"dataType": "Float",
|
||||||
|
"arrayType": "Scalar",
|
||||||
|
"direction": "Input",
|
||||||
|
"description": "",
|
||||||
|
"isSelected": true
|
||||||
|
},
|
||||||
|
{
|
||||||
|
"name": "EX",
|
||||||
|
"dataType": "Float",
|
||||||
|
"arrayType": "Scalar",
|
||||||
|
"direction": "Input",
|
||||||
|
"description": "",
|
||||||
|
"isSelected": true
|
||||||
|
},
|
||||||
|
{
|
||||||
|
"name": "KG",
|
||||||
|
"dataType": "Float",
|
||||||
|
"arrayType": "Scalar",
|
||||||
|
"direction": "Input",
|
||||||
|
"description": "",
|
||||||
|
"isSelected": true
|
||||||
|
},
|
||||||
|
{
|
||||||
|
"name": "KSS",
|
||||||
|
"dataType": "Float",
|
||||||
|
"arrayType": "Scalar",
|
||||||
|
"direction": "Input",
|
||||||
|
"description": "",
|
||||||
|
"isSelected": true
|
||||||
|
},
|
||||||
|
{
|
||||||
|
"name": "KKG",
|
||||||
|
"dataType": "Float",
|
||||||
|
"arrayType": "Scalar",
|
||||||
|
"direction": "Input",
|
||||||
|
"description": "",
|
||||||
|
"isSelected": true
|
||||||
|
},
|
||||||
|
{
|
||||||
|
"name": "KKSS",
|
||||||
|
"dataType": "Float",
|
||||||
|
"arrayType": "Scalar",
|
||||||
"direction": "Input",
|
"direction": "Input",
|
||||||
"description": "",
|
"description": "",
|
||||||
"isSelected": true
|
"isSelected": true
|
||||||
|
|||||||
BIN
Fortran/1.xlsx
BIN
Fortran/1.xlsx
Binary file not shown.
@ -1,60 +1,74 @@
|
|||||||
SUBROUTINE SUB_XAJMX( N ,& ! 单元出流数组大小 //输入变量
|
SUBROUTINE SUB_XAJMX( N ,& ! 单元出流数组大小 //输入变量
|
||||||
M ,& ! 无因次单位数组大小 //输入变量
|
M ,& ! 无因次单位数组大小 //输入变量
|
||||||
PAR ,& ! //输入变量
|
WUM ,& ! 上层张力水容量wum //输入变量
|
||||||
AREA ,& ! 单元面积 //输入变量
|
W1 ,& ! 下层张力水容量wl //输入变量
|
||||||
UH ,& ! 无因次单位线 //输入变量
|
WDM ,& ! 深层张力水容量wdm //输入变量
|
||||||
DT ,& ! 时段步长 //输入变量
|
KC ,& !蒸发能力折算系数KC //输入变量
|
||||||
P ,& ! 降雨系列 //输入变量
|
C ,& !深层蒸发系数c //输入变量
|
||||||
EP ,& ! 蒸发皿蒸发能力 //输入变量
|
B ,& !张力水蓄水容量系数b //输入变量
|
||||||
W ,& ! 土壤含水层 1.上层 2.下层 3.深层 //输入变量
|
IMP1 ,& !不透水面积比率imp1 //输入变量
|
||||||
FR ,& ! 初始产流面积 //输入变量
|
SM ,& !自由水蓄水容量sm //输入变量
|
||||||
S ,& ! 初始自由水深 //输入变量
|
EX ,& !自由水蓄水容量指数ex //输入变量
|
||||||
QRSS0 ,& ! 初始壤中流流量 //输入变量
|
KG ,& !地下水出流系数kg //输入变量
|
||||||
QRG0 ,& ! 初始地下水径流量 //输入变量
|
KSS ,& !壤中流出流系数kss //输入变量
|
||||||
QR ) ! 单元出流 //输出变量
|
KKG ,& !地下水出流系数kkg //输入变量
|
||||||
|
KKSS ,& !壤中流出流系数kkss //输入变量
|
||||||
|
AREA ,& ! 单元面积 //输入变量
|
||||||
|
UH ,& ! 无因次单位线 //输入变量
|
||||||
|
DT ,& ! 时段步长 //输入变量
|
||||||
|
P ,& ! 降雨系列 //输入变量
|
||||||
|
EP ,& ! 蒸发皿蒸发能力 //输入变量
|
||||||
|
W ,& ! 土壤含水层 1.上层 2.下层 3.深层 //输入变量
|
||||||
|
FR ,& ! 初始产流面积 //输入变量
|
||||||
|
S ,& ! 初始自由水深 //输入变量
|
||||||
|
QRSS0 ,& ! 初始壤中流流量 //输入变量
|
||||||
|
QRG0 ,& ! 初始地下水径流量 //输入变量
|
||||||
|
QR ) ! 单元出流 //输出变量
|
||||||
|
|
||||||
IMPLICIT NONE
|
IMPLICIT NONE
|
||||||
|
|
||||||
!///////////////////////////////////////变量声明//////////////////////////////////////////////
|
!///////////////////////////////////////变量声明//////////////////////////////////////////////
|
||||||
INTEGER::N !
|
INTEGER::N !
|
||||||
INTEGER::M !
|
INTEGER::M !
|
||||||
REAL::PAR(13) ! 1.上层张力水容量wum 2.下层张力水容量wl 3.深层张力水容量wdm
|
REAL::PAR(13) ! 1.上层张力水容量wum 2.下层张力水容量wl 3.深层张力水容量wdm
|
||||||
! 4.蒸发能力折算系数KC.深层蒸发系数c 6.张力水蓄水容量系数b
|
! 4.蒸发能力折算系数KC.深层蒸发系数c 6.张力水蓄水容量系数b
|
||||||
! 7.不透水面积比率imp1 8.自由水蓄水容量sm 9.自由水蓄水容量指数ex
|
! 7.不透水面积比率imp1 8.自由水蓄水容量sm 9.自由水蓄水容量指数ex
|
||||||
!10.地下水出流系数kg 11.壤中流出流系数kss 12.地下水出流系数kkg
|
!10.地下水出流系数kg 11.壤中流出流系数kss 12.地下水出流系数kkg
|
||||||
!13.壤中流出流系数kkss
|
!13.壤中流出流系数kkss
|
||||||
REAL::AREA ! 单元面积
|
REAL::AREA ! 单元面积
|
||||||
REAL::UH(M) ! 无因次单位线
|
REAL::UH(M) ! 无因次单位线
|
||||||
REAL::DT ! 时段步长
|
REAL::DT ! 时段步长
|
||||||
REAL::P(N) ! 降雨系列
|
REAL::P(N) ! 降雨系列
|
||||||
REAL::EP(N) ! 蒸发皿蒸发能力
|
REAL::EP(N) ! 蒸发皿蒸发能力
|
||||||
REAL::QR(N) ! 单元出流
|
REAL::QR(N) ! 单元出流
|
||||||
REAL::W(3) ! 土壤含水层 1.上层 2.下层 3.深层
|
REAL::W(3) ! 土壤含水层 1.上层 2.下层 3.深层
|
||||||
REAL::FR ! 初始产流面积
|
REAL::FR ! 初始产流面积
|
||||||
REAL::S ! 初始自由水深
|
REAL::S ! 初始自由水深
|
||||||
REAL::QRSS0 ! 初始壤中流流量
|
REAL::QRSS0 ! 初始壤中流流量
|
||||||
REAL::QRG0 ! 初始地下水径流量
|
REAL::QRG0 ! 初始地下水径流量
|
||||||
|
|
||||||
INTEGER::D
|
INTEGER::D
|
||||||
REAL::KSSD
|
REAL::KSSD
|
||||||
REAL::KGD
|
REAL::KGD
|
||||||
REAL::E(3)
|
REAL::E(3)
|
||||||
REAL::WM(3)
|
REAL::WUM !上层张力水容量wum
|
||||||
REAL::KC ! 蒸发能力折算系数
|
REAL::W1 !下层张力水容量wl
|
||||||
REAL::C ! 深层蒸发系数
|
REAL::WDM ! 深层张力水容量wdm
|
||||||
REAL::B ! 张力水蓄水容量系数
|
REAL::KC ! 蒸发能力折算系数
|
||||||
REAL::IMP1 ! 不透水面积比率
|
REAL::C ! 深层蒸发系数
|
||||||
REAL::SM ! 自由水蓄水容量
|
REAL::B ! 张力水蓄水容量系数
|
||||||
REAL::EX ! 自由水蓄水容量指数
|
REAL::IMP1 ! 不透水面积比率
|
||||||
REAL::KG ! 地下水出流系数
|
REAL::SM ! 自由水蓄水容量
|
||||||
REAL::KSS ! 壤中流出流系数
|
REAL::EX ! 自由水蓄水容量指数
|
||||||
REAL::KKG ! 地下水出流系数
|
REAL::KG ! 地下水出流系数
|
||||||
REAL::KKSS ! 壤中流出流系数
|
REAL::KSS ! 壤中流出流系数
|
||||||
|
REAL::KKG ! 地下水出流系数
|
||||||
|
REAL::KKSS ! 壤中流出流系数
|
||||||
|
|
||||||
!以下变量是原vb程序中未声明的变量
|
!以下变量是原vb程序中未声明的变量
|
||||||
INTEGER::I ! 计数器 //临时变量
|
INTEGER::I ! 计数器 //临时变量
|
||||||
INTEGER::J ! 计数器 //临时变量
|
INTEGER::J ! 计数器 //临时变量
|
||||||
INTEGER::ICHECK ! 判断计算时段长度是否合适的识别码
|
INTEGER::ICHECK ! 判断计算时段长度是否合适的识别码
|
||||||
INTEGER::NN
|
INTEGER::NN
|
||||||
|
|
||||||
REAL(KIND=8),PARAMETER::C5=5.000000000
|
REAL(KIND=8),PARAMETER::C5=5.000000000
|
||||||
@ -89,30 +103,13 @@ SUBROUTINE SUB_XAJMX( N ,& ! 单元出流数组大小 //输入变量
|
|||||||
REAL::QRG
|
REAL::QRG
|
||||||
REAL::QTR
|
REAL::QTR
|
||||||
|
|
||||||
!///////////////////////////////////////变量声明//////////////////////////////////////////////
|
!///////////////////////////////////////变量声明//////////////////////////////////////////////
|
||||||
|
|
||||||
!///////////////////////////////////////计算区域//////////////////////////////////////////////
|
!///////////////////////////////////////计算区域//////////////////////////////////////////////
|
||||||
|
|
||||||
! 赋值
|
! 赋值
|
||||||
ICHECK = 1
|
ICHECK = 1
|
||||||
|
|
||||||
DO I = 1, 3
|
|
||||||
|
|
||||||
WM(I) = PAR(I)
|
|
||||||
|
|
||||||
END DO
|
|
||||||
|
|
||||||
KC = PAR(4)
|
|
||||||
C = PAR(5)
|
|
||||||
B = PAR(6)
|
|
||||||
IMP1 = PAR(7)
|
|
||||||
SM = PAR(8)
|
|
||||||
EX = PAR(9)
|
|
||||||
KG = PAR(10)
|
|
||||||
KSS = PAR(11)
|
|
||||||
KKG = PAR(12)
|
|
||||||
KKSS = PAR(13)
|
|
||||||
|
|
||||||
DO I = 1, N
|
DO I = 1, N
|
||||||
|
|
||||||
QR(I) = 0.0
|
QR(I) = 0.0
|
||||||
@ -126,8 +123,7 @@ SUBROUTINE SUB_XAJMX( N ,& ! 单元出流数组大小 //输入变量
|
|||||||
D = 24 / DT
|
D = 24 / DT
|
||||||
CI = KKSS ** (1.0 / D)
|
CI = KKSS ** (1.0 / D)
|
||||||
CG = KKG ** (1.0 / D)
|
CG = KKG ** (1.0 / D)
|
||||||
KSSD = (1.0 - (1.0 - (KG + KSS)) ** (1.0 / D)) &
|
KSSD = (1.0 - (1.0 - (KG + KSS)) ** (1.0 / D)) / (1.0 + KG / KSS)
|
||||||
/ (1.0 + KG / KSS)
|
|
||||||
KGD = KSSD * KG / KSS
|
KGD = KSSD * KG / KSS
|
||||||
|
|
||||||
ELSE
|
ELSE
|
||||||
@ -145,11 +141,11 @@ SUBROUTINE SUB_XAJMX( N ,& ! 单元出流数组大小 //输入变量
|
|||||||
END IF
|
END IF
|
||||||
|
|
||||||
EP(I) = EP(I) * KC
|
EP(I) = EP(I) * KC
|
||||||
WM0 = WM(1) + WM(2) + WM(3)
|
WM0 = WUM + W1 + WDM
|
||||||
W0 = W(1) + W(2) + W(3)
|
W0 = W(1) + W(2) + W(3)
|
||||||
PE = P(I) - EP(I)
|
PE = P(I) - EP(I)
|
||||||
|
|
||||||
! 赋初值
|
! 赋初值
|
||||||
R = 0.0
|
R = 0.0
|
||||||
RIMP = 0.0
|
RIMP = 0.0
|
||||||
|
|
||||||
@ -163,15 +159,13 @@ SUBROUTINE SUB_XAJMX( N ,& ! 单元出流数组大小 //输入变量
|
|||||||
|
|
||||||
ELSE
|
ELSE
|
||||||
|
|
||||||
A = WMM * (1.0 - (1.0 - W0 / WM0) **&
|
A = WMM * (1.0 - (1.0 - W0 / WM0) ** (1.0 / (1.0 + B)))
|
||||||
(1.0 / (1.0 + B)))
|
|
||||||
|
|
||||||
END IF
|
END IF
|
||||||
|
|
||||||
IF ((PE + A) .LT. WMM) THEN
|
IF ((PE + A) .LT. WMM) THEN
|
||||||
|
|
||||||
R = PE - WM0 + W0 + WM0 * ((1.0 - (PE + A) /&
|
R = PE - WM0 + W0 + WM0 * ((1.0 - (PE + A) / WMM) ** (1.0 + B))
|
||||||
WMM) ** (1.0 + B))
|
|
||||||
|
|
||||||
ELSE
|
ELSE
|
||||||
|
|
||||||
@ -192,9 +186,9 @@ SUBROUTINE SUB_XAJMX( N ,& ! 单元出流数组大小 //输入变量
|
|||||||
ELSE
|
ELSE
|
||||||
|
|
||||||
E(1) = W(1) + P(I)
|
E(1) = W(1) + P(I)
|
||||||
E(2) = (EP(I) - E(1)) * W(2) / WM(2)
|
E(2) = (EP(I) - E(1)) * W(2) / W1
|
||||||
|
|
||||||
IF (W(2) .LE.( C * WM(2))) THEN
|
IF (W(2) .LE.( C * W1)) THEN
|
||||||
|
|
||||||
E(2) = C * (EP(I) - E(1))
|
E(2) = C * (EP(I) - E(1))
|
||||||
E(3) = 0.0
|
E(3) = 0.0
|
||||||
@ -218,15 +212,15 @@ SUBROUTINE SUB_XAJMX( N ,& ! 单元出流数组大小 //输入变量
|
|||||||
W(2) = W(2) - E(2)
|
W(2) = W(2) - E(2)
|
||||||
W(3) = W(3) - E(3)
|
W(3) = W(3) - E(3)
|
||||||
|
|
||||||
IF (W(1) .GT. WM(1)) THEN
|
IF (W(1) .GT. WUM) THEN
|
||||||
|
|
||||||
W(2) = W(1) - WM(1) + W(2)
|
W(2) = W(1) - WUM + W(2)
|
||||||
W(1) = WM(1)
|
W(1) = WUM
|
||||||
|
|
||||||
IF (W(2) .GT. WM(2)) THEN
|
IF (W(2) .GT. W1) THEN
|
||||||
|
|
||||||
W(3) = W(3) + W(2) - WM(2)
|
W(3) = W(3) + W(2) - W1
|
||||||
W(2) = WM(2)
|
W(2) = W1
|
||||||
|
|
||||||
END IF
|
END IF
|
||||||
|
|
||||||
@ -247,10 +241,9 @@ SUBROUTINE SUB_XAJMX( N ,& ! 单元出流数组大小 //输入变量
|
|||||||
S = X * S / FR
|
S = X * S / FR
|
||||||
SS = S
|
SS = S
|
||||||
Q = R / FR
|
Q = R / FR
|
||||||
NN = INT(Q / C5) + 1 ! 在vb程序中C5是双精度常数5
|
NN = INT(Q / C5) + 1 ! 在vb程序中C5是双精度常数5
|
||||||
Q = Q / NN
|
Q = Q / NN
|
||||||
KSSDD = (1.0 - (1.0 - (KGD + KSSD)) ** (1.0 / NN))&
|
KSSDD = (1.0 - (1.0 - (KGD + KSSD)) ** (1.0 / NN)) / (1.0 + KGD / KSSD)
|
||||||
/ (1.0 + KGD / KSSD)
|
|
||||||
KGDD = KSSDD * KGD / KSSD
|
KGDD = KSSDD * KGD / KSSD
|
||||||
RS = 0.0
|
RS = 0.0
|
||||||
RSS = 0.0
|
RSS = 0.0
|
||||||
@ -295,8 +288,7 @@ SUBROUTINE SUB_XAJMX( N ,& ! 单元出流数组大小 //输入变量
|
|||||||
|
|
||||||
ELSE IF ((Q + AU) .LT. SMMF) THEN
|
ELSE IF ((Q + AU) .LT. SMMF) THEN
|
||||||
|
|
||||||
RSD = (Q - SMF + S + SMF * (1.0 - (Q + AU)&
|
RSD = (Q - SMF + S + SMF * (1.0 - (Q + AU) / SMMF) ** (1.0 + EX)) * FR
|
||||||
/ SMMF) ** (1.0 + EX)) * FR
|
|
||||||
RSSD = (S + Q - RSD / FR) * KSSDD * FR
|
RSSD = (S + Q - RSD / FR) * KSSDD * FR
|
||||||
RGD = (S + Q - RSD / FR) * KGDD * FR
|
RGD = (S + Q - RSD / FR) * KGDD * FR
|
||||||
S = S + Q - (RSD + RSSD + RGD) / FR
|
S = S + Q - (RSD + RSSD + RGD) / FR
|
||||||
|
|||||||
@ -1,41 +1,67 @@
|
|||||||
subroutine XAJMX( N ,&
|
subroutine XAJMX( N ,& ! 单元出流数组大小 //输入变量
|
||||||
M ,&
|
M ,& ! 无因次单位数组大小 //输入变量
|
||||||
PAR ,&
|
WUM ,& ! 上层张力水容量wum //输入变量
|
||||||
AREA ,&
|
W1 ,& ! 下层张力水容量wl //输入变量
|
||||||
UH ,&
|
WDM ,& ! 深层张力水容量wdm //输入变量
|
||||||
DT ,&
|
KC ,& !蒸发能力折算系数KC //输入变量
|
||||||
P ,&
|
C ,& !深层蒸发系数c //输入变量
|
||||||
EP ,&
|
B ,& !张力水蓄水容量系数b //输入变量
|
||||||
W ,&
|
IMP1 ,& !不透水面积比率imp1 //输入变量
|
||||||
FR ,&
|
SM ,& !自由水蓄水容量sm //输入变量
|
||||||
S ,&
|
EX ,& !自由水蓄水容量指数ex //输入变量
|
||||||
QRSS0 ,&
|
KG ,& !地下水出流系数kg //输入变量
|
||||||
QRG0 ,&
|
KSS ,& !壤中流出流系数kss //输入变量
|
||||||
QR )&
|
KKG ,& !地下水出流系数kkg //输入变量
|
||||||
BIND(C, NAME="XAJMX")
|
KKSS ,& !壤中流出流系数kkss //输入变量
|
||||||
|
AREA ,& ! 单元面积 //输入变量
|
||||||
|
UH ,& ! 无因次单位线 //输入变量
|
||||||
|
DT ,& ! 时段步长 //输入变量
|
||||||
|
P ,& ! 降雨系列 //输入变量
|
||||||
|
EP ,& ! 蒸发皿蒸发能力 //输入变量
|
||||||
|
W ,& ! 土壤含水层 1.上层 2.下层 3.深层 //输入变量
|
||||||
|
FR ,& ! 初始产流面积 //输入变量
|
||||||
|
S ,& ! 初始自由水深 //输入变量
|
||||||
|
QRSS0 ,& ! 初始壤中流流量 //输入变量
|
||||||
|
QRG0 ,& ! 初始地下水径流量 //输入变量
|
||||||
|
QR )& ! 单元出流 //输出变量
|
||||||
|
BIND(C, NAME="XAJMX")
|
||||||
|
|
||||||
!DEC$ ATTRIBUTES DLLEXPORT::XAJMX
|
!DEC$ ATTRIBUTES DLLEXPORT::XAJMX
|
||||||
|
|
||||||
IMPLICIT NONE
|
IMPLICIT NONE
|
||||||
|
|
||||||
INTEGER::N ! 单元出流数组大小
|
INTEGER::N ! 单元出流数组大小
|
||||||
INTEGER::M ! 无因次单位数组大小
|
INTEGER::M ! 无因次单位数组大小
|
||||||
REAL::PAR(13) ! 1.上层张力水容量wum 2.下层张力水容量wl 3.深层张力水容量wdm
|
REAL::PAR(13) ! 1.上层张力水容量wum 2.下层张力水容量wl 3.深层张力水容量wdm
|
||||||
! 4.蒸发能力折算系数kc 5.深层蒸发系数c 6.张力水蓄水容量系数b
|
! 4.蒸发能力折算系数kc 5.深层蒸发系数c 6.张力水蓄水容量系数b
|
||||||
! 7.不透水面积比率imp1 8.自由水蓄水容量sm 9.自由水蓄水容量指数ex
|
! 7.不透水面积比率imp1 8.自由水蓄水容量sm 9.自由水蓄水容量指数ex
|
||||||
!10.地下水出流系数kg 11.壤中流出流系数kss 12.地下水出流系数kkg
|
!10.地下水出流系数kg 11.壤中流出流系数kss 12.地下水出流系数kkg
|
||||||
!13.壤中流出流系数kkss
|
!13.壤中流出流系数kkss
|
||||||
REAL::AREA ! 单元面积
|
REAL::AREA ! 单元面积
|
||||||
REAL::UH(M) ! 无因次单位线
|
REAL::UH(M) ! 无因次单位线
|
||||||
REAL::DT ! 时段步长
|
REAL::DT ! 时段步长
|
||||||
REAL::P(N) ! 降雨系列
|
REAL::P(N) ! 降雨系列
|
||||||
REAL::EP(N) ! 蒸发皿蒸发能力
|
REAL::EP(N) ! 蒸发皿蒸发能力
|
||||||
REAL::QR(N) ! 单元出流
|
REAL::QR(N) ! 单元出流
|
||||||
REAL::W(3) ! 土壤含水层 1.上层 2.下层 3.深层
|
REAL::W(3) ! 土壤含水层 1.上层 2.下层 3.深层
|
||||||
REAL::FR ! 初始产流面积
|
REAL::FR ! 初始产流面积
|
||||||
REAL::S ! 初始自由水深
|
REAL::S ! 初始自由水深
|
||||||
REAL::QRSS0 ! 初始壤中流流量
|
REAL::QRSS0 ! 初始壤中流流量
|
||||||
REAL::QRG0 ! 初始地下水径流量
|
REAL::QRG0 ! 初始地下水径流量
|
||||||
|
|
||||||
|
REAL::WUM !上层张力水容量wum
|
||||||
|
REAL::W1 !下层张力水容量wl
|
||||||
|
REAL::WDM ! 深层张力水容量wdm
|
||||||
|
REAL::KC ! 蒸发能力折算系数
|
||||||
|
REAL::C ! 深层蒸发系数
|
||||||
|
REAL::B ! 张力水蓄水容量系数
|
||||||
|
REAL::IMP1 ! 不透水面积比率
|
||||||
|
REAL::SM ! 自由水蓄水容量
|
||||||
|
REAL::EX ! 自由水蓄水容量指数
|
||||||
|
REAL::KG ! 地下水出流系数
|
||||||
|
REAL::KSS ! 壤中流出流系数
|
||||||
|
REAL::KKG ! 地下水出流系数
|
||||||
|
REAL::KKSS ! 壤中流出流系数
|
||||||
|
|
||||||
INTEGER::I
|
INTEGER::I
|
||||||
INTEGER::YEAR1
|
INTEGER::YEAR1
|
||||||
@ -46,19 +72,31 @@ subroutine XAJMX( N ,&
|
|||||||
INTEGER::ED
|
INTEGER::ED
|
||||||
INTEGER::EH
|
INTEGER::EH
|
||||||
|
|
||||||
CALL SUB_XAJMX( N ,& ! 单元出流数组大小 //输入变量
|
CALL SUB_XAJMX( N ,& ! 单元出流数组大小 //输入变量
|
||||||
M ,& ! 无因次单位数组大小 //输入变量
|
M ,& ! 无因次单位数组大小 //输入变量
|
||||||
PAR ,& ! //输入变量
|
WUM ,& ! 上层张力水容量wum //输入变量
|
||||||
AREA ,& ! 单元面积 //输入变量
|
W1 ,& ! 下层张力水容量wl //输入变量
|
||||||
UH ,& ! 无因次单位线 //输入变量
|
WDM ,& ! 深层张力水容量wdm //输入变量
|
||||||
DT ,& ! 时段步长 //输入变量
|
KC ,& !蒸发能力折算系数KC //输入变量
|
||||||
P ,& ! 降雨系列 //输入变量
|
C ,& !深层蒸发系数c //输入变量
|
||||||
EP ,& ! 蒸发皿蒸发能力 //输入变量
|
B ,& !张力水蓄水容量系数b //输入变量
|
||||||
W ,& ! 土壤含水层 1.上层 2.下层 3.深层 //输入变量
|
IMP1 ,& !不透水面积比率imp1 //输入变量
|
||||||
FR ,& ! 初始产流面积 //输入变量
|
SM ,& !自由水蓄水容量sm //输入变量
|
||||||
S ,& ! 初始自由水深 //输入变量
|
EX ,& !自由水蓄水容量指数ex //输入变量
|
||||||
QRSS0 ,& ! 初始壤中流流量 //输入变量
|
KG ,& !地下水出流系数kg //输入变量
|
||||||
QRG0 ,& ! 初始地下水径流量 //输入变量
|
KSS ,& !壤中流出流系数kss //输入变量
|
||||||
QR ) ! 单元出流 //输出变量
|
KKG ,& !地下水出流系数kkg //输入变量
|
||||||
|
KKSS ,& !壤中流出流系数kkss //输入变量
|
||||||
|
AREA ,& ! 单元面积 //输入变量
|
||||||
|
UH ,& ! 无因次单位线 //输入变量
|
||||||
|
DT ,& ! 时段步长 //输入变量
|
||||||
|
P ,& ! 降雨系列 //输入变量
|
||||||
|
EP ,& ! 蒸发皿蒸发能力 //输入变量
|
||||||
|
W ,& ! 土壤含水层 1.上层 2.下层 3.深层 //输入变量
|
||||||
|
FR ,& ! 初始产流面积 //输入变量
|
||||||
|
S ,& ! 初始自由水深 //输入变量
|
||||||
|
QRSS0 ,& ! 初始壤中流流量 //输入变量
|
||||||
|
QRG0 ,& ! 初始地下水径流量 //输入变量
|
||||||
|
QR ) ! 单元出流 //输出变量
|
||||||
|
|
||||||
end subroutine XAJMX
|
end subroutine XAJMX
|
||||||
|
|||||||
Loading…
Reference in New Issue
Block a user