Small advances in calibration of Halpha images
Simplification of program, removal of nlow, nhigh, output is now a single parameter...
Add WHERE sources = 0 to avoid the sources image when collecting the airmasses
Read_image_with_mask can now open parts of images, not just whole images. if_dir_not_exists_create is now recursive, so that you can have several levels missing and they will be created