Bug in gamrnd (prior_draw)

For some types of gamma samples, like

mean=1.5 and std=0.25,

which implies a=36, b= 0.0417

the gamrnd dynare routine provides wrong values (too high), always rejected since the values are always above the upper bound.