* Added plotting of total gas maps compared to AV maps

* Added calculation of angle theta on unprojected disk, in prep for calculating model of f_red. Implemented as switch on get_major_axis()

* added test code for f_red distribution in test_maj_axis()