tdt3d = 0.
tdt3d = tdt3d + cond_dt_tg
U tdt3d, olr_2d, tsr_2d,
diff_dt_tg = tdt3d
U udt3d, vdt3d, tdt3d,
O tdt3d, qdt3d,
diff_dt_tg = tdt3d - diff_dt_tg
atmPhys_dT(1:sNx,1:sNy,kc,bi,bj) = tdt3d(:,:,k)
_RL tdt3d (sNx,sNy,Nr)