32 |
|
|
33 |
|
int rev_orient = 0; /* shall we reverse surface orientation? */ |
34 |
|
|
35 |
+ |
double lim_graze = 0; /* limit scattering near grazing (deg above) */ |
36 |
+ |
|
37 |
|
/* Compare incident angles */ |
38 |
|
static int |
39 |
|
cmp_indir(const void *p1, const void *p2) |
180 |
|
fclose(fp); |
181 |
|
return(1); |
182 |
|
} |
183 |
< |
if (rev_orient) { /* reverse Z-axis to face outside */ |
183 |
> |
/* check if scatter angle is too low */ |
184 |
> |
if (fabs(theta_out - 90.) < lim_graze-FTINY) |
185 |
> |
continue; |
186 |
> |
if (rev_orient) { /* reverse Z-axis to face outside? */ |
187 |
|
theta_out = 180. - theta_out; |
188 |
|
phi_out = 360. - phi_out; |
189 |
|
} |
229 |
|
main(int argc, char *argv[]) |
230 |
|
{ |
231 |
|
extern int nprocs; |
232 |
+ |
int auto_grazing = 0; |
233 |
|
const char *symmetry = "0"; |
234 |
|
int ninpfiles, totinc; |
235 |
|
int a, i; |
246 |
|
case 's': |
247 |
|
symmetry = argv[++a]; |
248 |
|
break; |
249 |
+ |
case 'g': |
250 |
+ |
if (toupper(argv[a+1][0]) == 'A') |
251 |
+ |
auto_grazing = 1; |
252 |
+ |
else |
253 |
+ |
lim_graze = atof(argv[a+1]); |
254 |
+ |
++a; |
255 |
+ |
break; |
256 |
|
default: |
257 |
|
goto userr; |
258 |
|
} |
264 |
|
inpfile = (PGINPUT *)malloc(sizeof(PGINPUT)*totinc); |
265 |
|
if (inpfile == NULL) |
266 |
|
return(1); |
267 |
< |
for (i = 0; i < ninpfiles; i++) |
267 |
> |
if (auto_grazing) |
268 |
> |
lim_graze = 90.; |
269 |
> |
for (i = 0; i < ninpfiles; i++) { |
270 |
|
if (!init_pabopto_inp(i, argv[a+i])) |
271 |
|
return(1); |
272 |
< |
|
272 |
> |
if (auto_grazing && fabs(inpfile[i].theta - 90.) < lim_graze) |
273 |
> |
lim_graze = fabs(inpfile[i].theta - 90.); |
274 |
> |
} |
275 |
|
for (i = ninpfiles; i < totinc; i++) { /* copy for "up" symmetry */ |
276 |
|
inpfile[i] = inpfile[i-ninpfiles]; |
277 |
|
inpfile[i].phi += 180.; /* invert duplicate data */ |
325 |
|
save_bsdf_rep(stdout); /* complete header + data */ |
326 |
|
return(0); |
327 |
|
userr: |
328 |
< |
fprintf(stderr, "Usage: %s [-t][-n nproc][-s symmetry] meas1.dat meas2.dat .. > bsdf.sir\n", |
328 |
> |
fprintf(stderr, "Usage: %s [-t][-n nproc][-s symmetry][-g angle|'A'] meas1.dat meas2.dat .. > bsdf.sir\n", |
329 |
|
progname); |
330 |
|
return(1); |
331 |
|
} |