wx_pbl_dts_merge Subroutine

public subroutine wx_pbl_dts_merge(knon, dtime, ypplay, ypaprs, sigw, beta, wcstar, wdens, AT_x, AT_w, BT_x, BT_w, AcoefT0, AcoefQ0, BcoefT0, BcoefQ0, AcoefT, AcoefQ, BcoefT, BcoefQ, HTphiT_b, dd_HTphiT, HTphiQ_b, dd_HTphiQ, HTRn_b, dd_HTRn, phiT0_b, dphiT0, phiQ0_b, dphiQ0, Rn0_b, dRn0, g_T, g_Q, Gamma_phiT, Gamma_phiQ, dTs_ins, dqsatsrf_ins)

Compute Cp, Lv, qsat, dqsat_dT.

!! DO j = 1,knon !! IF (mod(j,20) .eq.0) THEN !! print , ' j dd_QQ QQ_b dd_HQphiQ dd_HQphiT dd_HQRn HQphiQ_b HQphiT_b HQRn_b ' !! ENDIF !! print 1789, j, dd_QQ(j), QQ_b(j), dd_HQphiQ(j), dd_HQphiT(j), dd_HQRn(j), HQphiQ_b(j), HQphiT_b(j), HQRn_b(j) !! 1789 FORMAT( I4, 10(1X,E10.2)) !! ENDDO g_Q(1:knon) = - (dd_QQ(1:knon)/QQ_b(1:knon)) * & (sigx(1:knon)-sigw(1:knon)-sigw(1:knon)sigx(1:knon)dd_HQphiQ(1:knon)/HQphiQ_b(1:knon)) & - sqrt(tau(1:knon))/(Inert(1:knon)QQ_b(1:knon)) * & ( HQphiT_b(1:knon) + L_v(1:knon)HQphiQ_b(1:knon) + HQRn_b(1:knon) + & (sigx(1:knon) - sigw(1:knon) - sigw(1:knon)sigx(1:knon)dd_HQphiQ(1:knon)/HQphiQ_b(1:knon)) * & (dd_HQphiT(1:knon) + L_v(1:knon)dd_HQphiQ(1:knon) + dd_HQRn(1:knon)) )

g_Q(1:knon) = - (sqrt(tau(1:knon))/(Inert(1:knon)QQ_b(1:knon))) * & ( HQphiT_b(1:knon) + L_v(1:knon)HQphiQ_b(1:knon) + HQRn_b(1:knon) ) & - (sigx(1:knon) - sigw(1:knon) - sigw(1:knon)sigx(1:knon)dd_HQphiQ(1:knon)/HQphiQ_b(1:knon)) * & ( dd_QQ(1:knon)/QQ_b(1:knon) & + (sqrt(tau(1:knon))/(Inert(1:knon)QQ_b(1:knon))) & * (dd_HQphiT(1:knon) + L_v(1:knon)dd_HQphiQ(1:knon) + dd_HQRn(1:knon)) )

Gamma_phiQ(1:knon) = (beta(1:knon)/(HQphiQ_b(1:knon)QQ_b(1:knon))) * & ( dd_QQ(1:knon) & + (sqrt(tau(1:knon))/(Inert(1:knon))) * & (dd_HQphiT(1:knon) + L_v(1:knon)dd_HQphiQ(1:knon) + dd_HQRn(1:knon)) )

Gamma_phiQ(1:knon) = (1/(HQphiQ_b(1:knon)QQ_b(1:knon))) & * ( dd_QQ(1:knon) & + (sqrt(tau(1:knon))/(Inert(1:knon))) & * (dd_HQphiT(1:knon) + L_v(1:knon)dd_HQphiQ(1:knon) + dd_HQRn(1:knon)) )

 AcoefQ(j) = AcoefQ0(j) - sigw(j)*sigx(j)*(dd_KQs(j)/Kech_Qs(j))*    &
             (dqsatsrf0(j) + (dqsatsrf_ins(j)-(beta(j)/QQ_b(j))*dTs0(j)-Gamma_phiQ(j)*phiQ0_b(j))/(1 - g_Q(j)))/ &
             beta(j)
 BcoefQ(j) = BcoefQ0(j) - sigw(j)*sigx(j)*(dd_KQs(j)/Kech_Qs(j))*Gamma_phiQ(j)/(1 - g_Q(j))/(beta(j)*dtime)

Arguments