Checking the stability regions of omegas for the parametrized Cumulant Just omega(2) is given in input. omega(2)=-1 means omega2=omegaBulk. Limiters read from input. lim(N)=10^10 means unlimited. lim(N) is for omega(N+2). Just omega(3:5) are limited as in the paper. omega(6:10) = 1._rk
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=rk), | intent(in) | :: | omegaVisc(:) |
vector of omegas in the level |
||
| real(kind=rk), | intent(in) | :: | omegaBulk |
omega bulk value in the level |
||
| real(kind=rk), | intent(in) | :: | omegaIn(:) |
vector of omegas as given in musubi.lua |
||
| integer, | intent(in) | :: | nSolve |
number of elements solved in kernel |
||
| integer, | intent(in) | :: | level |
current level |