有哪些信誉好的足球投注网站- 1、本文档共4页，可阅读全部内容。
- 2、有哪些信誉好的足球投注网站（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
Fortran 读取程序和求解潮汐各分潮的调和常数 范例.doc
!程序中数组的各元素依次为分潮MSf,Q1,O1,K1,P1,K2,N2,M2,S2,MK3,M4,MS6,M6
program tiaohechangshu
implicit none
integer,dimension(78)::n=(/0,2,-2,0,0,0,1,-2,0,1,0,0,1,-1,0,0,0,0,
1,1,0,0,0,0,1,1,-2,0,0,0,2,2,0,0,0,0,2,-1,0,1,0,0,2,0,0,0,0,0,2,2,-2,0,0,0,3,
1,0,0,0,0,4,0,0,0,0,0,4,2,-2,0,0,0,6,0,0,0,0,0/) !变量n中从左开始每隔六个数为一个分潮的杜德森数（n1,n2,n3,n4,n5,n6
real,dimension(6)::o=(/1400
000)
integer::y=1996,Y1=int((1996-1901)/4),D=198,t=3,i,j,ot,N2,g,k,g1 !经计算观测数据的中间时刻为1996年7月16日11时,转化成格林时间为t-8
real,dimension(13)::OO,V0,f,u,b,b1,R,C,H1,gg !其中变量OO、V0、f、u、H1、gg分别为各分潮的角速率、幅角、交点因子、交点订正角、调和函数（H1和gg）
real::T0,s0,h0,p0,N0,PP0,TT,s,h,p,NN,PP,N1,fO1,fK1,fM2,fK2,uO1,uK1,uM2,uK2
real,dimension(0:13,0:13)::AA,AA1
real,dimension(0:13,0:14)::AAA
real,dimension(13,13)::BB,BB1
real,dimension(13,14)::BBB
real,dimension(0:13)::a1,a
real,dimension(0:13,1)::F1,F11
real,dimension(13,1)::F2,F22
integer,dimension(62)::year,month,day
integer,dimension(1:24)::time1=(/00,01,02,03,04,05,06,07,08,09,
10,11,12,13,14,15,16,17,18,19,20,21,22,23/)
integer,dimension(62,12)::hr1
integer,dimension(-371:371)::hr !!!经计算N=(743-1)/2=371
!从(XM_July1996.dat)中读数据并创建.txt格式的文件
open(8,file=F:\Fortran90\XM_July1996.dat)
read(8,*)
do i=1,62
read(8,(11x,i4,i2,i2,x,12i5)),year(i),month(i),day(i),hr1(i,:)
end do
close(8)
open(4,file=F:\Fortran90\XM_July1996lie.txt)
do i=1,62
do j=1,12
if (mod(i,2)==0) then
write(4,(i4,i1,i2,i2,i5)),year(i),month(i),day(i),time1(12+j),hr1(i,j)
else
write(4,(i4,i1,i2,i2,i5)),year(i),month(i),day(i),time1(j),hr1(i,j)
end if
end do
end do
close(4)
open(5,file=F:\Fortran90\XM_July1996lie.txt)
do j=-371,371
read(5,(10x,i5))hr(j)
end do
close(5)
!!!各天文参量每日Greenwish零时的值
T0=180.0
s0=277.025+129.38481*(y-1900)+13.17640*(D+Y1)
h0=280.190-0.23872*(y-1900)+0.98565*(D+Y1)
p0=334.385+40.66249*(y-1900)+0.11140*(D+Y1)
N0=259.157-19.32818*(y-1900)-0.05295*(D+Y1)
PP0=281.221+0.01718*(y-1900)+0.000047*(D+Y1)
!!!计算某时刻的角速率
TT=o(1)*t+T0
s=o(2)*t+s0
h=o(3
文档评论（0）