def write_notes(system, Vnominal=240.0):
stationClass = 1
name, usaf = geo.closestUSAF( geo.zipToCoordinates(system.zipcode), stationClass)
mintemp = epw.minimum(usaf)
twopercentTemp = epw.twopercent(usaf)
fields = []
for i in set(system.shape):
print "PV Module Ratings @ STC"
print "Module Make:", i.array.make
fields.append(('Text1ModuleMake',i.array.make))
print "Module Model:", i.array.model
fields.append(('Text1ModuleModel',i.array.model))
print "Max Power-Point Current (Imp):",i.array.Impp
fields.append(('MAX POWERPOINT CURRENT IMP',i.array.Impp))
print "Max Power-Point Voltage (Vmp):",i.array.Vmpp
fields.append(('MAX POWERPOINT VOLTAGE VMP',i.array.Vmpp))
print "Open-Circuit Voltage (Voc):",i.array.Voc
fields.append(('OPENCIRCUIT VOLTAGE VOC',i.array.Voc))
print "Short-Circuit Current (Isc):",i.array.Isc
fields.append(('SHORTCIRCUIT CURRENT ISC',i.array.Isc))
fields.append(('MAX SERIES FUSE OCPD','15'))
print "Maximum Power (Pmax):",i.array.Pmax
fields.append(('MAXIMUM POWER PMAX',i.array.Pmax))
print "Module Rated Max Voltage:",i.array.Vrated
fields.append(('MAX VOLTAGE TYP 600VDC',i.array.Vrated))
fields.append(('VOC TEMP COEFF mVoC or oC',round(i.array.TkVoc,2)))
fields.append(('VOC TEMP COEFF mVoC','On'))
print "Inverter Make:",i.make
fields.append(('INVERTER MAKE',i.make))
print "Inverter Model:",i.model
fields.append(('INVERTER MODEL',i.model))
print "Max Power", i.Paco
fields.append(('MAX POWER 40oC',i.model))
fields.append(('NOMINAL AC VOLTAGE',240))
print "Max AC Current: %s" % round(i.Paco/Vnominal,2)
fields.append(('MAX AC CURRENT', round(i.Paco/Vnominal,2)))
fields.append(('MAX DC VOLT RATING',i.model))
print "Max AC OCPD Rating: %s" % ee.ocpSize(i.Paco/Vnominal*1.25)
print "Max System Voltage:",round(i.array.Vmax(mintemp),1)
print "AC Output Current: %s" % \
round(sum([i.Paco for i in system.shape])/Vnominal,2)
fields.append(('AC OUTPUT CURRENT', \
round(sum([i.Paco for i in system.shape])/Vnominal,2)))
print "Nominal AC Voltage: %s" % Vnominal
fields.append(('NOMINAL AC VOLTAGE_2',i.ac_voltage))
print "Minimum Temperature: %s C" % mintemp
print "2 Percent Max: %s C" % twopercentTemp
from fdfgen import forge_fdf
fdf = forge_fdf("",fields,[],[],[])
fdf_file = open("data.fdf","w")
fdf_file.write(fdf)
fdf_file.close()
import shlex
from subprocess import call
cmd = shlex.split("pdftk Example2-Micro-Inverter.pdf fill_form data.fdf output output.pdf flatten")
rc = call(cmd)
def fill(inverter, zipcode, acDcRatio = 1.2, mount="Roof", stationClass = 1, Vmax = 600, bipolar= True):
"""maximize array"""
import geo
import epw
tDerate = {"Roof":30,
"Ground":25,
"Pole":20}
name, usaf = geo.closestUSAF( geo.zipToCoordinates(zipcode), stationClass)
maxV = inverter.array.panel.Vmax(epw.minimum(usaf))
#NREL suggests that long term degradation is primarily current not voltage
derate20 = .97
minV = inverter.array.panel.Vmin(epw.twopercent(usaf),tDerate[mount]) * derate20
#print "MinV", minV
if inverter.vdcmax != 0:
Vmax = inverter.vdcmax
smax = int(Vmax/maxV)
#range to search
pTol = .30
inverterNominal = inverter.Paco
print inverter.make, inverter.model
print "Nominal AC Power:", inverterNominal
psize = inverter.array.panel.Pmax
print "Nominal Panel Power:", round(psize,1)
solutions = []
Imax = max(inverter.idcmax,inverter.Pdco*1.0/inverter.mppt_low)
stringMax = int(round(Imax/inverter.array.panel.Impp))+1
#Diophantine equation
for s in range(smax+1):
if (s*minV) >= inverter.mppt_low:
for p in range(stringMax):
pRatio = p*s*psize*1.0/inverterNominal
if pRatio < (acDcRatio*(1+pTol)) and \
pRatio > (acDcRatio*(1-pTol)):
sol ="%sW : %sS x %sP : ratio %s : %s - %s V" % (round(s*p*psize,1),s,p, round(pRatio,2), round(s*minV,0), round(s*maxV))
solutions.append(sol)
inverter.array.series = s
inverter.array.parallel = p
if len(solutions) ==0:
solutions.append("Error: Solution not found")
return solutions
def fill(inverter, zipcode, acDcRatio = 1.2, mount="Roof", stationClass = 1, \
Vmax = 600, bipolar= True):
import geo
"""String sizing"""
tDerate = {"Roof":30,
"Ground":25,
"Pole":20}
#csv is performance hit
name, usaf = geo.closestUSAF( geo.zipToCoordinates(zipcode), stationClass)
maxV = inverter.array.panel.Vmax(epw.minimum(usaf))
#NREL suggests that long term degradation is primarily current not voltage
derate20 = .97
minV = inverter.array.panel.Vmin(epw.twopercent(usaf),tDerate[mount]) * \
derate20
if inverter.vdcmax != 0:
Vmax = inverter.vdcmax
smax = int(Vmax/maxV)
#range to search
pTol = .30
inverterNominal = inverter.Paco
psize = inverter.array.panel.Pmax
solutions = []
Imax = max(inverter.idcmax,inverter.Pdco*1.0/inverter.mppt_low)
stringMax = int(round(Imax/inverter.array.panel.Impp))+1
#Diophantine equation
for s in range(smax+1):
if (s*minV) >= inverter.mppt_low:
for p in range(stringMax):
pRatio = p*s*psize*1.0/inverterNominal
if pRatio < (acDcRatio*(1+pTol)) and \
pRatio > (acDcRatio*(1-pTol)):
inverter.array.shape = [s]*p
t = copy.deepcopy(inverter)
t.minV = s*minV
t.maxV = s*maxV
solutions.append(t)
return solutions
#print header[1]
#print self.latitude, self.longitude
def __iter__(self):
return self
def next(self):
t = self.tmy_data.next()
sd = t['Date (MM/DD/YYYY)'] +' '+ t['Time (HH:MM)']
tz = -5
t['utc_datetime'] = strptime(sd,tz)
t['datetime'] = strptime(sd)
return t
def __del__(self):
self.csvfile.close()
if __name__ == "__main__":
import geo
tilt = 32.0
#import matplotlib.pyplot as plt
#place = zipToCoordinates(17601) #Lancaster
place = geo.zipToCoordinates(19113) #Philadelphia
name, usaf = geo.closestUSAF(place)
t = 0
for r in data(usaf):
output = irradiation.irradiation(r,place,tilt)
t += output
print t/1000
print t/(1000*365.0)
def setZipcode(self,zipcode):
self.zipcode = zipcode
self.place= geo.zipToCoordinates(self.zipcode)
self.tz = geo.zipToTZ(self.zipcode)
self.name, self.usaf = geo.closestUSAF(self.place)
def string_notes(system, run=0.0, stationClass = 1):
"""page 5"""
name, usaf = geo.closestUSAF( geo.zipToCoordinates(system.zipcode), stationClass)
mintemp = epw.minimum(usaf)
twopercentTemp = epw.twopercent(usaf)
ac_rated = 0.0
dc_rated = 0.0
for i in system.shape:
dc_rated += i.array.Pmax
try:
if i.phase == 1:
ac_rated += i.current * i.ac_voltage
else:
ac_rated += i.phase * i.current * i.ac_voltage/ 3**.5
except:
ac_rated += i.Paco
pass
notes = []
notes.append("%s KW AC RATED" % round(ac_rated/1000.0,2))
notes.append("%s KW DC RATED" % round(dc_rated/1000.0,2))
#BUG: This doesn't work for unbalanced 3 phase
if system.phase == 1:
Aac = round(ac_rated/i.ac_voltage,1)
else:
Aac = round(ac_rated/i.ac_voltage/3**.5,1)
notes.append( "System AC Output Current: %s A" % Aac)
notes.append("Nominal AC Voltage: %s V" % i.ac_voltage)
notes.append("")
notes.append("Minimum Temperature: %s C" % mintemp)
notes.append("2 Percent Max Temperature: %s C" % twopercentTemp)
notes.append("Weather Source: %s %s" % (name, usaf))
notes.append("")
di, dp = system.describe()
aMax = 0
for i in system.shape:
if dp.has_key(i.array.panel.model):
notes.append( "PV Module Ratings @ STC")
notes.append("Module Make: %s" % i.array.panel.make)
notes.append("Module Model: %s" % i.array.panel.model)
notes.append("Quantity: %s" % dp[i.array.panel.model])
notes.append("Max Power-Point Current (Imp): %s A" % i.array.panel.Impp)
notes.append("Max Power-Point Voltage (Vmp): %s V" % i.array.panel.Vmpp)
notes.append("Open-Circuit Voltage (Voc): %s V" % i.array.panel.Voc)
notes.append("Short-Circuit Current (Isc): %s A" % i.array.panel.Isc)
notes.append("Maximum Power (Pmax): %s W" % round(i.array.panel.Pmax,1))
#notes.append("Module Rated Max Voltage: %s V" % i.array.panel.Vrated)
notes.append("")
dp.pop(i.array.panel.model)
if di.has_key(i.model):
notes.append("Inverter Make: %s" % i.make)
notes.append("Inverter Model: %s" % i.model)
notes.append("Quantity: %s" % di[i.model])
notes.append("Max Power: %s KW" % round(i.Paco/1000.0,1))
#this is hack... This should be calculated based upon power cores
if hasattr(i,'current'):
notes.append("Max AC Current: %s A" % round(i.current,1))
elif i.ac_voltage == 480:
notes.append("Max AC Current: %s A" % round(i.Paco*1.0/i.ac_voltage/3**.5,1))
else:
notes.append("Max AC Current: %s A" % round(i.Paco*1.0/i.ac_voltage,1))
#greater than 1 in parallel
if len(i.array.shape) > 1:
pass
notes.append("DC Operating Current: %s A" % \
round(i.array.panel.Impp*len(i.array.shape),1))
notes.append("DC Short Circuit Current: %s A" % \
round(i.array.panel.Isc*len(i.array.shape),1))
#greater than 1 in series
if max(i.array.shape)> 1:
notes.append("DC Operating Voltage: %s V" % round(i.array.Vdc(),1))
notes.append("System Max DC Voltage: %s V" % round(i.array.Vmax(mintemp),1))
if i.array.Vmax(mintemp) > 600:
print "WARNING: Array exceeds 600V DC"
notes.append("Pnom Ratio: %s" % round((i.array.Pmax/i.Paco),2))
if (i.array.Vdc(twopercentTemp) *.9) < i.mppt_low:
print "WARNING: Array IV Knee drops out of Inverter range"
if (i.array.Pmax/i.Paco) < 1.1:
print "WARNING: Array potentially undersized"
notes.append("")
di.pop(i.model)
if i.array.Vmax(mintemp) > aMax:
aMax = i.array.Vmax(mintemp)
notes.append("Array Azimuth: %s Degrees" % system.azimuth)
notes.append("Array Tilt: %s Degrees" % system.tilt)
s9 = system.solstice(9)
s15 = system.solstice(15)
notes.append("December 21 9:00 AM Sun Azimuth: %s Degrees" % \
(round(degrees(s9[1]),1)))
notes.append("December 21 9:00 AM Sun Altitude: %s Degrees" % \
(round(degrees(s9[0]),1)))
notes.append("December 21 3:00 PM Sun Azimuth: %s Degrees" % \
(round(degrees(s15[1]),1)))
notes.append("December 21 3:00 PM Sun Altitude: %s Degrees" % \
(round(degrees(s9[0]),1)))
if 'geomag' in sys.modules:
notes.append("Magnetic declination: %s Degrees" % \
round(geomag.declination(dlat=system.place[0],dlon=system.place[1])))
notes.append("Minimum Row space ratio: %s" % \
round(system.minRowSpace(1.0),2))
print "\n".join(notes)
print ""
print "Minimum Bundle"
minC = vd.vd(Aac,5,verbose=False)
try:
ee.assemble(minC,Aac,conduit='STEEL')
if run > 0:
print "Long Run"
minC = vd.vd(Aac,run,v=i.ac_voltage,tAmb=15,pf=.95,material='AL',verbose=False)
ee.assemble(minC,Aac,conduit='PVC')
except:
print "Warning: Multiple sets of conductors"
return notes
print perez(458,84,457,0.558505360638,0.905073340885,1.30063355891)
print perez(365,52,283,0.558505360638,1.11854235428,1.44247248195)
print perez(83,22,75,0.558505360638,1.33812695863,1.59866999006)
print perez(0,0,0,0.558505360638,1.58526434614,1.7938116797)
print perez(0,0,0,0.558505360638,1.82430194526,1.98580885248)
print perez(0,0,0,0.558505360638,2.0648085416,2.18380405772)
print perez(0,0,0,0.558505360638,2.30391057936,2.38276512577)
print perez(0,0,0,0.558505360638,2.5366580297,2.57519603892)
print perez(0,0,0,0.558505360638,2.74856500007,2.7436087649)
print perez(0,0,0,0.558505360638,2.87979696349,2.83820540821)
#from scipy.interpolate import interp1d
#import numpy as np
import geo
#thorizon = interp1d(np.array([-180.0,180.0]),np.array([0.0,0.0]))
timestamp = datetime.datetime.now()
place = geo.zipToCoordinates('17603')
tilt = 1
azimuth = 180
print blave(timestamp,place,tilt,azimuth)
print bblave(timestamp,place,tilt,azimuth)
print "Bird",bird(0.478253620172,1000.0)
azimuth = 90
print blave(timestamp,place,tilt,azimuth)
azimuth = 270
print blave(timestamp,place,tilt,azimuth)
irradiance = 300.
utc_datetime = datetime.datetime.utcnow()
print utc_datetime
solarAz, solarAlt = ephemSun(place,utc_datetime)
surfaceTilt = radians(15)
surfaceAz = radians(180)
import geo
parser = argparse.ArgumentParser(description='Model a PV system. Currently displays annual output and graph')
#import sys
#opts, args = getopt.getopt(sys.argv[1:], 'f:h')
parser.add_argument('-z', '--zipcode',required=True)
parser.add_argument('-m', '--mname')
parser.add_argument('-v', '--voltage',type=int,default=600)
args = vars(parser.parse_args())
#print args
try:
#start program
zip = args['zipcode']
maxVoltage = args['voltage']
stationClass = 1
name, usaf = geo.closestUSAF( geo.zipToCoordinates(zip), stationClass)
print "%s USAF: %s" % (name, usaf)
print "Minimum Temperature: %s C" % minimum(usaf)
print "2%% Max: %s C" % twopercent(usaf)
print "Heating Degree days: %s" % hdd(usaf)
print "Cooling Degree days: %s" % cdd(usaf)
if args['mname']:
print ""
import modules
models = modules.model_search(args['mname'].split(' '))
m = None
if len(models) > 1:
for i in models:
print i
sys.exit(1)
elif len(models) == 1:
def string_notes(system, run=0.0):
"""page 5"""
stationClass = 3
name, usaf = geo.closestUSAF( geo.zipToCoordinates(system.zipcode), stationClass)
mintemp = epw.minimum(usaf)
twopercentTemp = epw.twopercent(usaf)
ac_rated = 0
dc_rated = 0
for i in system.shape:
dc_rated += i.array.Pmax
ac_rated += i.Paco
notes = []
notes.append("%s KW AC RATED" % round(ac_rated/1000.0,2))
notes.append("%s KW DC RATED" % round(dc_rated/1000.0,2))
#BUG: This doesn't work for 3 phase
if system.phase == 1:
Aac = round(sum([i.Paco for i in system.shape])/i.ac_voltage,1)
else:
Aac = round(sum([i.Paco for i in system.shape])/i.ac_voltage/3**.5,1)
notes.append( "System AC Output Current: %s A" % Aac)
notes.append("Nominal AC Voltage: %s V" % i.ac_voltage)
notes.append("")
notes.append("Minimum Temperature: %s C" % mintemp)
notes.append("2 Percent Max Temperature: %s C" % twopercentTemp)
notes.append("")
di, dp = system.describe()
aMax = 0
for i in system.shape:
if dp.has_key(i.array.panel.model):
notes.append( "PV Module Ratings @ STC")
notes.append("Module Make: %s" % i.array.panel.make)
notes.append("Module Model: %s" % i.array.panel.model)
notes.append("Quantity: %s" % dp[i.array.panel.model])
notes.append("Max Power-Point Current (Imp): %s A" % i.array.panel.Impp)
notes.append("Max Power-Point Voltage (Vmp): %s V" % i.array.panel.Vmpp)
notes.append("Open-Circuit Voltage (Voc): %s V" % i.array.panel.Voc)
notes.append("Short-Circuit Current (Isc): %s A" % i.array.panel.Isc)
notes.append("Maximum Power (Pmax): %s W" % round(i.array.panel.Pmax,1))
#notes.append("Module Rated Max Voltage: %s V" % i.array.panel.Vrated)
notes.append("")
dp.pop(i.array.panel.model)
if di.has_key(i.model):
notes.append("Inverter Make: %s" % i.make)
notes.append("Inverter Model: %s" % i.model)
notes.append("Quantity: %s" % di[i.model])
notes.append("Max Power: %s KW" % round(i.Paco/1000.0,1))
#this is hack... This should be calculated based upon power cores
if i.ac_voltage == 480:
notes.append("Max AC Current: %s A" % round(i.Paco*1.0/i.ac_voltage/3**.5,1))
else:
notes.append("Max AC Current: %s A" % round(i.Paco*1.0/i.ac_voltage,1))
#greater than 1 in parallel
if len(i.array.shape) > 1:
pass
notes.append("DC Operating Current: %s A" % \
round(i.array.panel.Impp*len(i.array.shape),1))
notes.append("DC Short Circuit Current: %s A" % \
round(i.array.panel.Isc*len(i.array.shape),1))
#greater than 1 in series
if max(i.array.shape)> 1:
notes.append("DC Operating Voltage: %s V" % round(i.array.Vdc(),1))
notes.append("System Max DC Voltage: %s V" % round(i.array.Vmax(mintemp),1))
notes.append("Pnom Ratio: %s" % round((i.array.Pmax/i.Paco),2))
notes.append("")
di.pop(i.model)
if i.array.Vmax(mintemp) > aMax:
aMax = i.array.Vmax(mintemp)
notes.append("Array Azimuth: %s Degrees" % system.azimuth)
notes.append("Array Tilt: %s Degrees" % system.tilt)
notes.append("December 21 9:00 AM Sun Azimuth: %s Degrees" % \
int(round(degrees(system.solstice(9)[1]),0)))
notes.append("December 21 3:00 PM Sun Azimuth: %s Degrees" % \
int(round(degrees(system.solstice(15)[1]),0)))
if sys.modules['geomag']:
notes.append("Magnetic declination: %s Degrees" % \
round(geomag.declination(dlat=system.place[0],dlon=system.place[1])))
notes.append("Minimum Row space ratio: %s" % \
round(system.minRowSpace(1.0),2))
print "\n".join(notes)
print ""
print "Minimum Bundle"
minC = vd.vd(Aac,5)
ee.assemble(minC,Aac,conduit='STEEL')
if run > 0:
print "Long Run"
minC = vd.vd(Aac,run,v=i.ac_voltage,tAmb=15,pf=.95,material='AL')
ee.assemble(minC,Aac,conduit='PVC')
return notes