中期观测资料十三个分潮的调和分析代码

一首简单的歌

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
!=============================================================================================

mofangbao

介绍最好能再详细些~

平流层的萝卜

20多k的代码，好高端的样子！

lqouc

平流层的萝卜 发表于 2014-4-9 20:37
20多k的代码，好高端的样子！

确实看起来很炫，光变量声明就拿一大堆，不过潮汐什么的貌似很多都是用matlab做的，使用某些函数包会极大地简化。

小阿呆

看起来挺不错的，下来看看

vickycomic

输公式和迭代简直是一大硬伤T.T谢谢楼主分享~~

two_Drifters

我的天，为什么潮汐调和分析的帖子都要贡献值啊，没有怎么破

two_Drifters

我的天，为什么潮汐调和分析的帖子都要贡献值啊，没有怎么破

li1015271178

分析完潮汐类型后的程序没运行 怎么回事

fuqingjun

怎样增加贡献值啊？