基于 EXCEL 2003时段雨强的快速计算

For j = 1 To lxkj - T If It ( j) > It( j + 1) Then It( j + 1) = It( j) End If Next j Cells(M , N ). Offset( , 3). Value = Str( T) + "分钟雨强 " Cells(M , N ). Offset(1, 3). Value = It( lxkj - T + 1) Cells(M , N ). Offset(1, 3). Select End If E lse M sgBox "没有符合条件的最大时段雨强 ,请检查数据选择 域 ! " , vbExclamation, "警告 ! " End If End Sub
0. 17 0. 12 0. 09 0. 08 0. 07 0. 06
0. 17 0. 12 0. 09 0. 08 0. 07 0. 06
0. 81 0. 42 0. 29 0. 24 0. 20 0. 17
0. 81 0. 42 0. 29 0. 23 0. 20 0. 17
1. 49 1. 30 1. 06 0. 87 0. 72 0. 61
新干县水土保持监督管理的 做法及存在问题与对策
邓珠根 1 ,傅绍娟 2
(1. 新干县水务局 ,江西 新干 331300; 2. 新干县黄泥埠水库管理局 ,江西 新干 331300)
[关键词 ] 水土流失 ;监督管理 ;开发建设项目 ;新干县 [摘 要 ] 介绍了新干县开发建设项目所造成水土流失的现状及危害 ,总结了该县水土保持监督管理的做法 :一是加大 宣传力度 ,增强全民的水保法制观念 ;二是建立和完善水保配套法规体系 ;三是建立执法体系 ,全面提升水保执法队伍素 质 ;四是紧抓“三权 ”,强化监管力度 。针对存在的主要问题提出了加强部门协作和宣传 、加大对违法行为的查处力度等ቤተ መጻሕፍቲ ባይዱ建议 。 [中图分类号 ] S157 [文献标识码 ] C [文章编号 ] 1000 - 0941 (2010) 06 - 0060 - 02
2 程序设计
基于以上数据记录特征 ,采取如下三步实现程序计算 : ① 将记录时间间隔 > 1 m in的数据拆分为间隔时间为 1 m in的数 据 。 ②遍历所有符合特定时段内的雨强 。 ③筛选出特定时段 内的雨强最大值 。这样设计有利于计算特定时间段 T的雨强 , 同时也可避免在遍历过程中遗漏某些特定时段条件的雨强值 , 从而快速计算出精确的时段雨强 。 2. 1 拆分算法
以 It ( ) 数组为基础 , 从数组第一个数开始 , 判断其下一个 数与其大小 ,如果前者大于后者则将该两数进行替换 , 然后再 把第二个数同第三个数进行比较 , 依此类推 , 最后数组中的最 大值即被传递到 It ( ) 数组的最后一位 。
崔普伟等 :基于 EXCEL 2003时段雨强的快速计算
I10
I20
I30
I40
I50
I60
1. 71 1. 48 1. 26 1. 05 0. 87 0. 74
1. 71 1. 48 1. 26 1. 05 0. 87 0. 74
1. 28 1. 00 0. 90 0. 77 0. 68 0. 63
1. 28 1. 00 0. 90 0. 77 0. 68 0. 63
3. 江西省水土保持科学研究所 ,江西 南昌 330029)
[关键词 ] 时段雨强 ;计算方法 ; EXCEL 2003; VBA 代码 [摘 要 ] 时段雨强的计算是土壤流失量计算中的重要内容 。在分析陕北团山沟径流场 1961—1969年实测降雨过程数 据的基础上 ,以时间间隔 1 m in为最小计算单元 ,将计算程序概化为三步 :将记录时间间隔 > 1 m in的数据拆分为间隔时 间为 1 m in的数据 ;遍历所有符合特定时段内的雨强 ;筛选出特定时段内雨强的最大值 。计算过程在 EXCEL 2003中用 VBA语言编码实现 ,并对其计算结果和手工计算结果作了对比 ,对比结果表明程序计算是准确的 ,可用于土壤流失量计 算中的雨强快速计算 。 [中图分类号 ] S157. 1 [文献标识码 ] A [文章编号 ] 1000 - 0941 (2010) 06 - 0058 - 03
(陈晓安同志为本研究提供了手工计算的雨强结果值 ,为 验证本算法的可靠性提供了帮助 ,特此致谢 ! )
[参考文献 ]
[ 1 ] 王万忠. 黄土地区降雨特性与土壤流失关系的研究 [ J ]. 水 土保持通报 , 1983, 3 (4) : 7 - 13.
·60·
中国水土保持 SW CC 2010年第 6期
·58·
中国水土保持 SW CC 2010年第 6期
基于 EXCEL 2003时段雨强的快速计算
崔普伟 1, 2 ,蔡强国 2 ,张利超 3 ,郑明国 2
(1. 华中农业大学 资源与环境学院 ,湖北 武汉 430070; 2. 中国科学院 地理科学与资源研究所 陆地水循环与地表过程重点实验室 ,北京 100101;
输入时间 T (m in) ,以数组 lxk j( ) 为基础 ,从数组中第一个 位置开始计算其 T时段内的平均雨强 , 保存到数组 It ( ) , 然后 再从第二个位置开始计算其 T时段内的平均雨强并保存到数组 It ( ) ,依此类推遍历数组中所有连续 T 时段的雨强 , 依次存储 到数组 It ( ) 中 。 2. 3 筛选算法
时段雨强的计算是土壤流失量计算中的重要内容 ,我国学 者 [1 - 5 ]在研究不同地区的降雨侵蚀力时 ,使用最大 5、10、20、 30、60 m in等不同时段雨强及它们与降雨总动能的乘积的组合 形式来研究降雨与土壤侵蚀的关系 ,这些研究体现了时段雨强 在土壤侵蚀预报模型中的重要地位 。不论是实际应用还是科 学研究 ,时段雨强的计算都是一个很重要的基础工作 ,传统的 手工计算往往会产生较大的误差且计算费时费工 。鉴于此 ,我 们以陕 北 团 山 沟 径 流 场 的 实 测 降 雨 过 程 资 料 为 基 础 (见 《1961—1969年黄河流域子洲径流实验站水文实验资料 》) ,考 虑到数据自身的记录情况 ,对时段雨强的快速计算进行了研 究。
1 降雨数据特征
时段雨强是指一场降雨过程中某连续时段内最大的平均 雨强 ,它一般发生在降雨过程的雨峰附近 。如最大 10 m in雨强 是指一场降雨中连续 10 m in累计雨量最大的那个时间段所对 应的平均雨强 。
对径流场 1961—1969年的自记雨量计数据进行分析 ,可知 记录的数据有如下一些特征 : ①降雨过程数据包括时间和时段 雨量两个数据项 ,一般记录得都比较详细 。 ②记录间隔时间不 等 ,雨强越大间隔时间越小 ,最短的时间间隔为 1 m in,最长的 时间间隔为 293 m in (间隔雨强为 0. 003 8 mm /m in) ,一般间隔 时间 > 30 m in的其间隔雨强 < 0. 01 mm /m in。 ③短历时降雨时 间间隔通常比长历时降雨短 ,如 > 0. 1 mm /m in的雨强 ,在短历 时降雨中记录的间隔时间为 1 m in,而在长历时降雨中却达到 十几至几十分钟 。 ④次降雨记录总时间从 28～2 530 m in不等 。 ⑤短历时降雨中 ,雨强 ≥0. 1 mm /m in时 ,记录的间隔时间通常 为 1～3 m in;雨强为 0. 1～0. 01 mm /m in时 ,间隔时间为 3 ～30
Op tion Base 1 Sub cui( ) D im M , N , Zongshu A s Integer M = Selection. Cells(1, 1). row N = Selection. Cells(1, 1). Column Zongshu = Selection. Rows. count D im lx1m in ( ) A s Single D im lxkj A s Integer For i = 1 To Zongshu D im quzhengI A s Integer quzhengI = Round (Cells( i + M - 1, N ) , 0) If quzhengI = 1 Then lxkj = lxkj + 1 ReD im Preserve lx1m in ( lxkj) lx1m in ( lxkj) = Cells( i + M - 1, N ). Offset( , 1). Value E lse lxkj = lxkj + quzhengI ReD im Preserve lx1m in ( lxkj) For j = 1 To quzhengI lx1m in ( lxkj - j + 1) = Cells ( i + M - 1, N ). Offset ( , 1). Value Next j End If Next i D im It( ) A s Single D im T A s Integer T = Val ( InputBox ( "请输入时段长 ,要求为正整数 ! " , "自 定义时段长 T" ) ) If lxkj - T > = 0 Then ReD im It( lxkj - T + 1) For i = 1 To lxkj - T + 1 puwei = 0 For j = 1 To T puwei = puwei + lx1m in ( i + j - 1) Next j It( i) = puwei / T Next i If lxkj - T = 0 Then Cells(M , N ). Offset(1, 3). Value = It(1) Cells(M , N ). Offset( , 3). Value = Str( T) + "分钟雨强 " E lse