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", |
1132 |
|
perez_param[0], perez_param[1], perez_param[2], perez_param[3], perez_param[4]); |
1133 |
|
|
1134 |
|
for (i = 1; i < nskypatch; i++) |
1145 |
|
|
1146 |
|
/* Calculate patch luminance */ |
1147 |
|
parr[3*i] = CalcRelLuminance(sspa, zsa); |
1150 |
– |
fprintf(stderr, "CalcRelLuminance(%f, %f) = %f\n", sspa, zsa, parr[3*i]); |
1148 |
|
if (parr[3*i] < 0) parr[3*i] = 0; |
1149 |
|
parr[3*i+2] = parr[3*i+1] = parr[3*i]; |
1150 |
|
} |