地球流体力学第三章各种波的图形.docxVIP

重力惯性波Fortran程序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(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,6 write(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,