!======================================= ! 2D cumulus model - kaminari ! - subroutine vel_u ! ! Author : TAKAHASHI Koko ! Date : 2003/11/21 最終更新 ! 2003/11/04 新規作成 ! Note : 水平速度 u ! !======================================= subroutine vel_u(i,k,im,km,cp,alp,dts,g13z,div_x_u,div_z_u, & & dx,dz,vptemp_bs,pi,u_adv,u_cori,u_trb,u) implicit none integer(8), intent(in) :: im integer(8), intent(in) :: km integer(8), intent(in) :: i, k real(8), intent(in) :: cp, alp, dts real(8), intent(in) :: dx, dz real(8), intent(in) :: vptemp_bs(-2:im+1,-2:km+1) real(8), intent(in) :: pi(-2:im+1,-2:km+1) real(8), intent(in) :: g13z(-2:im+1,-2:km+1) real(8), intent(in) :: div_x_u(-2:im+1,-2:km+1) real(8), intent(in) :: div_z_u(-2:im+1,-2:km+1) real(8), intent(in) :: u_adv(-2:im+1,-2:km+1) real(8), intent(in) :: u_cori(-2:im+1,-2:km+1) real(8), intent(in) :: u_trb(-2:im+1,-2:km+1) real(8), intent(inout) :: u(-2:im+1,-2:km+1) !--- 速度 u u(i,k) = u(i,k) & & + dts*( & & - cp*(vptemp_bs(i+1,k) + vptemp_bs(i,k)) & & /2.0D0 & & *( & & (pi(i,k)- pi(i-1,k))/dx & & - alp*div_x_u(i,k) & & + g13z(i,k)*( & & ( & & pi(i,k+1) + pi(i+1,k+1) & & - pi(i,k-1)- pi(i-1,k-1) & & )/(2.0D0*dz) & & - alp*div_z_u(i,k) & & ) & & ) & & - u_adv(i,k) + u_cori(i,k) + u_trb(i,k) & & ) end subroutine vel_u