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

bifacial_radiance/data/gencumsky/climateFile.cpp

@@ -0,0 +1,312 @@
+#include "climateFile.h"
+#include <stdio.h>
+#include <iostream>
+
+
+#include "paths.h"
+
+cClimateFile::cClimateFile(void)
+{
+	// set pointers to point to zero initially
+	m_ptIgh=0;
+	m_ptIdh=0;
+}
+
+cClimateFile::~cClimateFile(void)
+{
+	delete[] m_ptIgh;
+	delete[] m_ptIdh;
+}
+
+bool cClimateFile::ReadClimateFile(char *FileName, int HourConvention, eClimateFileFormat ClimateFileFormat,double StartTime, double EndTime, int StartDay, int EndDay, int StartMonth, int EndMonth)
+{
+	FILE *InputFile;
+	char *Line = new char[1000];
+
+	int hour,day;
+	int i;
+	float n1,n2,n3,n4,n5,n7,n8,n9,n10,n11,n12,n13,n14,n15,n16;
+	char line_string[1000]="";
+	char s7[1000]="";
+
+	// open file (checking that you can...)
+	if ((InputFile=LoadFile(FileName))==NULL)
+	{
+		fprintf(stderr,"Error opening: %s\n",FileName);
+		return false;
+	}
+	rewind(InputFile);
+	// TODO: Add error checking
+
+	// get rid of old climate data
+	delete[] m_ptIgh;
+	delete[] m_ptIdh;
+
+	//Assume we've got 8760 points
+	m_NumPoints=8760;
+
+	double *Col1 = new double[m_NumPoints];
+	double *Col2 = new double[m_NumPoints];
+
+	// read in points for each hour of each day
+	// points are assumed to be at 0:30, 1:30, ... from 1st Jan to 3st Dec
+	// The function assumes  a full hourly annual data file with 8760 lines.
+
+	//if EPW file format
+	if (ClimateFileFormat==GLOBAL_DIFFUSE_EPW)
+	{
+		//skip file header:
+		for (i=0; i<8; i++)
+		{
+			fscanf(InputFile,"%*[^\n]");
+			fscanf(InputFile,"%*[\n\r]");
+		}
+		for (day=0; day<365; day++)
+		{
+			for (hour=0; hour<24; hour++)
+			{
+				
+				fscanf(InputFile,"%f,%f,%f,%f,%f,%s",&n1,&n2,&n3,&n4,&n5,&line_string[0]);
+				//replace ',' with ' '
+				for (i=0;i<1000;i++){
+					if(line_string[i] == ',')
+						line_string[i] = ' ';
+				}
+				if (sscanf(line_string,"%s %f %f %f %f %f %f %f %f %f  %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f",
+					&s7[0],&n7,&n8,&n9,&n10,&n11,&n12,&n13,&n14,&n15,&n16,&n1,&n1,&n1,&n1,&n1,&n1,&n1,&n1,&n1,&n1,&n1,&n1,&n1,&n1,&n1)== EOF)
+				{
+					// ran out of file
+					fprintf(stderr,"Error processing climate file %s\n",FileName);
+					return false;
+				}
+				// StartDay, EndDay, StartMonth, EndMonth);
+				if(hour>= (StartTime) && hour< (EndTime-1))
+				{
+					if(JulianDate(StartMonth, StartDay)<=JulianDate(EndMonth, EndDay))
+					{
+						if(day>=JulianDate(StartMonth, StartDay) && day<=JulianDate(EndMonth, EndDay))
+						{
+							Col1[day*24+hour]=n14;
+							Col2[day*24+hour]=n16;
+						}
+					} else if(JulianDate(StartMonth, StartDay)>JulianDate(EndMonth, EndDay))
+					{
+						if(day>=JulianDate(StartMonth, StartDay) || day<=JulianDate(EndMonth, EndDay))
+						{
+							Col1[day*24+hour]=n14;
+							Col2[day*24+hour]=n16;
+						}
+					}
+				}
+				//printf("%.0f %.0f %.0f %.0f %.0f\n",n2,n3,n4,Col1[day*24+hour],Col2[day*24+hour]);
+			}
+		}
+
+
+
+	}else{
+		for (day=0; day<365; day++)
+		{
+			for (hour=0; hour<24; hour++)
+			{
+				if (fgets(Line,100,InputFile)==NULL)
+				{
+					// ran out of file
+					fprintf(stderr,"Error processing climate file %s\n",FileName);
+					return false;
+				}
+
+				sscanf(Line,"%lf %lf\n",&Col1[day*24+hour],&Col2[day*24+hour]);
+
+
+			}
+		}
+	}
+
+	// Now convert data into appropriate form
+	if (ClimateFileFormat==GLOBAL_DIFFUSE || ClimateFileFormat==GLOBAL_DIFFUSE_EPW)
+	{
+		m_ptIgh=Col1;
+		m_ptIdh=Col2;
+
+	}
+	else if (ClimateFileFormat == DIRECTHORIZONTAL_DIFFUSE)
+	{
+		m_ptIdh=Col2;
+		for (i=0; i<m_NumPoints; i++)
+		{
+			Col1[i]=Col1[i]+Col2[i];
+		}
+		m_ptIgh=Col1;
+	}
+	else
+	{
+		// ran out of file
+		fprintf(stderr,"Unknown climate file format!\n");
+		return false;
+	}
+
+
+	if (!ValidateData())
+	{
+		fprintf(stderr,"Invalid Data!\n");
+		return false;
+	}
+	fclose(InputFile);
+	return true;
+}
+
+
+bool cClimateFile::ValidateData()
+{
+	// TODO: Write this function!
+	return true;
+}
+
+int cClimateFile::JulianDate(int month, int day)
+{
+	// calculate julian date
+	int jd=0;
+	if(month==1)
+		jd=day;
+	if(month==2)
+		jd=31+day;
+	if(month==3)
+		jd=59+day;
+	if(month==4)
+		jd=90+day;
+	if(month==5)
+		jd=120+day;
+	if(month==6)
+		jd=151+day;
+	if(month==7)
+		jd=181+day;
+	if(month==8)
+		jd=212+day;
+	if(month==9)
+		jd=243+day;
+	if(month==10)
+		jd=273+day;
+	if(month==11)
+		jd=304+day;
+	if(month==12)
+		jd=334+day;
+	return (jd-1);
+}
+double cClimateFile::GetDirectRad(double hour, int day)
+{
+	// interpolate linearly between two closest hours
+	double hourdiff, hour1;
+	int pointer;
+
+	// find the half hour before current time (climate file has data for 0:30, 1:30,2:30, etc...)
+	hour1=(int)hour;
+	hourdiff=hour-hour1;
+	if (hourdiff >= 0.5) hour1+= 0.5;
+	else hour1-=0.5;
+
+	pointer=(day-1)*24+(int)(hour1-0.5);
+
+	// if pointer < 0, don't have data for specified day
+	if (pointer < 0) return -9999;
+
+	float ratio, diff;
+
+	ratio=hour-hour1;
+	// TODO: CHeck this (first/last hours of day)
+	diff=(m_ptIgh[pointer+1]-m_ptIdh[pointer+1])-(m_ptIgh[pointer]-m_ptIdh[pointer]);
+
+	return (m_ptIgh[pointer]-m_ptIdh[pointer]) + diff*ratio;
+}
+
+double cClimateFile::GetDiffuseRad(double hour, int day)
+{
+	// interpolate linearly between two closest hours
+	float hour1, hourdiff;
+	int pointer;
+
+	// find the half hour before current time (climate file has data for 0:30, 1:30,2:30, etc...)
+	hour1=int(hour);
+	hourdiff=hour-hour1;
+	if (hourdiff >= 0.5) hour1+= 0.5;
+	else hour1-=0.5;
+
+	pointer=(day-1)*24+(int)(hour1-0.5);
+
+	// if pointer < 0, don't have data for specified day
+	if (pointer < 0) return -9999;
+
+	float ratio, diff;
+
+	ratio=hour-hour1;
+	diff=m_ptIdh[pointer+1]-m_ptIdh[pointer];
+
+	return m_ptIdh[pointer] + diff*ratio;
+}
+
+double cClimateFile::GetGlobalRad(double hour, int day)
+{
+	// interpolate linearly between two closest hours
+	float hour1, hourdiff;
+	int pointer;
+
+	// find the half hour before current time (climate file has data for 0:30, 1:30,2:30, etc...)
+	hour1=int(hour);
+	hourdiff=hour-hour1;
+	if (hourdiff >= 0.5) hour1+= 0.5;
+	else hour1-=0.5;
+
+	pointer=(day-1)*24+(int)(hour1-0.5);
+
+	// if pointer < 0, don't have data for specified day
+	if (pointer < 0) return -9999;
+
+	float ratio, diff;
+
+	ratio=hour-hour1;
+	diff=m_ptIgh[pointer+1]-m_ptIgh[pointer];
+	return m_ptIgh[pointer] + diff*ratio;
+}
+
+
+// THIS ROUTINE 'BORROWED' FROM GENDAYLIT
+//FILE* cClimateFile::LoadFile( char *fname)			/* find file and open for reading */
+//{
+//	FILE  *fp;
+//	char  pname[MAXPATH];
+//	char *libpath=NULL;
+//	register char  *sp, *cp;
+
+//	if (fname == NULL)
+//		return(NULL);
+
+//	if (ISDIRSEP(fname[0]) || fname[0] == '.')	/* absolute path */
+//		return(fopen(fname, "r"));
+
+//	if (libpath == NULL) {			/* get search path */
+//		libpath = getenv(ULIBVAR);
+//		if (libpath == NULL)
+//			libpath = DEFPATH;
+//	}
+						/* check search path */
+//	sp = libpath;
+//	do {
+//		cp = pname;
+//		while (*sp && (*cp = *sp++) != PATHSEP)
+//			cp++;
+//		if (cp > pname && !ISDIRSEP(cp[-1]))
+//			*cp++ = DIRSEP;
+//		strcpy(cp, fname);
+//		if ((fp = fopen(pname, "r")) != NULL)
+//			return(fp);			/* got it! */
+//	} while (*sp);
+//						/* not found */
+//	return(NULL);
+//}
+
+FILE* cClimateFile::LoadFile( char *fname)	/*open file for reading*/
+{	FILE *Datei;
+	Datei = fopen(fname, "r");
+	if ( Datei == NULL) fprintf(stderr,		  "fatal warning FUNCTION LoadFile: open of %s for input failed.\n",fname);
+	return Datei;}
+

bifacial_radiance/data/gencumsky/climateFile.h

@@ -0,0 +1,37 @@
+#include <stdio.h>
+#include <stdlib.h>
+
+class cClimateFile
+{
+public:
+	enum eClimateFileFormat { GLOBAL_DIFFUSE, DIRECTHORIZONTAL_DIFFUSE,GLOBAL_DIFFUSE_EPW };
+
+	cClimateFile(void);
+	~cClimateFile(void);
+
+	// read in climate file
+	// must specify the hour numbering convention -1 = hour beginning, 0 = centered, 1 = hour ending
+	bool ReadClimateFile(char *FileName, int HourConvention, eClimateFileFormat ClimateFileFormat,double StartTime, double EndTime, int StartDay, int EndDay, int StartMonth, int EndMonth);
+
+	// check that the currently stored data is valid
+	bool ValidateData();
+
+	//caclulate julian date
+	int JulianDate(int month, int day);
+	//int cClimateFile::JulianDate(int month, int day);
+
+	// return direct radiation for a given day (Ibh=Igh-Idh)
+	// day 1 = 1st jan, hour is 0-23
+	double GetDirectRad(double hour, int day);
+	double GetDiffuseRad(double hour, int day);
+	double GetGlobalRad(double hour, int day);
+
+private:
+	// climate data
+	double *m_ptIgh, *m_ptIdh;
+	// number of data points (almost certainly = 24*365=8760)
+	int m_NumPoints;
+
+	FILE *LoadFile(char *fname);
+
+};

bifacial_radiance/data/gencumsky/cumulative.cal

@@ -0,0 +1,177 @@
+{ This .cal file was generated automatically by gencumulativesky }
+{ a.out +s2 -a 51.5 -h -0. -o 5.32 -G ./cardiff.txt  }
+
+skybright=row0+row1+row2+row3+row4+row5+row6+row7;
+
+row0=if(and(alt-0, 12-alt),select(floor(0.5+az/12.00)+1,
+	147651.992082,
+	148103.246847,
+	150502.848057,
+	155558.154335,
+	164752.750913,
+	180226.864316,
+	197570.397547,
+	212571.942436,
+	222715.516579,
+	229380.244970,
+	234364.512080,
+	238724.153999,
+	241687.185432,
+	243156.625884,
+	244030.858216,
+	245287.056661,
+	247567.996594,
+	250779.009968,
+	253345.247304,
+	253483.003437,
+	251441.577429,
+	248002.276951,
+	241899.370743,
+	231427.225726,
+	213790.778047,
+	191226.997091,
+	171722.489942,
+	159518.523684,
+	152621.012229,
+	149021.950490,
+	147651.992082),0);
+
+row1=if(and(alt-12, 24-alt),select(floor(0.5+az/12.00)+1,
+	134651.743021,
+	135053.083596,
+	137437.590856,
+	142420.217077,
+	151274.635553,
+	165867.004901,
+	187839.978362,
+	209363.153328,
+	225731.564198,
+	236851.521055,
+	243759.655986,
+	250020.623769,
+	255473.229880,
+	259620.185535,
+	262936.931112,
+	265882.298494,
+	267908.171928,
+	269448.819016,
+	270151.095768,
+	269505.390092,
+	266983.874338,
+	261333.811133,
+	252929.903479,
+	233597.041538,
+	206371.066468,
+	178675.901426,
+	159135.575806,
+	147001.511408,
+	139925.800854,
+	136142.123374,
+	134651.743021),0);
+
+row2=if(and(alt-24, 36-alt),select(floor(0.5+az/15.00)+1,
+	125460.067551,
+	126133.634995,
+	129825.423207,
+	137678.123424,
+	151862.426574,
+	174915.500053,
+	205419.015484,
+	231446.717882,
+	247028.536763,
+	255945.617304,
+	263118.023397,
+	268575.983401,
+	272514.849105,
+	275030.006794,
+	276372.309258,
+	275336.988349,
+	270847.601692,
+	260369.168109,
+	233101.528953,
+	193610.515939,
+	163540.284837,
+	144358.286209,
+	133324.367364,
+	127615.967316,
+	125460.067551),0);
+
+row3=if(and(alt-36, 48-alt),select(floor(0.5+az/15.00)+1,
+	120757.784075,
+	121324.947226,
+	124596.806475,
+	131288.409911,
+	142715.794089,
+	160550.960188,
+	185946.724621,
+	216780.394666,
+	240382.834313,
+	257157.830766,
+	264688.413196,
+	270191.252794,
+	275185.652566,
+	277586.860612,
+	276764.236991,
+	274128.718481,
+	266130.058101,
+	240671.534981,
+	206302.822953,
+	175137.206287,
+	152331.985627,
+	137185.051179,
+	127891.996560,
+	122785.088706,
+	120757.784075),0);
+
+row4=if(and(alt-48, 60-alt),select(floor(0.5+az/20.00)+1,
+	120553.214945,
+	121701.979622,
+	127264.207825,
+	138523.623777,
+	157362.196532,
+	184812.504009,
+	218174.248601,
+	244754.046815,
+	260098.694850,
+	268529.254096,
+	270470.791981,
+	260536.512623,
+	234959.274998,
+	200472.817437,
+	168938.370311,
+	145955.737209,
+	131451.317663,
+	123544.008800,
+	120553.214945),0);
+
+row5=if(and(alt-60, 72-alt),select(floor(0.5+az/30.00)+1,
+	125136.944026,
+	127825.032912,
+	138889.529604,
+	159856.698781,
+	189095.038539,
+	217166.299527,
+	231209.021293,
+	224177.613259,
+	198836.252107,
+	168210.150741,
+	144143.189720,
+	130187.870112,
+	125136.944026),0);
+
+row6=if(and(alt-72, 84-alt),select(floor(0.5+az/60.00)+1,
+	135639.142343,
+	143917.146301,
+	168017.896278,
+	185262.329733,
+	172003.700483,
+	146816.867657,
+	135639.142343),0);
+
+row7=if(alt-84,154421.194038,0);
+
+alt=asin(Dz)*180/PI;
+
+az=if(azi,azi,azi+360);
+azi=atan2(Dx,Dy)*180/PI;
+

bifacial_radiance/data/gencumsky/cumulative.rad

@@ -0,0 +1,14 @@
+void brightfunc skyfunc
+2 skybright cumulative.cal
+0
+0 
+skyfunc glow sky_glow
+0
+0
+4 1 1 1 0
+sky_glow source sky
+0
+0
+4 0 0 1 180
+
+

bifacial_radiance/data/gencumsky/cumulativesky_rotated.rad

@@ -0,0 +1,2 @@
+!xform -rz -97 cumulativesky.rad
+!xform -rz -97 SunFile.rad