| 1129 |
|
double sspa; /* Sun-sky point angle */ |
| 1130 |
|
double zsa; /* Zenithal sun angle */ |
| 1131 |
|
|
| 1132 |
– |
fprintf(stderr, "Alt, azi: %f %f\n", altitude, azimuth); |
| 1133 |
– |
fprintf(stderr, "Perez parameters: %f %f %f %f %f\n", |
| 1134 |
– |
perez_param[0], perez_param[1], perez_param[2], perez_param[3], perez_param[4]); |
| 1135 |
– |
|
| 1132 |
|
for (i = 1; i < nskypatch; i++) |
| 1133 |
|
{ |
| 1134 |
|
/* Calculate sun-sky point azimuthal angle */ |
| 1143 |
|
|
| 1144 |
|
/* Calculate patch luminance */ |
| 1145 |
|
parr[3*i] = CalcRelLuminance(sspa, zsa); |
| 1150 |
– |
fprintf(stderr, "CalcRelLuminance(%f, %f) = %f\n", sspa, zsa, parr[3*i]); |
| 1146 |
|
if (parr[3*i] < 0) parr[3*i] = 0; |
| 1147 |
|
parr[3*i+2] = parr[3*i+1] = parr[3*i]; |
| 1148 |
|
} |
