9 |
|
my $windoz = ($^O eq "MSWin32" or $^O eq "MSWin64"); |
10 |
|
use File::Temp qw/ :mktemp /; |
11 |
|
sub userror { |
12 |
< |
print STDERR "Usage: iso2klems [-t][-W][-s \"x=string;y=string\"] [input.dat]\n"; |
12 |
> |
print STDERR "Usage: iso2klems [-t][-f \"x=string;y=string\"][-u unit] [input.dat]\n"; |
13 |
|
exit 1; |
14 |
|
} |
15 |
|
my ($td,$rmtmp,$cmd); |
29 |
|
while ($#ARGV >= 0) { |
30 |
|
if ("$ARGV[0]" eq "-t") { |
31 |
|
$reverse = ! $reverse; |
32 |
< |
} elsif ("$ARGV[0]" eq "-s") { |
32 |
> |
} elsif ("$ARGV[0]" =~ /^-[fs]$/) { |
33 |
|
$wrapper .= " -f \"$ARGV[1]\""; |
34 |
|
shift @ARGV; |
35 |
+ |
} elsif ("$ARGV[0]" eq "-u") { |
36 |
+ |
$wrapper .= " -u $ARGV[1]"; |
37 |
+ |
shift @ARGV; |
38 |
|
} elsif ("$ARGV[0]" =~ /^-./) { |
39 |
|
userror(); |
40 |
|
} else { |
91 |
|
qq{-e "Rd1=if($h2h[1]-.001,Rdiff(tideg)*Rdiff(tsdeg)/$h2h[1],0)" } . |
92 |
|
qq{-e "Td2=if($h2h[2]-.001,Tdiff(180-tideg)*Tdiff(180-tsdeg)/$h2h[2],0)" } . |
93 |
|
qq{-e "Rd2=if($h2h[3]-.001,Rdiff(180-tideg)*Rdiff(180-tsdeg)/$h2h[3],0)" } . |
94 |
< |
qq{-e "\$1=if(Tspec(0),if(diag,Tspec(tideg)/om,Td1*corr),Td1)" } . |
95 |
< |
qq{-e "\$2=if(Rspec(0),if(diag,Rspec(tideg)/om,Rd1*corr),Rd1)" } . |
96 |
< |
qq{-e "\$3=if(Tspec(180),if(diag,Tspec(180-tideg)/om,Td2*corr),Td2)" } . |
97 |
< |
qq{-e "\$4=if(Rspec(180),if(diag,Rspec(180-tideg)/om,Rd2*corr),Rd2)" }; |
94 |
> |
q{-e "$1=if(Tspec(0),if(diag,Tspec(tideg)/om,Td1*corr),Td1)" } . |
95 |
> |
q{-e "$2=if(Rspec(0),if(diag,Rspec(tideg)/om,Rd1*corr),Rd1)" } . |
96 |
> |
q{-e "$3=if(Tspec(180),if(diag,Tspec(180-tideg)/om,Td2*corr),Td2)" } . |
97 |
> |
q{-e "$4=if(Rspec(180),if(diag,Rspec(180-tideg)/om,Rd2*corr),Rd2)" }; |
98 |
|
} else { |
99 |
|
$cmd = qq{cnt 145 145 } . |
100 |
|
qq{| rcalc -f $funcfile } . |
103 |
|
qq{-e 'Rd1=if($h2h[1]-.001,Rdiff(tideg)*Rdiff(tsdeg)/$h2h[1],0)' } . |
104 |
|
qq{-e 'Td2=if($h2h[2]-.001,Tdiff(180-tideg)*Tdiff(180-tsdeg)/$h2h[2],0)' } . |
105 |
|
qq{-e 'Rd2=if($h2h[3]-.001,Rdiff(180-tideg)*Rdiff(180-tsdeg)/$h2h[3],0)' } . |
106 |
< |
qq{-e '\$1=if(Tspec(0),if(diag,Tspec(tideg)/om,Td1*corr),Td1)' } . |
107 |
< |
qq{-e '\$2=if(Rspec(0),if(diag,Rspec(tideg)/om,Rd1*corr),Rd1)' } . |
108 |
< |
qq{-e '\$3=if(Tspec(180),if(diag,Tspec(180-tideg)/om,Td2*corr),Td2)' } . |
109 |
< |
qq{-e '\$4=if(Rspec(180),if(diag,Rspec(180-tideg)/om,Rd2*corr),Rd2)' }; |
106 |
> |
q{-e '$1=if(Tspec(0),if(diag,Tspec(tideg)/om,Td1*corr),Td1)' } . |
107 |
> |
q{-e '$2=if(Rspec(0),if(diag,Rspec(tideg)/om,Rd1*corr),Rd1)' } . |
108 |
> |
q{-e '$3=if(Tspec(180),if(diag,Tspec(180-tideg)/om,Td2*corr),Td2)' } . |
109 |
> |
q{-e '$4=if(Rspec(180),if(diag,Rspec(180-tideg)/om,Rd2*corr),Rd2)' }; |
110 |
|
} |
111 |
|
system qq{$cmd | rsplit "-t " $TK1dataf $RK1dataf $TK2dataf $RK2dataf}; |
112 |
|
die "Error running rcalc or rsplit" if ( $? ); |