正壓原始方程模式實習(xí)報告(文件)
【正文】 w(i,j)=a(i,j)+s*(a(i1,j)+a(i+1,j)+a(i,j1)+a(i,j+1)4*a(i,j))/ 執(zhí)行正平滑enddoenddodo i=2,m1 do j=2,n1 a(i,j)=w(i,j)enddoEnddoreturnElsedo i=2,m1 do j=2,n1 w(i,j)=a(i,j)+s*(a(i1,j)+a(i+1,j)+a(i,j1)+a(i,j+1)4*a(i,j))/ 執(zhí)行正平滑 enddoenddodo i=2,m1do j=2,n1a(i,j)=w(i,j)enddoenddodo i=2,m1do j=2,n1 w(i,j)=a(i,j)+(s)*(a(i1,j)+a(i+1,j)+a(i,j1)+a(i,j+1)4*a(i,j))/ 執(zhí)行正逆平滑enddoenddodo i=2,m1 do j=2,n1 a(i,j)=w(i,j) enddoenddoendifreturnend 結(jié)束子程序：integer 定義數(shù)組p=m1i=1,mj=2,qdo 利用公式（）y方向風(fēng)速va(m,j)=rm(m,j)**(za(m,j)za(m1,j))/(f(m,j)*d) 循環(huán)x方向上格點數(shù)va(i,j)=rm(i,j)**(za(i+1,j)za(i1,j))/(*f(i,j)*d)endcgw
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號-1