计算程序

计算程序,第1张

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”麻烦解释为什么。等相关内容解答,如果想了解更多相关内容,可以关注我们,你们的支持是我们更新的动力!

欢迎分享,转载请注明来源:内存溢出

原文地址:https://54852.com/zz/9320109.html

(0)
打赏 微信扫一扫微信扫一扫 支付宝扫一扫支付宝扫一扫
上一篇 2023-04-27
下一篇2023-04-27

发表评论

登录后才能评论

评论列表(0条)

    保存