GDOU-B-11-112廣東海洋大學學生實驗報告書（學生用表）實驗名稱海洋數(shù)值模擬實驗 課程名稱海洋數(shù)值模式課程號學院（系）海洋數(shù)值模擬 專業(yè)物理海洋 班級海洋1101學生姓名學號— 實驗地點科技樓實驗日期實驗目的:對海洋數(shù)值模擬的過程有一個較為初步的了解。熟悉模式編程語言以及操作系統(tǒng)。熟悉常用的數(shù)值模式的特點。實驗原理（差分格式）：1.科氏力作用下的海水運動原理：=utl+△”伊andvfl+l=驢一AtfunImprovedScheme2:TheLocal*RotationApproachTheCoriolisforceoperates制arightiingletovelocityandduesnotchangethe

tif[noiion,onlythedirection,Fieneeh fejHitrecanKesinmlatedbyalocn]notationofthevelocitysector^〔hati%ufl+1=cos(a)uri+sin(ff)irt,成十’=coy值｝『一sin(tv)wriJtuitigcoineiricconsiderations,thelot&tionanglecanbedetetnuaedaiw=2UJTsiaifO^Aff).FsAjf\<<1,ihistynbeapproxlmatedbyw&Af/,三.實驗方法(程序代碼)：1科氏力作用下的海水運動：(1)fortran程序代碼：⑵1)數(shù)值計算程序：programcoriolisforceimplicitnoneinteger,parameter::n=1000real::u(n)real::v(n)real::dt(n) !choiceonedaydt=24*3600/1000real::f(n) !f=2*omig*sin(fai)fai取0~60real::alf(n)real::beta(n)real::fai(n)integeri,j,ku(1)=0.51u(2:n)=0.0v(1)=0.499v(2:n)=0.0dt(1)=0.0dt(2:n)=0.0fai(1)=2.0fai(2:n)=0.0alf(1:n)=0.0f(1:n)=0.0open(unit=10,file='no4.txt')write(10,*)u(1),v(1)!***************數(shù)據計算********************************************dok=2,ndt(k)=dt(k-1)+24.0/1000.0fai(k)=fai(k-1)+60.0/1000.0enddodoi=1,nf(i)=0.0042*2.0*sin(fai(i))alf(i)=2.0*asin(0.5*dt(i)*f(i))enddodoj=2,nu(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!===========================================end2)繪圖程序:clear,clcV=load('no4.txt');u=V(:,1);v=V(:,2);plot(u,v)title('北緯0到60度科氏力作