This module contains functions for initializing Cumulant relaxation paramaters author: Gregorio Gerardo Spinelli
Type | Visibility | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|---|
integer, | private, | parameter | :: | = | 27 | ||
integer, | private, | parameter | :: | q000 | = | 27 |
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 |