Harmonic mappings onto parallel slit domains

We consider typically real harmonic univalent functions in the unit disk ID whose range is the complex plane slit along infinite intervals on each of the lines x +/- ib, b > 0. They are obtained via the shear construction of conformal mappings of D onto the plane without two or four half-lines symmetric with respect to the real axis.