ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/doc/man/man1/gendaylit.1
(Generate patch)

Comparing ray/doc/man/man1/gendaylit.1 (file contents):
Revision 1.4 by greg, Fri Aug 9 16:44:19 2013 UTC vs.
Revision 1.5 by greg, Tue Oct 9 22:33:31 2018 UTC

# Line 1 | Line 1
1 < .\" RCSid "$Id$"
2 < .TH GENDAYLIT 1 4/12/94 "RADIANCE ISE/ADEME EXTENSIONS"
3 < .SH NAME
1 > .TH "GENDAYLIT" "1" "8/31/18" "RADIANCE ISE/ADEME EXTENSIONS" ""
2 > .SH "NAME"
3   gendaylit - generates a RADIANCE description of the daylight sources using Perez models for direct and diffuse components
5 .SH SYNOPSIS
6 .B "gendaylit month day hour [-P|-W|-L|-G|-E] input_value(s) "
7 [
8 .B options
9 ]
4   .br
11 .B "gendaylit -ang altitude azimuth [-P|-W|-L|-G|-E] input_value(s) "
12 [
13 .B options
14 ]
15 .SH DESCRIPTION
16 .I Gendaylit
17 produces a RADIANCE scene description based on an angular distribution of the
18 daylight sources (direct+diffuse) for the given atmospheric conditions
19 (direct and diffuse component of the solar radiation), date and
20 local standard time. The default output is the radiance of the sun (direct) and the sky (diffuse)
21 integrated over the visible spectral range (380-780 nm). We have used the
22 calculation of the sun's position and the ground brightness models which
23 were programmed in
24 .I gensky.                                              
5  
6 < The diffuse angular distribution is calculated using the Perez et al.
7 < sky luminance distribution model (see Solar Energy Vol. 50, No. 3, pp. 235-245, 1993) which, quoting Perez,
8 < describes "the mean instantaneous sky luminance angular distribution patterns for all sky
9 < conditions from overcast to clear, through partly cloudy, skies". The correctness of the
10 < resulting sky radiance/luminance values in this simulation is ensured through the normalization of the modelled
31 < sky diffuse to the measured sky diffuse irradiances/illuminances.
6 > .SH "SYNOPSIS"
7 > \fBgendaylit month day hour [-P|-W|-L|-G|-E] input_value(s) \fR[ \fBoptions \fR]
8 > .br
9 > \fBgendaylit -ang altitude azimuth [-P|-W|-L|-G|-E] input_value(s) \fR[ \fBoptions \fR]
10 > .br
11  
12 < As described below, the radiation can be defined with the pairs direct-normal and diffuse-horizontal irradiance
13 < (-W option), direct-horizontal and diffuse-horizontal irradiance (-G option), direct-normal and diffuse-horizontal
14 < illuminance (-L option) or global-horizontal irradiation alone (-E option). The direct-normal radiation
36 < is understood here as the radiant flux coming from the sun and an area of approximately 3 degrees around the sun
37 < (World Meteorological Organisation specifications for measuring the direct radiation.
38 < The aperture angle of a pyrheliometer is approximately 6 degrees).
39 < To simplify the calculations for the direct radiation, the sun is represented as a disk and no  
40 < circumsolar radiation is modelled in the 3 degrees around the sun. This means that
41 < all the measured/evaluated direct radiation is added to the 0.5 degree sun source.
12 > .SH "DESCRIPTION"
13 > \fIGendaylit\fR produces a RADIANCE scene description based on an angular distribution of the daylight sources (direct+diffuse) for the given atmospheric conditions (direct and diffuse component of the solar radiation), date and local standard time. The default output is the radiance of the sun (direct) and the sky (diffuse) integrated over the visible spectral range (380-780 nm). We have used the calculation of the sun's position and the ground brightness models which were programmed in \fIgensky.\fR
14 > .br
15  
16 < The direct and diffuse solar irradiances/illuminances are the inputs needed for the calculation.
17 < These quantities are the commonly accessible data from radiometric measurement centres or from the
45 < Test Reference Year. The use of such data is the recommended method for achieving the most accurate
46 < simulation results.
16 > The diffuse angular distribution is calculated using the Perez et al.  sky luminance distribution model (see Solar Energy Vol. 50, No. 3, pp. 235-245, 1993) which, quoting Perez, describes "the mean instantaneous sky luminance angular distribution patterns for all sky conditions from overcast to clear, through partly cloudy, skies". The correctness of the resulting sky radiance/luminance values in this simulation is ensured through the normalization of the modelled sky diffuse to the measured sky diffuse irradiances/illuminances.
17 > .br
18  
19 < The atmospheric conditions are modelled with the Perez et al. parametrization  
20 < (see Solar Energy Vol. 44, No 5, pp. 271-289, 1990), which is dependent on the values for
50 < the direct-normal and the diffuse-horizontal irradiances. The three parameters
51 < are epsilon, delta and the solar zenith angle. "Epsilon variations express the transition from
52 < a totally overcast sky (epsilon=1) to a low turbidity clear sky (epsilon>6); delta
53 < variations reflect the opacity/thickness of the clouds". Delta can vary from 0.05
54 < representing a dark sky to 0.5 for a very bright sky. Not every combination of
55 < epsilon, delta and solar zenith angle is possible. For a clear day, if
56 < epsilon and the solar zenith angle are known, then delta can be determined. For intermediate or overcast
57 < days, the sky can be dark or bright, giving a range of possible values for delta
58 < when epsilon and the solar zenith are fixed. The relation between epsilon and delta
59 < is represented in a figure on page 393 in Solar Energy Vol.42, No 5, 1989. Note that the
60 < epsilon parameter is a function of the solar zenith angle. It means that a clear day
61 < will not be defined by fixed values of epsilon and delta. Consequently the input
62 < parameters, epsilon, delta and the solar zenith angle, have to be determined on a graph.
63 < It might be easier to work with the measured direct and diffuse components (direct normal irradiance/illuminance
64 < and diffuse horizontal irradiance/illuminance) than with the epsilon and delta parameters.
19 > As described below, the radiation can be defined with the pairs direct-normal and diffuse-horizontal irradiance (-W option), direct-horizontal and diffuse-horizontal irradiance (-G option), direct-normal and diffuse-horizontal illuminance (-L option) or global-horizontal irradiation alone (-E option). The direct-normal radiation is understood here as the radiant flux coming from the sun and an area of approximately 3 degrees around the sun (World Meteorological Organisation specifications for measuring the direct radiation.  The aperture angle of a pyrheliometer is approximately 6 degrees).  To simplify the calculations for the direct radiation, the sun is represented as a disk and no circumsolar radiation is modelled in the 3 degrees around the sun. This means that all the measured/evaluated direct radiation is added to the 0.5 degree sun source.
20 > .br
21  
22 < The conversion of irradiance into illuminance for the direct and the diffuse
23 < components is determined by the luminous efficacy models of Perez et al. (see
68 < Solar Energy Vol. 44, No 5, pp. 271-289, 1990). To convert the luminance values
69 < into radiance integrated over the visible range of the spectrum,
70 < we devide the luminance by the white light efficacy factor of  
71 < 179 lm/W. This is consistent with the RADIANCE calculation because the luminance
72 < will be recalculated from the radiance integrated over the visible range by:
22 > The direct and diffuse solar irradiances/illuminances are the inputs needed for the calculation.  These quantities are the commonly accessible data from radiometric measurement centres or from the Test Reference Year. The use of such data is the recommended method for achieving the most accurate simulation results.
23 > .br
24  
25 < luminance = radiance_integrated_over_visible_range * 179   or
26 < luminance = (RED*.263 + GREEN*.655 + BLUE*.082) * 179    with the capability
76 < to model colour (where radiance_integrated_over_visible_range == (RED + GREEN + BLUE)/3).
25 > The atmospheric conditions are modelled with the Perez et al. parametrization (see Solar Energy Vol. 44, No 5, pp. 271-289, 1990), which is dependent on the values for the direct-normal and the diffuse-horizontal irradiances. The three parameters are epsilon, delta and the solar zenith angle. "Epsilon variations express the transition from a totally overcast sky (epsilon=1) to a low turbidity clear sky (epsilon>6); delta variations reflect the opacity/thickness of the clouds". Delta can vary from 0.05 representing a dark sky to 0.5 for a very bright sky. Not every combination of epsilon, delta and solar zenith angle is possible. For a clear day, if epsilon and the solar zenith angle are known, then delta can be determined. For intermediate or overcast days, the sky can be dark or bright, giving a range of possible values for delta when epsilon and the solar zenith are fixed. The relation between epsilon and delta is represented in a figure on page 393 in Solar Energy Vol.42, No 5, 1989. Note that the epsilon parameter is a function of the solar zenith angle. It means that a clear day will not be defined by fixed values of epsilon and delta. Consequently the input parameters, epsilon, delta and the solar zenith angle, have to be determined on a graph.  It might be easier to work with the measured direct and diffuse components (direct normal irradiance/illuminance and diffuse horizontal irradiance/illuminance) than with the epsilon and delta parameters.
26 > .br
27  
28 < From
29 < .I gensky, if the hour is preceded by a plus sign ('+'), then it is interpreted as local solar time instead of standard time.
30 < The second form gives the solar angles explicitly. The altitude is measured in degrees above the horizon, and the
31 < azimuth is measured in degrees west of South. The x axis points east, the y axis points north, and the z axis
32 < corresponds to the zenith. The actual material and surface(s) used for the sky is left up to the user.
33 < In addition to the specification of a sky distribution function,
34 < .I gendaylit
35 < suggests an ambient value in a comment at the beginning of the description to use with the
36 < .I \-av
37 < option of the RADIANCE rendering programs. (See rview(1) and rpict(1).) This value is the cosine-weighted
38 < radiance of the sky in W/sr/m^2.
39 < .PP
40 < .I Gendaylit
41 < can be used with the following input parameters. They offer three possibilities
42 < to run it: with the Perez parametrization, with irradiance values and with illuminance values.
43 < .TP 10n
44 < .BR \-P
45 < .I epsilon
46 < .I delta  
47 < (these are the Perez parameters)
48 < .TP
49 < .BR \-W
50 < .I direct-normal-irradiance
51 < (W/m^2),
52 < .I diffuse-horizontal-irradiance
53 < (W/m^2)
54 < .TP
55 < .BR \-G
56 < .I direct-horizontal-irradiance
57 < (W/m^2),
108 < .I diffuse-horizontal-irradiance
109 < (W/m^2)
110 < .TP
111 < .BR \-L
112 < .I direct-normal-illuminance
113 < (lm/m^2),
114 < .I diffuse-horizontal-illuminance
115 < (lm/m^2)
116 < .TP
117 < .BR \-E
118 < .I global-horizontal-irradiance
119 < (W/m^2)
120 < .PP
121 < The -E option calculates the diffuse irradiance fraction with the model of Erbs, Klein and Duffie (Solar Energy 28/4, 1982),
122 < being followed by the calculation of the -G option. Due to the high uncertainty of the model, the results have to be handled
123 < with care. A second irradiance value, if available, is definitely recommended.
124 < .PP
28 > The conversion of irradiance into illuminance for the direct and the diffuse components is determined by the luminous efficacy models of Perez et al. (see Solar Energy Vol. 44, No 5, pp. 271-289, 1990). To convert the luminance values into radiance integrated over the visible range of the spectrum, we devide the luminance by the white light efficacy factor of 179 lm/W. This is consistent with the RADIANCE calculation because the luminance will be recalculated from the radiance integrated over the visible range by:
29 > .br
30 >
31 > luminance = radiance_integrated_over_visible_range * 179   or luminance = (RED*.263 + GREEN*.655 + BLUE*.082) * 179    with the capability to model colour (where radiance_integrated_over_visible_range == (RED + GREEN + BLUE)/3).
32 > .br
33 >
34 > From \fIgensky,\fR \fIif\fR \fIthe\fR \fIhour\fR \fIis\fR \fIpreceded\fR \fIby\fR \fIa\fR \fIplus\fR \fIsign\fR \fI('+'),\fR \fIthen\fR \fIit\fR \fIis\fR \fIinterpreted\fR \fIas\fR \fIlocal\fR \fIsolar\fR \fItime\fR \fIinstead\fR \fIof\fR \fIstandard\fR \fItime.\fR  The second form gives the solar angles explicitly. The altitude is measured in degrees above the horizon, and the azimuth is measured in degrees west of South. The x axis points east, the y axis points north, and the z axis corresponds to the zenith. The actual material and surface(s) used for the sky is left up to the user.  In addition to the specification of a sky distribution function, \fIgendaylit\fR suggests an ambient value in a comment at the beginning of the description to use with the \fI-av\fR option of the RADIANCE rendering programs. (See rview(1) and rpict(1).) This value is the cosine-weighted radiance of the sky in W/sr/m^2.
35 > .br
36 >
37 > \fIGendaylit\fR can be used with the following input parameters. They offer three possibilities to run it: with the Perez parametrization, with irradiance values and with illuminance values.
38 > .br
39 >
40 > \fB-P        \fR\fIepsilon\fR \fIdelta\fR (these are the Perez parameters)
41 > .br
42 >
43 > \fB-W        \fR\fIdirect-normal-irradiance\fR (W/m^2), \fIdiffuse-horizontal-irradiance\fR (W/m^2)
44 > .br
45 >
46 > \fB-G        \fR\fIdirect-horizontal-irradiance\fR (W/m^2), \fIdiffuse-horizontal-irradiance\fR (W/m^2)
47 > .br
48 >
49 > \fB-L        \fR\fIdirect-normal-illuminance\fR (lm/m^2), \fIdiffuse-horizontal-illuminance\fR (lm/m^2)
50 > .br
51 >
52 > \fB-E        \fR\fIglobal-horizontal-irradiance\fR (W/m^2)
53 > .br
54 >
55 > The -E option calculates the diffuse irradiance fraction with the model of Erbs, Klein and Duffie (Solar Energy 28/4, 1982), being followed by the calculation of the -G option. Due to the high uncertainty of the model, the results have to be handled with care. A second irradiance value, if available, is definitely recommended.
56 > .br
57 >
58   The output can be set to either the radiance of the visible radiation, the solar radiance (full spectrum) or the luminance.
59 < .TP 10n
127 < .BR \-O [0|1|2]  
128 < (0=output in W/m^2/sr visible radiation (default), 1=output in W/m^2/sr solar radiation, 2=output in lm/m^2/sr luminance).
129 < .PP
130 < .I Gendaylit
131 < supports the following options.
132 < .TP 10n
133 < .BR \-s
134 < The source description of the sun is not generated.
135 < .TP
136 < .BR \-w
137 < Suppress warning messages
138 < .TP
139 < .BI -g \ rfl
140 < Average ground reflectance is
141 < .I rfl.
142 < This value is used to compute
143 < .I skyfunc
144 < when Dz is negative.
145 < .PP
146 < The following options do not apply when the solar
147 < altitude and azimuth are given explicitly.
148 < .TP
149 < .BI -a \ lat
150 < The site latitude is
151 < .I lat
152 < degrees north.
153 < (Use negative angle for south latitude.)
154 < This is used in the calculation of sun angle.
155 < .TP
156 < .BI -o \ lon
157 < The site longitude is
158 < .I lon
159 < degrees west.
160 < (Use negative angle for east longitude.)
161 < This is used in the calculation of solar time and sun angle.
162 < Be sure to give the corresponding standard meridian also!
163 < If solar time is given directly, then this option has no effect.
164 < .TP
165 < .BI -m \ mer
166 < The site standard meridian is
167 < .I mer
168 < degrees west of Greenwich.
169 < (Use negative angle for east.)
170 < This is used in the calculation of solar time.
171 < Be sure to give the correct longitude also!
172 < If solar time is given directly, then this option has no effect.
173 < .TP
174 < .BI -i \ time_interval [min]
175 < If gendaylit is used with weather files, the specified instantaneous points of time may be incorrect. This error occurs
176 < due to the fact that measurement results are frequently defined for time intervals, not for specific points of time.
177 < Although gendaylit is working correctly, this may lead to wrong outputs especially at low sun altitudes.
178 < The -i option allows to specify the time interval of the measurements in minutes, causing the solar position to be corrected for low sun
179 < altitudes. A warning message is returned if a correction has been performed.
59 > .br
60  
61 < .SH EXAMPLES
61 > \fB-O\fR[0|1|2] (0=output in W/m^2/sr visible radiation (default), 1=output in W/m^2/sr solar radiation, 2=output in lm/m^2/sr luminance).
62 > .br
63 >
64 > \fIGendaylit\fR can generate a colored sky description based on the model from A. Diakite, TU-Berlin. It will be invoked by the \fB-C \fRoption.
65 > .br
66 >
67 > \fB-C\fR                        generate colored sky description
68 > .br
69 >
70 > At the moment the color information is generated only for overcast skies (epsilon < 1.065) or clear skies (epsilon >4.5). For intermediate skies there is no color output generated for now until the model has been developed also for these sky types. For the daylight locus following options have been implemented:
71 > .br
72 >
73 >                        daylight locus determined by the measured skies of Berlin (\fBdefault\fR)
74 > .br
75 > \fB-CIE\fR                      CIE daylight locus
76 > .br
77 > \fB-l\fR \fIl1\fR \fIl2\fR \fIl3\fR \fIl4\fR \fIl5\fR \fIl6\fR \fIl7\fR user defined daylight locus, according to following equations:
78 > .br
79 >
80 >                        (xd=l1/T^3+l2/T^2+l3/T+l4 ; yd=l5*xd^2+l6*xd+l7)
81 >
82 > .br
83 > \fIGendaylit\fR supports the following options.
84 > .br
85 >
86 > \fB-s        \fRThe source description of the sun is not generated.
87 > .br
88 >
89 > \fB-w        \fRSuppress warning messages
90 > .br
91 >
92 > \fB-g \fR\fIrfl\fR    Average ground reflectance is \fIrfl.\fR  This value is used to compute \fIskyfunc\fR when Dz is negative.
93 > .br
94 >
95 > The following options do not apply when the solar altitude and azimuth are given explicitly.
96 > .br
97 >
98 > \fB-a \fR\fIlat\fR The site latitude is \fIlat\fR degrees north.  (Use negative angle for south latitude.)  This is used in the calculation of sun angle.
99 > .br
100 >
101 > \fB-o \fR\fIlon\fR The site longitude is \fIlon\fR degrees west.  (Use negative angle for east longitude.)  This is used in the calculation of solar time and sun angle.  Be sure to give the corresponding standard meridian also!  If solar time is given directly, then this option has no effect.
102 > .br
103 >
104 > \fB-m \fR\fImer\fR The site standard meridian is \fImer\fR degrees west of Greenwich.  (Use negative angle for east.)  This is used in the calculation of solar time.  Be sure to give the correct longitude also!  If solar time is given directly, then this option has no effect.
105 > .br
106 >
107 > \fB-i \fR\fItime_interval\fR\fB[min]\fR
108 > .br
109 >       If gendaylit is used with weather files, the specified instantaneous points of time may be incorrect. This error occurs due to the fact that measurement results are frequently defined for time intervals, not for specific points of time.  Although gendaylit is working correctly, this may lead to wrong outputs especially at low sun altitudes.  The -i option allows to specify the time interval of the measurements in minutes, causing the solar position to be corrected for low sun altitudes. A warning message is returned if a correction has been performed.
110 > .br
111 >
112 > .SH "EXAMPLES"
113   A clear non-turbid sky for a solar altitude of 60 degrees and an azimut of 0 degree might be defined by:
114 < .IP "" .2i
115 < gendaylit -ang 60 0 -P 6.3 0.12 or gendaylit -ang 60 0 -W 840 135
116 < This sky description corresponds to the clear sky standard of the CIE.
117 < .PP
114 > .br
115 >
116 >  gendaylit -ang 60 0 -P 6.3 0.12 or gendaylit -ang 60 0 -W 840 135 This sky description corresponds to the clear sky standard of the CIE.
117 > .br
118 >
119   The corresponding sky with a high turbidity is:
120 < .IP "" .2i
121 < gendaylit -ang 60 0 -P 3.2 0.24 or gendaylit -ang 60 0 -W 720 280
122 < .PP
123 < The dark overcast sky (corresponding to the CIE overcast standard, see CIE draft standard,
124 < Pub. No. CIE DS 003, 1st Edition, 1994) is obtained by:
125 < .IP "" .2i
126 < gendaylit -ang 60 0 -P 1 0.08
127 < .PP
120 > .br
121 >
122 >  gendaylit -ang 60 0 -P 3.2 0.24 or gendaylit -ang 60 0 -W 720 280
123 > .br
124 >
125 > The dark overcast sky (corresponding to the CIE overcast standard, see CIE draft standard, Pub. No. CIE DS 003, 1st Edition, 1994) is obtained by:
126 > .br
127 >
128 >  gendaylit -ang 60 0 -P 1 0.08
129 > .br
130 >
131   A bright overcast sky is modelled with a larger value of delta, for example:
132 < .IP "" .2i
133 < gendaylit -ang 60 0 -P 1 0.35
134 < .PP
135 < To generate the same bright overcast sky for March 2th at 3:15pm standard time at a site
136 < latitude of 42 degrees, 108 degrees west longitude, and a 110 degrees standard meridian:
137 < .IP "" .2i
138 < gendaylit 3 2 15.25 -a 42 -o 108 -m 110 -P 1 0.35
139 < .PP
140 < .SH FILES
132 > .br
133 >
134 >  gendaylit -ang 60 0 -P 1 0.35
135 > .br
136 >
137 > To generate the same bright overcast sky for March 2th at 3:15pm standard time at a site latitude of 42 degrees, 108 degrees west longitude, and a 110 degrees standard meridian:
138 > .br
139 >
140 >  gendaylit 3 2 15.25 -a 42 -o 108 -m 110 -P 1 0.35
141 > .br
142 >
143 > .SH "FILES"
144   /usr/local/lib/ray/perezlum.cal
145 < .SH AUTHOR
146 < Jean-Jacques Delaunay, Jan Wienold, Wendelin Sprenger, Fraunhofer ISE (Freiburg i.B., Germany) ([email protected])
147 < .SH ACKNOWLEDGEMENTS
148 < The first work on this program was supported by the German Federal Ministry for Research
149 < and Technology (BMFT) under contract No. 0329294A, and a scholarship from
150 < the French Environment and Energy Agency (ADEME) which was co-funded by Bouygues.
151 < Many thanks to Peter Apian-Bennewitz, Arndt Berger, Christian Reetz, Ann Kovach, R. Perez, C. Gueymard and G. Ward for their help.
145 > .br
146 > /usr/local/lib/ray/perezlum_c.cal
147 > .br
148 >
149 > .SH "AUTHOR"
150 > Jean-Jacques Delaunay, Jan Wienold, Wendelin Sprenger, Fraunhofer ISE (Freiburg i.B., Germany)
151 > .br
152 > Jan Wienold, EPFL ([email protected])
153 > .br
154 >
155 > .SH "ACKNOWLEDGEMENTS"
156 > The first work on this program was supported by the German Federal Ministry for Research and Technology (BMFT) under contract No. 0329294A, and a scholarship from the French Environment and Energy Agency (ADEME) which was co-funded by Bouygues.  Many thanks to A. Diakite, Peter Apian-Bennewitz, Arndt Berger, Christian Reetz, Ann Kovach, R. Perez, C. Gueymard and G. Ward for their help.
157 > .br
158 >
159   .SH "SEE ALSO"
160   gensky(1), rpict(1), rview(1), xform(1)
161 + .br

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines