# -*- coding: utf-8 -*-

"""

Created on Wed Dec 27 05:36:53 2017

@author: shubaraghavan

"""

from scipy.stats import weibull_min

# from HDD_CDD_UValues import *

from Inputs_Energy import *

from Inputs_ElecRate import *

from RefrigerantCalc import Refrigerant

from AnnualHeatDemand import *

import numpy as np

from numpy import *

HPCond = {} # indexed by year, stock number

HPCondHeatEng = {} # Total Heat Energy usage by stock in a year

HPCondCoolEng = {} # Total cooling energy usage by stock in a year

class FuelType:

self.UnitEmissions = UnitEmissions

# Time = 2020

UltimYr0 = 15 # MeanLife #vintage <2000..normal life till ThisYear+3..and then killed

UltimYr1 = 15 # MeanLife #vintage between 2000 and THisYear =2020..early demise for these ineff adoptions

UltimYr2 = 15 # MeanLife # after 2020 adoptions are assumed to be efficient ..so will live normal life

# print "UNITNG",UnitNG

NG = FuelType("NG", UnitNG, NGCostYrly, NGEmisYrly)

Elec = FuelType("Elec", UnitElec, ElecCostYrly, ElecEmisYrly)

Ref1 = Refrigerant(2088, 1, 0.005, 0.05, 0.2) # R-410A

Ref2 = Refrigerant(675, 1, 0.005, 0.05,

0.2) # HFC-32 https://www.epa.gov/sites/production/files/2016-12/documents/international_transitioning_to_low-gwp_alternatives_in_res_and_com_ac_chillers.pdf

Ref3 = Refrigerant(3, 1, 0.005, 0.05, 0.2) # 1

Ref11 = Refrigerant(1725, 1, 0.005, 0.1, 0.3) # ..

Ref12 = Refrigerant(650, 1, 0.005, 0.1, 0.3) # ...

Ref13 = Refrigerant(1, 1, 0.005, 0.1, 0.3)

class Device:

return N

class NGWH1(Device, object):

return super(NGWH1, self).AnnualWaterEngUsage(self.dailyVol, self.IncTemp)

class NGWH2(Device, object):

return super(NGWH2, self).AnnualWaterEngUsage(self.dailyVol, self.IncTemp)

class ERWH1(Device, object):

return super(ERWH1, self).AnnualWaterEngUsage(self.dailyVol, self.IncTemp)

class HPWH1(Device, object):

return super(HPWH1, self).AnnualWaterEngUsage(self.dailyVol, self.IncTemp)

class HPWH2(Device, object):

return super(HPWH2, self).AnnualWaterEngUsage(self.dailyVol, self.IncTemp)

class HPWH3(Device, object):

return super(HPWH3, self).AnnualWaterEngUsage(self.dailyVol, self.IncTemp)

class NGHeater0(Device, object):

"only works for heater"

class NGHeater1(Device, object):

"only works for heater"

class NGHeater2(Device, object):

"only works for heater"

class NGHeater3(Device, object):

"only works for heater"

class ERHeater0(Device, object):

"only works for heater"

class ERHeater1(Device, object):

"only works for heater"

class HP1(Device, object):

"only works for heater"

class HP2(Device, object):

"only works for heater"

class HP3(Device, object):

"only works for heater"

class Cooler0(Device, object):

return super(Cooler0, self).AnnualCoolEngUsage_BTU()

class Cooler1(Device, object):

return super(Cooler1, self).AnnualCoolEngUsage_BTU()

class Cooler2(Device, object):

return super(Cooler2, self).AnnualCoolEngUsage_BTU()

class Cooler3(Device, object):

return super(Cooler3, self).AnnualCoolEngUsage_BTU()

# ========================

# =====================

class Cond1(Device, object): # assuming this is going to be all "ELectric" with no NG in it

"only works for heater"

class Cond2(Device, object): # assuming this is going to be all "ELectric" with no NG in it

"only works for heater"

class Cond3(Device, object): # assuming this is going to be all "ELectric" with no NG in it

"only works for heater"

def NGHeater(C, s1, s2, yr, cost, is_retrofit=False, is_new=False):

return NGHeater3(C, s1, s2, r0, r1, r2, yr, totalcost, is_retrofit, is_new)

def ERHeater(C, s1, s2, yr, cost, is_new=False):

return ERHeater1(C, s1, s2, r0, r1, r2, yr, totalcost, is_new)

def HP(C, s1, s2, yr):

return HP3(C, s1, s2, yr)

def Cooler(C, s1, s2, yr, cost, is_new=False): # is_new == True means new home, false means retrofit

return Cooler3(C, s1, s2, r0, r1, r2, yr, totalcost, is_new)

def Cond(C, s1, s2, yr, cost, is_retrofit=False, is_fuel_switch=False, is_new=False):

return Cond3(C, s1, s2, r0, r1, r2, yr, totalCost, is_retrofit, is_fuel_switch)

def HPWH(C, yr):

return HPWH3(C, yr)

def NGWH(C, yr):

return NGWH3(C, yr)

def ERWH(C, yr):

return ERWH1(C, yr)

# #==============================================================================

# import matplotlib.pyplot as plt

# from matplotlib.pyplot import *

# import matplotlib.patches as mpatches

# fig1 = plt.figure(figsize=(10.0, 8.0))

# a1 = fig1.add_subplot(1,1, 1)

# fig2 = plt.figure(figsize=(10.0, 8.0))

# a2 = fig2.add_subplot(1,1, 1)

# colors = [ 'orange', 'purple', 'g','blue', 'violet','k']

# cnt = 0

# this = 1980

# for yr0 in range(this,2050,10):

# s=0

# k= 3

# # Heat = NGHeater(k,size1,size2 , yr0) # a heater created every 10 years

# Heat1 = Cooler(k,size1, size2, yr0)

# # Heat2 = NGHeater(k,size1,size2, yr0)

# dev1 = [Heat1]*1

# for years in range(yr0,EndYear+5):

# num = 0

# for d in dev1:

# num += 1

# if d.vintage > years: # heater not created

# plot.hold(True)

# plot.hold(True)

# p = 0.0

# dead1 = 1-d.deadsofar(years) #actually live

# a1.scatter(years , p, s = 5, color = 'white')

# a2.scatter(years , dead, s = 5, color = 'red')

# else:

# p = d.death_prob(years) # prob of death ..once heater is created..

# s = d.dead_alive(years)

# dead = 1-d.deadsofar(years)

# a1.scatter(years , p, s = 5, color = 'red')

# a2.scatter(years , dead, s = 5, color ='green')

# print "p", years, d.vintage, d.death_prob(years), p, s

# cnt +=1

# fig1.tight_layout()

# plt.show()

# fig2.tight_layout()

# plt.show()

###yr = ThisYear

# k = 3

##Stck = 100

## #

# devices = []

# devices += [ NGHeater(k,size1,size2, yr)]

##devices += [HPWH1(50,70,k)]

##devices += [HP1(k,size1,size2,yr)]

##

##waterEng = devices[1].AnnualHeatEngUsage_BTU()

##WHCost = devices[1].AnnualEngCost(yr,waterEng)

##WHEmis = devices[1].AnnEmissions(yr,waterEng)

##WHNPVEmis = devices[1].NPVEmissions(yr,waterEng)

###WHAnnEmis = devices[1].annualizedEmissions(yr,waterEng)

##

# SHEng = devices[0].AnnualHeatEngUsage_BTU()

# print "SH Stats", k, ThisYear, devices[0].lt,SHEng/Therm_BTU #, SHCost, SHEmis, SHNPVEmis

##SHCost =devices[0].AnnualEngCost(yr,SHEng)

# SHEmis = devices[0].AnnEmissions(yr,SHEng)

# SHNPVEmis = devices[0].NPVEmissions(yr,SHEng)

#

# SHEng1 = devices[2].AnnualHeatEngUsage_BTU()

# SHCost1 =devices[2].AnnualEngCost(yr,SHEng1)

# SHEmis1 = devices[2].AnnEmissions(yr,SHEng1)

# SHNPVEmis1 = devices[2].NPVEmissions(yr,SHEng1)

##SHAnnEmis = devices[0].annualizedEmissions(yr,SHEng)