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

bifacial_radiance/HPCScripts/BasicSimulations/addNewModule.py

@@ -0,0 +1,15 @@
+import bifacial_radiance
+import os
+
+testfolder = 'TEMP'
+
+if not os.path.exists(testfolder):
+    os.makedirs(testfolder)
+
+rad_obj = bifacial_radiance.RadianceObj('makemod', testfolder)
+
+moduletype='tutorial-module'
+x = 2
+y = 1
+
+rad_obj.makeModule(name=moduletype, x=x, y=y)

bifacial_radiance/HPCScripts/BasicSimulations/dask_on_node.sh

@@ -0,0 +1,11 @@
+#!/usr/bin/bash
+
+PATH=$PATH:$HOME/Radiance/bin
+LD_LIBRARY_PATH=$LD_LIBRARY_PATH:$HOME/Radiance/lib
+RAYPATH=$RAYPATH:$HOME/Radiance/lib
+export PATH
+export LD_LIBRARY_PATH
+export RAYPATH
+
+dask-worker --interface=ib0 --nprocs 18 --nthreads 1 \
+    --scheduler-file=/scratch/sayala/dask_testing/scheduler.json

bifacial_radiance/HPCScripts/BasicSimulations/run_sbatch.sbatch

@@ -0,0 +1,51 @@
+#!/usr/bin/bash
+
+#SBATCH --job-name="demo_run"
+#SBATCH --nodes=2
+#SBATCH --ntasks-per-node=36
+#SBATCH --time=15
+#SBATCH --export=ALL
+#SBATCH --account=pvsoiling
+# --partition=debug
+#SBATCH --mail-type=ALL
+#SBATCH --mail-user=silvana.ovaitt@nrel.gov
+
+# Save info
+python3 /home/sayala/BasicSimulations/addNewModule.py
+
+#-----------------------
+export BASE=/scratch/$USER
+mkdir -p $BASE/$SLURM_JOB_ID
+cd $BASE/$SLURM_JOB_ID
+# Record starting time
+date
+
+# Save info
+cat $0 > $SLURM_JOB_ID.script
+printenv > $SLURM_JOB_ID.env
+cp /home/sayala/BasicSimulations/simulate_tracking_gendaylit.py .
+
+# Start up dask scheduler
+dask-scheduler --interface ib0 \
+    --scheduler-file=/scratch/sayala/dask_testing/scheduler.json &
+
+# Wait for scheduler to start
+sleep 5
+
+# Start up dask worker on all nodes (Note, script is used to also set
+# environment variables on all the nodes. If these were set by default
+# (using bash_profile for example), the commented command below could
+# be used to start up workers.
+srun /home/sayala/BasicSimulations/dask_on_node.sh & 
+
+mkdir RUNS
+cd RUNS
+
+# Wait for workers to start
+sleep 5
+
+# Run script to submit tasks
+python3 /home/sayala/BasicSimulations/simulate_tracking_gendaylit.py
+
+# Record ending time
+date

bifacial_radiance/HPCScripts/BasicSimulations/simulate_fixedtilt_gencumsky.py

@@ -0,0 +1,110 @@
+import numpy as np
+import os
+import pandas as pd
+import time
+import math
+from itertools import chain
+from itertools import product
+import bifacial_radiance
+from dask.distributed import Client
+import math
+
+# Generate spectra for DNI, DHI and albedo using smarts
+
+# Run simulation using the given timestamp and wavelength
+def simulate_single(tilt=None, results_folder_fmt=None, weather_file=None):    
+    
+     # Verify test_folder exists 
+    test_folder = results_folder_fmt.format(f'{tilt:02}')      
+    if not os.path.exists(test_folder):
+        os.makedirs(test_folder)
+
+    # Variables that stay the same
+    #Main Variables needed throughout
+    albedo = 0.6
+    sim_general_name = 'bifacial_example'
+    lat = 37.5
+    lon = -77.6
+    moduletype = 'tutorial-module'
+    pitch = 3
+    clearance_height = 0.2
+    azimuth = 180
+    nMods = 20
+    nRows = 7
+    hpc = True
+
+    pitch = 3
+    clearance_height = 0.2
+    azimuth = 90
+    nMods = 20
+    nRows = 7
+    hpc = True
+
+    sim_name = sim_general_name+'_'+str(tilt)
+    demo = bifacial_radiance.RadianceObj(sim_name,str(test_folder))  
+    demo.setGround(albedo)
+    metdata = demo.readWeatherFile(weather_file) 
+    demo.genCumSky(savefile = sim_name)
+    sceneDict = {'tilt':tilt,'pitch':pitch,'clearance_height':clearance_height,'azimuth':azimuth, 'nMods': nMods, 'nRows': nRows} 
+    scene = demo.makeScene(moduletype=moduletype,sceneDict=sceneDict, hpc=hpc, radname = sim_name)
+    octfile = demo.makeOct(filelist= demo.getfilelist(), octname = demo.basename , hpc=hpc)  
+    analysis = bifacial_radiance.AnalysisObj(octfile=octfile, name=sim_name)
+    frontscan, backscan = analysis.moduleAnalysis(scene=scene)
+    frontdict, backdict = analysis.analysis(octfile, name=sim_name, frontscan=frontscan, backscan=backscan)
+
+
+    results = 1
+
+    return results
+
+
+def run_simulations_dask(tilts, kwargs):
+    # Create client
+    
+    scheduler_file = '/scratch/sayala/dask_testing/scheduler.json'
+    client = Client(scheduler_file=scheduler_file)
+    
+    # Iterate over inputs
+    futures = []
+    
+    # Add Iterations HERE
+
+    for tilt in tilts:
+        futures.append(client.submit(simulate_single, tilt=tilt, **kwargs)) 
+
+    # Get results for all simulations
+    res = client.gather(futures)
+    
+    # Close all dask workers and scheduler
+    try:
+    	client.shutdown()
+    except:
+        pass
+
+    # Close client
+    client.close()
+
+    res = 'FINISHED!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'
+    return res
+
+
+if __name__ == "__main__":
+
+    print(" I AM HERE")
+    # Define locations within file system
+    weather_file = '/home/sayala/WeatherFiles/USA_CO_Boulder-Broomfield-Jefferson.County.AP.724699_TMY3.epw'  
+    results_folder_fmt = '/scratch/sayala/RadianceScenes/BasicSimulations/FixedTilt_Gencumsky/Tilt_{}' 
+
+    # Define inputs    
+    kwargs = {
+        'weather_file': weather_file,
+        'results_folder_fmt': results_folder_fmt
+    }
+    
+    # Array for all hours in the year
+    tilts = [25, 30]
+
+    # Pass variables being looped on, and kwargs
+    run_simulations_dask(tilts, kwargs)
+
+    print("*********** DONE ************")

bifacial_radiance/HPCScripts/BasicSimulations/simulate_fixedtilt_gendaylit.py

@@ -0,0 +1,102 @@
+import numpy as np
+import os
+import pandas as pd
+import time
+import math
+from itertools import chain
+from itertools import product
+import bifacial_radiance
+from dask.distributed import Client
+import math
+
+# Generate spectra for DNI, DHI and albedo using smarts
+
+# Run simulation using the given timestamp and wavelength
+def simulate_single(idx=None, results_folder_fmt=None, weather_file=None):    
+    
+     # Verify test_folder exists 
+    test_folder = results_folder_fmt.format(f'{idx:04}')      
+    if not os.path.exists(test_folder):
+        os.makedirs(test_folder)
+
+    # Variables that stay the same
+    #Main Variables needed throughout
+    albedo = 0.6
+    sim_general_name = 'bifacial_example'
+    lat = 37.5
+    lon = -77.6
+    moduletype = 'tutorial-module'
+    tilt = 10
+    pitch = 3
+    clearance_height = 0.2
+    azimuth = 180
+    nMods = 20
+    nRows = 7
+    hpc = True
+
+    sim_name = sim_general_name+'_'+str(idx)
+    demo = bifacial_radiance.RadianceObj(sim_name,str(test_folder), hpc=True)  
+    demo.setGround(albedo)
+    metdata = demo.readWeatherFile(weather_file) 
+    demo.gendaylit(idx)
+    sceneDict = {'tilt':tilt,'pitch':pitch,'clearance_height':clearance_height,'azimuth':azimuth, 'nMods': nMods, 'nRows': nRows} 
+    scene = demo.makeScene(module=moduletype,sceneDict=sceneDict,radname = sim_name)
+    octfile = demo.makeOct(octname = demo.basename)  
+    analysis = bifacial_radiance.AnalysisObj(octfile=octfile, name=sim_name, hpc=True)
+    frontscan, backscan = analysis.moduleAnalysis(scene=scene)
+    frontdict, backdict = analysis.analysis(octfile, name=sim_name, frontscan=frontscan, backscan=backscan)
+    results = 1
+
+    return results
+
+
+def run_simulations_dask(indices, kwargs):
+    # Create client
+    
+    scheduler_file = '/scratch/sayala/dask_testing/scheduler.json'
+    client = Client(scheduler_file=scheduler_file)
+    
+    # Iterate over inputs
+    futures = []
+    
+    # Add Iterations HERE
+
+    for idx in indices:
+        futures.append(client.submit(simulate_single, idx=idx, **kwargs)) 
+
+    # Get results for all simulations
+    res = client.gather(futures)
+    
+    # Close all dask workers and scheduler
+    try:
+    	client.shutdown()
+    except:
+        pass
+
+    # Close client
+    client.close()
+
+    res = 'FINISHED!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'
+    return res
+
+
+if __name__ == "__main__":
+
+    # Define locations within file system
+    weather_file = '/home/sayala/WeatherFiles/USA_CO_Boulder-Broomfield-Jefferson.County.AP.724699_TMY3.epw'  
+    results_folder_fmt = '/scratch/sayala/RadianceScenes/BasicSimulations/FixedTilt_Gendaylit/Timestamp_{}' 
+
+    # Define inputs    
+    kwargs = {
+        'weather_file': weather_file,
+        'results_folder_fmt': results_folder_fmt
+    }
+    
+    # Array for all hours in the year
+    indices = np.array(list(range(0, 8760)))    
+    indices = np.array(list(range(4020, 4022)))
+
+    # Pass variables being looped on, and kwargs
+    run_simulations_dask(indices, kwargs)
+
+    print("*********** DONE ************")

bifacial_radiance/HPCScripts/BasicSimulations/simulate_tracking_gendaylit.py

@@ -0,0 +1,126 @@
+import numpy as np
+import os
+import pandas as pd
+import time
+import math
+from itertools import chain
+from itertools import product
+import bifacial_radiance
+from dask.distributed import Client
+import math
+
+# Generate spectra for DNI, DHI and albedo using smarts
+
+# Run simulation using the given timestamp and wavelength
+def simulate_single(daydate=None, results_folder_fmt=None, weather_file=None):    
+    
+     # Verify test_folder exists 
+    test_folder = results_folder_fmt.format(f'{daydate}')      
+    if not os.path.exists(test_folder):
+        os.makedirs(test_folder)
+
+    # Variables that stay the same
+    #Main Variables needed throughout
+    albedo = 0.6
+    sim_general_name = 'bifacial_example'
+    lat = 37.5
+    lon = -77.6
+    moduletype = 'tutorial-module'
+    gcr = 0.35
+    hub_height = 0.2
+
+    nMods = 20
+    nRows = 7
+    hpc = True
+    cumulativesky = False
+
+    limit_angle = 60
+    backtrack = True
+
+
+    sim_name = sim_general_name+'_'+str(daydate)
+    demo = bifacial_radiance.RadianceObj(sim_name,str(test_folder))  
+    demo.setGround(albedo)
+    metdata = demo.readWeatherFile(weather_file, coerce_year=2021)
+    sceneDict = {'gcr':gcr,'hub_height':hub_height, 'nMods': nMods, 'nRows': nRows} 
+    trackerdict = demo.set1axis(limit_angle = limit_angle, backtrack = backtrack, gcr = gcr, cumulativesky = cumulativesky)
+
+    # Restrict trackerdict here
+    #foodict = {k: v for k, v in trackerdict.items() if k.startswith('21'+'_'+day_date)}
+    #trackerdict = demo.gendaylit1axis(trackerdict = foodict)
+    trackerdict = demo.gendaylit1axis(startdate=daydate, enddate=daydate)
+    trackerdict = demo.makeScene1axis(moduletype=moduletype,sceneDict=sceneDict, cumulativesky=cumulativesky, hpc=hpc) #makeScene creates a .rad file with 20 modules per row, 7 rows.
+    trackerdict = demo.makeOct1axis(customname = sim_name, hpc=hpc)
+    demo.analysis1axis(customname = sim_name, hpc=hpc)
+
+    results = 1
+
+    return results
+
+
+def run_simulations_dask(daylist, kwargs):
+    # Create client
+    
+    scheduler_file = '/scratch/sayala/dask_testing/scheduler.json'
+    client = Client(scheduler_file=scheduler_file)
+    
+    # Iterate over inputs
+    futures = []
+    
+    # Add Iterations HERE
+
+    for daydate in daylist:
+        futures.append(client.submit(simulate_single, daydate=daydate, **kwargs)) 
+
+    # Get results for all simulations
+    res = client.gather(futures)
+    
+    # Close all dask workers and scheduler
+    try:
+    	client.shutdown()
+    except:
+        pass
+
+    # Close client
+    client.close()
+
+    res = 'FINISHED!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'
+    return res
+
+
+if __name__ == "__main__":
+    print(">>>>>>>>>>>>>>>>>> STARTING !")
+    # Define locations within file system
+    weather_file = '/home/sayala/WeatherFiles/USA_CO_Boulder-Broomfield-Jefferson.County.AP.724699_TMY3.epw'  
+    results_folder_fmt = '/scratch/sayala/RadianceScenes/BasicSimulations/Gendaylit1axis/Day_{}' 
+
+    # Define inputs    
+    kwargs = {
+        'weather_file': weather_file,
+        'results_folder_fmt': results_folder_fmt
+    }
+    
+    import datetime
+
+    FullYear = False
+
+    if FullYear:
+        start = datetime.datetime.strptime("01-01-2021", "%d-%m-%Y")
+        end = datetime.datetime.strptime("31-12-2021", "%d-%m-%Y") # 2014 not a leap year.
+    else:
+        start = datetime.datetime.strptime("24-03-2021", "%d-%m-%Y")
+        end = datetime.datetime.strptime("27-03-2021", "%d-%m-%Y") # 2014 not a leap year.
+
+    date_generated = [start + datetime.timedelta(days=x) for x in range(0, (end-start).days)]
+
+    daylist = []
+    for date in date_generated:
+        daylist.append(date.strftime("%y_%m_%d"))
+    # loop doesn't add last day :
+    daylist.append('21_12_31')
+
+
+    # Pass variables being looped on, and kwargs
+    run_simulations_dask(daylist, kwargs)
+
+    print("*********** DONE ************")

bifacial_radiance/HPCScripts/Other Examples (unorganized)/PuertoRico.py

@@ -0,0 +1,168 @@
+import numpy as np
+import os
+import pandas as pd
+import time
+import math
+from itertools import chain
+import bifacial_radiance
+from dask.distributed import Client
+from math import sin, cos, radians
+
+
+# Run simulation using the given timestamp 
+def simulate_single(xgap=None, numpanels = None, sensorx=None,
+                   test_folder_fmt=None, weather_file=None):    
+
+    ft2m = 0.3048
+    hub_height = 8.0 * ft2m
+    y = 1
+    pitch = 0.001 # y * np.cos(np.radians(tilt))+D
+    ygap = 0.15
+    tilt = 18
+
+    sim_name = ('Coffee_'+str(numpanels)+'up_'+
+                str(round(xgap,1))+'_xgap_'+str(sensorx)+'posx')
+
+    # Verify test_folder exists before creating radiance obj
+    test_folder = test_folder_fmt.format(f'{numpanels}', 
+                                        f'{round(xgap,1)}',f'{sensorx:03}')
+
+    if not os.path.exists(test_folder):
+        os.makedirs(test_folder)
+
+    lat = 18.202142
+    lon = -66.759187
+    albedo = 0.25 # Grass value from Torres Molina, "Measuring UHI in Puerto Rico" 18th LACCEI 
+                # International Multi-Conference for Engineering, Education, and Technology
+    x = 1.64      
+
+    azimuth = 180
+    if numpanels == 3:
+        nMods = 9
+    if numpanels == 4:
+        nMods = 7
+    nRows = 1
+    moduletype = 'PR_'+str(numpanels)+'up_'+str(round(xgap,1))+'xgap'
+
+    hpc = False
+
+    demo = bifacial_radiance.RadianceObj(sim_name,str(test_folder))  
+    demo.setGround(albedo)
+    demo.readWeatherFile(weather_file)
+
+    sceneDict = {'tilt':tilt,'pitch':pitch,'hub_height':hub_height,'azimuth':azimuth, 'nMods': nMods, 'nRows': nRows} 
+    scene = demo.makeScene(moduletype=moduletype,sceneDict=sceneDict, hpc=hpc, radname = sim_name)
+
+    demo.genCumSky()
+
+    octfile = demo.makeOct(filelist = demo.getfilelist(), octname = demo.basename, hpc=hpc)  
+    analysis = bifacial_radiance.AnalysisObj(octfile=octfile, name=sim_name)
+
+    ii = 1
+    jj = 1
+    
+    simplesim=False
+    if simplesim:
+        sensorsy_front = 20  
+        sensorsy_back = 1
+        sensorsx_front = 1
+        sensorsx_back = 1
+    else:
+        sensorsy_front = 201
+        sensorsy_back = 1
+        sensorsx_front = 1
+        sensorsx_back = 1
+
+
+    frontscan, backscan = analysis.moduleAnalysis(scene, sensorsy_front=sensorsy_front,
+                                                        sensorsx_front=sensorsx_front,
+                                                        sensorsy_back=sensorsy_back,
+                                                        sensorsx_back=sensorsx_back)
+
+ 
+    xinc = 0.10 # cm . Sampling every 10 cm. 
+    yinc = 0.10 # cm  . Sampling every 10 cm
+
+    frontscan['zstart'] = 0.01
+    frontscan['xstart'] = -20+sensorx*xinc
+    frontscan['orient'] = '0 0 -1'
+    frontscan['zinc'] = 0
+    frontscan['yinc'] = yinc
+    frontscan['ystart'] = -10 # n
+    frontdict, backdict = analysis.analysis(octfile = octfile, name = 'xloc_'+str(sensorx), 
+                                            frontscan=frontscan, backscan=backscan)
+
+    results = 1
+
+    print("***** Finished simulation for "+ str(sim_name))
+
+    results=1
+    return results
+
+
+
+def run_simulations_dask(xgaps, numpanelss, sensorsxs, kwargs):
+    # Create client
+    
+    scheduler_file = '/scratch/sayala/dask_testing/scheduler.json'
+    client = Client(scheduler_file=scheduler_file)
+    
+    # Iterate over inputs
+    futures = []
+    
+    for nn in range (0, len(numpanelss)):
+        numpanels = numpanelss[nn]
+        for xx in range (0, len(xgaps)):
+            xgap = xgaps[xx]
+            for ii in sensorsxs:
+                futures.append(client.submit(simulate_single, xgap=xgap, numpanels=numpanels, sensorx=ii,**kwargs))
+
+    # Get results for all simulations
+    res = client.gather(futures)
+    
+    # Close all dask workers and scheduler
+    try:
+    	client.shutdown()
+    except:
+        pass
+
+    # Close client
+    client.close()
+
+    res = 'FINISHED!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'
+    return res
+
+if __name__ == "__main__":
+    # Define locations within file system
+
+    weather_file = '/scratch/sayala/JORDAN/PRI_Mercedita.AP.785203_TMY3.epw'
+    test_folder_fmt = '/scratch/sayala/JORDAN/PUERTO_RICO_NEW/numpanels_{}_xgap_{}_Posx_{}'
+    
+   
+    # Define inputs    
+    kwargs = {
+        'weather_file': weather_file,
+        'test_folder_fmt': test_folder_fmt
+    }
+    
+#    indices = np.array(list(range(2881, 6552)))
+    ft2m = 0.3048
+    lat = 18.202142
+
+    # Loops
+
+    # simplesim
+    simplesim=False
+    if simplesim:
+        xgaps = np.array([3, 4])*ft2m
+        numpanelss = np.array([3, 4])
+        sensorsxs = 2
+    else:
+        xgaps = np.array([3, 4, 6, 9, 12, 15, 18, 21]) * ft2m
+        numpanelss = np.array([3, 4])
+        sensorsxs = np.array(list(range(0, 401)))   
+
+    # Specify method for running simulation
+    use_dask = True
+    if use_dask:
+        run_simulations_dask(xgaps, numpanelss, sensorsxs, kwargs)