中期水位资料对潮汐进行调和分析

《中期水位资料对潮汐进行调和分析》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、!利用1996年7月厦门站的潮汐观测数据计算调和常数，并利用主要分潮和浅水分潮进行潮汐预报programworkimplicitnonecharacter*80::alcharacter(len=5)，dimension(62,16)::aainteger::bb(62,12),c(62,2),caita(-371:37l)，i，i1，i2,j,t1real::N0，n(13,6)，a(0:13,0:13)，b(l:13，l:13),s，s0，sl，s2，s3,sa，hh!n代表Doodson代码；a,b为

2、系数矩阵real::xiaoa(0:13),xiaob(l3)，ggl，gg2，pjchaocha，t，ma，mi!计算法方程所需的参数real，dimension(1:13)::w，u，f，V0，f1(0:13)，f2!fl和f2为法方程右边系数real，dimension(13)::sita，h，g，r，hO(13)，g0(13)，h1(13)，g1(13)!调和常数参数real,dimension(-371:371)::caital,caita3,caita4,caita8,caita9,caita5

3、,caitall!主分潮、浅水分潮的潮高数值real，dimension(:)，allocatable::hightide，lowtide，chaocha!高低潮数值integer,dimension(:),allocatable::hightrq,lowtrq,hight,hightt,lowt,lowtt!读取数据，把潮位数据赋值给bb，把年月份数据赋值给copen(unit=2，file=’XM_July1996.dat')read(2；(a)*)alprint*，*数据文件的第一行信息:'，aldo

4、i=l，62read(2，’(16a5)’)aa(i，:)enddodoi=l，62read(aa(i，5:16)，*)bb(i，:)read(aa(i，3:4)，*)c(i，:)enddodoi=l，62c(i，2)=int(real(c(i，2))/10.0)enddoclose(2)!计算分潮角速率ww=(/0.002822,0.037219,0.038731,0.041781,0.163845,0.241534,0.078999&0.080511,0.083333,0.122292,0.161023

5、,0.041553,0.083561/)w=360*wprint*print*,*角速率w:w!计算NO(middletime:1996-7-16;datasum:744,middlenumber:372)N0=259.157-19.32818*(1996-1900)-0.05295*(3P3+30*2+29+15+int((95.0)/4.0))!初始升交点平均黄经N0=-(0.00220641*3+N0)print*!转换成格林威治时间print*，’NO:’,NO!数字序号对应选取的分潮，但将5、6

6、(P1、K2)分别与12、13(MS4、M6)对调，其中Pl、K2为随从分潮!计算交点订正角uu(3)=10.8*sind(N0)-1.34*sind(2*N0)+0.19*sind(3*N0)u(4)=-8.86*sind(N0)+0.68*sind(2*N0)-0.07*sind(3*N0)u(8)=-2.14*sind(N0)u(13)=-17.74*sind(N0)+0.68*sind(2*N0)-0.04*sind(3*N0)u(l)=-u(8)u(2)=u(3)u(7)=u(8)u(10)=u(

7、8)+u(4)u(ll)=2*u(8)u(5)=u(8)u(6)=3*u(8)u(12)=0!print*print*/交点订正角u:u!计算交点因子ff(3)=1.0089+0.1871*cosd(N0)-0.147*cosd(2*N0)+0.0014*cosd(3*N0)f(4)=1.006+0J15*cosd(N0)-0.0088*cosd(2*N0)+0.0006*cosd(3*N0)f(8)=1.0004-0.0373*cosd(N0)+0.0003*cosd(2*N0)f(13)=1.0241

8、+0.2863*cosd(N0)+0.0083*cosd(2*N0)-0.0015*cosd(3*N0)f(l)=f(8)f(2)=f(3)f(7)=f(8)f(9)=lf(10)=f(8)*f(4)f(ll)=f(8)**2f(6)=f(8)**3f（12）=l!print*print*/交点因子ff，f!查表得到的Doodson代码n（l，:）=（/0,2，-2,0,0,0/）n（2，:）=（/l，-2,0，