- 积分
- 45
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-4-2
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
program zz123
implicit none
character aaa !读取数据文件中文字
integer i,j !循环变量
integer::N=743 !N为数据量
real,parameter::pi=3.1415926
integer::y=1997,Y1=int((1997-1901)/4),D=227,N2,g1,t=0,f0 !D为从1月1日起算到中间时刻经过的日期
integer::um(6)=(/3,4,6,7,9,10/) !回报中用到的主要分潮和浅水分潮的位置
integer,parameter::c=11,s=2 !c=主分潮数 s=随从分潮数
real,dimension(13)::OO,sigema,V0,VV0,f,ff,u,H1,g !变量OO、V0、f、u、H1、gg分别为各分潮的角速率、幅角、交点因子、交点订正角和调和函数(H1和g)
real,dimension(13)::VV !13个分潮的理论初相角
integer::na=23,L !na为随从分潮移到方程右边之后的方程阶数,L为子程序返回值(高斯赛德尔迭代)
dimension A(0:13,0:13),B(1:13,1:13),FA(0:11),FB(11),& !正规法方程的元素,A,B,F=FA,G=FB
EQUA(0:22,0:22),EQUB(0:22),EQUC(0:22),last_EQUB(0:22) !正规法方程的系数矩阵,解向量和右端项
double precision A,B,FA,FB,EQUA,EQUB,EQUC,last_EQUB
real::a12,a13,b12,b13,k12,k13,Q12,Q13 !定义随从分潮的(a,b)以及(差比关系)
real::EPS=1e-9 ! 误差限
real::T0,t00,s0,h0,p0,N0,PP0,tttt,s_,h,p,NN,PP,N1 !t=0时刻的天文参量 及t=中间时刻的天文参量
real,dimension(13,2)::fu !分别存放13个分潮的f和u,顺序依次为 MSf、Q1、O1、K1、N2、M2、S2、MK2、M4、MS4、M6、P1、K2
real,dimension(:)::V(13,3),fv(13,2)
integer,dimension(743)::tt
integer,dimension(:),allocatable::pred,tideleave
integer,dimension(1:743)::tide0 !除去最后一个的水位数据
real,dimension(-371:371)::tide
real::sss=0 !平均水位
real::typ !潮型数
integer,dimension(78)::dood=(/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/) !十三个分潮的杜德森数,依次为MSf,Q1,O1,K1,P1,K2,N2,M2,S2,MK3,M4,MS4,M6
!=====================================该段数据为各平太阳时====================================
real :: St=15
real :: Sq=14.49205212
real :: Ss=0.54901653
real :: Sh=0.04106864
real :: Sp=0.00464183
real :: Sn=0.00220641
real :: Sppie=0.00000196
!=============================================================================================
|
-
-
tide.f90
20.97 KB, 阅读权限: 10, 下载次数: 42, 下载积分: 金钱 -5
售价: 1 贡献 [记录]
加入差比关系的潮汐调和分析
评分
-
查看全部评分
|