| 123 |
|
# Output MGF description if requested |
| 124 |
|
if ( $geout ) { |
| 125 |
|
print "\t<Geometry format=\"MGF\">\n"; |
| 126 |
< |
print "<MGFblock unit=\"$gunit\">\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</MGFblock>\n"; |
| 132 |
> |
print "</MGFblock>\n"; |
| 133 |
|
print "\t</Geometry>\n"; |
| 134 |
|
} |
| 135 |
|
# Set up surface sampling |
| 277 |
|
q{-e '$1=(0.265*$1+0.670*$2+0.065*$3)/Omega' }; |
| 278 |
|
if ($pctcull >= 0) { |
| 279 |
|
$cmd .= "-of $td/" . ($bmodnm,$fmodnm)[$forw] . ".flt " . |
| 280 |
< |
"| rttree_reduce -a -h -ff -t $pctcull -r $tensortree -g $ttlog2"; |
| 280 |
> |
"| rttree_reduce -h -ff -t $pctcull -r $tensortree -g $ttlog2"; |
| 281 |
|
system "$cmd" || die "Failure running rttree_reduce"; |
| 282 |
|
} else { |
| 283 |
|
$cmd .= "$td/" . ($bmodnm,$fmodnm)[$forw] . ".flt"; |
