bifacial-radiance 0.5.1__py2.py3-none-any.whl

bifacial_radiance/data/default.ini

@@ -0,0 +1,65 @@
+#Version 0.3.0
+[simulationParamsDict]
+testfolder: C:\Users\cdeline\Documents\Python Scripts\Test2
+weatherFile: C:\Users\cdeline\Documents\Python Scripts\Test2\EPWs\USA_VA_Richmond.Intl.AP.724010_TMY.epw
+getEPW: True
+simulationname: Demo1
+moduletype: Longi
+rewriteModule: True
+cellLevelModule: False
+axisofrotationTorqueTube: False
+torqueTube: True
+tracking: True
+cumulativeSky: True
+selectTimes: True
+latitude: 37.5
+longitude: -77.6st
+
+[timeControlParamsDict]
+starttime: 06_21_12_00
+endtime: 06_21_12_00
+
+[moduleParamsDict]
+numpanels: 2
+x: 0.98
+y: 1.980
+bifi: 0.90
+xgap: 0.020
+ygap: 0.150
+zgap: 0.100
+
+[sceneParamsDict]
+gcrorpitch: pitch
+gcr: 0.350
+pitch: 10.0
+albedo: 0.30
+nMods: 20
+nRows: 7
+azimuth: 180
+tilt: 30
+clearance_height: 0.7
+hub_height: 2.35
+axis_azimuth: 180
+
+[trackingParamsDict]
+backtrack: True
+limit_angle: 60
+angle_delta: 30
+
+[torquetubeParamsDict]
+diameter: 0.10
+tubetype: Hex
+torqueTubeMaterial: Metal_Grey
+
+[analysisParamsDict]
+sensorsy: 9
+modWanted: 10
+rowWanted: 3
+
+[cellLevelModuleParamsDict]
+numcellsx: 12
+numcellsy: 6
+xcell: 0.150
+ycell: 0.150
+xcellgap: 0.100
+ycellgap: 0.100

bifacial_radiance/data/gencumsky/License.txt

@@ -0,0 +1,54 @@
+The Radiance Software License, Version 1.0
+
+Copyright (c) 1990 - 2018 The Regents of the University of California,
+through Lawrence Berkeley National Laboratory.   All rights reserved.
+
+Redistribution and use in source and binary forms, with or without
+modification, are permitted provided that the following conditions
+are met:
+
+1. Redistributions of source code must retain the above copyright
+        notice, this list of conditions and the following disclaimer.
+
+2. Redistributions in binary form must reproduce the above copyright
+      notice, this list of conditions and the following disclaimer in
+      the documentation and/or other materials provided with the
+      distribution.
+
+3. The end-user documentation included with the redistribution,
+          if any, must include the following acknowledgment:
+            "This product includes Radiance software
+                (http://radsite.lbl.gov/)
+                developed by the Lawrence Berkeley National Laboratory
+              (http://www.lbl.gov/)."
+      Alternately, this acknowledgment may appear in the software itself,
+      if and wherever such third-party acknowledgments normally appear.
+
+4. The names "Radiance," "Lawrence Berkeley National Laboratory"
+      and "The Regents of the University of California" must
+      not be used to endorse or promote products derived from this
+      software without prior written permission. For written
+      permission, please contact radiance@radsite.lbl.gov.
+
+5. Products derived from this software may not be called "Radiance",
+      nor may "Radiance" appear in their name, without prior written
+      permission of Lawrence Berkeley National Laboratory.
+
+THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESSED OR IMPLIED
+WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
+OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
+DISCLAIMED.   IN NO EVENT SHALL Lawrence Berkeley National Laboratory OR
+ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF
+USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
+ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
+OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT
+OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
+SUCH DAMAGE.
+====================================================================
+
+This software consists of voluntary contributions made by many
+individuals on behalf of Lawrence Berkeley National Laboratory.   For more
+information on Lawrence Berkeley National Laboratory, please see
+<http://www.lbl.gov/>.

bifacial_radiance/data/gencumsky/Makefile

@@ -0,0 +1,17 @@
+CC = g++
+CFLAGS = -D_XOPEN_SOURCE
+INCLUDES = -lm
+TARGET = gencumulativesky
+SOURCES = *.cpp 
+
+all: $(TARGET)
+
+$(TARGET): $(SOURCES)
+	$(CC) $(CFLAGS) $(SOURCES) $(INCLUDES)  -o $(TARGET)
+
+.PHONY : clean
+
+clean:
+	rm -f $(TARGET) *.o
+
+

bifacial_radiance/data/gencumsky/README.txt

@@ -0,0 +1,9 @@
+SOLAR RADIATION MAPPING
+========================
+This is directory contains the source for the gencumulativesky program.  To compile it using gcc:
+
+g++ -D_XOPEN_SOURCE *.cpp -lm -o gencumulativesky
+
+(The _XOPEN_SOURCE definition is required so that M_PI is defined when math.h is included - this will need changing depending on your compiler)
+
+bifacial_radiance update Oct 2019:  run python make_gencumskyexe.py from the command line to compile a platform-specific instance of gencumulativesky.exe

bifacial_radiance/data/gencumsky/Sun.cpp

@@ -0,0 +1,118 @@
+#include "Sun.h"
+
+#define _USE_MATH_DEFINES
+#include <math.h>
+#include <stdio.h>
+
+cSun::cSun (double latitude, int day, double hourangle, double meridian)
+{
+  _latitude=latitude;
+  if (!this->SetDay(day))
+      this->SetDay(0);
+
+  if (!this->SetHourAngle(hourangle))
+      this->SetHourAngle(M_PI/2);
+}
+
+cSun::~cSun(void)
+{
+}
+
+bool cSun::SetDay(int day)
+{
+  if (day >=1 && day <=365)
+  {
+    _day=2.0*M_PI*(day-1.0)/365.0;
+    _declination = 0.006918 - 0.399912*cos(_day) + 0.070257*sin(_day) 
+               - 0.006758*cos(2*_day) + 0.000907*sin(2*_day) - 0.002697*cos(3*_day) 
+               + 0.00148*sin(3*_day);
+
+    CalculateSunrise();
+	return true;
+  }
+  else
+    return false;
+}
+
+
+bool cSun::SetLatitude(double latitude)
+{ 
+  if (latitude >= -M_PI/2 && latitude <=M_PI/2)
+  {
+    _latitude = latitude; 
+	CalculateSunrise();
+    return true;
+  }
+  else
+    return false;
+}
+
+bool cSun::SetHourAngle(double hourangle)
+{
+  _hourangle=hourangle;
+  if (hourangle >= _sunrise && hourangle <= (2*M_PI - _sunrise))
+  {
+    return true;
+  }
+  else
+    return false;
+}
+
+void cSun::calculateAltitude()
+{
+  _altitude=asin(sin(_latitude)*sin(_declination) - cos(_latitude)*cos(_declination)*cos(_hourangle));
+
+}
+
+
+void cSun::calculateAzimuth()
+{
+  double temp;
+
+  temp = (-sin(_latitude)*sin(_altitude) + sin(_declination))/(cos(_latitude)*cos(_altitude));
+
+  if (temp > 1)
+	 _azimuth=0;
+  else if (temp < -1)
+	 _azimuth=M_PI;  
+  else if (_hourangle < M_PI)
+     _azimuth=acos(temp);
+  else
+     _azimuth=2*M_PI - acos(temp);
+
+
+
+}
+
+void cSun::GetPosition(double &Alt, double &Az)
+{
+  calculateAltitude();
+  calculateAzimuth();
+
+ 
+
+  
+  Alt=_altitude;
+  Az=_azimuth;
+}
+
+void cSun::CalculateSunrise()
+{
+	if (tan(_latitude)*tan(_declination)>=1)
+		_sunrise=0;
+	else if (tan(_latitude)*tan(_declination)<=-1)
+		_sunrise=M_PI;
+	else
+		_sunrise=acos(tan(_latitude)*tan(_declination));
+
+}
+
+double cSun::TimeDiff(double Longitude,double Meridian )
+{
+	double Et;
+
+	Et = 229.2 * (0.000075 + 0.001868*cos(_day) - 0.032077*sin(_day) - 0.014615*cos(2*_day) - 0.04089*sin(2*_day));
+	return (-4*(Longitude-Meridian)*180./M_PI + Et)/60;
+
+
+}

bifacial_radiance/data/gencumsky/Sun.h

@@ -0,0 +1,45 @@
+#define _USE_MATH_DEFINES
+#include <math.h>
+class cSun
+{
+public:
+	cSun (double latitude=54*2*M_PI/360, int day=1, double _hourangle=M_PI/2, double _meridian=0.0); // latitude in radians, day is integer 1-365
+	~cSun(void);
+
+	bool SetDay(int day);
+	int GetDay() { return int((_day*365/(2*M_PI)) + 1); }
+
+	bool SetLatitude(double latitude);
+	double GetLatitude() { return _latitude; }
+  
+	void SetMeridian(double meridian) { _meridian=meridian; }
+	double GetMeridian() { return _meridian; }
+
+	bool SetHourAngle(double hourangle);
+	double GetHourAngle() { return _hourangle; }
+
+	void GetPosition (double &Alt, double &Az);  
+	double GetSunrise () { return _sunrise; }
+
+	// Difference between solar and clock time (add result to clock time to get solar)
+	// return value in hours, longitude in degrees(?)
+	double TimeDiff (double Longitude,double Meridian);
+
+private:
+// all of these angles stored in radians
+	double _altitude;
+	double _azimuth;
+	double _latitude;
+	double _hourangle;
+	double _declination;
+	double _sunrise;
+
+	double _meridian;
+
+	double _day;
+
+	void calculateAltitude();
+	void calculateAzimuth();
+	void CalculateSunrise();
+
+};

bifacial_radiance/data/gencumsky/average_val.awk

@@ -0,0 +1,3 @@
+BEGIN { sum=0; count =0; }
+      { sum=sum+$3; count=count+1 }
+END   {printf " %f\n",sum/count; }

bifacial_radiance/data/gencumsky/cPerezSkyModel.cpp

@@ -0,0 +1,238 @@
+#include "cPerezSkyModel.h"
+
+#define _USE_MATH_DEFINES
+#include <math.h>
+#include <stdio.h>
+
+// Perez all weather sky model coefficients
+double cPerezSkyModel::m_a1[8]= {1.3525,-1.2219,-1.1000,-0.5484,-0.6000,-1.0156,-1.0000,-1.0500};
+double cPerezSkyModel::m_a2[8]= {-0.2576,-0.7730,-0.2515,-0.6654,-0.3566,-0.3670,0.0211,0.0289};
+double cPerezSkyModel::m_a3[8]= {-0.2690,1.4148,0.8952,-0.2672,-2.5000,1.0078,0.5025,0.4260};
+double cPerezSkyModel::m_a4[8]= {-1.4366,1.1016,0.0156,0.7117,2.3250,1.4051,-0.5119,0.3590};
+double cPerezSkyModel::m_b1[8]= {-0.7670,-0.2054,0.2782,0.7234,0.2937,0.2875,-0.3000,-0.3250};
+double cPerezSkyModel::m_b2[8]= {0.0007,0.0367,-0.1812,-0.6219,0.0496,-0.5328,0.1922,0.1156};
+double cPerezSkyModel::m_b3[8]= {1.2734,-3.9128,-4.5000,-5.6812,-5.6812,-3.8500,0.7023,0.7781};
+double cPerezSkyModel::m_b4[8]= {-0.1233,0.9156,1.1766,2.6297,1.8415,3.3750,-1.6317,0.0025};
+double cPerezSkyModel::m_c1[8]= {2.8000,6.9750,24.7219,33.3389,21.0000,14.0000,19.0000,31.0625};
+double cPerezSkyModel::m_c2[8]= {0.6004,0.1774,-13.0812,-18.3000,-4.7656,-0.9999,-5.0000,-14.5000};
+double cPerezSkyModel::m_c3[8]= {1.2375,6.4477,-37.7000,-62.2500,-21.5906,-7.1406,1.2438,-46.1148};
+double cPerezSkyModel::m_c4[8]= {1.0000,-0.1239,34.8438,52.0781,7.2492,7.5469,-1.9094,55.3750};
+double cPerezSkyModel::m_d1[8]= {1.8734,-1.5798,-5.0000,-3.5000,-3.5000,-3.4000,-4.0000,-7.2312};
+double cPerezSkyModel::m_d2[8]= {0.6297,-0.5081,1.5218,0.0016,-0.1554,-0.1078,0.0250,0.4050};
+double cPerezSkyModel::m_d3[8]= {0.9738,-1.7812,3.9229,1.1477,1.4062,-1.0750,0.3844,13.3500};
+double cPerezSkyModel::m_d4[8]= {0.2809,0.1080,-2.6204,0.1062,0.3988,1.5702,0.2656,0.6234};
+double cPerezSkyModel::m_e1[8]= {0.0356,0.2624,-0.0156,0.4659,0.0032,-0.0672,1.0468,1.5000};
+double cPerezSkyModel::m_e2[8]= {-0.1246,0.0672,0.1597,-0.3296,0.0766,0.4016,-0.3788,-0.6426};
+double cPerezSkyModel::m_e3[8]= {-0.5718,-0.2190,0.4199,-0.0876,-0.0656,0.3017,-2.4517,1.8564};
+double cPerezSkyModel::m_e4[8]= {0.9938,-0.4285,-0.5562,-0.0329,-0.1294,-0.4844,1.4656,0.5636};
+double cPerezSkyModel::m_PerezClearnessBin[8]={1.065,1.23,1.5,1.95,2.8,4.5,6.2,999999};
+
+// Perez global luminous efficacy coefficients
+double cPerezSkyModel::m_GlobLumEffya[8]={96.63,107.54,98.73,92.72,86.73,88.34,78.63,99.65};
+double cPerezSkyModel::m_GlobLumEffyb[8]={-.47,.79,.7,.56,.98,1.39,1.47,1.86};
+double cPerezSkyModel::m_GlobLumEffyc[8]={11.5,1.79,4.4,8.36,7.1,6.06,4.93,-4.46};
+double cPerezSkyModel::m_GlobLumEffyd[8]={-9.16,-1.19,-6.95,-8.31,-10.94,-7.6,-11.37,-3.15};
+
+// Perez diffuse luminous efficacy coefficients
+double cPerezSkyModel::m_DiffLumEffya[8]={97.24,107.22,104.97,102.39,100.71,106.42,141.88,152.23};
+double cPerezSkyModel::m_DiffLumEffyb[8]={-.46,1.15,2.96,5.59,5.94,3.83,1.90,.35};
+double cPerezSkyModel::m_DiffLumEffyc[8]={12.0,.59,-5.53,-13.95,-22.75,-36.15,-53.24,-45.27};
+double cPerezSkyModel::m_DiffLumEffyd[8]={-8.91,-3.95,-8.77,-13.9,-23.74,-28.83,-14.03,-7.98};
+
+// Perez direct luminous efficacy coefficients
+double cPerezSkyModel::m_BeamLumEffya[8]={57.2,98.99,109.83,110.34,106.36,107.19,105.75,101.18};
+double cPerezSkyModel::m_BeamLumEffyb[8]={-4.55,-3.46,-4.90,-5.84,-3.97,-1.25,.77,1.58};
+double cPerezSkyModel::m_BeamLumEffyc[8]={-2.98,-1.21,-1.71,-1.99,-1.75,-1.51,-1.26,-1.10};
+double cPerezSkyModel::m_BeamLumEffyd[8]={117.12,12.38,-8.81,-4.56,-6.16,-26.73,-34.44,-8.29};
+
+// TODO:  check sunrise/sunset
+
+cPerezSkyModel::cPerezSkyModel(void)
+{
+	m_coefficientsset=false;
+}
+
+cPerezSkyModel::~cPerezSkyModel(void)
+{
+}
+
+bool cPerezSkyModel::SetSkyConditions(double Idh, double Ibh, cSun *Sun)
+{
+	double SolarZenith;
+	double Ibn;
+
+	double PerezBrightness, PerezClearness;
+	double E0,day_angle, AirMass;
+
+	int i, intClearness;
+
+	// if no sun, return no luminance
+	if (Idh <= 0)
+	{
+		m_a=0;
+		m_b=0;
+		m_c=0;
+		m_d=0;
+		m_e=0;
+		m_coefficientsset=true;
+		return false;
+	}
+
+	Sun->GetPosition(m_SolarAlt,m_SolarAz);
+	
+	// store solar zenith
+	SolarZenith = M_PI/2 - m_SolarAlt;
+
+	// calculate clearness
+	if (m_SolarAlt > 0)
+		Ibn=Ibh/sin(m_SolarAlt);
+	else if (m_SolarAlt <=0 && Ibh>0)
+	{
+		// if there's direct horizontal radiation specified but sun is below horizon,
+		// lump it in with the diffuse
+		Idh=Idh+Ibh;
+		Ibn=0;
+	}
+	else
+		Ibn=0;
+
+	PerezClearness = ((Idh+Ibn)/Idh + 1.041*pow(SolarZenith,3))/(1+1.041*pow(SolarZenith,3));
+
+	// calculate brightness
+	// extra terrestrial radiation
+
+	// TODO: check IextraT and AirMass eqns with task 3 microclimate modelling paper
+	day_angle=Sun->GetDay()*2*M_PI/365;
+
+    E0 = 1367 * (1.00011+0.034221*cos(day_angle)+0.00128*sin(day_angle)
+								+0.000719*cos(2*day_angle)+0.000077*sin(2*day_angle));
+
+	// air optical mass
+	if (m_SolarAlt >= 10 * M_PI/180)
+		AirMass=1/sin(m_SolarAlt);
+	else
+//		AirMass=1/(sin(m_SolarAlt) + 0.15*pow(m_SolarAlt*180/M_PI + 3.885,-1.253));
+		AirMass=1/(sin(m_SolarAlt) + 0.50572*pow(180*m_SolarAlt/M_PI+6.07995,-1.6364));
+
+	// fix in case a very negative solar altitude is input
+//	if (m_SolarAlt*180/M_PI + 3.885 >=0)
+	if (m_SolarAlt*180/M_PI + 6.07995 >=0)
+		PerezBrightness=AirMass*Idh/E0;
+	else
+	{
+		// Idh is not zero, but sun altitude < 6 degrees, if Idh is very small (<10) neglect it
+        if (Idh <= 10)
+		{
+			m_a=0;
+			m_b=0;
+			m_c=0;
+			m_d=0;
+			m_e=0;
+			m_coefficientsset=true;
+			return false;
+		}
+
+		// Idh is > 10 and sun altitude v. low, flag up an error and blunder on anyway
+		fprintf(stderr,"Error!  Solar altitude is %.0f < -6 degrees and Idh = %.0f > 10 W/m^2 on day %d !Ibn is %.0f.  Attempting to continue!\n",m_SolarAlt*180/M_PI,Idh,Sun->GetDay(),Ibn);
+		PerezBrightness=0;
+	}
+
+	// TODO: Temporary bit!!!
+	if (PerezBrightness < 0.2 && (PerezClearness > 1.065 && PerezClearness < 2.8)) PerezBrightness=0.2;
+
+	// Now determine the model coefficients
+	// TODO: Error checking
+	if (PerezClearness <1)
+	{
+		//fprintf(stderr,"ERROR! CLEARNESS < 1\n");
+		return false;
+	}
+
+	// find which 'clearness bin' to use (note intClearness is set to one lower than the
+	// tradiational bin numbers (i.e. for clearness bin 1, intClearness=0)
+	for (i=7; i>=0; i--)
+		if (PerezClearness < m_PerezClearnessBin[i]) intClearness=i;
+
+	m_a = m_a1[intClearness] + m_a2[intClearness]*SolarZenith
+		+ PerezBrightness*(m_a3[intClearness] + m_a4[intClearness]*SolarZenith);
+	m_b = m_b1[intClearness] + m_b2[intClearness]*SolarZenith
+		+ PerezBrightness*(m_b3[intClearness] + m_b4[intClearness]*SolarZenith);
+	m_e = m_e1[intClearness] + m_e2[intClearness]*SolarZenith
+		+ PerezBrightness*(m_e3[intClearness] + m_e4[intClearness]*SolarZenith);
+
+	if (intClearness > 0)
+	{
+		m_c = m_c1[intClearness] + m_c2[intClearness]*SolarZenith
+			+ PerezBrightness*(m_c3[intClearness] + m_c4[intClearness]*SolarZenith);
+		m_d = m_d1[intClearness] + m_d2[intClearness]*SolarZenith
+			+ PerezBrightness*(m_d3[intClearness] + m_d4[intClearness]*SolarZenith);
+	}
+	else
+	{
+		// different equations for c & d in clearness bin no. 1
+		m_c=exp(pow(PerezBrightness*(m_c1[intClearness] + m_c2[intClearness]*SolarZenith),m_c3[intClearness])) - 1;
+		m_d=-exp(PerezBrightness*(m_d1[intClearness] + m_d2[intClearness]*SolarZenith)) + m_d3[intClearness] +   m_d4[intClearness]*PerezBrightness;
+	}
+
+	m_coefficientsset=true;
+	m_PerezClearness=PerezClearness;
+	m_IntPerezClearness=intClearness;
+	m_PerezBrightness=PerezBrightness;
+	return true;
+}
+
+double cPerezSkyModel::GetRelativeLuminance(double Alt, double Az)
+{
+	double cosSkySunAngle;
+	double lv;
+
+	if (!m_coefficientsset)
+	{
+		// trying to use model without setting it up
+		printf("Attempt to use model before coefficients are set!\n");
+		return -1;
+	}
+
+	cosSkySunAngle= sin(Alt)*sin(m_SolarAlt) + cos(m_SolarAlt)*cos(Alt)*cos(fabs(Az-m_SolarAz));
+
+	lv=(1 + m_a*exp(m_b/sin(Alt))) * (1 + m_c*exp(m_d*acos(cosSkySunAngle)) + m_e*cosSkySunAngle*cosSkySunAngle);
+	if (lv < 0) lv=0;
+	return lv;
+}
+
+// TODO: Work out W properly!
+// Td - three hourly surface dew point temp (degC)
+double cPerezSkyModel::GetDiffuseLumEffy(double SolarAlt, double Td)
+{
+	double W = 2.0;
+
+	return m_DiffLumEffya[m_IntPerezClearness] + m_DiffLumEffyb[m_IntPerezClearness]*W + m_DiffLumEffyc[m_IntPerezClearness]*sin(SolarAlt)
+		+ m_DiffLumEffyd[m_IntPerezClearness]*log(m_PerezBrightness);
+}
+
+// Td - three hourly surface dew point temp (degC)
+double cPerezSkyModel::GetBeamLumEffy(double SolarAlt, double Td)
+{
+	double W = 2.0;
+	double BeamLumEffy=m_BeamLumEffya[m_IntPerezClearness] + m_BeamLumEffyb[m_IntPerezClearness]*W + m_BeamLumEffyc[m_IntPerezClearness]*exp(5.73*(M_PI/2-SolarAlt)-5)
+		+ m_BeamLumEffyd[m_IntPerezClearness]*m_PerezBrightness;
+
+	if (BeamLumEffy>0)
+		return BeamLumEffy;
+	else
+		return 0;
+}
+
+// Td - three hourly surface dew point temp (degC)
+double cPerezSkyModel::GetGlobalLumEffy(double SolarAlt, double Td)
+{
+	double W = 2.0;
+	double GlobLumEffy=m_GlobLumEffya[m_IntPerezClearness] + m_GlobLumEffyb[m_IntPerezClearness]*W + m_GlobLumEffyc[m_IntPerezClearness]*sin(SolarAlt)
+		+ m_GlobLumEffyd[m_IntPerezClearness]*log(m_PerezBrightness);
+	// careful we don't return a negative efficacy!
+
+	if (GlobLumEffy>0)
+		return GlobLumEffy;
+	else
+		return 0;
+}

bifacial_radiance/data/gencumsky/cPerezSkyModel.h

@@ -0,0 +1,57 @@
+#include "Sun.h"
+class cPerezSkyModel 
+{
+
+public:
+	cPerezSkyModel(void);
+	virtual ~cPerezSkyModel(void);
+
+	// return whether or not it was succesful
+	// cSun object should be set up in the correct position before calling
+	virtual bool SetSkyConditions(double Idh, double Ibh, cSun *Sun);
+	virtual double GetRelativeLuminance(double Alt, double Az);
+	virtual double GetDiffuseLumEffy(double SolarAlt, double Td);
+	virtual double GetGlobalLumEffy(double SolarAlt, double Td);
+	virtual double GetBeamLumEffy(double SolarAlt, double Td);
+
+	double GetSkyClearness() { return m_PerezClearness; }
+
+private:
+	double m_SolarAlt, m_SolarAz; // in radians
+
+	// used to make sure we don't return a luminance before model is set up
+	bool m_coefficientsset;
+
+	// model coefficients
+	// TODO: probably shouldn't declare these as static, but need to find a way
+	// to initialise them
+	static double m_a1[8], m_a2[8], m_a3[8], m_a4[8];
+	static double m_b1[8], m_b2[8], m_b3[8], m_b4[8];
+	static double m_c1[8], m_c2[8], m_c3[8], m_c4[8];
+	static double m_d1[8], m_d2[8], m_d3[8], m_d4[8];
+	static double m_e1[8], m_e2[8], m_e3[8], m_e4[8];
+	static double m_PerezClearnessBin[8];
+
+	// Perez global luminous efficacy coefficients
+	static double m_GlobLumEffya[8];
+	static double m_GlobLumEffyb[8];
+	static double m_GlobLumEffyc[8];
+	static double m_GlobLumEffyd[8];
+
+	// Perez diffuse luminous efficacy coefficients
+	static double m_DiffLumEffya[8];
+	static double m_DiffLumEffyb[8];
+	static double m_DiffLumEffyc[8];
+	static double m_DiffLumEffyd[8];
+
+	// Perez direct luminous efficacy coefficients
+	static double m_BeamLumEffya[8];
+	static double m_BeamLumEffyb[8];
+	static double m_BeamLumEffyc[8];
+	static double m_BeamLumEffyd[8];
+
+	double m_a,m_b,m_c,m_d,m_e;
+	double m_PerezClearness;
+	int m_IntPerezClearness;
+	double m_PerezBrightness;
+};