A degree to radian conversion should be done

[Fengzi21]

In sampler_hi.f90

...
 !generate inclination 
 cos_i=rand()
 inclinations(i)=acos(cos_i)*180./pi


 !generate ellipticity starting from inclination - assumption on the 3D morphology. all HI galaxies have spiral morphology. start from an intrinsic axis ratio of 0.2
 q=0.2+random_normal()*0.05 
 q2=q**2 ! square alpha
 q=sqrt(q2+(cos(inclinations(i)))**2.*(1.-q2)) ! observed axis ratio, linked to inclination
...

inclinations are generated in degrees but later used in cos without being converted to radians when generating the observed axis ratio.