海洋要素计算作业之二——潮汐（威海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）')