资源描述:
《《河流模拟课设》word版》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、programmainparameter(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)dimensionx(mm,2,nn),rough(nn),dx(nn),zlevel(nn),Q(nn),Npoint(nn),b(nn),Q0(nn),QQ(nd)dimensiona(nn),
2、xw(nn),dxa(nn),s(nn),alow0(nn),alow(nn),h(nn)dimensionu(nn),sx(nn),gb(nn),dy(nn),nday(ny)dataNday/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)91format(3x,'断面号',6x'距坝里程(千米)',3x,'
3、深泓(米)')open(13,file='水面线.txt',status='unknown')write(13,92)92format(3x,'断面号',6x'距坝里程(千米)',3x,'水面线高程(米)')open(14,file='初始水位库容.txt',status='unknown')write(14,93)93format(3x,'年',6x'水位(米)',3x,'库容(亿立方米)')open(15,file='年淤积总量.txt',status='unknown')write(15,94)94format(3x'年'
4、,3x'淤积总量(万立方米)',3x'年径流量(万方)',3x'累计输沙量(万吨)',3x'年均流量(万方)')open(16,file='深泓(十年).txt',status='unknown')write(16,95)95format(3x'年',3x'断面号',3x'深泓(米)')open(17,file='坝前断面(十年).txt',status='unknown')write(17,96)96format(3x'年',3x'起点距(米)',3x'高程(米)')open(18,file='水位库容(十年).txt',st
5、atus='unknown')write(18,97)97format(3x,'年',6x'水位(米)',3x,'库容(亿立方米)')!===========================================================!======================读入地形数据=========================!===========================================================doi=1,npxtread(10,*)read(1
6、0,*)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))enddo!===========================================================!======================读入糙率=============================!====================================
7、=======================doi=1,npxtrough(i)=0.06enddo!===========================================================!====================断面间距计算===========================!===========================================================doi=1,30dx(i)=dxa(i)-dxa(i+1)enddo!=======
8、====================================================!====================判断读入数据的准确性===================!============================