原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
河流模拟课设
program main
parameter (nn=31,mm=80,nd=3653,ny=10,fai=1.0,ndisp=0,npxt=31,fai1=0.5)
parameter (rs=2650*9.8,r=9800,w=0.0012,rou0=1325.0,dt=86400.0,d=0.002,rou=1000,rous=2650)
dimension x(mm,2,nn),rough(nn),dx(nn),zlevel(nn),Q(nn),Npoint(nn),b(nn),Q0(nn),QQ(nd)
dimension a(nn),xw(nn),dxa(nn),s(nn),alow0(nn),alow(nn),h(nn)
dimension u(nn),sx(nn),gb(nn),dy(nn),nday(ny)
data Nday /365,731,1096,1461,1826,2192,2557,2922,3287,3653/
open(10,file=地形.txt,status=old)
open(12,file=深泓.txt,status=unknown)
write(12,91)
91 format(3x,断面号,6x距坝里程(千米),3x,深泓(米))
open(13,file=水面线.txt,status=unknown)
write(13,92)
92 format(3x,断面号,6x距坝里程(千米),3x,水面线高程(米))
open(14,file=初始水位库容.txt,status=unknown)
write(14,93)
93 format(3x,年,6x水位(米),3x,库容(亿立方米))
open(15,file=年淤积总量.txt,status=unknown)
write(15,94)
94 format(3x年,3x淤积总量(万立方米),3x年径流量（万方）,3x累计输沙量（万吨）,3x年均流量（万方）)
open(16,file=深泓（十年）.txt,status=unknown)
write(16,95)
95 format(3x年,3x断面号,3x深泓(米))
open(17,file=坝前断面（十年）.txt,status=unknown)
write(17,96)
96 format(3x年,3x起点距(米),3x高程(米))
open(18,file=水位库容（十年）.txt,status=unknown)
write(18,97)
97 format(3x,年,6x水位(米),3x,库容(亿立方米))
!===========================================================
!======================读入地形数据=========================
!===========================================================
do i=1,npxt
read(10,*)
read(10,*)n1,npoint(i),dx(i),d1,dxa(i),rough(i)
read(10,*)(x(j,1,i),j=1,npoint(i))
read(10,*)(x(j,2,i),j=1,npoint(i))
end do
!===========================================================
!======================读入糙率=============================
!===========================================================
do i=1,npxt
rough(i)=0.06
enddo
!===========================================================
!====================断面间距计算===========================
!===========================================================
do i=1,30
dx(i)=dxa(i)-dxa(i+1)
end do
!===========================================================
!====================判断读入数据的准确性===================
!===========
文档评论（0）