!======================================= ! 2D cumulus model - kaminari ! - subroutine init_value ! ! Author : TAKAHASHI Koko ! Date : 2003/11/11 最終更新 ! 2003/11/11 新規作成 ! Note : 初期値設定 ! !======================================= subroutine init_value(i,k,im,km,u,w,prss_bs,temp_bs,prss, & & temp,qv,qc,qr,km_sub) implicit none integer(8), intent(in) :: i, k integer(8), intent(in) :: im, km real(8), intent(in) :: prss_bs(-2:im+1,-2:km+1) real(8), intent(in) :: temp_bs(-2:im+1,-2:km+1) ! real(8), intent(in) :: x(-2:im+1,-2:km+1) ! real(8), intent(in) :: z(-2:im+1,-2:km+1) real(8), intent(out) :: u(-2:im+1,-2:km+1) real(8), intent(out) :: w(-2:im+1,-2:km+1) real(8), intent(out) :: prss(-2:im+1,-2:km+1) real(8), intent(out) :: temp(-2:im+1,-2:km+1) real(8), intent(out) :: qv(-2:im+1,-2:km+1) real(8), intent(out) :: qc(-2:im+1,-2:km+1) real(8), intent(out) :: qr(-2:im+1,-2:km+1) real(8), intent(out) :: km_sub(-2:im+1,-2:km+1) !--- 水平速度 u u(i,k) = 0.0d0 !--- 鉛直速度 w w(i,k) = 0.0d0 !--- 圧力 prss(i,k) = prss_bs(i,k) !--- 温度 temp(i,k) = temp_bs(i,k) ! temp(i,k) = temp_bs(i,k) & ! & + 1.0d2*exp( - (x(i,k) - 1.0d4/2.0d0)**2.0d0/(5.0d2**2.0d0)& ! & - z(i,k)**2.0d0/(2.0d2**2.0d0)) !--- 水蒸気の混合比 qv(i,k) = 0.00d0 !--- 雲水の混合比 qc(i,k) = 0.00d0 !--- 雨水の混合比 qr(i,k) = 0.00d0 !--- 渦混合係数 km_sub(i,k) = 10.0d0 end subroutine init_value