calculate subscription levels for 5 LECs with only load and no PV (first step)in file subscription_cal.py