This module contains functions for initializing Cumulant relaxation paramaters author: Gregorio Gerardo Spinelli

Type | Visibility | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|---|

integer, | private, | parameter | :: | QQ | = | 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 |