Let $Z=(Z^{1}, \ldots, Z^{d})$ be the d-dimensional Lévy {process} where {$Z^i$'s} are independent 1-dimensional Lévy {processes} with identical jumping kernel $ ν^1(r) =r^{-1}φ(r)^{-1}$. Here $φ$ is {an} increasing function with weakly scaling condition of order $\underline α, \overline α\in (0, 2)$. We consider a symmetric function $J(x,y)$ comparable to \begin{align*} \begin{cases} ν^1(|x^i - y^i|)\qquad&\text{ if $x^i \ne y^i$ for some $i$ and $x^j = y^j$ for all $j \ne i$}\\ 0\qquad&\text{ if $x^i \ne y^i$ for more than one index $i$}. \end{cases} \end{align*} Corresponding to the jumping kernel $J$, there exists an anisotropic Markov process $X$, see \cite{KW22}. In this article, we establish sharp two-sided Dirichlet heat kernel estimates for $X$ in $C^{1,1}$ open set, under certain regularity conditions. As an application of the main results, we derive the Green function estimates.