86 |
|
#include <string.h> |
87 |
|
#include <ctype.h> |
88 |
|
#include "rtmath.h" |
89 |
+ |
#include "rtio.h" |
90 |
|
#include "resolu.h" |
91 |
|
#include "platform.h" |
92 |
|
#include "color.h" |
556 |
|
break; |
557 |
|
case 'f': |
558 |
|
for (j = 0; j < ntsteps; j++) { |
559 |
< |
fwrite(mtx_data+mtx_offset, sizeof(float), 3, |
559 |
> |
putbinary(mtx_data+mtx_offset, sizeof(float), 3, |
560 |
|
stdout); |
561 |
|
mtx_offset += 3*nskypatch; |
562 |
|
} |
567 |
|
ment[0] = mtx_data[mtx_offset]; |
568 |
|
ment[1] = mtx_data[mtx_offset+1]; |
569 |
|
ment[2] = mtx_data[mtx_offset+2]; |
570 |
< |
fwrite(ment, sizeof(double), 3, stdout); |
570 |
> |
putbinary(ment, sizeof(double), 3, stdout); |
571 |
|
mtx_offset += 3*nskypatch; |
572 |
|
} |
573 |
|
break; |
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 |
+ |
|
1136 |
|
for (i = 1; i < nskypatch; i++) |
1137 |
|
{ |
1138 |
|
/* Calculate sun-sky point azimuthal angle */ |
1147 |
|
|
1148 |
|
/* Calculate patch luminance */ |
1149 |
|
parr[3*i] = CalcRelLuminance(sspa, zsa); |
1150 |
+ |
fprintf(stderr, "CalcRelLuminance(%f, %f) = %f\n", sspa, zsa, parr[3*i]); |
1151 |
|
if (parr[3*i] < 0) parr[3*i] = 0; |
1152 |
|
parr[3*i+2] = parr[3*i+1] = parr[3*i]; |
1153 |
|
} |