| 307 |
|
double rotation = 0; /* site rotation (degrees) */ |
| 308 |
|
double elevation; /* site elevation (meters) */ |
| 309 |
|
int dir_is_horiz; /* direct is meas. on horizontal? */ |
| 310 |
+ |
FILE *sunsfp = NULL; /* output file for individual suns */ |
| 311 |
|
float *mtx_data = NULL; /* our matrix data */ |
| 312 |
|
int avgSky = 0; /* compute average sky r.t. matrix? */ |
| 313 |
|
int ntsteps = 0; /* number of time steps */ |
| 368 |
|
skycolor[1] = atof(argv[++i]); |
| 369 |
|
skycolor[2] = atof(argv[++i]); |
| 370 |
|
break; |
| 371 |
+ |
case 'n': /* no matrix output */ |
| 372 |
+ |
avgSky = -1; |
| 373 |
+ |
rhsubdiv = 1; |
| 374 |
+ |
/* fall through */ |
| 375 |
|
case 'd': /* solar (direct) only */ |
| 376 |
|
skycolor[0] = skycolor[1] = skycolor[2] = 0; |
| 377 |
|
grefl[0] = grefl[1] = grefl[2] = 0; |
| 378 |
|
break; |
| 379 |
+ |
case 'D': /* output suns to file */ |
| 380 |
+ |
sunsfp = fopen(argv[++i], "w"); |
| 381 |
+ |
if (!sunsfp) { |
| 382 |
+ |
fprintf(stderr, "%s: cannot open '%s' for output\n", |
| 383 |
+ |
argv[0], argv[i]); |
| 384 |
+ |
exit(1); |
| 385 |
+ |
} |
| 386 |
+ |
fixed_sun_sa = PI/360.*0.533; |
| 387 |
+ |
fixed_sun_sa *= PI*fixed_sun_sa; |
| 388 |
+ |
break; |
| 389 |
|
case 's': /* sky only (no direct) */ |
| 390 |
|
suncolor[0] = suncolor[1] = suncolor[2] = 0; |
| 391 |
|
break; |
| 488 |
|
memset(mtx_data+mtx_offset, 0, sizeof(float)*3*nskypatch); |
| 489 |
|
continue; |
| 490 |
|
} |
| 476 |
– |
if (verbose && mo != last_monthly) |
| 477 |
– |
fprintf(stderr, "%s: stepping through month %d...\n", |
| 478 |
– |
progname, last_monthly=mo); |
| 491 |
|
/* compute solar position */ |
| 492 |
|
julian_date = jdate(mo, da); |
| 493 |
|
sda = sdec(julian_date); |
| 506 |
|
} |
| 507 |
|
/* compute sky patch values */ |
| 508 |
|
ComputeSky(mtx_data+mtx_offset); |
| 509 |
+ |
/* output sun if indicated */ |
| 510 |
+ |
if (sunsfp && (altitude > 0) & (dir_illum > 1e-4)) { |
| 511 |
+ |
double srad = dir_illum/(WHTEFFICACY * fixed_sun_sa); |
| 512 |
+ |
FVECT sv; |
| 513 |
+ |
vector(sv, altitude, azimuth); |
| 514 |
+ |
fprintf(sunsfp, "\nvoid light solar%d\n0\n0\n", ntsteps); |
| 515 |
+ |
fprintf(sunsfp, "3 %.3e %.3e %.3e\n", srad*suncolor[0], |
| 516 |
+ |
srad*suncolor[1], srad*suncolor[2]); |
| 517 |
+ |
fprintf(sunsfp, "\nsolar%d source sun%d\n0\n0\n", ntsteps, ntsteps); |
| 518 |
+ |
fprintf(sunsfp, "4 %.6f %.6f %.6f 0.533\n", sv[0], sv[1], sv[2]); |
| 519 |
+ |
} |
| 520 |
+ |
if (avgSky < 0) /* no matrix? */ |
| 521 |
+ |
continue; |
| 522 |
+ |
|
| 523 |
|
AddDirect(mtx_data+mtx_offset); |
| 524 |
|
/* update cumulative sky? */ |
| 525 |
|
for (i = 3*nskypatch*(avgSky&(ntsteps>1)); i--; ) |
| 526 |
|
mtx_data[i] += mtx_data[mtx_offset+i]; |
| 527 |
+ |
/* monthly reporting */ |
| 528 |
+ |
if (verbose && mo != last_monthly) |
| 529 |
+ |
fprintf(stderr, "%s: stepping through month %d...\n", |
| 530 |
+ |
progname, last_monthly=mo); |
| 531 |
|
} |
| 532 |
+ |
if (!ntsteps) { |
| 533 |
+ |
fprintf(stderr, "%s: no valid time steps on input\n", progname); |
| 534 |
+ |
exit(1); |
| 535 |
+ |
} |
| 536 |
|
/* check for junk at end */ |
| 537 |
|
while ((i = fgetc(stdin)) != EOF) |
| 538 |
|
if (!isspace(i)) { |
| 543 |
|
fputs(buf, stderr); fputc('\n', stderr); |
| 544 |
|
break; |
| 545 |
|
} |
| 546 |
< |
if (!ntsteps) { |
| 547 |
< |
fprintf(stderr, "%s: no valid time steps on input\n", progname); |
| 548 |
< |
exit(1); |
| 549 |
< |
} |
| 546 |
> |
|
| 547 |
> |
if (avgSky < 0) /* no matrix output? */ |
| 548 |
> |
goto alldone; |
| 549 |
> |
|
| 550 |
|
dif = 1./(double)ntsteps; /* average sky? */ |
| 551 |
|
for (i = 3*nskypatch*(avgSky&(ntsteps>1)); i--; ) |
| 552 |
|
mtx_data[i] *= dif; |
| 607 |
|
if (ferror(stdout)) |
| 608 |
|
goto writerr; |
| 609 |
|
} |
| 610 |
< |
if (fflush(stdout) == EOF) |
| 610 |
> |
alldone: |
| 611 |
> |
if (fflush(NULL) == EOF) |
| 612 |
|
goto writerr; |
| 613 |
|
if (verbose) |
| 614 |
|
fprintf(stderr, "%s: done.\n", progname); |
| 615 |
|
exit(0); |
| 616 |
|
userr: |
| 617 |
< |
fprintf(stderr, "Usage: %s [-v][-h][-A][-d|-s][-r deg][-m N][-g r g b][-c r g b][-o{f|d}][-O{0|1}] [tape.wea]\n", |
| 617 |
> |
fprintf(stderr, "Usage: %s [-v][-h][-A][-d|-s|-n][-D file][-r deg][-m N][-g r g b][-c r g b][-o{f|d}][-O{0|1}] [tape.wea]\n", |
| 618 |
|
progname); |
| 619 |
|
exit(1); |
| 620 |
|
fmterr: |
| 621 |
< |
fprintf(stderr, "%s: input weather tape format error\n", progname); |
| 621 |
> |
fprintf(stderr, "%s: weather tape format error in header\n", progname); |
| 622 |
|
exit(1); |
| 623 |
|
writerr: |
| 624 |
|
fprintf(stderr, "%s: write error on output\n", progname); |
