In this paper we continue our investigation of the potential theory of Markov processes with jump kernels decaying at the boundary. To be more precise, we consider processes in ${\mathbb R}^d_+$ with jump kernels of the form ${\mathcal B}(x,y) |x-y|^{-d-α}$ and killing potentials $κ(x)=cx_d^{-α}$, $0<α<2$. The boundary part ${\mathcal B}(x,y)$ is comparable to the product of three terms with parameters $β_1, β_2$, $β_3$ and $β_4$ appearing as exponents in these terms. The constant $c$ in the killing term can be written as a function of $α$, ${\mathcal B}$ and a parameter $p\in ((α-1)_+, α+β_1)$, which is strictly increasing in $p,$ decreasing to $0$ as $p\downarrow (α-1)_+$ and increasing to $\infty$ as $p\uparrowα+β_1$. We establish sharp two-sided estimates on the Green functions of these processes for all $p\in ((α-1)_+, α+β_1)$ and all admissible values of $β_1, β_2$, $β_3$ and $β_4$. Depending on the regions where $β_1$, $β_2$ and $p$ belong, the estimates on the Green functions are different. In fact, the estimates have three different forms depending on the regions the parameters belong to. As applications, we prove that the boundary Harnack principle holds in certain region of the parameters and fails in some other region of the parameters. Combined with the main results of \cite{KSV},we completely determine the region of the parameters where the boundary Harnack principle holds.