Fortran 读取程序和求解潮汐各分潮的调和常数 范例

发布时间：2023-01-10 20:45:30


!程序中数组的各元素依次为分潮
MSf,Q1,O1,K1,P1,K2,N2,M2,S2,MK3,M4,MS6,M6
programtiaohechangshu
implicitnone
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=(/14.49205212,0.54901653,0.04106864,&
0.00464183,0.00220641,0.00000196/
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格式的文件

Fortran 读取程序和求解潮汐各分潮的调和常数 范例

推荐内容

相关推荐