The basic model in this project (scm3Elm.py) solves the following equation for En:
∂t En = (g2m - g2l)tEn+l + (1 - g2m)tEn + (g2l - 1)tEn-m+l
where:
tEn = g-lknsin[(m-l)α]En3/2
with the two dimensional wavenumber defined as a complex number using kn = k0 (g eiα)n
in general different values of l and m are possible as described in [https://arxiv.org/abs/1903.09494]. And the case that combines l=2,m=3 and l=1, m=5 and considering four such spiral chains (which has the best angular coverage) is included as fcm6E.py
Finally the model for the case l=1,m=3 for the complex amplitude Φn is also given as scm3_13.py