qdt3d = 0.
qdt3d = qdt3d + cond_dt_qg
diff_dt_qg = qdt3d
I qdt3d,
O tdt3d, qdt3d,
diff_dt_qg = qdt3d - diff_dt_qg
atmPhys_dQ(1:sNx,1:sNy,kc,bi,bj) = qdt3d(:,:,k)
_RL qdt3d (sNx,sNy,Nr)