
c name:qun-yuefor,无界承压含水层干扰井群法预测开采降深场,有潜水越流进入;当有界时,则把虚井同实井一样代入即可(但其Q却有正负之分)
c 计算公式:第一越流系统计算公式——汉图士-雅可布公式
c 即:s(j)=∑QiF(ui,r/B)/(4πT),u(i)=Sr(i,j)2/(4Tt)
c B=sqrt(Tm’/k’)r/B=r/sqrt(Tm’/k’)越流系数=k’/m’
c 其中越流井函数F(u,r/B),本程序自动查表求其值;
c 均质,各向同性s:降深(m),S:储水系数;m:预测点个数n:抽水井个数(包括实虚)ihv:预测时段个数;T:承压水全区T值(mm/d),W:承压水全区S值,bk1m:弱透水层越流系数,即K’/m’(1/d)
parameter(m=4,n=2,ihv=40,ihv1=34)
double precision u1
c 抽水井的x、y坐标(单位:m)预测点的x、y坐标(单位:m),预测点可以是网格状剖分节点加非节点的预测点
dimension xn(n),yn(n),r(n,m)
dimension xm(m),ym(m)
dimension q(n),s(m,ihv),tt(ihv)
dimension bu(40),brb(30),bf(40,30)
character file3(m)8,fileb8,file2(m)4
cmm:越流井函数F(u,r/B)表列数,nn:越流井函数F(u,r/B)表行数
mm=28
nn=38
open(1,file=’biaoftxt’)
read(1,)(brb(j),j=1,mm)
do 210 i=1,nn
read(1,)bu(i),(bf(i,j),j=1,mm)
210 continue
close(1)
c 抽水井点(n眼)坐标
open(1,file=’xy-jingtxt’)
do 9 i=1,n
read(1,)xn(i),yn(i)
9 continue
close(1)
c 预测点(m个)坐标
open(1,file=’xy-yucetxt’)
do 10 i=1,m
read(1,)xm(i),ym(i)
10 continue
close(1)
open(1,file=’kmwbtxt’)
read(1,)ak,am,w,bk,bm
close(1)
T=akam
bk1m=bk/bm
b=sqrt(T/bk1m)
open(1,file=’qtxt’)
read(1,)(q(i),i=1,n)
close(1)
c’yucehaotxt’存放m个字符串,每串4字符,为预测点输出预测的历时水位曲线标注字符用
open(1,file=’yucehaotxt’)
read(1,42)(file2(i),i=1,m)
close(1)
42 format(5a4)
c’yucewentxt’存放m个字符串,为"预测号dat",每串8字符,为预测点输出预测的历时水位曲线准备文件名
open(1,file=’yucewentxt’)
read(1,44)(file3(i),i=1,m)
close(1)
44 format(5a8)
c’tl1’-预测时间系列:累计时间(d)
open(1,file=’tl1txt’)
read(1,)(tt(i),i=1,ihv)
close(1)
tt0=0
c 各预测点j至抽水井点i之距r(m)
do 30 j=1,m
s(j,0)=0
do 20 i=1,n
r(i,j)=sqrt((xn(i)-xm(j))2+(yn(i)-ym(j))2)
if(r(i,j)eq0)r(i,j)=015
20 continue
30 continue
close(1)
open(1,file=’sidstxt’)
do 60 k=1,ihv1
write(,)’k=’,k,’t=’,tt(k)
write(1,70)k,tt(k)
70 format(/2x,’K=’,i3,’t=’,f156,’日’)
c 下一大段对m个预测点逐点计算1~n眼抽水井在k时刻对每个预测点水位降之和s(j,k)
do 50 j=1,m
s(j,k)=0
c 计算r/B,及k时刻i~j泰斯井函数自变量u1
do 40 i=1,n
u1=r(i,j)2w/(40Ttt(k))
rb=r(i,j)/b
c 限制u1超出8
if(u1gt8)u1=8
c 限制r/B超出9
if(rbgt9)rb=9
c 越流井函数F(u,r/b),程序自查表求之
c 程序自查表,先查列,查计算值rb位于F表的哪两列brb(jj-1)~brb(jj)之间
do 220 jj=2,mm
jx=jj
if(rbgebrb(jj-1)andrblebrb(jj))goto 230
220 continue
c 计算出F表的brb(jj-1)~brb(jj)两列之间的列距dx
230 dx=brb(jx)-brb(jx-1)
c 计算值rb到F表brb(jj-1)列之距dx1、到F表brb(jj)列之距dx2
dx1=rb-brb(jx-1)
dx2=brb(jx)-rb
c 程序自查表,再查行,查计算值u1位于F表的哪两行bu(ii-1)~bu(ii)之间
do 240 ii=2,nn
iy=ii
if(u1gebu(ii-1)andu1lebu(ii))goto 250
240 continue
c 计算出F表的bu(ii-1)~bu(ii)两行之间的行距dy
250 dy=bu(iy)-bu(iy-1)
c 计算值u1到F表bu(ii-1)行之距dy1、到F表bu(ii)行之距dy2
dy1=u1-bu(iy-1)
dy2=bu(iy)-u1
c 找出计算值(u1,rb)在F表中最近四点的F值(即bf),其中,f11为左上角,f12右上角,f21左下角,f22右下角
f11=bf(iy-1,jx-1)
f12=bf(iy-1,jx)
f21=bf(iy,jx-1)
f22=bf(iy,jx)
f1=f11+dx1/dx(f12-f11)
f2=f21+dx1/dx(f22-f21)
f=f1+dy1/dy(f2-f1)
c 用越流泰斯公式计算k时刻第i眼抽水井对第j个预测点的水位降ss
ss=q(i)/(4031415926T)f
c 把第i眼(i=1~n)抽水井在k时刻对第j预测点水位降ss累加到j点的迭加降深s(j,k)中
s(j,k)=s(j,k)+ss
40 continue
50 continue
do 100 j=1,m 100 write(1,80)xm(j),ym(j),-s(j,k),j
60 continue
close(1)
80 format(2f101,f93,i4)
do 300 j=1,m
fileb=file3(j)
open(1,file=fileb)
write(1,310)tt0,s(j,0),file2(j)
do 320 k=1,ihv1
write(1,380)tt(k),-s(j,k)
320 continue
300 close(1)
380 format(f101,f82)
310 format(f101,f82,1x,a4)
cq=0
do 90 i=1,n
90 cq=cq+q(i)
open(1,file=’qqtxt’)
write(1,140)(i,q(i),i=1,n)
write(1,)
write(1,)’总采量=’,cq,’(m3/d)’
close(1)
140 format(5(i4,f80))
open(1,file=’okkdat’)
do 290 j=1,m
290 write(1,80)xm(j),ym(j),-s(j,ihv1),j
close(1)
stop
end
不需要。在ff14买了rb之后,插件已经齐备,就不需买单独插件了。插件是一种遵循一定规范的应用程序接口编写出来的程序。其只能运行在程序规定的系统平台下,而不能脱离指定的平台单独运行。
C语言中以追加方式打开一个文件应选择“a”参数。
打开文件之后,程序可以得到文件的相关信息,例如大小、类型、权限、创建者、更新时间等。在后续读写文件的过程中,程序还可以记录当前读写到了哪个位置,下次可以在此基础上继续 *** 作。
扩展资料:
C语言中以追加方式打开一个文件,不同的 *** 作需要不同的文件权限。例如,只想读取文件中的数据的话,“只读”权限就够了;既想读取又想写入数据的话,“读写”权限就是必须的了。
另外,文件也有不同的类型,按照数据的存储方式可以分为二进制文件和文本文件,它们的 *** 作细节是不同的。
C选项。
“文件指针名”必须是被说明为FILE 类型的指针变量;
“文件名”是被打开文件的文件名;
“使用文件方式”是指文件的类型和 *** 作要求。
“文件名”是字符串常量或字符串数组。
例如:
FILE fp;
fp=("file a","r");
其意义是在当前目录下打开文件file a,只允许进行“读” *** 作,并使fp指向该文件。
又如:
FILE fphzk
fphzk=("c:\\hzk16","rb")
其意义是打开C驱动器磁盘的根目录下的文件hzk16,这是一个二进制文件,只允许按二进制方式进行读 *** 作。两个反斜线“\\ ”中的第一个表示转义字符,第二个表示根目录。
扩展资料:
C语言包含的各种控制语句仅有9种,关键字也只有32 个,程序的编写要求不严格且以小写字母为主,对许多不必要的部分进行了精简。实际上,语句构成与硬件有关联的较少,且C语言本身不提供与硬件相关的输入输出、文件管理等功能,如需此类功能,需要通过配合编译系统所支持的各类库进行编程,故c语言拥有非常简洁的编译系统。
c语言的处理和表现能力都非常的强大,依靠非常全面的运算符和多样的数据类型,可以轻易完成各种数据结构的构建,通过指针类型更可对内存直接寻址以及对硬件进行直接 *** 作,因此既能够用于开发系统程序,也可用于开发应用软件。
参考资料来源:百度百科-C语言
以上就是关于计算程序全部的内容,包括:计算程序、ff14买了rb还要买单独插件吗、C语言中以追加方式打开一个文件应选择_()A、“r” B、“w” C、“rb” D、“a”麻烦解释为什么。等相关内容解答,如果想了解更多相关内容,可以关注我们,你们的支持是我们更新的动力!
欢迎分享,转载请注明来源:内存溢出
微信扫一扫
支付宝扫一扫
评论列表(0条)