#VIC汇流程序C++版 本程序用于VIC(Variable Infiltration Capacity)水文模型运行结果的汇流演算,从fortran版程序改写,改了一点东西。未来可能会给它搞搞并行化。 相比原程序去掉了XMASK文件。 修改后的输入参数文件样例:
# 流向文件路径
/项目位置/direc.v
# 流速文件路径,可设为F后带上均一值。
F
1.50
# 水力扩散文件路径,可设为F后带上均一值。
F
800
# 产流面积比例文件路径,可设为F后带上均一值。
/项目位置/fract.v
# 站点参数文件路径。
/项目位置/stnloc.v
# VIC输出文件路径,其文件名经纬度保留小数位数,以及VIC输出文件中产流数据所在列数。
/VIC输出路径/fluxes_
4
6
# 汇流文件输出路径。结尾记得带上“/”号。
/汇流文件输出路径/
# 进行汇流计算的开始和结束时间,以及实际输出到文件的数据开始和结束时间。
1949 01 1951 12
1949 01 1951 12
# 坡面汇流文件路径,可使用sethekin样例自带的。
/坡面汇流文件路径/UH.all
其他格式不变。
##使用方式
对于流向文件里方向的编码方式有两种,一种为用数字18编码的原版方式,另一种则为ArcInfo的用2的07次幂编码的方式,直接从ArcGIS中导出的流向栅格ASCII文件便是这种方式。
- 原版的方式为: 北-1 东北-2 东-3 东南-4 南-5 西南-6 西-7 西北-8
- ArcInfo的方式为: 东-1 东南-2 南-4 西南-8 西-16 西北-32 北-64 东北-128
本程序经行默认编译得到的方式为原版方式。
若你要使用的流向文件的编码方式不同则在使用make
的时候应输入相应命令。
在终端中打开源代码目录:cd .../CRout/src
编译:make
若你的流向文件为ArcInfo方式编码,则输入make arcinfo
经行编译。
运行:./crout
。直接输入可查看编译所得的程序使用的流向编码方式。
- 注意:流向文件应用不大于0的整数标示无数据的网格。