| 1 |
|
.\" RCSid "$Id$" |
| 2 |
+ |
.TH GENDAYLIT 1 4/12/94 "RADIANCE ISE/ADEME EXTENSIONS" |
| 3 |
|
.SH NAME |
| 4 |
< |
gendaylit - generates a RADIANCE description of the daylit sources using Perez models for diffuse and direct components |
| 4 |
> |
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] direct_value diffuse_value " |
| 6 |
> |
.B "gendaylit month day hour [-P|-W|-L|-G|-E] input_value(s) " |
| 7 |
|
[ |
| 8 |
|
.B options |
| 9 |
|
] |
| 10 |
|
.br |
| 11 |
< |
.B "gendaylit -ang altitude azimuth [-P|-W|-L] direct_value diffuse_value " |
| 11 |
> |
.B "gendaylit -ang altitude azimuth [-P|-W|-L|-G|-E] input_value(s) " |
| 12 |
|
[ |
| 13 |
|
.B options |
| 14 |
|
] |
| 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 (diffus) |
| 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 |
| 30 |
|
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. |
| 32 |
|
|
| 33 |
< |
The direct radiation is understood here as the radiant flux coming from the sun |
| 34 |
< |
and an area of approximately 3 degrees around the sun (World Meteorological Organisation specifications |
| 35 |
< |
for measuring the direct radiation. The aperture angle of a pyrheliometer is approximately 6 degrees). |
| 33 |
> |
As described below, the radiation can be defined with the pairs direct-normal and diffuse-horizontal irradiance |
| 34 |
> |
(-W option), direct-horizontal and diffuse-horizontal irradiance (-G option), direct-normal and diffuse-horizontal |
| 35 |
> |
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. |
| 42 |
|
|
| 43 |
< |
.I The direct and diffuse |
| 44 |
< |
.I solar irradiances/illuminances |
| 45 |
< |
.I are the inputs needed |
| 46 |
< |
.I for the calculation. |
| 43 |
< |
These quantities are the commonly accessible data from radiometric measurement centres, conversion models |
| 44 |
< |
(e.g. global irradiance to direct irradiance), or from the Test Reference Year. The use of such |
| 45 |
< |
data is the recommended method for achieving the most accurate simulation results. |
| 43 |
> |
The direct and diffuse solar irradiances/illuminances are the inputs needed for the calculation. |
| 44 |
> |
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. |
| 47 |
|
|
| 47 |
– |
|
| 48 |
|
The atmospheric conditions are modelled with the Perez et al. parametrization |
| 49 |
|
(see Solar Energy Vol. 44, No 5, pp. 271-289, 1990), which is dependent on the values for |
| 50 |
< |
the direct and the diffuse irradiances. The three parameters |
| 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 |
| 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, |
| 60 |
< |
or can be obtained from the author of this RADIANCE extension upon request. Note that the |
| 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. |
| 65 |
|
|
| 67 |
– |
|
| 66 |
|
The conversion of irradiance into illuminance for the direct and the diffuse |
| 67 |
|
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 : |
| 72 |
> |
will be recalculated from the radiance integrated over the visible range by: |
| 73 |
|
|
| 74 |
|
luminance = radiance_integrated_over_visible_range * 179 or |
| 77 |
– |
|
| 75 |
|
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). |
| 77 |
|
|
| 78 |
|
From |
| 79 |
< |
.I gensky |
| 80 |
< |
, if the hour is preceded by a plus sign ('+'), then it is interpreted |
| 81 |
< |
as local solar time instead of standard time. |
| 82 |
< |
The second form gives the solar angles explicitly. |
| 83 |
< |
The altitude is measured in degrees above the horizon, and the |
| 87 |
< |
azimuth is measured in degrees west of South. |
| 88 |
< |
.PP |
| 89 |
< |
The x axis points east, |
| 90 |
< |
the y axis points north, and the z axis corresponds to the zenith. |
| 91 |
< |
The actual material and surface(s) used for the sky is left |
| 92 |
< |
up to the user. |
| 93 |
< |
.PP |
| 94 |
< |
In addition to the specification of |
| 95 |
< |
a sky distribution function, |
| 79 |
> |
.I gensky, if the hour is preceded by a plus sign ('+'), then it is interpreted as local solar time instead of standard time. |
| 80 |
> |
The second form gives the solar angles explicitly. The altitude is measured in degrees above the horizon, and the |
| 81 |
> |
azimuth is measured in degrees west of South. The x axis points east, the y axis points north, and the z axis |
| 82 |
> |
corresponds to the zenith. The actual material and surface(s) used for the sky is left up to the user. |
| 83 |
> |
In addition to the specification of a sky distribution function, |
| 84 |
|
.I gendaylit |
| 85 |
< |
suggests an ambient value in a comment at the beginning of the |
| 98 |
< |
description to use with the |
| 85 |
> |
suggests an ambient value in a comment at the beginning of the description to use with the |
| 86 |
|
.I \-av |
| 87 |
< |
option of the RADIANCE rendering programs. |
| 88 |
< |
(See rview(1) and rpict(1).) |
| 102 |
< |
This value is the cosine-weighted radiance of the sky in |
| 103 |
< |
W/sr/m^2. |
| 87 |
> |
option of the RADIANCE rendering programs. (See rview(1) and rpict(1).) This value is the cosine-weighted |
| 88 |
> |
radiance of the sky in W/sr/m^2. |
| 89 |
|
.PP |
| 90 |
|
.I Gendaylit |
| 91 |
|
can be used with the following input parameters. They offer three possibilities |
| 92 |
< |
to run it: with the Perez parametrization, with the irradiance values |
| 108 |
< |
and with the illuminance values. |
| 92 |
> |
to run it: with the Perez parametrization, with irradiance values and with illuminance values. |
| 93 |
|
.TP 10n |
| 94 |
|
.BR \-P |
| 95 |
|
.I epsilon |
| 102 |
|
.I diffuse-horizontal-irradiance |
| 103 |
|
(W/m^2) |
| 104 |
|
.TP |
| 105 |
+ |
.BR \-G |
| 106 |
+ |
.I direct-horizontal-irradiance |
| 107 |
+ |
(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 output can be set to either the radiance of the visible radiation (default), the solar radiance (full spectrum) or the luminance. |
| 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 |
| 125 |
> |
The output can be set to either the radiance of the visible radiation, the solar radiance (full spectrum) or the luminance. |
| 126 |
|
.TP 10n |
| 127 |
|
.BR \-O [0|1|2] |
| 128 |
< |
(0=output in W/m^2/sr visible radiation, 0=output in W/m^2/sr solar radiation, 2=output in lm/m^2/sr luminance) |
| 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. |
| 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. |
| 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. |
| 180 |
+ |
|
| 181 |
|
.SH EXAMPLES |
| 182 |
|
A clear non-turbid sky for a solar altitude of 60 degrees and an azimut of 0 degree might be defined by: |
| 183 |
|
.IP "" .2i |
| 186 |
|
.PP |
| 187 |
|
The corresponding sky with a high turbidity is: |
| 188 |
|
.IP "" .2i |
| 189 |
< |
gendaylit -ang 60 0 -P 3.2 0.24 or gendaylit -ang 60 0 -W 720 280 |
| 189 |
> |
gendaylit -ang 60 0 -P 3.2 0.24 or gendaylit -ang 60 0 -W 720 280 |
| 190 |
|
.PP |
| 191 |
|
The dark overcast sky (corresponding to the CIE overcast standard, see CIE draft standard, |
| 192 |
|
Pub. No. CIE DS 003, 1st Edition, 1994) is obtained by: |
| 193 |
|
.IP "" .2i |
| 194 |
< |
gendaylit -ang 60 0 -P 1 0.08 |
| 194 |
> |
gendaylit -ang 60 0 -P 1 0.08 |
| 195 |
|
.PP |
| 196 |
|
A bright overcast sky is modelled with a larger value of delta, for example: |
| 197 |
|
.IP "" .2i |
| 205 |
|
.SH FILES |
| 206 |
|
/usr/local/lib/ray/perezlum.cal |
| 207 |
|
.SH AUTHOR |
| 208 |
< |
Jean-Jacques Delaunay, FhG-ISE Freiburg, ([email protected]) |
| 208 |
> |
Jean-Jacques Delaunay, Jan Wienold, Wendelin Sprenger, Fraunhofer ISE (Freiburg i.B., Germany) ([email protected]) |
| 209 |
|
.SH ACKNOWLEDGEMENTS |
| 210 |
< |
The work on this program was supported by the German Federal Ministry for Research |
| 210 |
> |
The first work on this program was supported by the German Federal Ministry for Research |
| 211 |
|
and Technology (BMFT) under contract No. 0329294A, and a scholarship from |
| 212 |
|
the French Environment and Energy Agency (ADEME) which was co-funded by Bouygues. |
| 213 |
< |
Many thanks to Peter Apian-Bennewitz, Arndt Berger, Ann Kovach, R. Perez, C. Gueymard and G. Ward for their help. |
| 213 |
> |
Many thanks to Peter Apian-Bennewitz, Arndt Berger, Christian Reetz, Ann Kovach, R. Perez, C. Gueymard and G. Ward for their help. |
| 214 |
|
.SH "SEE ALSO" |
| 215 |
|
gensky(1), rpict(1), rview(1), xform(1) |
