原创力文档- 1、本文档共3页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
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°的海水运
动，时间为一天。
文档评论（0）