海洋要素计算(潮汐)

海洋要素计算作业之二——潮汐（威海2013年五月份）

一．本次潮汐调和分析共选取了十三个分潮：

MSf，Q1，O1，K1，P1，K2，N2，M2，S2，MK3，M4，MS4，M6

为使您查看方便，将本次大作业的放在本文件夹各文件内，具体参考如下：

1.原数据为：qd.dat;

2.Fortran编程见该文件夹内：tide.f90文件；

3.求各分潮调和常数H、g的值及其中间过程得到的各值见：qd_tide.dat文件；二．对比回报值和实测值：

1. 回报1968年一月份的水位值见：huibao.dat；

2. 用matlab绘制的潮汐过程曲线见：潮汐过程曲线.bmp

3. 用给定的六个分潮求得的高潮和低潮发生的时刻及潮位值见—：gaodichao.dat;

运行tide.f90后求得威海地区2013年5月份的平均潮差。

由图可知：由于只计算了一个月的潮汐数据，所以回报值和实测值相符的不是很好，如果计算一年的数据，应该会取得比较良好的结果。

三．程序

%% 潮汐过程曲线图

clear,clc

%%

huibao=load('G:\chaoxi\huibao.dat');

% huibao=fread(fhuibao);

shice=load('G:\chaoxi\qd.dat');

% shice=fread(fshice);

%huibao_y=zeros(1,12*62);

%shice_y=zeros(1,12*62);

huibao=double(huibao');

huibao_y=double(huibao(:));

%shice_y=reshape(shice',1,[])

%for i=1:12;

% for j=1:62

% huibao_y(i)=huibao(i,j)

% shice_y(i)=shice(i,j)

%end

%end

shice=double(shice');

shice_y=double(shice(:));

x=linspace(1,31,length(huibao_y));

plot(x,huibao_y,'r-')

hold on

plot(x,shice_y,'b-')

title('威海(37°31′N ,122°08′E)2013年五月潮汐调和分析图') legend('回报值','实测值')

xlabel('时间（2013年五月份）')

ylabel('水位（m）')