加载中…加载中…加载中…20120204，周六
零点时分，睡前站在阳台窗边，抬望夜空，山坡外中信的灯光已熄灭，夜空灰暗，忽见一条淡淡的带状无序的云白物体在闪动，东头从正面山坡，西端没入侨源阁后中信以远，不敢相信，等了眼细看，图案晃来晃去，翻来覆去，上下左右飞旋，纷乱又不断重复，略有点像急速流动的庞大的无穷大符号，惊奇中在阳台不同位置约凝视了三、五分钟，匆忙回房间看时钟的时间是后半夜零点30分，再出来时好像还晃了几个来回，忽然就不动了，静止的稀疏云团正是方才纷乱的原型，再转眼间就消失了，连天空的薄云都消失了。如果刚才看到的不明飞旋体是团状并更亮一些，估计会怀疑看到飞碟了，但从最后静止的图案看应该是白云而已，可是一片稀疏云团怎么可能纷飞旋辗转转得如此之快，难道是幻觉是视觉错乱，貌似梦境，其实非也，不可思议。
这个世界有多少不可思议，在平淡的夜空，在平淡的人生。
流连周六在家休习，明日全校开学典礼。
PS：后来，纽约时间2月4日，就像这片无穷旋转的流云，林书豪横空出世；不同之处在于，前者，无以证实，没有后来；后者，刹那间惊艳了整个NBA，一而再，再而三（2月10日PS）。
| 分类： 2012 |
20120121
傍晚，弟弟一家从南昌驾车至山居；
20120122 廿九 除夕
清晨，小雨，步行至东站，7点母亲和安江到；上午，姐和然乘机来，于燕塘站接；
20120123 初一
二老与弟、弟媳去香港游；
20120126 初四
与姐、弟、弟媳及黄阿姨等中午到汕尾，午餐后到百仕达入住，下午游遮浪；
20120127 初五
弟弟一家早餐后返程；与姐乘中午乘加班车返；
20120129 初七
下午，与母亲、姐游广州塔沿岸；
20120130 初八
上班，补一红包(200.0)；姐和然乘下午列车返程；
20120131 初九
上班，又补一红包(200.0)；
下了N日（5-6天）的雨，湿气膨胀；
1月16日 周一 阴晴 下午，参加高三年级家长会，为一锁头，浑然迟到了:(打的赶到礼堂，只听到最后一位老师的讲话；这个寒假高三学生都要在校学习，中间也只有7天休息，学校要求不可放松；从区一摸到市调研，流连连续3次在班上进入年级理科前50名，但也只是在重本的边缘；邢老师安排流连在班级家长会上介绍了经验；
1月20日 周五 阴
下午来办公室，已贴封条，揭开一半进来；应该是兔年最后一次上网；山居的第一个春节，没有网络，没有电视，有的是家人团聚，母亲，姐姐，弟弟，，，十余亲人，其乐融融；
龙年如水而至，恭祝各位新春快乐！
一觉连梦醒来，听到滴滴答答的声音，天光朦胧，在下雨，新年第一场雨；
柏木餐桌+4张椅子+鞋柜（5100），雨中送到，桌面椅面布满水滴，中午；
和W&Y到省厅2A，见莫总；出来时下小雨，下午；
20120106 阴有小雨 公司年饭于佰好酒楼；
20120108 阴间晴 全院年饭于一方渔家（100）；
20120109 晴 蓝&甘来山居安装4扇黑胡桃+竹绘铝合金玻璃门（2550）；
今天，是一个属于十八岁的成人礼，
十七岁的你也站在这里，果然这样！
当你的小伙伴迈着蹒跚的脚步注视你读书时候，
你已学会了特立独行。
你的名字就是你的座右铭，
汉字的简洁和博大，原本如此！
刘连——流连当下，执着现在，
鲁诣——纵然愚鲁亦有造诣，无畏前行直达目的。
你的脚印从小到大，背着考试的书包，
有失落，也有沉迷，尽管这般！
你喜欢独自思考，不断进步，
大概你已明白：流连其所在，鲁诣必争先。
今天，同样属于十七岁的你，
我早就说过你是神奇的鲁诣，注定如此！
你母亲的祖国令你独一无二，
我们的传奇才刚刚开始。
花城新中轴
燕岭公园畔
地铁沉双轨
绿道出海心
水是生命之源、生产之要、生态之基。除水害、兴水利，治水、理水，历来是富民安邦的大事，事关国计民生、经济发展、社会进步。
建国60年来，特别是改革开放以来，汕尾市的水利建设取得了辉煌成就。1988年（建市）~2009年间，全市以年均超过1亿元的投入，兴建、维修加固了大批水利工程，区域内逐渐形成防洪（潮）、治涝、节水、灌溉、供水、发电、环保、水土保持、旅游、渔业及航运等多功能的水利体系，在保障人民生命财产安全和国民经济的发展中发挥了巨大的作用。
洪潮涝旱灾害频繁仍然是汕尾市的心腹大患，农田水利建设滞后一直是影响农业稳定发展和粮食安全的严重硬伤，水利设施薄弱是基础设施的明显短板，水资源供需矛盾深化是可持续发展的主要瓶颈。
新形势下，水利是保障城乡居民饮水安全的基本生存体系，是现代农业建设不可或缺的首要条件，是经济社会发展不可替代的基础支撑，是生态环境改善不可分割的保障系统，具有很强的公益性、基础性、战略性。加快汕尾市流域水利改革发展并进行规划，不仅关系到防洪安全、供水安全、粮食安全，而且关系到经济安全、生态安全、社会安全；不仅事关农业农村发展，而且事关经济社会发展全局。
本次规划修编工作是贯彻落实国务院水利部《关于开展流域综合规划修编工作的意见》（国办发[2007]44号），并根据省水利厅有关要求，针对汕尾市当前和今后水利建设中的紧要问题，进行深入系统的综合调查研究，坚持以人为本、人与自然和谐相处，致力全面规划、统筹协调、开发与保护并重，突出综合治理、因地制宜、远近结合，提出了《广东省汕尾市流域综合规划修编报告（附图、附表）》。这一修编成果，覆盖近期、远期以及可以展望的未来，对于汕尾市的流域治理，水资源开发、利用、配置、节约和保护，促进水利可持续发展，努力走出一条具有汕尾特色的民生水利道路，应当起到重要的指导作用。
本次规划方针概括为十六个字：兴利除害，疏蓄结合，调丰保枯，循环利用。我们不仅要关注防洪安全，同时要关心供水安全、粮食安全，实现经济社会及生态安全。疏、蓄的出发点，应力求符合水的自然属性和人对于水的依存性，不应该简单地求泄，将宝贵的水资源泄于无形。在疏的过程中找到蓄的所在，疏而不漏，是疏蓄结合的要义。水有丰枯，丰时能够减少弃泄，才会有枯时不竭的源头，疏蓄结合的主要目的之一是供水，即为调丰保枯服务。要大力改善蓄的能力，而大型水库的有效利用和开发建设是提高蓄水能力的重要举措和保障。水来水去，绿色的水资源和水能蕴藏其间，人类、自然和环境造化其间。水是可持续发展的重要因素，在水的良性循环过程中利用水，让水资源和水力资源发挥更大的效益，实现人水和谐，是流域综合规划应尽的职责。
规划目标概括为四丰：河丰，海丰，陆丰，人丰。体现了让汕尾市江河丰沛、湖海丰饶、陆地丰茂、人民丰盛的理念。
流域综合规划修编涉及面广，十分复杂。市水利局认真组织并委托广东华南水电高新技术开发有限公司承担，有关技术人员夜以继日地工作，付出了大量的心血，凝结了广泛的智慧，四易其稿，历时一年半完成。
规划成果经过市水利局组织的内部评审，经过省水利厅特邀7人专家组的审查，并根据有关专家以及相关部门提出的意见与建议，反复进行了认真细致的修改。修编工作得到省、市水利有关领导、专家和技术人员的关心和支持，在此深表感谢。特别感谢张火金、范锦春、黎开志、刘伟成、李学展、汪利民、叶云晋、吴小明、林振勋、刘霞、骆志雄（排名不分先后）等专家的指导。特别感谢罗水庚、蔡家和、陈军强、林乔生等合作双方的领导及专家创造了优良的工作环境并提供了大量的基础资料。还要感谢网络时代的信息援助，特别感谢十二棵橡树之梦、葛维亚先生、科学中国人等网友及家人超越专业和专家界线的理解和支持。正是由于上述帮助，使修编工作充满挑战的快乐。
编制这样的综合性规划报告书，是一件浩繁的工作。能否于浩繁的数字，洞见无限可能的未来？修编人员以务实性与前瞻性兼修的规划思路从事这次综合规划修编，为了力求全面、准确、严肃、科学、合理、可行，有限的修编人员已尽了莫大的努力，但由于技术水平、经费、时间等条件限制，难免有疏漏和谬误之处。今后应根据汕尾市社会经济的迅猛发展对水利提出的更高需求及时进行补充和调整，使之更臻完善。
修编者
2011年3月1日
标签：
杂谈 |
分类： DigitWater.For |
The water modeling，a brief history of time
abstract:
One waterhole,like a water,like water clock,like
black holes,will swallow me,could be another world,for
myself,haha,and for you or not for you.
题记：
一个水洞，像一滴水，像水时钟，像黑洞，会吞没我，成为另一个世界，为我自己，哈哈，也为你或者不为你。
keywords: waterhole,digitwater,LSModel,理水诣
!
!
!
!
WYH_Mod.F90为舞阳河项目增加 20081026 1:19
!
变量、常数和数组封装模块
module
WYH_Mod
implicit none
!
!
!
因某源水调度项目增加,但无必然关系 20110120
character(len=201) :: Filename104 ! 测站关系信息表
!
输入文件名，...\GN_SM_SUB\104ST_STRL_B.CSV ! 20110214
21:35
character(len=201) :: Filename137 ! 水厂及取水口信息表
!
输入文件名，...\GN_SM_SUB\137ST_WFCT_B.CSV ! 20110209
22:21
character(len=201) :: Filename206 ! 水库水情表
!
输入文件名，...\GN_SM_SUB\206ST_RSVR_R.CSV ! 20110121
15:42
character(len=201) :: Filename207 ! 堰闸水情表
!
输入文件名，...\GN_SM_SUB\207ST_WASG_R.CSV ! 20110121
15:41
!character(len=201) :: Filename208 ! 闸门启闭情况表
!
输入文件名，...\GN_SM_SUB\208ST_GATE_R.CSV ! 20110120
9:50
character(len=201) :: Filename209 ! 泵站水情表 ! 踏遍莲花少一人 20110121
!
输入文件名，...\GN_SM_SUB\209ST_PUMP_R.CSV
character(len=201) :: Filename304 ! 堰闸水情预报表 20110120 13:39
! 输出文件名，...\GN_FD_SUB\304ST_WAS_F.CSV
character(len=201) :: Filename305 ! 闸门启闭情况预报表 20110121 13:42
! 输出文件名，...\GN_FD_SUB\305ST_GATE_F.CSV
character(len=201) :: Filename313 ! 泵站水情预报表(未输出) 20110121
14:30
! 输出文件名，...\GN_FD_SUB\305ST_PUMP_F.CSV
!
!character(len=2),allocatable :: PDCHCD(:,:) !
表209中的引排特征码，1引水,0或O排水；:,: M_JH,NT ! 20110121 11:01
!
调节点指定流量也存储在CHLEN(:) 20110125 12:32
!
integer ::
M_DP ! 节点图中的泵站总数 ! 20110123 22:16
integer ::
M_DD ! 节点图中的堰闸总数 ! 20110123 23:11
integer ::
iDP,iDD ! 泵站、堰闸节点序号 20110125 15:07
integer,allocatable :: ID_DP(:) ! 泵站节点ID号；: M_DP ! 2110125
15:33
integer,allocatable :: IDPhash(:) ! 泵站哈稀数；: N_ND
integer,allocatable :: ID_DD(:) ! 堰闸节点ID号；: M_DD ! 未用
20110125
integer,allocatable :: IDDhash(:) ! 泵站哈稀数；: N_ND
integer ::
jND,jJH,jDP,jNR ! 节点序号,与iND,iJH等区别 20110125 15:45
integer,allocatable :: L_DP(:) ! 泵站优先级；: M_DP ! 2110125 16:07
!
泵站运行优先级代码,1优先级1,2优先级2,3优先级3,,,0无优先级，针对局部工程，实为就近工程优先 20101125
11:33
! 泵站优先级首先存储在CHLEN(:) 20110125 12:24
integer ::
L_DPmax ! 泵站最大优先级数 20110125 16:03
integer ::
iTJ,M_TJ ! 调节点序号,调节点总数 20110125 18:01
integer,allocatable :: ID_TJ(:) ! 调节节点ID号；: M_TJ ! 2110125
18:08
integer,allocatable :: ITJhash(:) ! 调节点哈稀数；: N_ND
real,allocatable :: Q_TJ(:) ! 调节点指定流量；: M_TJ ! 20100125 17:05
!
real,allocatable :: PMPQ(:) !
抽水流量；: M_DP ! 20110207 15:29
real,allocatable :: PMCL(:) !
抽水含氯度；: M_DP ! 20110207 16:01
! 注销以下3行 20110216 0:16
!real,allocatable ::
PPUPZ(:),PPDWZ(:) ! 泵站上游、下游控制水位(不变)；: M_DP ! 20110203 14:50
!character(len=2),allocatable
:: PDCHCD(:) ! 表209中的引排特征码(不变)，1引水,0或O排水；:,: M_DP ! 20110207
13:55
!integer,allocatable ::
IPRCD(:) ! 泵站控制水库编码；: M_DP ! 2110213 1:38
!
! 表130 泵站(水库电站、径流水电站)发电用 20110126 10:29
integer :: iTP ! 类型,程序中有3类：1水库电站,2泵站,3径流电站
20110126 11:46
integer
:: mTP ! 站数
integer
:: nMGN ! 机组号数 20110126 12:26
integer,allocatable
:: M_TP(:) ! 分类站数；: iTP 20110126
18:49
integer,allocatable
:: M_MGN(:,:)
! 机组号数,一个机组号对应表130的一行数据；:,: iTP,mTP 20110126 11:41
real,allocatable
:: PTDSQ(:,:) ! 机组总设计流量；:,: iTP,mTP 20110126 13:05
integer,allocatable
:: PMCHCNT(:,:,:) ! (泵站)装机台数；:,:,: iTP,mTP,nMGN
real,allocatable
:: PSGRTPWR(:,:,:) ! 单机额定功率；:,:,: iTP,mTP,nMGN
real,allocatable
:: PSGDSQ(:,:,:) ! 单机设计流量；:,:,: iTP,mTP,nMGN
!
表130 原为水库发电用,可兼容；老的程序应修改 20110126 11:11
integer
:: M_PS ! 节点图中的径流式水电站总数；也隐含在M_TP数组中；暂无 ! 20110126 12:05
!
! 表131 堰闸补充 20110214
22:37
real,allocatable :: SWDSQ(:) !
水库堰闸单孔设计流量；: M_NR
real,allocatable :: TWDSQ(:) !
水库堰闸总设计流量；: M_NR
integer :: M_DD !
节点图中的堰闸总数(水库隐含的堰闸不在此中) ! 20110123 23:11
integer,allocatable :: ID_DD(:)
! 堰闸节点ID号；: M_DD ! 未用 20110125
integer,allocatable ::
IDDhash(:) ! 水闸哈稀数；: N_ND
character(len=20),allocatable
:: DGTTP(:) ! 闸门型式；: M_DD 20110215 18:21
real,allocatable :: DWRCREL(:)
! 堰顶高程；: M_DD
integer,allocatable ::
DTGTNUM(:) ! 总孔数；: M_DD
real,allocatable :: DSGWD(:) !
单孔宽；: M_DD
real,allocatable :: DAVGTWD(:)
! 平均孔宽；: M_DD
real,allocatable :: DTGTWD(:) !
总孔宽；: M_DD
real,allocatable :: DGTHD(:) !
闸门高度；: M_DD
real,allocatable :: DSWDSQ(:) !
水库堰闸单孔设计流量；: M_DD
real,allocatable :: DTWDSQ(:) !
水库堰闸总设计流量；: M_DD
!
! 表104 测站关系及起调水位 20110215
14:22
integer :: iRL,M_RL !
关系站序号；关系站总数,最多为泵站,水闸及用水户的总和 20110215 14:38
integer,allocatable :: RSTCD(:)
! 关系站编码；:,: M_WU ! 20110215 19:11
character(len=2),allocatable ::
RLMRK(:) ! 关系标志；:,: M_WU ! 20110215 13:55
integer,allocatable ::
UPSTCD(:),DWSTCD(:) ! 上游关联站码,下游关联站码；: M_WU
real,allocatable ::
UPBGZ(:),DWBGZ(:) ! 上游起调水位,下游起调水位；: M_WU
！
real dQ38 ! 调整节点流量与其指定流量的差值 20110128 22:43
real dW38 !
调整用水户缺水流量 20110128 23:12
!
以下实施比较繁，当下无需采用（根据：水洞收缩原理） 20110130 18:13
!integer ::
M_ZF,iZF ! 节点图中的分流节点 ! 20110129 3:56
!integer,allocatable :: M_ZFO(:) ! 分流节点相邻的分流河段数；: M_ZF ! 20110130
11:58
!integer,allocatable :: M_ZFI(:) ! 分流节点相邻的入流河段数；: M_ZF
!integer,allocatable :: ID_ZFO(:,:) ! 分流河段ID号；:,:
M_ZF,M_NJLmax
!integer,allocatable :: ID_ZFI(:,:) ! 入流河段ID号；:,:
M_ZF,M_NJLmax
!real,allocatable :: QO_ZFO(:) ! 总分流；: M_ZF
!real,allocatable :: QI_ZFI(:) ! 总入流；: M_ZF
!real,allocatable :: F_ZFO(:,:) ! 分流比；: M_ZF,M_ZFO(:) 20110130
12:38
!integer,allocatable :: IZFhash(:) ! 分流点哈稀数；: N_ND 20110130
15:23
!integer,allocatable :: ID_ZF(:) ! 分流点ID号；: M_ZF
!
! 表137
水厂(用水户)及取水口信息(隐含泵站) 20110210 10:41
integer,allocatable :: LSTN(:) ! 相邻取水口数(以测站编码统计,无序)；: M_WU
integer,allocatable :: WTTP(:) ! 取水类型；: M_WU
! 1按相邻节点顺序取水以满足所需(省缺值),2按设计流量比例分配以满足所需,
! 3先河水后库水且按相邻节点顺序,4先河水后库水且库容大者优先
real,allocatable :: CAPACITY(:) ! 日生产能力(万m3/d)；: M_WU ! 20110210
11:36
!
integer ::
LSSTNmax ! 最大相关取水口数 ! 20110210 18:04
integer ::
LSNDNmax ! 最大相邻取水节点数 ! 20110210 20:53
integer,allocatable :: LSCD(:) ! 临时相关取水口编码；:
LSTN(iWU)<=LSSTNmax
real,allocatable :: LDNQ(:) ! 临时相关取水口设计流量；: LSTN(iWU) ! 20110210
22:50
integer,allocatable :: LSTCD(:,:) ! 相关取水口编码；:,:
M_WU,LSTN(iWU)
real,allocatable :: LDSNQ(:,:)
! 相关取水口设计流量(m3/s,按表137相关取水口编码顺序)；:,:
M_WU,LSTN(iWU) ! 20110210 12:10
!
real,allocatable :: QDSN(:,:)
! 取水节点设计流量(m3/s,按表333相邻取水节点ID号顺序)；:,:
M_WU,LNDN(iND) ! 20110210 16:14
! QDSN(:,:)的iND与LDSNQ(:,:)的iWU相匹配(即为同一节点)
20110210 18:57
!real,allocatable :: QWR(:) ! 需水量(隐含泵站抽水流量)；: M_WU 20110208 15:32
不可用,与警戒水位流量重名
real,allocatable :: QWUR(:,:) ! 需水流量(隐含泵站抽水流量)；:,: M_WU,LNDN(iND)
20110210 14:42
! 规定退水点唯一,最后一个相关节点表示的需水量即为总需水量20110208
16:22
! QWUR(:,:)与QDSN(:,:)相对应,QWUR为实时所需,QDSN为设计值
20110210 15:56
end
module
!
!
!
!
GNread333.F90 20080602 22:00
! 读入数据
subroutine
GNread333 ! 20080602
use MSK_Mod
use GN_Mod
use WYH_Mod ! 20081028 10:49
use GT_Mod ! 20090419 21:38
implicit none
!
!
!
if(STTP(iND).eq.'DP') then ! 为珠海新增 泵站取水点 20110125 19:12
JH_TYPE(iJH) = 13 ! 泵站与用水户性质相近,共用一个数值 20110130
23:26
BJ_type(iJH) = 1 ! 所有泵站均为边界,或直接或间接 20110130
23:27
!JH_TYPE(iJH) = 14 ! 未用
!!BJ_type(iJH) = 1 ! 暂时做不到直接为边界 20110124 1:30 可以了
20110124 19:13
endif
!
!
DRNA(N_ND) = 1.0 !
用于计算方程系数A(i,k) 20110211 12:45
NDTP(N_ND) = 51 ! 可省 20081126
21:47 ! 改为需要 20110211 12:46
STTP(N_ND) = 'ZC' !
是否计算方程系数A(i,k)时有用，勿缺 20110123 9:06 ! 改为可省 20110211 12:46
LNDID(iND,j) = NDID(N_ND) !
替换相邻节点(对应表333) 20110130 16:27
! GN_Read.F90 20080603
! 读入数据
subroutine GN_Read
!
! 从表137中移出 20110211 23:51
do iWU = 1,M_WU
jND = ID_WU(iWU)
!
jJH = IJHhash(jND) ! 20110211 23:56
if(NT_JH(jJH).eq.1) then ! 边界条件只有1个数据 20110211
23:59
!
LS = LNDN(jND)
if(STTP(jND).eq.'DP') then !
泵站
QWUR(iWU,LS)
= 0. ! 需水总流量赋初值 20110211 17:55
jDP =
IDPhash(jND)
if(PMCL(jDP).gt.250.0) then ! 咸度超标
QWUR(iWU,1) = 0.0 ! 20110211 18:41
else
if(PMPQ(jDP).gt.0.0) then
QWUR(iWU,1) = PMPQ(jDP) ! PMPQ(jDP)从表209读入 20110211 18:13
else
QWUR(iWU,1) = QDSN(iWU,1) ! 按设计流量 20110211
18:52
endif
endif
!
if(LS.eq.2)
then ! 按程序设定这时只有1个取水口,量入为出
QWUR(iWU,LS) = QWUR(iWU,1)
!
elseif(LS.eq.3) then ! 有2个进水口,非此即彼 （!洪湾泵站）
if(QWUR(iWU,1).gt.0.0) then
QWUR(iWU,2) = 0.0
else
if(PMPQ(jDP).gt.0.0) then
QWUR(iWU,2) = PMPQ(jDP) ! 20110211 19:02
else
QWUR(iWU,2) = QDSN(iWU,1) !
20110211 18:52
endif
endif
QWUR(iWU,LS) = QWUR(iWU,1) + QWUR(iWU,2)
!
elseif(LS.ge.4) then ! 有3个以上进水口,争先恐后（!广昌泵站） 20110211 19:26
if(QWUR(iWU,1).gt.0.0) then
do i = 2,LS - 1
QWUR(iWU,i) = 0.0
enddo
QWUR(iWU,LS) = QWUR(iWU,1)
!
else
QWUR(iWU,1) = 0.0
if(PMPQ(jDP).gt.0.0) then
QWUR(iWU,LS) = PMPQ(jDP)
else
QWUR(iWU,LS) = QDSN(iWU,1) !
20110211 20:18
endif
!
do i = 2,LS - 1
QWUR(iWU,i) = 0.0
enddo
RLS = 0.
do i = 2,LS - 1
iND =
INDhash(LNDID(jND,i))
!
if(STTP(iND).eq.'DP') then !
经相邻泵站取水
iDP =
IDPhash(iND)
if(PMCL(iDP).gt.250.0) then ! 咸度超标
QWUR(iWU,i) = 0.0 ! 20110211 19:48
else
if(PMPQ(iDP).gt.0.0) then
QWUR(iWU,i) = PMPQ(iDP) ! 20110211 19:02
else
QWUR(iWU,i) = QDSN(iWU,i) ! 20110211 19:56
endif
RLS = RLS + QWUR(iWU,i)
if(RLS.ge.QWUR(iWU,LS)) then
QWUR(iWU,i) = QWUR(iWU,i) - RLS +
QWUR(iWU,LS)
goto 120
endif
endif
!
else ! 最后通过其它节点取水 20110211
20:46
QWUR(iWU,i)
= QWUR(iWU,LS) - RLS
goto
120
endif
enddo
120 continue
endif
endif
!
else ! 用水户 20110211
20:50
if(LS.eq.2)
then ! 按程序设定这时只有1个取水口,量入为出（!香洲水厂）
QWUR(iWU,1) = QWUR(iWU,2) ! QWUR(iWU,2)从表261读入
!
elseif(LS.ge.3) then ! 有2个以上取水口,分类协调（!拱北水厂）
do i = 1,LS - 1
QWUR(iWU,i) = 0.0
enddo
RLS = 0.0
if(WTTP(iWU).eq.1) then ! 排队索取
do i = 1, LS - 1
QWUR(iWU,i) =
QDSN(iWU,i)
RLS = RLS + QWUR(iWU,i)
if(RLS.ge.QWUR(iWU,LS))
then
QWUR(iWU,i)
= QWUR(iWU,i) - RLS + QWUR(iWU,LS)
goto 250 !
20110211 21:24
endif
enddo
250 continue
! 哈哈
elseif(WTTP(iWU).eq.2)
then ! 比例分配 20110211 21:48
RLS1 = 0.0
do i = 1, LS - 1
RLS1 = RLS1 +
QDSN(iWU,i)
enddo
do i = 1, LS - 1
QWUR(iWU,i) = QWUR(iWU,LS) *
QDSN(iWU,i) / RLS1
enddo
elseif(WTTP(iWU).eq.3) then ! 流水争先
! 咖啡麦当劳(21.0) 20110211 22:25
elseif(WTTP(iWU).eq.4) then ! 遇强则强
!
endif
endif
endif
!
do i2 = 1,NT
Q_JH(jJH,i2)
= QWUR(iWU,LS) ! 20110212 0:14
QQ_JH(jJH,i2) = Q_JH(jJH,i2)
enddo
endif
!
enddo
!
!
!
do i = 1, M_NR !
读表130前
TDSQ(i) = 0.0 !
避免水库空缺时无值(下午,由于随意删去一个水库编码,结果一塌糊涂) 20110215 23:59
enddo
!
!
!
! ZHread209.F90 20110121 9:23 ！ (莲花山外山，水网洞中洞 20110214 0:53)
! 读入表209,泵站水情
subroutine ZHread209
use MSK_Mod
use GN_Mod
use WYH_Mod ! 20081201 14:58
implicit none
!
inquire(file=Filename209,exist=alive)
if(alive) then
open(209,file=Filename209,status='old')
read(209,*)
read(209,*)
iDP = 0 !
统计泵站数 20110128 12:13
do
while(.true.)
do i = 3,14
CLS(i) = '-100' ! 下面读入有空缺时赋值为-100 20110207
16:16
enddo
read(209,*,iostat=status) k,TM_ymd,TM_hms,CLS(3),CLS(4),
& ! 20080914 11:08
CLS(5),CLS(6),CLS(7),CLS(8),CLS(9),CLS(10),CLS(11),CLS(12),CLS(13),
& ! 20110121 9:25
CLS(14) ! 20110213 1:08
if(status/=0) exit
if(any(ISTCD(:).eq.k)) then
if(ISThash(k).le.O_ND) then ! 该节点为河网节点(避免无关节点)
20081205 18:19
iND = ISThash(k)
iJH = IJHhash(iND)
if(STTP(iND).eq.'DP') then !
泵站节点 20110121 9:27
iDP = iDP +
1 ! 20110128 12:14
!PDCHCD(iDP)
= CLS(12) ! 20110207 14:01
read(CLS(7),*) PMPQ(iDP) ! 20110207 15:57
read(CLS(8),*) RLS3 ! 调整位置 20110213 1:19i
!read(CLS(13),*) RLS3 ! 20110122 18:39
PMCL(iDP) =
RLS3 ! 20110207 16:02
if(CLS(13).eq.'1') then ! 20110207 13:48 (happy two year!)
read(CLS(3),*)
PPUPZ(iDP) ! 20110203 15:00 (水洞方除夕)
elseif(CLS(13).eq.'0'.or.CLS(12).eq.'O') then
read(CLS(4),*) PPDWZ(iDP)
endif
PDCHCD(iDP)
= CLS(13) ! 20110213 18:32
read(CLS(14),*) IPRCD(iDP) ! 泵站引排目标水库 20110213 1:33
! (!南沙湾泵站——大镜山水库,洪湾泵站——竹仙洞水库,OK）
! (!广昌泵站——南屏水库,增加一个虚拟泵站节点,以避免影响南沙湾及洪湾泵站 20110214 15:01）
! (!蛇地坑,银坑泵站从属洪湾泵站)
call TM_mdhm ! 20080619 9:17
if(iSCHTP_now.ne.2) then ! 20081201 14:41
if(TM_year.ge.BGTM_year.and.TM_year.le.EDTM_year)
then ! 增加if/end行，过滤掉计算时间以外的数据 20081012 16:32
goto 12
endif
elseif(iSCHTP_now.eq.2) then !
if(TM_year.ge.REFTM_year.and.TM_year.le.REDTM_year)
then
if(TM_mon.eq.2.and.RMON_TM(TM_year,2).eq.29) then ! 合并为一句 20081215
23:50
TM_year = BGTM_year + TM_year
- REFTM_year
if(MON_TM(TM_year,2).eq.29)
then
goto
12
endif
else
TM_year = BGTM_year + TM_year
- REFTM_year
goto 12
endif ! 时间转换到当前计算年份 20081201 23:16
endif
endif
goto 24 !
20110120 12:34
12
continue
!if(STTP(iND).eq.'DP') then ! 20110201 2:40
NT_JH(iJH) = NT_JH(iJH) + 1
i = NT_JH(iJH)
OBT(iJH,i) = (MONDAY_TTM(TM_year,TM_mon) + TM_day - 1)*24. + TM_hh
+ TM_mm/60. ! 20081013 10:55
!if(NDTP(iND).eq.22) then ! 泵站为直接边界 20110120
14:29
if(RLS3.le.250.0) then ! 与引排特征值无关 20110125
14:48
read(CLS(7),*) OBQ(iJH,i) ! 给定边界流量
elseif(RLS3.gt.250.0) then
OBQ(iJH,i) = 0.0 ! 咸度超标,停止外江抽水 20110124
18:38
endif
if(NDTP(iND).eq.69) then ! 泵站为间接边界 20110124
19:04
do j = 1,LNDN(iND)
iNDL = INDhash(LNDID(iND,j))
if(NDTP(iNDL).eq.22) then !
相邻节点为边界
iJHL = IJHhash(iNDL)
NT_JH(iJHL)
= i ! 注意定义 20110201 3:50
OBT(iJHL,i) = (MONDAY_TTM(TM_year,TM_mon) + TM_day - 1)*24. + TM_hh
+ TM_mm/60. ! 20110201 2:17
if(RLS3.le.250.0) then ! 外江取水口引水,另设外边界
OBQ(iJHL,i) = OBQ(iJH,i) ! 20110201 2:47
elseif(RLS3.gt.250.0) then ! 内河道或管网排水,内边界 20110121 18:29
OBQ(iJHL,i) = 0.0 ! 咸度超标 20110124 18:37
endif
endif
enddo
endif
24
continue
!
!
if(iDP.ne.M_DP) then ! 2110128 11:15
I_error = I_error + 1
EWSHOW = '表209的泵站数与表333中的不匹配'
write(14,"(i3,',',a1,',',i3,',',2(a1,','),a64,',')")
&
I_error,Blank,DBID(209),Blank,Blank,EWSHOW
endif
!
close(209)
else
I_warning = I_warning + 1 !
20090419 3:18
EWSHOW = '缺少表209文件'
write(14,"(a1,',',i6,',',i3,',',2(a1,','),a33,',')")
&
Blank,I_warning,DBID(209),Blank,Blank,EWSHOW ! 20080616 20:56
endif
!
end subroutine ZHread209
!
!
!
!
GNread137.F90 20110210 10:58
!
读入表137数据水厂及取水口信息
subroutine
GNread137
!
use MSK_Mod
use GN_Mod
use WYH_Mod ! 20081218 15:48
implicit none
!
! 水厂(用水户)与取水口(隐含泵站)
allocate (LSTN(M_WU),WTTP(M_WU)) ! 表137 20110210
16:50
allocate (CAPACITY(M_WU))
allocate (QDSN(M_WU,LNDNmax)) ! 20110210
16:53
allocate (QWUR(M_WU,LNDNmax)) ! 20110208
16:08
!
LSNDNmax = 0 ! 统计用水户最大相邻取水节点数
do iWU = 1,M_WU
iND = ID_WU(iWU)
LSTN(iWU) = 0 ! 赋初值 20110210 23:36
if(STTP(iND)(1:1).eq.'W') then ! 排除泵站 20110210
22:17
if((LNDN(iND)-1).gt.LSSTNmax)
then
LSNDNmax =
LNDN(iND) - 1 ! 20110210 21:03
endif
endif
enddo
!
LSSTNmax = 0 ! 统计最大相关取水口数
inquire(file=Filename137,exist=alive)
if(alive) then
open(137,file=Filename137,status='unknown')
read(137,*)
read(137,*)
LS = 0 ! 统计用水户数
do while(.true.)
do j =
1,8
CLS(j) = '0' ! 20110210 17:17
enddo
CLS(9) = '1'
! 20110210 17:19
read(137,*,iostat=status)
ILS(1),CLS(2),CLS(3),CLS(4),CLS(5),CLS(6),CLS(7),CLS(8),CLS(9) !
20110210 17:26
if(status/=0) exit
!
read(CLS(8),*) ILS(8)
if(ILS(8).gt.LSSTNmax) LSSTNmax = ILS(8) ! 不管是否有用,都加入统计,以免溢出
20110211 3:55
!
if(any(ISTCD(:).eq.ILS(1))) then
iND = ISThash(ILS(1))
if(STTP(iND)(1:1).eq.'W') then
LS = LS +1 !
iWU = IWUhash(iND) ! 20110210 17:33
read(CLS(6),*) CAPACITY(iWU)
!read(CLS(8),*) LSTN(iWU)
LSTN(iWU) = ILS(8) ! 20110211 3:57
read(CLS(9),*) WTTP(iWU)
!if(LSTN(iWU).gt.LSSTNmax) LSSTNmax = LSTN(iWU)
! 可能产生溢出 20110211 3:56
endif
else
if(any(ISTCD(:).ne.ILS(1))) then
I_warning = I_warning + 1
EWSHOW = '没有用到的水厂，测站编码：' ! 20110210 17:41
write(14,"(a1,',',i6,',',i3,',',a1,',',a1,',',a26,i8,',')")
&
Blank,I_warning,DBID(137),Blank,Blank,EWSHOW,ILS(1)
endif
enddo
if(LSNDNmax.le.LSSTNmax)
then
!
elseif(LSNDNmax.gt.LSSTNmax)
then ! 20110210 21:17
if(LSNDNmax.ge.2) then
I_error = I_error + 1
EWSHOW = '表137水厂取水口不够多' ! 宁多勿少
write(14,"(i3,',',a1,',',i3,',',2(a1,','),a33,',')")
&
I_error,Blank,DBID(137),Blank,Blank,EWSHOW
endif
endif
if((LS+M_DP).ne.M_WU) then !
2110210 20:41
I_warning = I_warning + 1
EWSHOW = '表137的水厂数与表333中的不匹配'
write(14,"(a1,',',i3,',',i3,',',2(a1,','),a36,',')")
&
Blank,I_warning,DBID(137),Blank,Blank,EWSHOW
endif
close(137)
else
if(LSNDNmax.ge.2) then ! 20110210 22:21
I_error = I_error + 1
EWSHOW = '缺少表137文件'
write(14,"(i3,',',a1,',',i3,',',2(a1,','),a33,',')")
&
I_error,Blank,DBID(137),Blank,Blank,EWSHOW
endif
endif
!
allocate (LSCD(LSSTNmax),LSTCD(M_WU,LSSTNmax)) !
20110210 19:03
allocate (LDNQ(LSSTNmax),LDSNQ(M_WU,LSSTNmax)) !
20110210 22:55
!
inquire(file=Filename137,exist=alive)
if(alive) then
open(137,file=Filename137,status='unknown')
read(137,*)
read(137,*)
do while(.true.)
read(137,*,iostat=status)
ILS(1),CLS(2),CLS(3),CLS(4),CLS(5),CLS(6),CLS(7),ILS(8),CLS(9)
&
,(LSCD(j),LDNQ(j),j=1,ILS(8))
if(status/=0) exit
if(any(ISTCD(:).eq.ILS(1))) then
iND = ISThash(ILS(1))
if(STTP(iND)(1:1).eq.'W') then
iWU = IWUhash(iND) ! 20110210 17:33
if(LSTN(iWU).ge.1) then
do j = 1,LSTN(iWU)
LSTCD(iWU,j)
= LSCD(j) ! 20110210 19:47
LDSNQ(iWU,j)
= LDNQ(j) ! 20110210 23:00
enddo
endif
endif
endif
enddo
close(137)
!
do iWU = 1,M_WU
jND = ID_WU(iWU)
LS = LNDN(jND)
QDSN(iWU,LS) = 0. !
取水总设计流量赋初值
if(LSTN(iWU).ge.1) then
do i = 1,LS
- 1 ! 20110210 0:57
iND = INDhash(LNDID(jND,i))
do j = 1,LSTN(iWU)
if(LSTCD(iWU,j).eq.ISTCD(iND)) then
QDSN(iWU,i) = LDSNQ(iWU,j)
QDSN(iWU,LS) = QDSN(iWU,LS) +
QDSN(iWU,i)
endif
enddo
enddo
elseif(LSTN(iWU).eq.0)
then
if(STTP(jND)(1:1).eq.'W') then ! 用水户只有一个取水口时LSTN可省缺为0 20110210
1:04
! 不用考虑QDSN
endif
if(STTP(jND).eq.'DP') then ! 泵站不统计LSTN省缺为0
jDP = IDPhash(jND)
if(LS.ge.3) then ! 有2个以上取水口 20110211 17:07
!do i = 1,LNDN(jND) - 1 ! 20110210 0:57
do i = 1,LS - 1 ! 20110211 17:08
iND =
INDhash(LNDID(jND,i))
if(STTP(iND).eq.'DP') then !
经相邻泵站取水
iDP =
IDPhash(iND)
QDSN(iWU,i) = PTDSQ(2,iDP)
else ! 经泵站本身取水口(可以有多处选择)
20110211 17:38
QDSN(iWU,i) = PTDSQ(2,jDP)
endif
QDSN(iWU,LS) = QDSN(iWU,LS) +
QDSN(iWU,i) ! 意义不大 20110210 1:55
enddo
endif
endif
endif
!
enddo
!
endif
!
end
subroutine GNread137 ! 20110210 10:57
!
! GNread104.F90 20110215
0:07
! 读入表104数据 测站关系信息
subroutine ZHread104
!
use MSK_Mod
use GN_Mod
use WYH_Mod ! 20081218 15:48
implicit none
!
allocate (UPSTCD(M_WU),UPBGZ(M_WU),RLMRK(M_WU))
! 20110215 14:59
allocate
(DWSTCD(M_WU),DWBGZ(M_WU),RSTCD(M_WU))
!
inquire(file=Filename104,exist=alive)
if(alive) then
open(104,file=Filename104,status='unknown')
read(104,*)
read(104,*)
M_RL = 0 ! 统计关系站数
do while(.true.)
do j = 1,6
CLS(j) = '0'
enddo
read(104,*,iostat=status) CLS(1),CLS(2),CLS(3),CLS(4),CLS(5),CLS(6)
! 20110215 17:56
if(status/=0) exit
read(CLS(1),*) ILS(1) ! 部分表没这样做 20110215 15:11
if(any(ISTCD(:).eq.ILS(1))) then
iND
= ISThash(ILS(1)) ! 20110215 17:36
if(STTP(iND).eq.'DP'.or.STTP(iND).eq.'DD'.or.STTP(iND)(1:1).eq.'W')
then
M_RL = M_RL +1 !
RSTCD(M_RL) = ILS(1) ! 20110215 19:14
read(CLS(2),*) RLMRK(M_RL)
if(RLMRK(M_RL).eq.'U') then
read(CLS(3),*) UPSTCD(M_RL)
read(CLS(4),*) UPBGZ(M_RL)
elseif(RLMRK(M_RL).eq.'D') then
read(CLS(5),*) DWSTCD(M_RL)
read(CLS(6),*) DWBGZ(M_RL)
elseif(RLMRK(M_RL).eq.'B') then
read(CLS(3),*) UPSTCD(M_RL)
read(CLS(4),*) UPBGZ(M_RL)
read(CLS(5),*) DWSTCD(M_RL)
read(CLS(6),*) DWBGZ(M_RL)
endif
endif
else if(any(ISTCD(:).ne.ILS(1))) then
I_warning = I_warning + 1
EWSHOW = '没有用到的关系站，测站编码：' ! 20110215 15:25
write(14,"(a1,',',i6,',',i3,',',a1,',',a1,',',a26,i8,',')")
&
Blank,I_warning,DBID(104),Blank,Blank,EWSHOW,ILS(1)
endif
enddo
if(M_RL.lt.M_WU) then ! 2110215 15:29
I_warning = I_warning + 1
EWSHOW = '表104的关系站数没有覆盖表333中的全部'
write(14,"(a1,',',i3,',',i3,',',2(a1,','),a36,',')")
&
Blank,I_warning,DBID(104),Blank,Blank,EWSHOW
endif
close(104)
else
if(M_RL.eq.0) then ! 20110215 15:31
I_warning = I_warning + 1
EWSHOW = '缺少表104文件或不存在关系站'
write(14,"(a1,',',i3,',',i3,',',2(a1,','),a33,',')")
&
Blank,I_warning,DBID(104),Blank,Blank,EWSHOW
endif
endif
end subroutine ZHread104 !
20110215 0:08
!
!
!
! GN_VFWater.F90 20080703
10:49
!PROGRAM
GN_VFWater
subroutine
GN_Main
!
use MSK_Mod
use GN_Mod
use dfcom
use WYH_Mod ! 20081119 10:20
implicit none
!
!
!
! 关系节点 20110215 19:04
do iRL = 1,M_RL
iND = ISThash(RSTCD(iRL)) !
20110215 19:15
iJH = IJHhash(iND)
iWU = IWUhash(iND)
if(STTP(iND).eq.'DP')
then
if(RLMRK(iRL).eq.'U') jND = ISThash(UPSTCD(iRL))
if(RLMRK(iRL).eq.'D') jND = ISThash(DWSTCD(iRL))
if(STTP(jND).eq.'RR') then ! 20110215 19:49
jNR = INRhash(jND) ! 目标水库
if(RLMRK(iRL).eq.'U'.and.H1_NR(jNR).ge.DWBGZ(iRL)) then ! 20110215
19:31
Q_JH(iJH,i2) = 0.0 ! 无需动用调咸泵站 20110213 19:43
QQ_JH(iJH,i2) = Q_JH(iJH,i2)
elseif(RLMRK(iRL).eq.'D'.and.H1_NR(jNR).ge.DWBGZ(iRL)) then !
20110215 19:29
RLS = FSLTDZ_tm(jNR,i2) - 0.01 ! 避免有数字接近但有微小误差而产生的问题 29119213
22:49
if(H1_NR(jNR).ge.RLS) then
Q_JH(iJH,i2) = 0.0 ! 关闭泵站
QQ_JH(iJH,i2) = Q_JH(iJH,i2)
if(LNDN(iND).ge.3) then ! 20110213 20:46
do k = 1, LNDN(iND)
QWUR(iWU,k) = 0.0 ! 目前与时间无关 20110213 23:26
enddo
endif
else
if(i2.ge.2.and.Q_JH(iJH,(i2-1)).eq.0.0) then
Q_JH(iJH,i2) = 0.0 ! 继续关闭泵站
QQ_JH(iJH,i2) = Q_JH(iJH,i2)
if(LNDN(iND).ge.3) then ! 20110213 20:46
do k = 1, LNDN(iND)
QWUR(iWU,k) = 0.0
enddo
endif
endif
endif
if(NDTP(iND).eq.69) then ! 泵站为间接边界 20110213 21:18
do j = 1, LNDN(iND) - 1
iNDL = INDhash(LNDID(iND,j))
if(NDTP(iNDL).eq.22) then ! 相邻节点为边界
iJHL
= IJHhash(iNDL)
Q_JH(iJHL,i2) = Q_JH(iJH,i2)
QQ_JH(iJHL,i2) = Q_JH(iJHL,i2)
endif
enddo
endif
!
endif
!
endif
endif
enddo
!
!
!
if(M_TJ.ge.1) then !
存在调节点
do iTJ = 1,M_TJ ! 都与泵站相关
20110201 4:54 ! 与水闸等也有关 20110218 15:39
jND =
ID_TJ(iTJ)
jJH =
IJHhash(jND)
dQ38 =
QQ_JH(jJH,i2) - Q_TJ(iTJ) ! 20110128 22:45
if(dQ38.gt.0.00001) then ! 避免数据接近0 20110218 11:34
do j = 2,LNDN(jND) ! 调节点从第2个节点开始
iND = INDhash(LNDID(jND,j))
iJH = IJHhash(iND)
LS = LNDN(iND) ! 20110218 11:54
if(LS.eq.1) then ! 外边界 20110218 11:55
if(dQ38.gt.0.0) then ! 20110218 13:47
RLS = Q_JH(iJH,i2) - dQ38 ! 20110218 13:50
if(RLS.ge.0.0) then
Q_JH(iJH,i2) = Q_JH(iJH,i2) - dQ38 ! OK,此为新的边界流量 20110128
22:37
dQ38 = 0.
else
dQ38 = dQ38 - Q_JH(iJH,i2) ! 还要继续
Q_JH(iJH,i2) = 0.0
endif
endif
elseif(LS.eq.2) then ! 传统用水户边界(一进一出) 20110218 12:00
if(dQ38.gt.0.0) then ! 20110218 13:47
RLS = Q_JH(iJH,i2) - dQ38 ! 20110218 13:50
if(RLS.ge.0.0) then
Q_JH(iJH,i2) = Q_JH(iJH,i2) - dQ38 ! OK,此为新的边界流量 20110128
22:37
dQ38 = 0.
else
dQ38 = dQ38 - Q_JH(iJH,i2) ! 还要继续
Q_JH(iJH,i2) = 0.0 ! 新的边界流量为0
endif
iNJ = INJhash(iND)
Q_NJ(iNJ,i2)
= Q_JH(iJH,i2) ! 给内节点赋值
iNDL
= INDhash(LNDID(iND,1))
iJHL
= IJHhash(iNDL)
Q_JH(iJHL,i2) = Q_JH(iJH,i2) ! 20110218 12:19
QQ_JH(iJHL,i2) = Q_JH(iJHL,i2)
endif
endif
elseif(LS.ge.3) then ! 水洞型用水户边界(多进一出) 20110218 12:25
iWU
= IWUhash(iND) ! 20110218 13:16
QQWUR(iWU,LS)
= 0. ! 需水总流量赋初值 20110218 12:59
do k = 1, LS - 1
if(dQ38.gt.0.0) then ! 20110218 13:47
RLS = QQWUR(iWU,k) - dQ38 ! 20110218 13:50
if(RLS.ge.0.0) then
QQWUR(iWU,k) = QQWUR(iWU,k) - dQ38 ! OK,此为新的边界流量之一 20110218
12:49
dQ38 = 0.
else
dQ38 = dQ38 - QQWUR(iWU,k) ! 还要继续,直到结束k 20110218 13:31
QQWUR(iWU,k) = 0.0
iNDL
= INDhash(LNDID(iND,k))
if(NDTP(iND).eq.69.and.NDTP(iNDL).eq.22) then ! 泵站为间接边界 20110218
15:46
iJHL
= IJHhash(iNDL)
Q_JH(iJHL,i2) = QQWUR(iWU,k) ! 20110218 12:19
QQ_JH(iJHL,i2) = Q_JH(iJHL,i2)
endif
endif
QQWUR(iWU,LS)
= QQWUR(iWU,LS) + QQWUR(iWU,k)
endif
enddo
Q_JH(iJH,i2)
= QQWUR(iWU,LS) ! 20110218 13:03
iNJ =
INJhash(iND)
Q_NJ(iNJ,i2)
= Q_JH(iJH,i2)
endif
!
QQ_JH(iJH,i2)
= Q_JH(iJH,i2)
iWU =
IWUhash(iND)
W_Q(iWU,i2) =
Q_JH(iJH,i2) ! 20110201 7:47
!
enddo ! j
!
elseif(dQ38.eq.0.0) then ! 20110128 23:08
dQ38 = 0.0 ! 含义不同了，统计用水户缺水流量 20110128 23:24
do j = 2,LNDN(jND)
iND = INDhash(LNDID(jND,j))
if(NDTP(iND).ge.70.and.NDTP(iND).le.79) then
iWU = IWUhash(iND)
if(W_QQ(iWU,i2).lt.W_Q(iWU,i2)) then ! 用水小于需水 20110128 21:38
dQ38 = dQ38 + W_Q(iWU,i2) - W_QQ(iWU,i2)
endif
endif
enddo
if(dQ38.gt.0.0) then
do j = 2,LNDN(jND)
iND = INDhash(LNDID(jND,j))
iJH = IJHhash(iND) ! 20110202 21:39
if(dQ38.gt.0.0.and.NDTP(iND).eq.22.and.STTP(iND).eq.'DP') then !
指定泵站 20110128 23:49
iDP = IDPhash(iND)
RLS1 = PTDSQ(2,iDP) - Q_JH(iJH,i2)
if(RLS1.ge.dQ38) then ! 允许流量足够大 20110128 23:38
Q_JH(iJH,i2)= Q_JH(iJH,i2) + dQ38 ! OK,此为新的边界流量 20110128
23:29
dQ38 = 0.
else
Q_JH(iJH,i2) = Q_JH(iJH,i2) + RLS1
dQ38 = dQ38 - RLS1 ! 还要继续
endif
QQ_JH(iJH,i2) = Q_JH(iJH,i2)
iWU = IWUhash(iND) ! 20110201 9:28
W_Q(iWU,i2) = Q_JH(iJH,i2) ! 20110201 7:47
if(LNDN(iND).ge.2) then ! 内边界 20110125 20:58
iNJ = INJhash(iND)
Q_NJ(iNJ,i2)
= Q_JH(iJH,i2) ! 给内节点赋值
endif
elsif(dQ38.gt.0.0.and.NDTP(iND).eq.69) then ! 泵站为间接边界 20110124
19:04
do k = 1,LNDN(iND)
iNDL = INDhash(LNDID(iND,k))
iJHL
= IJHhash(iNDL)
if(dQ38.gt.0.0.and.NDTP(iNDL).eq.22) then ! 相邻节点为边界
iDP = IDPhash(iND)
RLS1 = PTDSQ(2,iDP) - Q_JH(iJHL,i2)
if(RLS1.ge.dQ38) then ! 允许流量足够大 20110128 23:38
Q_JH(iJHL,i2)= Q_JH(iJHL,i2) + dQ38 ! OK,此为新的边界流量 20110128
23:29
Q_JH(iJH,i2) = Q_JH(iJHL,i2) ! 20110202 21:20
dQ38 = 0.
else
Q_JH(iJHL,i2) = Q_JH(iJHL,i2) + RLS1
Q_JH(iJH,i2) = Q_JH(iJHL,i2) ! 20110202 21:20
dQ38 = dQ38 - RLS1 ! 还要继续
endif
QQ_JH(iJHL,i2) = Q_JH(iJHL,i2)
QQ_JH(iJH,i2) = Q_JH(iJH,i2) ! 20110202 21:24
iWU = IWUhash(iND) ! 20110201 9:28
W_Q(iWU,i2) = Q_JH(iJH,i2) ! 20110201 7:47
if(LNDN(iND).ge.2) then ! 20110202 22:09
iNJ = INJhash(iND) ! 20110201 22:05
Q_NJ(iNJ,i2) = Q_JH(iJH,i2) ! 20110202 21:59
endif
endif
enddo ! k
endif
enddo
! j
call
Model_MSK_Main ! 20110129 0:13
call
Model_MSK_Main
elseif(dQ38.lt.0.0) then ! 20110129
0:11
I_warning = I_warning + 1
EWSHOW = '供水破坏'
write(14,"(i3,',',(i4,'-',i2.2,'-',i2.2,'
',i2,':00:00'),a1,',',',',2(a1,','),a33,',')")
&
I_warning,ITM_year(i),ITM_mon(i),ITM_day(i),ITM_hh(i),Blank,Blank,Blank,EWSHOW
! 20110125 22:10
endif
endif
！
！
！
! Model_MSK.F90 20081118 19:03
! 建立河段方程、求解、计算节点流量的主程序
subroutine Model_MSK_Main
use MSK_Mod
use GN_Mod
use WYH_Mod ! 20081106 17:56
use IMSL
implicit none
!
!
!
建立河段方程
do i =
1,M_NH ! 对河段循环，计算河段方程的系数和常数项
!
!
do j = 1,LNDN(iND) !
水网恢恢,疏蓄结合,来去无常 20110212 17:56
iNDL =
INDhash(LNDID(iND,j))
iJHL = IJHhash(iNDL) ! 首节点的相邻节点在全部节点中的序号 20081127 16:33
if(BJ_type(iJHL).eq.1) then ! 20081115 23:04
if(JH_type(iJHL).eq.13) then ! 相邻用水户节点
iNJL = INJhash(iNDL)
if(ID_NJ(iNJ).eq.ID_NJL(iNJL,M_NJL(iNJL)))
then
! 相邻退水点
!
如果增加简单节点以自动增加用水户退水虚拟河段,则不会出现该情况 20110211 15:37
!B(i) = B(i) + C0_ZHD(i)**M_ZHD(i) * Q_NJ(iNJL,i2) ! 20081218
10:21
! iND无论是否某个用水户的退水点都不用继续计算B(i)
20110212 17:06
! 相邻取水点 20081124 2:40
! 1、水洞收缩原理：原则上不允许从用水户边界节点取水形成分叉口 20110212
16:55
! 2、即：出口唯一,不允许节点iND下游存在多个相邻取水点及退水点 20110212
17:15
else ! 3、为安全计,应同时符合膨胀原理
20110211 15:41
!
问题：目前水流QWUR存储在形式上采用收缩原理,边界复杂时可能崩塌 20110212 17:41
if(LNDN(iNDL).le.2) then
B(i) = B(i) - C0_ZHD(i)**M_ZHD(i) * Q_JH(iJHL,i2) ! 20081218
10:21
!B(i) = B(i) - C0_ZHD(i)**M_ZHD(i) * Q_JH(iJHL,i1) ! 20081124
2:43
elseif(LNDN(iNDL).ge.3) then !
允许节点iNDL存在多个相邻取水点 20110212 15:14
do k = 1,
LNDN(iNDL) - 1
jND = INDhash(LNDID(iNDL,k))
if(jND.eq.iND) then
B(i) = B(i) - C0_ZHD(i)**M_ZHD(i) * QWUR(iWU,k) ! 20110212
17:21
endif
enddo
endif
endif
!
endif
endif
enddo ! 20081127 16:34
!
!
elseif(BJ_type(iJH).eq.0) then ! 首节点为平衡内节点 ! 20081022 18:00
!
!
if(LNDN(iNDL).eq.2) then
B(i) = B(i) - C0_ZHD(i)**M_ZHD(i) * Q_JH(iJHL,i2) ! 20081218
10:23
elseif(LNDN(iNDL).ge.3) then !
允许存在多个相邻取水点 20110212 1:32
iWU =
IWUhash(iNDL) ! 20110212 1:35
do k = 1,
LNDN(iNDL) - 1
jND = INDhash(LNDID(iNDL,k))
if(jND.eq.iND) then
B(i) = B(i) - C0_ZHD(i)**M_ZHD(i) * QWUR(iWU,k)
! 子夜乱加,黄昏再减 20110212 18:20
!B(i) = B(i) + C0_ZHD(i)**M_ZHD(i) * QWUR(iWU,k)
! 20110212 1:48
endif
!
!
else
if(BJ_type(iJHL).eq.0) then ! 20081115 23:51
if(JH_type(iJHL).eq.5) then ! 相邻节点为河段节点 ! 包括虚拟河段 20081126
k = INHhash(ID_NJL(iNJ,j)) ! 相邻河段在河段中的序号
if(ID_NJ(iNJ).eq.ID_MJ(k)) then ! 首节点为相邻河段的末节点 20081115 21:58
!iND = ID_NH(i) ! 20110118 !
水洞依稀从此诣 20110214 1:12
jND = ID_NH(i) !
不要和iND混淆,下面两行相应修改 20110212 1:19
if(NDTP(jND).eq.51) then !
与下面相当 20110123 9:12
A(i,k) = A(i,k) - DRNA(jND)*C0_ZHD(i)**M_ZHD(i) ! 取代下面的 20110211
12:49
! 这里的DRNA(iND)为出流系数或消退系数,
! 正常河段为1.0,分流河段为分流比,用水户虚拟退水河段为回归水系数 20110211
11:43
endif
endif
else if(JH_type(iJHL).ne.5) then ! ,,,太难处理了,还是只用河段来连接的好 20081116
2:41
!
!
! 解方程
call
LSARG(M_NH,A,M_NH,B,1,X) ! 20071130
!
!
!
if(LNDN(iNDL).ge.3) then
do k = 1, LNDN(iNDL) - 1
jND = INDhash(LNDID(iNDL,k))
if(jND.eq.iND) then
Q_NJ(iNJ,i2) = Q_NJ(iNJ,i2) - QWUR(iWU,k) ! 暂不考虑缺水 20110212
2:24
if(Q_NJ(iNJ,i2).lt.ZLMNQ(iJH)) then ! 不能让节点取水后的流量小于生态最小值ZLMNQ
dQ_LS = ZLMNQ(iJH) - Q_NJ(iNJ,i2) ! 20081106 11:09
Q_NJ(iNJ,i2) = ZLMNQ(iJH)
QWUR(iWU,k) = QWUR(iWU,k) - dQ_LS
endif
endif
enddo
else
! 由于增加了退水虚拟河段,以下暂不修改 20110212 2:53
if(Q_NJ(iNJ,i2).gt.ZLMNQ(iJH)) then ! 增加此if判别 20081112 17:51
Q_NJ(iNJ,i2)
= Q_NJ(iNJ,i2) - W_Q(iWU,i2)
W_QQ(iWU,i2) = W_Q(iWU,i2) ! 20081202 16:00
!
W_QR(iWU,i2) = W_QQ(iWU,i2) * F_REFL(iWU,i2) ! 20110201 16:11
Q_JH(iJHL,i2)
= W_QQ(iWU,i2)
Q_NJ(iNJL,i2)
= W_QR(iWU,i2)
! 漏得太久了,秒杀了N多时间 20110202 17:44 (虫洞欲除夕)
!
!
do iWU =
1,M_WU
iND = ID_WU(iWU)
iJH = IJHhash(iND) ! 20110211 23:56
if(NT_JH(iJH).eq.1) then ! 边界条件只有1个数据 20110211
23:59
LS = LNDN(iND)
if(STTP(iND)(1:1).eq.'W') then ! 用水户 20110212
23:49
if(LS.ge.3) then ! 为了拱北 20110212 0:04
QWUR(iWU,LS)
= 0. ! 需水总流量赋初值 20110211 17:55
do k = 1, LS
- 1
QWUR(iWU,LS) = QWUR(iWU,LS) + QWUR(iWU,k)
enddo
W_QQ(iWU,i2)
= QWUR(iWU,LS)
W_QR(iWU,i2)
= W_QQ(iWU,i2) * F_REFL(iWU,i2) ! 用水户回归水 20110209 1:06
iNJ = INJhash(iND)
Q_JH(iJH,i2)
= W_QQ(iWU,i2) ! 20110213 0:02
Q_NJ(iNJ,i2)
= W_QR(iWU,i2)
QQ_JH(iJH,i2) = W_QQ(iWU,i2)
endif
endif
endif
enddo
!
!
!
!write(5,"('OK:
广宁水模型计算结束',',')")
write(5,"('OK:
理水模型计算结束',',')") ! 20110215 18:59
!
!
!
! owali 20110218 17:15
水是生命之源、生产之要、生态之基。
兴水利、除水害，事关人类生存、经济发展、社会进步，历来是治国安邦的大事。
促进经济长期平稳较快发展和社会和谐稳定，夺取全面建设小康社会新胜利，必须下决心加快水利发展，切实增强水利支撑保障能力，实现水资源可持续利用。
近年来我国频繁发生的严重水旱灾害，造成重大生命财产损失，暴露出农田水利等基础设施十分薄弱，必须大力加强水利建设。
水利面临的新形势(水情)：
新中国成立以来，特别是改革开放以来，党和国家始终高度重视水利工作，领导人民开展了气壮山河的水利建设，取得了举世瞩目的巨大成就，为经济社会发展、人民安居乐业作出了突出贡献。但必须看到，
人多水少、水资源时空分布不均是我国的基本国情水情。
洪涝灾害频繁是中华民族的心腹大患，水资源供需矛盾突出是可持续发展的主要瓶颈，农田水利建设滞后是影响农业稳定发展和国家粮食安全的最大硬伤，水利设施薄弱是国家基础设施的明显短板。
水利面临的形势更趋严峻，增强防灾减灾能力要求越来越迫切，强化水资源节约保护工作越来越繁重，加快扭转农业主要“靠天吃饭”局面任务越来越艰巨。
２０１０年西南地区发生特大干旱、多数省区市遭受洪涝灾害、部分地方突发严重山洪泥石流，再次警示我们加快水利建设刻不容缓。
新形势下水利的地位和作用(水位)：
水利是现代农业建设不可或缺的首要条件，是经济社会发展不可替代的基础支撑，是生态环境改善不可分割的保障系统，具有很强的公益性、基础性、战略性。
加快水利改革发展，不仅事关农业农村发展，而且事关经济社会发展全局；
不仅关系到防洪安全、供水安全、粮食安全，而且关系到经济安全、生态安全、国家安全。
要把水利工作摆上党和国家事业发展更加突出的位置，着力加快农田水利建设，推动水利实现跨越式发展。
指导思想：
全面贯彻党的十七大和十七届三中、四中、五中全会精神，
以邓小平理论和“三个代表”重要思想为指导，
深入贯彻落实科学发展观，
把水利作为国家基础设施建设的优先领域，
把农田水利作为农村基础设施建设的重点任务，
把严格水资源管理作为加快转变经济发展方式的战略举措，
注重科学治水、依法治水，
突出加强薄弱环节建设，
大力发展民生水利，
不断深化水利改革，
加快建设节水型社会，
促进水利可持续发展，
努力走出一条中国特色水利现代化道路。
目标任务：
力争通过５年到１０年努力，从根本上扭转水利建设明显滞后的局面。
到２０２０年，基本建成防洪抗旱减灾体系，重点城市和防洪保护区防洪能力明显提高，抗旱能力显著增强。
“十二五”期间，基本完成重点中小河流（包括大江大河支流、独流入海河流和内陆河流）重要河段治理、全面完成小型水库除险加固和山洪灾害易发区预警预报系统建设；
基本建成水资源合理配置和高效利用体系，全国年用水总量力争控制在６７００亿立方米以内，城乡供水保证率显著提高，城乡居民饮水安全得到全面保障，万元国内生产总值和万元工业增加值用水量明显降低，农田灌溉水有效利用系数提高到０.５５以上，“十二五”期间新增农田有效灌溉面积４０００万亩；
基本建成水资源保护和河湖健康保障体系，主要江河湖泊水功能区水质明显改善，城镇供水水源地水质全面达标，重点区域水土流失得到有效治理，地下水超采基本遏制；
基本建成有利于水利科学发展的制度体系，最严格的水资源管理制度基本建立，水利投入稳定增长机制进一步完善，有利于水资源节约和合理配置的水价形成机制基本建立，水利工程良性运行机制基本形成。
基本原则：
一要坚持民生优先。着力解决群众最关心最直接最现实的水利问题，推动民生水利新发展。
二要坚持统筹兼顾。注重兴利除害结合、防灾减灾并重、治标治本兼顾，促进流域与区域、城市与农村、东中西部地区水利协调发展。
三要坚持人水和谐。顺应自然规律和社会发展规律，合理开发、优化配置、全面节约、有效保护水资源。
四要坚持政府主导。发挥公共财政对水利发展的保障作用，形成政府社会协同治水兴水合力。
五要坚持改革创新。加快水利重点领域和关键环节改革攻坚，破解制约水利发展的体制机制障碍。
大兴农田水利建设：
到２０２０年，基本完成大型灌区、重点中型灌区续建配套和节水改造任务。
结合全国新增千亿斤粮食生产能力规划实施，在水土资源条件具备的地区，新建一批灌区，增加农田有效灌溉面积。
实施大中型灌溉排水泵站更新改造，加强重点涝区治理，完善灌排体系。
健全农田水利建设新机制，中央和省级财政要大幅增加专项补助资金，市、县两级政府也要切实增加农田水利建设投入，引导农民自愿投工投劳。
加快推进小型农田水利重点县建设，优先安排产粮大县，加强灌区末级渠系建设和田间工程配套，促进旱涝保收高标准农田建设。
因地制宜兴建中小型水利设施，支持山丘区小水窖、小水池、小塘坝、小泵站、小水渠等“五小水利”工程建设，重点向革命老区、民族地区、边疆地区、贫困地区倾斜。
大力发展节水灌溉，推广渠道防渗、管道输水、喷灌滴灌等技术，扩大节水、抗旱设备补贴范围。
积极发展旱作农业，采用地膜覆盖、深松深耕、保护性耕作等技术。
稳步发展牧区水利，建设节水高效灌溉饲草料地。
加快中小河流治理和小型水库除险加固：
中小河流治理要优先安排洪涝灾害易发、保护区人口密集、保护对象重要的河流及河段，加固堤岸，清淤疏浚，使治理河段基本达到国家防洪标准。
巩固大中型病险水库除险加固成果，加快小型病险水库除险加固步伐，尽快消除水库安全隐患，恢复防洪库容，增强水资源调控能力。
推进大中型病险水闸除险加固。
山洪地质灾害防治要坚持工程措施和非工程措施相结合，抓紧完善专群结合的监测预警体系，加快实施防灾避让和重点治理。
抓紧解决工程性缺水问题：
加快推进西南等工程性缺水地区重点水源工程建设，坚持蓄引提与合理开采地下水相结合，以县域为单元，尽快建设一批中小型水库、引提水和连通工程，支持农民兴建小微型水利设施，显著提高雨洪资源利用和供水保障能力，基本解决缺水城镇、人口较集中乡村的供水问题。
提高防汛抗旱应急能力：
尽快健全防汛抗旱统一指挥、分级负责、部门协作、反应迅速、协调有序、运转高效的应急管理机制。
加强监测预警能力建设，加大投入，整合资源，提高雨情汛情旱情预报水平。
建立专业化与社会化相结合的应急抢险救援队伍，着力推进县乡两级防汛抗旱服务组织建设，健全应急抢险物资储备体系，完善应急预案。
建设一批规模合理、标准适度的抗旱应急水源工程，建立应对特大干旱和突发水安全事件的水源储备制度。
加强人工增雨（雪）作业示范区建设，科学开发利用空中云水资源。
继续推进农村饮水安全建设：
到２０１３年解决规划内农村饮水安全问题，“十二五”期间基本解决新增农村饮水不安全人口的饮水问题。
积极推进集中供水工程建设，提高农村自来水普及率。有条件的地方延伸集中供水管网，发展城乡一体化供水。
加强农村饮水安全工程运行管理，落实管护主体，加强水源保护和水质监测，确保工程长期发挥效益。
制定支持农村饮水安全工程建设的用地政策，确保土地供应，对建设、运行给予税收优惠，供水用电执行居民生活或农业排灌用电价格。
继续实施大江大河治理：
进一步治理淮河，搞好黄河下游治理和长江中下游河势控制，继续推进主要江河河道整治和堤防建设，加强太湖、洞庭湖、鄱阳湖综合治理，全面加快蓄滞洪区建设，合理安排居民迁建。搞好黄河下游滩区安全建设。
“十二五”期间抓紧建设一批流域防洪控制性水利枢纽工程，不断提高调蓄洪水能力。
加强城市防洪排涝工程建设，提高城市排涝标准。
推进海堤建设和跨界河流整治。
加强水资源配置工程建设：
完善优化水资源战略配置格局，在保护生态前提下，尽快建设一批骨干水源工程和河湖水系连通工程，提高水资源调控水平和供水保障能力。
加快推进南水北调东中线一期工程及配套工程建设，确保工程质量，适时开展南水北调西线工程前期研究。
积极推进一批跨流域、区域调水工程建设。
着力解决西北等地区资源性缺水问题。
大力推进污水处理回用，积极开展海水淡化和综合利用，高度重视雨水、微咸水利用。
搞好水土保持和水生态保护：
实施国家水土保持重点工程，采取小流域综合治理、淤地坝建设、坡耕地整治、造林绿化、生态修复等措施，有效防治水土流失。
进一步加强长江上中游、黄河上中游、西南石漠化地区、东北黑土区等重点区域及山洪地质灾害易发区的水土流失防治。
继续推进生态脆弱河流和地区水生态修复，加快污染严重江河湖泊水环境治理。
加强重要生态保护区、水源涵养区、江河源头区、湿地的保护。
实施农村河道综合整治，大力开展生态清洁型小流域建设。
强化生产建设项目水土保持监督管理。
建立健全水土保持、建设项目占用水利设施和水域等补偿制度。
合理开发水能资源：
在保护生态和农民利益前提下，加快水能资源开发利用。
统筹兼顾防洪、灌溉、供水、发电、航运等功能，科学制定规划，积极发展水电，加强水能资源管理，规范开发许可，强化水电安全监管。
大力发展农村水电，积极开展水电新农村电气化县建设和小水电代燃料生态保护工程建设，搞好农村水电配套电网改造工程建设。
强化水文气象和水利科技支撑：
加强水文气象基础设施建设，扩大覆盖范围，优化站网布局，着力增强重点地区、重要城市、地下水超采区水文测报能力，加快应急机动监测能力建设，实现资料共享，全面提高服务水平。
健全水利科技创新体系，强化基础条件平台建设，加强基础研究和技术研发，力争在水利重点领域、关键环节和核心技术上实现新突破，获得一批具有重大实用价值的研究成果，加大技术引进和推广应用力度。
提高水利技术装备水平。
建立健全水利行业技术标准。
推进水利信息化建设，全面实施“金水工程”，加快建设国家防汛抗旱指挥系统和水资源管理信息系统，提高水资源调控、水利管理和工程运行的信息化水平，以水利信息化带动水利现代化。
加强水利国际交流与合作。
加大公共财政对水利的投入：
多渠道筹集资金，力争今后１０年全社会水利年平均投入比２０１０年高出一倍。
发挥政府在水利建设中的主导作用，将水利作为公共财政投入的重点领域。
各级财政对水利投入的总量和增幅要有明显提高。
进一步提高水利建设资金在国家固定资产投资中的比重。
大幅度增加中央和地方财政专项水利资金。
从土地出让收益中提取１０％用于农田水利建设，充分发挥新增建设用地土地有偿使用费等土地整治资金的综合效益。
进一步完善水利建设基金政策，延长征收年限，拓宽来源渠道，增加收入规模。
完善水资源有偿使用制度，合理调整水资源费征收标准，扩大征收范围，严格征收、使用和管理。
有重点防洪任务和水资源严重短缺的城市要从城市建设维护税中划出一定比例用于城市防洪排涝和水源工程建设。
切实加强水利投资项目和资金监督管理。
加强对水利建设的金融支持：
综合运用财政和货币政策，引导金融机构增加水利信贷资金。
有条件的地方根据不同水利工程的建设特点和项目性质，确定财政贴息的规模、期限和贴息率。
在风险可控的前提下，支持农业发展银行积极开展水利建设中长期政策性贷款业务。
鼓励国家开发银行、农业银行、农村信用社、邮政储蓄银行等银行业金融机构进一步增加农田水利建设的信贷资金。
支持符合条件的水利企业上市和发行债券，探索发展大型水利设备设施的融资租赁业务，积极开展水利项目收益权质押贷款等多种形式融资。
鼓励和支持发展洪水保险。
提高水利利用外资的规模和质量。
广泛吸引社会资金投资水利：
鼓励符合条件的地方政府融资平台公司通过直接、间接融资方式，拓宽水利投融资渠道，吸引社会资金参与水利建设。
鼓励农民自力更生、艰苦奋斗，在统一规划基础上，按照多筹多补、多干多补原则，加大一事一议财政奖补力度，充分调动农民兴修农田水利的积极性。
结合增值税改革和立法进程，完善农村水电增值税政策。
完善水利工程耕地占用税政策。
积极稳妥推进经营性水利项目进行市场融资。
发博文
添加到我的博客