地球流体力学第三章各种波的图形.pdf
( 4.5 )《地球流体力学第三章各种波的图形.pdf》由会员分享，可在线阅读，更多相关《地球流体力学第三章各种波的图形.pdf（10页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、重力惯性波重力惯性波FortranFortran 程序程序Program my_001IMPLICIT NoneINTEGER,PARAMETER:M=16,n=2,Lx=16,Ly=16,T=6REAL,PARAMETER:pi=3.1415926REAL,PARAMETER:f=pi/6REAL x(M),y(M),sigmaREAL eta(M,M,6),u(M,M,6),v(M,M,6)INTEGER i,j,kdo i=1,Mx(i)=(i-1)*pi*2/Lxy(i)=(i-1)*pi*2/Lyenddosigma=2*pi/Tdo k=0,5do j=1,Mdo i=1,Meta
2、(i,j,k+1)=cos(x(i)+y(j)-sigma*k)u(i,j,k+1)=sigma*cos(x(i)+y(j)-pi*k)-f*sin(x(i)+y(j)-sigma*k)v(i,j,k+1)=sigma*cos(x(i)+y(j)-pi*k)+f*sin(x(i)+y(j)-sigma*k)enddo;enddo;enddoopen(11,file=zlgxb.grd,form=binary)do k=1,6write(11) (eta(i,j,k),i=1,M),j=1,M)write(11) (u(i,j,k),i=1,M),j=1,M)write(11) (v(i,j,k
3、),i=1,M),j=1,M)enddoclose(11)endgsgs 文件文件REINITopen D:myggradszlgxb.ctlenable print D:myggradszlgxb.pngset x 1 16set y 1 16set z 1set mproj offset grads offk=1while(k=6)set t kd etaset gxout vectorset arrscl 0.5 1.5d u;vdraw title t= kprintck=k+1endwhiledisable printreinit图图邦加莱波邦加莱波FortranFortran 程序
4、程序Program my_002IMPLICIT NoneINTEGER,PARAMETER:M=16,n=2,Lx=16,Ly=16,T=6REAL,PARAMETER:pi=3.1415926REAL,PARAMETER:f=pi/6REAL x(M),y(M),sigmaREAL eta(M,M,6),u(M,M,6),v(M,M,6)INTEGER i,j,kdo i=1,My(i)=(i-1)*pi/(M-1)x(i)=(i-1)*pi*2/Lxenddosigma=2*pi/Tdo k=0,5do j=1,Mdo i=1,Meta(i,j,k+1)=cos(y(j)*cos(x(i
5、)-sigma*k)u(i,j,k+1)=cos(y(j)*cos(x(i)-sigma*k)v(i,j,k+1)=-sin(y(j)*sin(x(i)-sigma*k)enddo;enddo;enddoopen(11,file=bjlb.grd,form=binary)do k=1,6write(11) (eta(i,j,k),i=1,M),j=1,M)write(11) (u(i,j,k),i=1,M),j=1,M)write(11) (v(i,j,k),i=1,M),j=1,M)enddoclose(11)endgsgs 文件文件和重力惯性波一样，此处略。图图开尔文波南边界开尔文波南边界
- 配套讲稿：
如PPT文件的首页显示word图标，表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制：
部分文档作品中含有的国旗、国徽等图片，仅作为作品整体效果示例展示，禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词：
- 地球 流体力学 第三 各种 图形
淘文阁 - 分享文档赚钱的网站所有资源均是用户自行上传分享，仅供网友学习交流，未经上传用户书面授权，请勿作他用。
限制150内