物理海洋数值模拟 .pdfVIP

GDOU-B-11-112

广东海洋大学学生实验报告书（学生用表）

实验名称海洋数值模拟实验课程名称海洋数值模式课程号

学院(系)海洋数值模拟专业物理海洋班级海洋1101

学生姓名学号实验地点科技楼实验日期

一.实验目的：

1.对海洋数值模拟的过程有一个较为初步的了解。

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

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

二．实验原理（差分格式）：

1.科氏力作用下的海水运动原理：

三．实验方法（程序代码）

：1科氏力作用下的海水运动：

（1）fortran程序代码：

（2）1）数值计算程序：

programcoriolisforce

implicitnone

integer,parameter::n=1000

real::u(n)

real::v(n)

real::dt(n)!choiceonedaydt=24*3600/1000

real::f(n)!f=2*omig*sin(fai)fai取0~60

real::alf(n)

real::beta(n)

real::fai(n)

integeri,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)

!***************数据计算********************************************

dok=2,n

dt(k)=dt(k-1)+24.0/1000.0

fai(k)=fai(k-1)+60.0/1000.0

enddo

doi=1,n

f(i)=0.0042*2.0*sin(fai(i))

alf(i)=2.0*asin(0.5*dt(i)*f(i))

enddo

doj=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)

enddo

!===========================================

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°的海水运

动，时间为一天。