`xa` and `xb` can be included in the struct `GRT_MODEL1D`

[Dengda98]

For fixed $\omega$ and $k$，parameters like xa, xb, caca and cbcb for each layer are determined.

[Dengda98]

fixed in #127