  We present a numerical code that solves the forward and inverse problem of
the polarized radiative transfer equation in geometrical scale under the Zeeman
regime. The code is fully parallelized, making it able to easily handle large
observational and simulated datasets. We checked the reliability of the forward
and inverse modules through different examples. In particular, we show that
even when properly inferring various physical parameters (temperature, magnetic
field components, and line-of-sight velocity) in optical depth, their
reliability in height-scale depends on the accuracy with which the gas-pressure
or density are known. The code is made publicly available as a tool to solve
the radiative transfer equation and perform the inverse solution treating each
pixel independently. An important feature of this code, that will be exploited
in the future, is that working in geometrical-scale allows for the direct
calculation of spatial derivatives, which are usually required in order to
estimate the gas pressure and/or density via the momentum equation in a
three-dimensional volume, in particular the three-dimensional Lorenz force.
