O t3d, q3d, u3d, v3d,
I deltaT, t3d, q3d, pFull3d, pHalf3d, coldT,
t3d_tmp = t3d + cond_dt_tg
t3d_tmp = t3d
I sNx,sNy, myTime, lat2d, pHalf3d, t3d, q3d,
I t3d(:,:,Nr), q3d(:,:,Nr), u3d(:,:,Nr), v3d(:,:,Nr),
I sNx,sNy, myTime, lat2d, pHalf3d, t_surf, t3d,
I u3d, v3d, t3d, q3d,
_RL t3d (sNx,sNy,Nr)
t3d(:,:,k) = theta(1:sNx,1:sNy,kc,bi,bj)*conv_theta2T
t3d(:,:,k) = t3d(:,:,k)*pStarFacK(1:sNx,1:sNy,bi,bj)