wx_pbl0_merge Subroutine

public subroutine wx_pbl0_merge(knon, ypplay, ypaprs, sigw, dTs_forcing, dqs_forcing, yt_x, yt_w, yq_x, yq_w, yu_x, yu_w, yv_x, yv_w, ycdragh_x, ycdragh_w, ycdragq_x, ycdragq_w, ycdragm_x, ycdragm_w, AcoefT_x, AcoefT_w, AcoefQ_x, AcoefQ_w, AcoefU_x, AcoefU_w, AcoefV_x, AcoefV_w, BcoefT_x, BcoefT_w, BcoefQ_x, BcoefQ_w, BcoefU_x, BcoefU_w, BcoefV_x, BcoefV_w, AcoefT, AcoefQ, AcoefU, AcoefV, BcoefT, BcoefQ, BcoefU, BcoefV, ycdragh, ycdragq, ycdragm, yt1, yq1, yu1, yv1)

   BcoefT(j) = (sigw(j)*Kech_h_w(j)*Kech_T_pw(j)*BcoefT_w(j) + &
            sigx(j)*Kech_h_x(j)*Kech_T_px(j)*BcoefT_x(j) )/(Kech_h(j)*Kech_Tp(j))
   BcoefQ(j) = (sigw(j)*Kech_q_w(j)*Kech_Q_pw(j)*BcoefQ_w(j) + &
            sigx(j)*Kech_q_x(j)*Kech_Q_px(j)*BcoefQ_x(j) )/(Kech_q(j)*Kech_Qp(j))
   BcoefU(j) = (sigw(j)*Kech_m_w(j)*Kech_U_pw(j)*BcoefU_w(j) + &
            sigx(j)*Kech_m_x(j)*Kech_U_px(j)*BcoefU_x(j) )/(Kech_m(j)*Kech_Up(j))
   BcoefV(j) = (sigw(j)*Kech_m_w(j)*Kech_V_pw(j)*BcoefV_w(j) + &
            sigx(j)*Kech_m_x(j)*Kech_V_px(j)*BcoefV_x(j) )/(Kech_m(j)*Kech_Vp(j))

Print ,'YYYYpbl0: BcoefT_x, sigw, sigx, dd_Kh, dd_KTp, Kech_h_w ', & BcoefT_x, sigw, sigx, dd_Kh, dd_KTp, Kech_h_w Print ,'YYYYpbl0: Kech_T_pw, dd_BT, Kech_h, Kech_Tp ', & Kech_T_pw, dd_BT, Kech_h, Kech_Tp jyg yt1(j) = yt_x(j,1) + sigw(j)dd_t(1.+sigx(j)dd_Kh/KCT) yq1(j) = yq_x(j,1) + sigw(j)dd_q(1.+sigx(j)dd_Kh/KCQ) yu1(j) = yu_x(j,1) + sigw(j)dd_u(1.+sigx(j)dd_Km/KCU) yv1(j) = yv_x(j,1) + sigw(j)dd_v(1.+sigx(j)dd_Km/KCV)

Arguments