Category: DHSVM

DHSVM结果文件处理程序

DHSVM的一些水文成分(比如降雨蒸散发等)是输出在Mass.Balance这个文件中,按照模拟的时间步长输出的。

如下:

但是实际我们可能关系的是在daily,monthly or yearly的尺度上的结果,方便分析比较。为此我写了一个MATLAB程序能把这个文件输出到三个不同的文件,分别代表daily, montly, yearly的结果。

下面是代码的头文件,需要用户自定义一些:

%% This Code is used to sumarized the default output of DHSVM (Mass.Balance)

%  I[……]

Read more

VIC模型对气象数据进行时间降尺度

DHSVM模型需要subdaily的气象数据,而一般气象站只有daily的数据。为此,我将VIC模型的自带的气象数据的Meteorological Forcing Disaggregation模块,比较粗略的抠出来了,使其可以方便的使用于气象数据时间降尺度,从而满足DHSVM模型模型的数据要求。否则,直接采用原来的VIC模型,对于新手,比较复杂。 抠出来的代码我已经编译了,然后放在这里共享。程序输入为daily的气象台站数据,然后输出subdaily数据,格式完全和DHSVM对应。因此可以直接驱动DHSVM数据。 关于这个程序,我写了个简单说明: #**********************[……]

Read more

DHSVM Subdailiy to Daily streamflow

Update 2016/11/30

今天我又重新写了subdaily to monthly的程序,现在放在这里一并和大家分享。

截个图:

DHSVM_stream.rar

DHSVM模型流域出口的河川径流的格式(3小时为例):

现在要将Subdaily的streamflow转化为daily,我写了一个C的程序来实现,先贴结果.

源代码如下:

#include <stdio.h>

#include <stdlib.h>

#include <string.h>

#define M[……]

Read more

圣维南方程组

1. 简介

描述水道和其他具有自由表面的浅水体中渐变不恒定水流运动规律的偏微分方程组。由反映质量守恒律的连续方程和反映动量守恒律的运动方程组成。1871年由法国科学家A.J.C.B.de圣维南提出,故名。
一百多年来,虽然为了考虑更多的因素和实际应用方便对它的基本假定作了某些简化或改进,产生出多种不同的表达形式,但其实质没有变化。主要进展表现在求解方法的改进和创新。1877年法国工程师克莱茨提出了瞬态法。1938年苏联С.А.赫里斯季安诺维奇提出另一类解法──特征线法。但均因计算量较大,不得不进行各种简化处理,使实际应用受到限制。自50年代以来,随着电子计算机的普及,研究和提出了一[……]

Read more

河道汇流连续方程推导

现在任意取长为Δx的微小河段,图中Q(x)为河段入口处(x)处的流量,Q(x+Δx)为出口(x+Δx)处的流量,q为单位长度的侧向径流量,包括该单元净雨产生的地表径流和从其它单位汇入的壤中流。

则根据质量守恒定律,可以得到下面的推导:
这就是著名的圣维南方程组的中的连续方程。该式表示河道非恒定流中,流量沿程的变化率与过水断面随时间变化率之和等于单位长度的侧向汇流流量。
在明渠恒定流时,由于水面的坡度不断随时间变化,故断面上的水文-流量关系成单值关系,即一个水位对应一个固定的流量。而对于非恒定流,则同一水位下明渠断面上的水文与流量之间呈多值关系。例如,涨水过程中[……]

Read more

曼宁公式

曼宁公式是一个估测液体在开放管道(即明渠流)或非满管流中平均速度的经验公式。开放管道中的液体是因重力而流动。曼宁公式最早是由法国工程师Philippe Gauckler在1867年提出,在1890年爱尔兰工程师罗伯特·曼宁也提出相同的公式。

Manning formula:
其中:
  • V为截面的平均速度(L/T,ft/s或m/s)
  • k是转换系数,因次L1/3/T,国际标准制为1 m1/3/s,若用美制,则是1.4859 ft1/3/s(注:(1 m)1/3/s = (3.2808399 ft) 1/3/s = 1.4859 ft1/3/s)();
  • n为Gauckler[……]

    Read more

单位线法

1. 概述:

在给定的流域下,单位时段内均匀分布的单位地面(直接)净雨量,在流域出口断面形成的地面(直接)径流过程线,称为单位线。单位线主要用于推求洪峰流量的过程。单位净雨量通常取雨深10毫米;单位历时可以是一个时段如1小时、3小时、6小时等,也可以是瞬时,即净雨历时趋于无限小的情况。相应于前者的单位线称为时段单位线,相应于后者的称为瞬时单位线。[百度百科]http://baike.baidu.com/view/1572309.htm
2. 基本假设
将流域看线性系统,即流域满足倍加性和叠加性
(1)倍比性
 如果单位时段内的净雨不是一个单位而是k个单位,则形成的流量过程是单[……]

Read more