module parameters
integer, parameter :: rho=2
end module parameters

