Functions to compare ISDs

- Modify function weights (renamed as AICc_ISD) in module fit_other_dist
to compute the AICc values for the four forms of ISD given data from one
site and MLE parameters.
- Add function AICc_ISD_to_file in module working_functions, which calls
AICc_ISD and save the AICcs for multiple sites to an output file.