Arcswat运行
学习arcswat几个月,现对处理数据和运行模型过程中碰到的一些问题和对应的解决方法作记录。
1.hru矢量数据的生成
在用landuse soil slope划分hru的界面,选择create hrus,就能生成矢量的hru。
2.hru划分
从文献中得, 日尺度的swat模拟对hru划分的个数敏感,所以要多次尝试得到平衡计算效率和结果准确度的hru个数,日尺度的比例时5%/5%/5%。
月尺度的模拟对于hru不时特别敏感,所以如果后续在swatcup中矫正的结果不好的话(针对P-factor, R-factor, NS,PIBAs的值),大多数情况不是hru划分的问题,而是矫正参数的选择以及参数范围的划定问题。
对于月尺度的hru划分来说,25%/25%/25%或者20%/20%/20%都是可以接受的划分。
3.swat running
在这个过程中报错的话,一半是在txtin文件夹里边wnd或者hmd文件中,有些行多出几个字符,解决方法是把数值前有-符号的值,后边紧跟的0删掉,保证每一行的字符一致。
4.swatcup 矫正
对于软件参数的设置和说明,可以上网查找manual和Google group的swat cup小组,有详细的说明和报错讨论以及回答。
对于矫正结果的说明,最重要其实是子流域的划分,尽可能详细,然后是hru的划分,主要是提高计算效率,过高的hru会大大增加计算负担。在初期测试矫正结果和参数的时候,可以设定较小的iteration次数,50次左右就能得到有一定参考价值的statistic summary for each observed station.
swatcup在新建的时候,是激活了observe.rch module. 但是如果需要矫正在hru or sub,则需要额外在executable files 里边勾选对应的module,否则swat时不会执行到hru or sub.