ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/Development/ray/src/gen/genssky.c
(Generate patch)

Comparing ray/src/gen/genssky.c (file contents):
Revision 2.3 by greg, Fri Aug 2 18:47:25 2024 UTC vs.
Revision 2.4 by greg, Thu Aug 8 02:00:48 2024 UTC

# Line 179 | Line 179 | static double get_overcast_brightness(const double dz,
179  
180   static void write_header(const int argc, char **argv, const double cloud_cover,
181                           const double grefl, const int res) {
182 +  int i;
183    printf("# ");
184 <  for (int i = 0; i < argc; i++) {
184 >  for (i = 0; i < argc; i++) {
185      printf("%s ", argv[i]);
186    }
187    printf("\n");
# Line 195 | Line 196 | static void write_rad(const double *sun_radiance, cons
196      printf("void spectrum sunrad\n0\n0\n22 380 780 ");
197      /* Normalize to one */
198      double sum = 0.0;
199 <    for (int i = 0; i < NSSAMP; ++i) {
199 >    int i;
200 >    for (i = 0; i < NSSAMP; ++i) {
201        sum += sun_radiance[i];
202      }
203      double mean = sum / NSSAMP;
204 <    for (int i = 0; i < NSSAMP; ++i) {
204 >    for (i = 0; i < NSSAMP; ++i) {
205        printf("%.3f ", sun_radiance[i] / mean);
206      }
207      double intensity = mean * WVLSPAN;
# Line 248 | Line 250 | int gen_spect_sky(DATARRAY *tau_clear, DATARRAY *scat_
250    FVECT view_point = {0, 0, ER + 10};
251    const double radius = VLEN(view_point);
252    const double sun_ct = fdot(view_point, sundir) / radius;
253 <  for (int j = 0; j < yres; ++j) {
254 <    for (int i = 0; i < xres; ++i) {
253 >  int i, j, k;
254 >  for (j = 0; j < yres; ++j) {
255 >    for (i = 0; i < xres; ++i) {
256        SCOLOR radiance = {0};
257        SCOLR sky_sclr = {0};
258  
# Line 274 | Line 277 | int gen_spect_sky(DATARRAY *tau_clear, DATARRAY *scat_
277                           radiance);
278        }
279  
280 <      for (int k = 0; k < NSSAMP; ++k) {
280 >      for (k = 0; k < NSSAMP; ++k) {
281          radiance[k] *= WVLSPAN;
282        }
283  
# Line 283 | Line 286 | int gen_spect_sky(DATARRAY *tau_clear, DATARRAY *scat_
286          double grndbr = zenithbr * GNORM;
287          double skybr = get_overcast_brightness(rdir[2], zenithbr);
288          if (rdir[2] < 0) {
289 <          for (int k = 0; k < NSSAMP; ++k) {
289 >          for (k = 0; k < NSSAMP; ++k) {
290              radiance[k] = wmean2(radiance[k], grndbr * D6415[k], cloud_cover);
291            }
292          } else {
293 <          for (int k = 0; k < NSSAMP; ++k) {
293 >          for (k = 0; k < NSSAMP; ++k) {
294              radiance[k] = wmean2(radiance[k], skybr * D6415[k], cloud_cover);
295            }
296          }
# Line 306 | Line 309 | int gen_spect_sky(DATARRAY *tau_clear, DATARRAY *scat_
309    if (cloud_cover > 0) {
310      double zenithbr = get_zenith_brightness(sundir);
311      double skybr = get_overcast_brightness(sundir[2], zenithbr);
312 <    for (int i = 0; i < NSSAMP; ++i) {
312 >    int i;
313 >    for (i = 0; i < NSSAMP; ++i) {
314        sun_radiance[i] =
315            wmean2(sun_radiance[i], D6415[i] * skybr / WVLSPAN, cloud_cover);
316      }
# Line 414 | Line 418 | int main(int argc, char *argv[]) {
418    char mie_name[20] = "mie_ca";
419    char lstag[3];
420    char *ddir = ".";
421 +  int i;
422  
423    if (argc == 2 && !strcmp(argv[1], "-defaults")) {
424      printf("-i %d\t\t\t\t#scattering order\n", sorder);
# Line 451 | Line 456 | int main(int argc, char *argv[]) {
456      exit(1);
457    }
458  
459 <  for (int i = 4; i < argc; i++) {
459 >  for (i = 4; i < argc; i++) {
460      if (argv[i][0] == '-') {
461        switch (argv[i][1]) {
462        case 'a':

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines