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

bifacial_radiance/data/gencumsky/gencumulativesky.cpp

@@ -0,0 +1,269 @@
+// GenCumulativeSky was conceived, developed and validated by Darren Robinson and Andrew Stone for efficient solar irradiation modelling using RADIANCE
+// When using GenCumulativeSky they would be pleased if you would ackowledge their work by referring to the following article: "Robinson, D., Stone, A.,
+// Irradiation modeling made simple the cumulative sky approach and its applications, Proc. PLEA 2004, Eindhoven 2004."
+//
+// GenCumulativeSky.cpp : Defines the entry point for the console application.
+//
+
+#include "cSkyVault.h"
+
+
+int main(int argc, char* argv[])
+{
+	cSkyVault sky;
+	double (*patchData)[5] = new double[145][5];
+	int i,j, counter;
+
+	double *CumSky;
+	double dx,dy,dz;
+	double alt,az,deltaalt,deltaaz;
+	double ScalingFactor=1.0;
+
+	double rowdeltaaz[7]={12,12,15,15,20,30,60};
+	//double rowdeltaaz[14]={6,6,6,6,7.5,7.5,7.5,7.5,10,10,15,15,30,30};
+	int rowdeltaalt=12;
+
+	char *filename=argv[argc-1];
+	double hourshift=0;
+
+	bool DoIlluminance=false;
+	bool DoRadiance_179=false;
+
+	cSkyVault::eSunType SunType;
+	cClimateFile::eClimateFileFormat ClimateFileFormat;
+
+	bool DoDiffuse;
+
+	if (!(argc>1))
+	{
+		// User didn't give any command line arguments
+		fprintf(stderr,"GenCumulativeSky: Error - invalid input parameters\n");
+		goto USAGEERROR;
+	}
+
+	counter=1;
+	sky.StartTime=0.0;
+	sky.EndTime=24.0;
+	sky.StartDay=1;
+	sky.EndDay=31;
+	sky.StartMonth=1;
+	sky.EndMonth=12;
+
+	// Set the default parameters
+	SunType=cSkyVault::NO_SUN;
+	DoDiffuse=true;
+	sky.SetLatitude(51.7*M_PI/180);
+	sky.SetLongitude(0*M_PI/180);
+	ClimateFileFormat=cClimateFile::GLOBAL_DIFFUSE;
+
+	while (counter<argc-1)
+	{
+		if (argv[counter][0]=='+' && argv[counter][1]=='s' && argv[counter][2]=='1')
+		{
+			SunType=cSkyVault::CUMULATIVE_SUN;
+			counter++;
+		}
+		else if (argv[counter][0]=='+' && argv[counter][1]=='s' && argv[counter][2]=='2')
+		{
+			SunType=cSkyVault::MANY_SUNS;
+			counter++;
+		}
+		else if (argv[counter][0]=='-' && argv[counter][1]=='d' && argv[counter][2]==' ' )
+		{
+			DoDiffuse=false;
+			counter++;
+		}
+		else if (argv[counter][0]=='-' && argv[counter][1]=='G' )
+		{
+			ClimateFileFormat=cClimateFile::GLOBAL_DIFFUSE;
+			counter++;
+		}
+		else if (argv[counter][0]=='-' && argv[counter][1]=='B' )
+		{
+			ClimateFileFormat=cClimateFile::DIRECTHORIZONTAL_DIFFUSE;
+			counter++;
+		}
+		else if (argv[counter][0]=='-' && argv[counter][1]=='E' )
+		{
+			ClimateFileFormat=cClimateFile::GLOBAL_DIFFUSE_EPW;
+			counter++;
+		}
+		else if (argv[counter][0]=='-' && argv[counter][1]=='l' )
+		{
+			DoIlluminance=true;
+			counter++;
+		}
+		else if (argv[counter][0]=='-' && argv[counter][1]=='r' )
+		{
+			ScalingFactor=1.0/179000;;
+			counter++;
+		}
+		else if (argv[counter][0]=='-' && argv[counter][1]=='p' )
+		{
+			ScalingFactor=1.0/1000;;
+			counter++;
+		}
+		else if (argv[counter][0]=='-' && argv[counter][1]=='a')
+		{
+			if ((argc-counter)>2)
+			{
+				sky.SetLatitude(atof(argv[counter+1])*M_PI/180.);
+				counter+=2;
+			}
+			else
+			{
+				goto USAGEERROR;
+			}
+		}
+		else if (argv[counter][0]=='-' && argv[counter][1]=='h')
+		{
+			if ((argc-counter)>2)
+			{
+				// hourshift - used to alter the default time convention
+				hourshift=atof(argv[counter+1]);
+				counter+=2;
+			}
+			else
+			{
+				goto USAGEERROR;
+			}
+		}
+		else if (argv[counter][0]=='-' && argv[counter][1]=='o' && (argc-counter)>2)
+		{
+			if ((argc-counter)>2)
+			{
+				sky.SetLongitude(atof(argv[counter+1])*M_PI/180.);
+				counter+=2;
+			}
+			else
+			{
+				goto USAGEERROR;
+			}
+		}
+		else if (argv[counter][0]=='-' && argv[counter][1]=='m' && (argc-counter)>2)
+		{
+			if ((argc-counter)>2)
+			{
+				sky.SetMeridian(atof(argv[counter+1])*M_PI/180.);
+				counter+=2;
+			}
+			else
+			{
+				goto USAGEERROR;
+			}
+		}
+		else if (argv[counter][0]=='-' && argv[counter][1]=='t' && argv[counter][2]=='i' && argv[counter][3]=='m' && argv[counter][4]=='e'  )
+		{
+			if ((argc-counter)>3 )
+			{
+				sky.StartTime=atof(argv[counter+1]);
+				sky.EndTime=atof(argv[counter+2]);
+				counter+=3;
+				if(sky.StartTime<0 || sky.StartTime >24 || sky.EndTime<0 || sky.EndTime >24 || sky.EndTime <= sky.StartTime)
+				{
+					fprintf(stderr,"GenCumulativeSky: Error - invalid start time (%.3f) or end time (%.3f)\n",sky.StartTime,sky.EndTime );
+					goto USAGEERROR;
+				}
+			}
+			else
+			{
+				goto USAGEERROR;
+			}
+
+		}
+		else if (argv[counter][0]=='-' && argv[counter][1]=='d' && argv[counter][2]=='a' && argv[counter][3]=='t' && argv[counter][4]=='e'  )
+		{
+			if ((argc-counter)>5 )
+			{
+				sky.StartMonth=atoi(argv[counter+1]);
+				sky.StartDay=atoi(argv[counter+2]);
+				sky.EndMonth=atoi(argv[counter+3]);
+				sky.EndDay=atoi(argv[counter+4]);
+				counter+=5;
+				if(sky.StartDay<0 || sky.StartDay >31.0 ||sky.EndDay<0.0 || sky.EndDay >31.0 ||sky.StartMonth<0.0 || sky.StartMonth >12.0 ||sky.EndMonth<0.0 || sky.EndMonth >12.0)
+				{
+					fprintf(stderr,"GenCumulativeSky: Error - invalid start date (%2.0f//%2.0f) or end date (%2.0f//%2.0f)\n",(double)sky.StartMonth,(double)sky.StartDay,(double)sky.EndMonth,(double)sky.EndDay );
+					goto USAGEERROR;
+				}
+			}
+			else
+			{
+				goto USAGEERROR;
+			}
+
+		}
+		else
+		{
+			fprintf(stderr,"GenCumulativeSky: Error - invalid input parameter '%s'\n\n",argv[counter]);
+			goto USAGEERROR;
+		}
+	}
+
+
+	
+	if (!sky.LoadClimateFile(filename, ClimateFileFormat, sky.StartTime, sky.EndTime, sky.StartDay, sky.EndDay, sky.StartMonth, sky.EndMonth))
+	{
+		fprintf(stderr,"GenCumulativeSky: Error reading climate file %s\n\n",filename);
+		goto USAGEERROR;
+	}
+
+	sky.CalculateSky(SunType, DoDiffuse, DoIlluminance, hourshift);
+
+	CumSky=sky.GetCumulativeSky();
+	sky.GetPatchDetails(patchData);
+
+
+
+	printf("{ This .cal file was generated automatically by GenCumulativeSky }\n");
+	printf("{ ");
+	for (j=0; j<argc; j++)
+		printf("%s ",argv[j]);
+	printf(" }\n\n");
+	printf("skybright=");
+	for (j=0; j<7; j++)
+	{
+		printf("row%d+",j);
+	}
+	printf("row7;\n\n");
+
+	counter=0;
+	for (j=0; j<7; j++)
+	{
+		// note first patch split into two parts - first part (> 0 deg) and last patch (<360)
+		printf("row%d=if(and(alt-%d, %d-alt),select(floor(0.5+az/%5.2f)+1,\n",j,j*rowdeltaalt,(j+1)*rowdeltaalt,rowdeltaaz[j]);
+		for (i=0+counter; i<counter+360/int(rowdeltaaz[j]); i++)
+		{
+			printf("\t%f,\n",CumSky[i]*ScalingFactor);
+		}
+		printf("\t%f),0);\n\n",CumSky[counter]*ScalingFactor);
+		counter+=(int)(360/rowdeltaaz[j]);
+	}
+
+	printf("row7=if(alt-84,%f,0);\n\n",CumSky[144]*ScalingFactor);
+
+	printf("alt=asin(Dz)*180/PI;\n\n");
+
+	printf("az=if(azi,azi,azi+360);\n");
+	printf("azi=atan2(Dx,Dy)*180/PI;\n\n");
+
+	return 0;
+
+USAGEERROR:
+	fprintf(stderr,"Usage: GenCumulativeSky [-d] [+s1|+s2] [-a latitude] [-o longitude] [-l] [-m standard meridian] [-h hourshift] [-G|-B|-E] <climate file>\n");
+	fprintf(stderr,"(Note: longitude +ve East of Greenwich)\n\n");
+	fprintf(stderr,"\t-d\tIgnore diffuse irradiance\n");
+	fprintf(stderr,"\t+s1\tUse \"smeared sun\" approach (default)\n");
+	fprintf(stderr,"\t+s2\tUse \"binned sun\" approach\n");
+	fprintf(stderr,"\t-l\tOutput luminance instead of radiance\n");
+	fprintf(stderr,"\t-r\tOutput radiance/179000 (ensures that units in the Radiance Image Viewer are in kWhm-2)\n");
+	fprintf(stderr,"\t-p\tOutput radiance/1000 (ensures that units in the Radiance RGB data file  are in kWhm-2)\n");
+	fprintf(stderr,"\t-G\tFile format is col1=global irradiance (W/m2), col2=diffuse irradiance (W/m2)\n");
+	fprintf(stderr,"\t-B\tFile format is col1=direct horizontal irradiance (W/m2), col2=diffuse irradiance (W/m2)\n");
+	fprintf(stderr,"\t-E\tFile format is an energyplus weather file (*.epw) The gprogram uses the global irradiance (W/m2) and diffuse irradiance (W/m2) data columns.\n");
+	fprintf(stderr,"\t\tIn combination with \'-E' the considered time interval can be specified:\n");
+	fprintf(stderr,"\t\t-time <start time of day> <end time of day>\n");
+	fprintf(stderr,"\t\t-date mm_start dd_start mm_end dd_end (if start-date after end-date then the winter interval is considered)\n");
+	fprintf(stderr,"\n");
+	return -1;
+}
+

bifacial_radiance/data/gencumsky/make_gencumskyexe.py

@@ -0,0 +1,107 @@
+#! /usr/bin/env python
+"""
+Multiplatform make gencumulativesky executable for Windows, Linux, & Mac OS X
+"""
+
+import logging
+import os
+import sys
+try:
+    from setuptools import distutils
+except ImportError:
+    sys.exit('setuptools was not detected - please install setuptools and pip')
+
+# setup basic logging
+logging.basicConfig()
+LOGGER = logging.getLogger(__name__)
+
+# compiler options
+CC = distutils.ccompiler.new_compiler()  # initialize compiler object
+LIBRARIES = ['m']  # use math library
+# include extra definitions for functions in X/Open Portability Guide
+# https://www.gnu.org/software/libc/manual/html_node/XPG.html
+MACROS = [('_XOPEN_SOURCE', None)]
+# full path to this file
+GENCUMSKY_DIR = os.path.abspath(os.path.dirname(__file__))
+EXE_FILE = 'gencumulativesky'  # name of compiled executable file
+# use g++ executables instead of gcc
+LINKER_EXE = COMPILER = COMPILER_CXX = ['g++']
+
+# platform specific constants
+PLATFORM = sys.platform
+if PLATFORM == 'win32':
+    # math library and XOPG macro are not used on Windows
+    LIBRARIES = None
+    MACROS = [('WIN32', None)]
+elif PLATFORM == 'darwin':
+    pass
+elif PLATFORM in ['linux', 'linux2']:
+    PLATFORM = 'linux'  # consolidate linux names
+else:
+    sys.exit('Platform "%s" is unknown or unsupported.' % PLATFORM)
+# on linux and mac set the compiler and linker to use g++ instead of gcc
+if PLATFORM in ['linux', 'darwin']:
+    CC.set_executables(
+        linker_exe=LINKER_EXE, compiler=COMPILER, compiler_cxx=COMPILER_CXX)
+
+
+def compile_gencumsky(compiler=CC, output_dir=GENCUMSKY_DIR, macros=None,
+                      output_progname=EXE_FILE, libraries=None, logger=LOGGER):
+    """
+    Create platform specific gencumulativesky executable from download.
+    
+    Parameters
+    ----------
+    compiler : `distutils.ccompiler.CCompiler`
+        defaults has linker and cc set to g++ for linux and mac os x
+    output_dir : str
+        path to output folder, default is absolute path to this folder
+    macros : list of 2-item tuples
+        if ``None`` then defaults to ``[('WIN32', None)]`` on Windows, and
+        ``[('_XOPEN_SOURCE', None)]`` for linux and mac os x
+    output_progname : str
+        default is ``'gencumulativesky'``, windows automatically adds ``.exe``
+    libraries : list
+        if ``None`` then defualts to ``['m']`` for linux and mac os x
+    logger : `logging.Logger`
+        default is only critical logging
+    
+    Raises
+    ------
+    `distutils.errors.UnknownFileError`, `distutils.errors.CompileError`,
+    `distutils.errors.LinkError`
+    """
+    logger.debug('compiler: %r', compiler)
+    logger.debug('output directory: %s', output_dir)
+    logger.debug('output program name: %s', output_progname)
+    if macros is None:
+        macros = MACROS  # use default macros
+    logger.debug('macros: %r', macros)
+    if libraries is None:
+        libraries = LIBRARIES  # use default libraries
+    logger.debug('libraries: %r', libraries)
+    # a list of *.cpp source code to compile
+    src = [os.path.join(output_dir, f) for f in os.listdir(output_dir)
+           if f.endswith('.cpp')]
+    # log *.cpp source to compile
+    logger.debug('*.cpp source files:\n\t%s', '\n\t'.join(src))
+    # root of this file, used for compiler output directory
+    root, _ = os.path.splitdrive(output_dir)
+    root += os.sep  # append platform specific path separator
+    logger.debug('root: %s', root)
+    # compile source to objects in build directory
+    objs = compiler.compile(src, output_dir=root, macros=macros)
+    # log objects to compile
+    logger.debug('compiled objects:\n\t%s', '\n\t'.join(objs))
+    # link objects to executable in build directory
+    compiler.link_executable(
+        objs, output_progname, output_dir=output_dir, libraries=libraries)
+
+
+if __name__ == '__main__':
+    LOGGER.setLevel(logging.DEBUG)  # set logging level to debug
+    LOGGER.debug('compiling gencumsky ...\n')
+    compile_gencumsky(logger=LOGGER)  # compile
+    LOGGER.debug(' ... %s successfully compiled!\n', EXE_FILE)
+    os.system(os.path.join(GENCUMSKY_DIR, EXE_FILE))
+    sys.exit(0)  # successful exit!

bifacial_radiance/data/gencumsky/paths.h

@@ -0,0 +1,62 @@
+/* Copyright (c) 1992 Regents of the University of California */
+
+/* SCCSid "@(#)paths.h 2.7 7/21/93 LBL" */
+/* ANSI-C apian@ise.fhg.de */
+
+/*
+ * Definitions for paths on different machines
+ */
+
+#ifdef MSDOS
+
+#define DIRSEP		'/'
+#define ISDIRSEP(c)	((c)=='/' || (c)=='\\')
+#define CASEDIRSEP	case '/': case '\\'
+#define PATHSEP		';'
+#define MAXPATH		128
+#define TEMPLATE	"rtXXXXXX"
+#define TEMPLEN		8
+#define ULIBVAR		"RAYPATH"
+#ifndef DEFPATH
+#define DEFPATH		"c:\\Documents and Settings\\stonea\\My Documents"
+//#define DEFPATH		";c:/ray/lib"
+#endif
+extern char  *fixargv0();
+
+#else
+#ifdef AMIGA
+
+#define DIRSEP		'/'
+#define PATHSEP		';'
+#define MAXPATH		128
+#define TEMPLATE	"/tmp/rtXXXXXX"
+#define TEMPLEN		13
+#define ULIBVAR		"RAYPATH"
+#ifndef DEFPATH
+#define DEFPATH		";/ray/lib"
+#endif
+#define	 fixargv0(a0)	(a0)
+
+#else
+
+#define DIRSEP		'/'
+#define PATHSEP		':'
+#define MAXPATH		256
+#define TEMPLATE	"/usr/tmp/rtXXXXXX"
+#define TEMPLEN		17
+#define ULIBVAR		"RAYPATH"
+#ifndef DEFPATH
+#define DEFPATH		"/usr/local/lib/ray"
+#endif
+#define	 fixargv0(a0)	(a0)
+
+#endif
+#endif
+
+#ifndef ISDIRSEP
+#define ISDIRSEP(c)	((c)==DIRSEP)
+#endif
+#ifndef CASEDIRSEP
+#define CASEDIRSEP	case DIRSEP
+#endif
+

bifacial_radiance/data/ground.rad

@@ -0,0 +1,83 @@
+#ground materials for use in defining ground reflectance etc.
+#experimental reflectance values taken at R(650nm) G(510nm) B(470nm)
+#normalized reflectance: (0.265*R+0.670*G+0.065*B)
+
+#no reflections 0%R
+void plastic black
+0
+0
+5 0.01 0.01 0.01 0 0    
+
+#G173 light soil .21R
+void plastic litesoil
+0
+0
+5 .29 .187 .163 0 0
+
+#Concrete (from SMARTS) albedo 0.28
+void plastic concrete
+0
+0
+5 .364 .257 .2235 0 0
+
+#White EPDM (measured) albedo .807
+void plastic white_EPDM
+0
+0
+5 .863 .791 .745 0 0
+
+#CU Boulder Beige roof (clean) albedo 0.43
+void plastic beigeroof
+0
+0
+5 .547 .395 .339 0 0
+
+#CU Boulder Beige roof (lite soil Case 3) albedo 0.30 case 2: 5 .43 .30 .258 0 0
+void plastic beigeroof_lite
+0
+0
+5 .388 .273 .234 0 0
+
+#CU Boulder Beige roof (heavy soil Case 4) albedo 0.25
+void plastic beigeroof_heavy
+0
+0
+5 .321 .228 .197 0 0
+
+#Asphalt albedo 0.1
+void plastic asphalt
+0
+0
+5 .1 .1 .1 0 .2
+
+#Greyroof albedo 0.62 from SolarWorld calculator
+void plastic greyroof
+0
+0
+5 .62 .62 .62 0 0
+
+void plastic Metal_Aluminum_Anodized
+0
+0
+5 0.792157 0.819608 0.858824 0.2 0.1
+#map_Kd MonoPanel/Metal_Aluminum_Anodized.jpg
+
+void plastic Color_I11
+0
+0
+5 0.000000 0.000000 0.600000 0.1 0.1
+
+void plastic M_0010_Snow
+0
+0
+5 1.000000 0.980392 0.980392 0.1 0.1
+
+void plastic Metal_Grey
+0
+0
+5  0.43 0.448 0.465  0.5 0.5
+
+void glass stock_glass
+0
+0
+3 0.96 0.96 0.96

bifacial_radiance/data/module.json

@@ -0,0 +1,103 @@
+{
+    "PVmod": {
+        "bifi": 1,
+        "glass": false,
+        "modulefile": null,
+        "modulematerial": "black",
+        "numpanels": 1,
+        "offsetfromaxis": 0,
+        "scenex": 2.01,
+        "sceney": 1.0,
+        "scenez": 0.1,
+        "text": "! genbox black PVmod 2 1 0.02 | xform -t -1.0 -0.5 0 -a 1 -t 0 1.0 0",
+        "x": 2,
+        "xgap": 0.01,
+        "y": 1,
+        "ygap": 0.0,
+        "z": 0.02,
+        "zgap": 0.1
+    },
+    "PrismSolar-Bi60": {
+        "bifi": 1,
+        "cellModule": {
+            "centerJB": null,
+            "numcellsx": 6,
+            "numcellsy": 10,
+            "xcell": 0.152,
+            "xcellgap": 0.01,
+            "ycell": 0.152,
+            "ycellgap": 0.01
+        },
+        "glass": false,
+        "modulefile": null,
+        "modulematerial": "black",
+        "numpanels": 1,
+        "offsetfromaxis": 0,
+        "scenex": 0.972,
+        "sceney": 1.61,
+        "scenez": 0.1,
+        "text": "! genbox black cellPVmodule 0.152 0.152 0.02 | xform -t -0.40499999999999997 -0.805 0 -a 6 -t 0.162 0 0 -a 10 -t 0 0.162 0 -a 1 -t 0 1.61 0",
+        "x": 0.962,
+        "xgap": 0.01,
+        "y": 1.61,
+        "ygap": 0.0,
+        "z": 0.02,
+        "zgap": 0.1
+    },
+    "basic-module": {
+        "bifi": 1,
+        "glass": false,
+        "modulefile": null,
+        "modulematerial": "black",
+        "numpanels": 1,
+        "offsetfromaxis": 0,
+        "scenex": 2.01,
+        "sceney": 1.0,
+        "scenez": 0.1,
+        "text": "! genbox black basic-module 2.0 1.0 0.02 | xform -t -1.0 -0.5 0 -a 1 -t 0 1.0 0",
+        "x": 2.0,
+        "xgap": 0.01,
+        "y": 1.0,
+        "ygap": 0.0,
+        "z": 0.02,
+        "zgap": 0.1
+    },
+    "test": {
+        "bifi": 1,
+        "glass": false,
+        "glassEdge": 0.01,
+        "modulefile": null,
+        "modulematerial": "black",
+        "numpanels": 1,
+        "offsetfromaxis": 0,
+        "scenex": 1.6,
+        "sceney": 0.95,
+        "scenez": 0.1,
+        "text": "! genbox black test 1.59 0.95 0.02 | xform -t -0.795 -0.475 0 -a 1 -t 0 0.95 0",
+        "x": 1.59,
+        "xgap": 0.01,
+        "y": 0.95,
+        "ygap": 0.0,
+        "z": 0.02,
+        "zgap": 0.1
+    },
+    "test-module": {
+        "bifi": 1,
+        "glass": false,
+        "glassEdge": 0.01,
+        "modulefile": null,
+        "modulematerial": "black",
+        "numpanels": 1,
+        "offsetfromaxis": 0,
+        "scenex": 1.01,
+        "sceney": 2.0,
+        "scenez": 0.1,
+        "text": "! genbox black test-module 1 2 0.02 | xform -t -0.5 -1.0 0 -a 1 -t 0 2.0 0",
+        "x": 1,
+        "xgap": 0.01,
+        "y": 2,
+        "ygap": 0.0,
+        "z": 0.02,
+        "zgap": 0.1
+    }
+}