gendaylit warning statement if negative elevation angle passed.
@@ -719,7 +719,8 @@ class RadianceObj:
# We should already be filtering for elevation >0. But just in case...
if sunalt <= 0:
sunalt = np.arcsin((ghi-dhi)/dni)*180/np.pi # reverse engineer elevation from ghi, dhi, dni
-
+ print('Warning: negative sun elevation passed: {:0.2} with positive ghi. '.format(solpos.elevation)+
+ 'Re-calculated sun elevation: {:0.2}'.format(sunalt))
#" -L %s %s -g %s \n" %(dni/.0079, dhi/.0079, self.ground.ReflAvg) + \
skyStr = ("# start of sky definition for daylighting studies\n"
"# location name: " + str(locName) + " LAT: " + str(lat)
