| 41 |
|
char *kfullfn = "klems_full.cal"; |
| 42 |
|
char *khalffn = "klems_half.cal"; |
| 43 |
|
char *kquarterfn = "klems_quarter.cal"; |
| 44 |
+ |
char *ciefn = "cieskyscan.cal"; |
| 45 |
|
|
| 46 |
|
/* string indicating parameters */ |
| 47 |
|
const char PARAMSTART[] = "@rfluxmtx"; |
| 500 |
|
calfn = kquarterfn; kquarterfn = NULL; |
| 501 |
|
binf = "kqbin"; |
| 502 |
|
nbins = "Nkqbins"; |
| 503 |
+ |
} else if (!strcasecmp(curparams.hemis, "cie")) { |
| 504 |
+ |
calfn = ciefn; ciefn = NULL; |
| 505 |
+ |
sprintf(sbuf, "rNx=%g,rNy=%g,rNz=%g,Ux=%g,Uy=%g,Uz=%g,RHS=%c1", |
| 506 |
+ |
curparams.nrm[0], curparams.nrm[1], curparams.nrm[2], |
| 507 |
+ |
curparams.vup[0], curparams.vup[1], curparams.vup[2], |
| 508 |
+ |
curparams.sign); |
| 509 |
+ |
binv = "cbin"; |
| 510 |
+ |
nbins = "Ncbins"; |
| 511 |
|
} else { |
| 512 |
|
fprintf(stderr, "%s: unrecognized hemisphere sampling: h=%s\n", |
| 513 |
|
progname, curparams.hemis); |
| 1292 |
|
yrs = argv[++a]; |
| 1293 |
|
na = 0; |
| 1294 |
|
continue; |
| 1295 |
< |
case 'c': /* number of samples */ |
| 1296 |
< |
sampcnt = atoi(argv[++a]); |
| 1297 |
< |
if (sampcnt <= 0) |
| 1298 |
< |
goto userr; |
| 1299 |
< |
na = 0; /* we re-add this later */ |
| 1300 |
< |
continue; |
| 1295 |
> |
case 'c': /* spectral sampling or count */ |
| 1296 |
> |
switch (argv[a][2]) { |
| 1297 |
> |
case 's': |
| 1298 |
> |
na = 2; |
| 1299 |
> |
break; |
| 1300 |
> |
case 'w': |
| 1301 |
> |
na = 3; |
| 1302 |
> |
break; |
| 1303 |
> |
case '\0': |
| 1304 |
> |
sampcnt = atoi(argv[++a]); |
| 1305 |
> |
if (sampcnt <= 0) |
| 1306 |
> |
goto userr; |
| 1307 |
> |
na = 0; /* we re-add this later */ |
| 1308 |
> |
continue; |
| 1309 |
> |
} |
| 1310 |
> |
break; |
| 1311 |
|
case 'I': /* only for pass-through mode */ |
| 1312 |
|
case 'i': |
| 1313 |
|
iropt = argv[a]; |
