133 |
|
print "\t</Geometry>\n"; |
134 |
|
} |
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); |
136 |
> |
my $nx = int(sqrt($nsamp*($dim[1]-$dim[0])/($dim[3]-$dim[2])) + 1); |
137 |
> |
my $ny = int($nsamp/$nx + 1); |
138 |
|
$nsamp = $nx * $ny; |
139 |
|
my $ns = 2**$ttlog2; |
140 |
|
my (@pdiv, $disk2sq, $sq2disk, $tcal, $kcal); |
278 |
|
if ($pctcull >= 0) { |
279 |
|
$cmd .= "-of $td/" . ($bmodnm,$fmodnm)[$forw] . ".flt " . |
280 |
|
"| rttree_reduce -h -ff -t $pctcull -r $tensortree -g $ttlog2"; |
281 |
+ |
$cmd .= " -a" if ($tensortree == 3); |
282 |
|
system "$cmd" || die "Failure running rttree_reduce"; |
283 |
|
} else { |
284 |
|
$cmd .= "$td/" . ($bmodnm,$fmodnm)[$forw] . ".flt"; |