Closed DWesl closed 5 months ago
I saw this after creating a pull request to fix the same problem. See: https://github.com/jswhit/pyspharm/pull/12
It probably doesn't make a big difference which one is merged, but the comment lines suggest that the theta
array should in fact be of length l
rather than 1
, hence the choice in the above PR.
https://github.com/jswhit/pyspharm/blob/ba80d44636bb0d0df1e4f2bdb74465b08c71c724/src/alf.f#L203-L211
Superceeded by #12
My compiler complained about a rank mismatch, so I copied the declarations from the called functions to the calling functions.
Edit: Refocused to be just the Fortran updates, since I can do a lot of packaging things. I'll put those changes in a new PR soon.