二维潮波数值模拟.docVIP

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! 二维潮波数值模拟与预报 ! 计算海域包括渤海和黄海北部（117.50 -126.67 E,34.00 - 41.00 N） ! 空间分辨率是10x10 ! 开边界处的h和g由TMD TOOLBOX求得 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! program hw4 implicit none interface subroutine h_get(h,dx,dy,dt,u,v,d,cd) integer::i,j real,intent(in)::dt,dy integer,dimension(43,56),intent(in)::cd real,dimension(43,56),intent(in)::d,u,v real,dimension(43,56),intent(inout)::h real,dimension(43),intent(in)::dx end subroutine h_get subroutine v_get(v,u,cu,cv,dx,dy,dt,d,h,f) integer::i,j real::kk=1.5e-3,gg=9.8 real,intent(in)::dt,dy integer,dimension(43,56),intent(in)::cu,cv real,dimension(43,56),intent(in)::d,h,u real,dimension(43,56),intent(inout)::v real,dimension(43),intent(in)::dx,f real,dimension(43,56)::uu,ss end subroutine v_get subroutine u_get(u,v,cu,cv,dx,dy,dt,d,h,f) integer::i,j real::kk=1.5e-3,gg=9.8 real,intent(in)::dt,dy integer,dimension(43,56),intent(in)::cu,cv real,dimension(43,56),intent(in)::d,h,v real,dimension(43,56),intent(inout)::u real,dimension(43),intent(in)::dx,f real,dimension(43,56)::vv,rr end subroutine u_get end interface !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! PART 1 ! dx, dy : 网格宽度，单位（m） ! dt : 时间步长，取为M2分潮周期的1/240，单位（s） ! d : 各个网格点的水深 ! h : 海水相对静止海面的高度 ! u, v : 海水水平流速 ! uu, vv : u、v的平均值 ! f : 科氏参量 ! cd : 水深控制场 ! cu, cv : 速度控制场 ! h0, g0 : 开边界处的振幅和迟角 ! h_check : 用于检验迭代是否达到要求的精度 ! eps : 迭代精度 ! PART 2 ! h_3d : 迭代完成后得到的水位 ! a,b,g,hh : 调和常数 ! l1,l2,l3,l4,y1,y2 :最小二乘法正规方程组的系数 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! real,parameter::pi=3dy=6371004*2*pi/360/6 real,parameter::dt=24.8412/2*3600/240,w=360/(240*dt),eps=0.01 !M2 !real,parameter::dt=24.8412*3600/240,w=360/(240*dt),eps=0.01 !K1 integer::i,j,k,m,n,t integer,dimension(43,56)::cd,cu,cv real,dimension(43,56)::d,h,u,v,uu,vv,h_check real,dimension(43)::dx,f real,dimension(19:54)::h0,