广东海洋大学学生实验报告书（学生用表） 实验名称 海洋数值模拟实验 课程名称
海洋数值模式 课程号 学院(系) 海洋数值模拟 专业
物理海洋 班级 海洋1101 学生姓名 学号 实验地点 科技楼 实验日期
一. 实验目的：
1.对海洋数值模拟的过程有一个较为初步的了解。

2.熟悉模式编程语言以及操作系统。

3.熟悉常用的数值模式的特点。

二．实验原理（差分格式）：
1. 科氏力作用下的海水运动原理：
三．实验方法（程序代码）：1科氏力作用下的海水运动：
（1）fortran 程序代码：
（2）1）数值计算程序：
program coriolisforce GDOU-B-11-112
implicit none
integer,parameter:: n=1000
real:: u(n)
real:: v(n)
real:: dt(n) !choice one day dt=24*3600/1000 real:: f(n) !f=2*omig*sin(fai) fai取0~60 real:: alf(n)
real:: beta(n)
real:: fai(n)
integer i,j,k
u(1)=0.51
u(2:n)=0.0
v(1)=0.499
v(2:n)=0.0
dt(1)=0.0
dt(2:n)=0.0
fai(1)=2.0
fai(2:n)=0.0
alf(1:n)=0.0
f(1:n)=0.0
open(unit=10,file='no4.txt')
write(10,*) u(1),v(1)
!***************数据计算******************************************** do k=2,n
dt(k)=dt(k-1)+24.0/1000.0
fai(k)=fai(k-1)+60.0/1000.0
end do
do i=1,n
f(i)=0.0042*2.0*sin(fai(i))
alf(i)=2.0*asin(0.5*dt(i)*f(i))
end do
do j=2,n
u(j)=cos(alf(j-1))*u(j-1)+sin(alf(j-1))*v(j-1)
v(j)=cos(alf(j-1))*v(j-1)-sin(alf(j-1))*u(j-1)
write(10,*) u(j),v(j)
end do
!===========================================
end
2）绘图程序：
clear,clc
V=load('no4.txt');
u=V(:,1);
v=V(:,2);
plot(u,v)
title('北纬0到60度科氏力作用下的海水运动')
xlabel('速度u')
ylabel('速度v')
四．实验结果（作图）
1.科氏力作用下海水从北纬0°到北纬60°的海水运动，时间为一天。