| 24 |
|
my $doforw = 0; |
| 25 |
|
my $doback = 1; |
| 26 |
|
my $pctcull = 90; |
| 27 |
< |
my $gunit = "Meter"; |
| 27 |
> |
my $gunit = "meter"; |
| 28 |
|
my @dim; |
| 29 |
|
# Get options |
| 30 |
|
while ($#ARGV >= 0) { |
| 108 |
|
print "<!-- File produced by: genBSDF @savedARGV -->\n"; |
| 109 |
|
print |
| 110 |
|
'<WindowElementType>System</WindowElementType> |
| 111 |
+ |
<FileType>BSDF</FileType> |
| 112 |
|
<Optical> |
| 113 |
|
<Layer> |
| 114 |
|
<Material> |
| 118 |
|
printf "\t\t<Thickness unit=\"$gunit\">%.3f</Thickness>\n", $dim[5] - $dim[4]; |
| 119 |
|
printf "\t\t<Width unit=\"$gunit\">%.3f</Width>\n", $dim[1] - $dim[0]; |
| 120 |
|
printf "\t\t<Height unit=\"$gunit\">%.3f</Height>\n", $dim[3] - $dim[2]; |
| 121 |
< |
print "\t\t<DeviceType>Integral</DeviceType>\n"; |
| 121 |
> |
print "\t\t<DeviceType>Other</DeviceType>\n"; |
| 122 |
> |
print " </Material>\n"; |
| 123 |
|
# Output MGF description if requested |
| 124 |
|
if ( $geout ) { |
| 125 |
< |
print "\t\t<Geometry format=\"MGF\" unit=\"$gunit\">\n"; |
| 125 |
> |
print "\t<Geometry format=\"MGF\">\n"; |
| 126 |
> |
print "\t\t<MGFblock unit=\"$gunit\">\n"; |
| 127 |
|
printf "xf -t %.6f %.6f 0\n", -($dim[0]+$dim[1])/2, -($dim[2]+$dim[3])/2; |
| 128 |
|
open(MGFSCN, "< $mgfscn"); |
| 129 |
|
while (<MGFSCN>) { print $_; } |
| 130 |
|
close MGFSCN; |
| 131 |
|
print "xf\n"; |
| 132 |
< |
print "\t\t</Geometry>\n"; |
| 132 |
> |
print "</MGFblock>\n"; |
| 133 |
> |
print "\t</Geometry>\n"; |
| 134 |
|
} |
| 131 |
– |
print " </Material>\n"; |
| 135 |
|
# Set up surface sampling |
| 136 |
|
my $nx = int(sqrt($nsamp*($dim[1]-$dim[0])/($dim[3]-$dim[2])) + .5); |
| 137 |
|
my $ny = int($nsamp/$nx + .5); |
| 306 |
|
print "\t\t\t<WavelengthDataDirection>Reflection $side</WavelengthDataDirection>\n"; |
| 307 |
|
print |
| 308 |
|
' <AngleBasis>LBNL/Shirley-Chiu</AngleBasis> |
| 309 |
< |
<ScatteringDataType>BRDF</ScatteringDataType> |
| 309 |
> |
<ScatteringDataType>BTDF</ScatteringDataType> |
| 310 |
|
<ScatteringData> |
| 311 |
|
'; |
| 312 |
|
$cmd = "rcalc -if3 -e 'Omega:PI/($ns*$ns)' " . |
| 407 |
|
system "$cmd" || die "Failure running: $cmd\n"; |
| 408 |
|
@tfarr = `$rccmd $td/$fmodnm.flt`; |
| 409 |
|
die "Failure running: $rccmd $td/$fmodnm.flt\n" if ( $? ); |
| 410 |
+ |
chomp(@tfarr); |
| 411 |
|
@rfarr = `$rccmd $td/$bmodnm.flt`; |
| 412 |
|
die "Failure running: $rccmd $td/$bmodnm.flt\n" if ( $? ); |
| 413 |
+ |
chomp(@rfarr); |
| 414 |
|
} |
| 415 |
|
if ( $doback ) { |
| 416 |
|
$cmd = "cnt $ndiv $ny $nx | rcalc -of -e '$tcal' " . |
| 543 |
|
<WavelengthDataDirection>Reflection Front</WavelengthDataDirection> |
| 544 |
|
<ColumnAngleBasis>LBNL/Klems Full</ColumnAngleBasis> |
| 545 |
|
<RowAngleBasis>LBNL/Klems Full</RowAngleBasis> |
| 546 |
< |
<ScatteringDataType>BRDF</ScatteringDataType> |
| 546 |
> |
<ScatteringDataType>BTDF</ScatteringDataType> |
| 547 |
|
<ScatteringData> |
| 548 |
|
'; |
| 549 |
|
# Output front reflection (transposed order) |
| 593 |
|
<WavelengthDataDirection>Reflection Back</WavelengthDataDirection> |
| 594 |
|
<ColumnAngleBasis>LBNL/Klems Full</ColumnAngleBasis> |
| 595 |
|
<RowAngleBasis>LBNL/Klems Full</RowAngleBasis> |
| 596 |
< |
<ScatteringDataType>BRDF</ScatteringDataType> |
| 596 |
> |
<ScatteringDataType>BTDF</ScatteringDataType> |
| 597 |
|
<ScatteringData> |
| 598 |
|
'; |
| 599 |
|
# Output back reflection (transposed order) |
