t3d(:,:,k) = theta(1:sNx,1:sNy,kc,bi,bj)*conv_theta2T

          t3d(:,:,k) = t3d(:,:,k)*pStarFacK(1:sNx,1:sNy,bi,bj)

     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,

     I           u3d, v3d, t3d, q3d,

      _RL t3d     (sNx,sNy,Nr)