|
54
|
1 |
#======================================================================
|
|
|
2 |
# W R I T E _ R E S I D U A L _ P R O F S . P L
|
|
|
3 |
# doc: Fri May 15 20:22:54 2020
|
|
|
4 |
# dlm: Fri May 15 20:56:06 2020
|
|
|
5 |
# (c) 2020 A.M. Thurnherr
|
|
|
6 |
# uE-Info: 21 53 NIL 0 0 72 2 2 4 NIL ofnI
|
|
|
7 |
#======================================================================
|
|
|
8 |
|
|
|
9 |
# HISTORY:
|
|
|
10 |
# May 25, 2020: - created from [plot_residual_profs.pl]
|
|
|
11 |
|
|
|
12 |
sub write_residual_profs($)
|
|
|
13 |
{
|
|
|
14 |
my($ofn) = @_;
|
|
|
15 |
|
|
|
16 |
@antsNewLayout = ('depth','dc_residual12','dc_residual34','uc_residual12','uc_residual34');
|
|
|
17 |
|
|
|
18 |
open(STDOUT,">$ofn") || error("$ofn: $!\n");
|
|
|
19 |
undef($antsActiveHeader) unless ($ANTS_TOOLS_AVAILABLE);
|
|
|
20 |
|
|
|
21 |
for (my($bi)=0; $bi<=$#{$DNCAST{MEAN_RESIDUAL12}}; $bi++) {
|
|
|
22 |
my(@out);
|
|
|
23 |
push(@out,($bi+0.5) * $opt_o); # depth
|
|
|
24 |
push(@out,(numberp($DNCAST{MEAN_RESIDUAL12}[$bi]) && $DNCAST{N_SAMP}[$bi]>=$minsamp) ?
|
|
|
25 |
$DNCAST{MEAN_RESIDUAL12}[$bi] : nan);
|
|
|
26 |
push(@out,(numberp($DNCAST{MEAN_RESIDUAL34}[$bi]) && $DNCAST{N_SAMP}[$bi]>=$minsamp) ?
|
|
|
27 |
$DNCAST{MEAN_RESIDUAL34}[$bi] : nan);
|
|
|
28 |
push(@out,(numberp($UPCAST{MEAN_RESIDUAL12}[$bi]) && $UPCAST{N_SAMP}[$bi]>=$minsamp) ?
|
|
|
29 |
$UPCAST{MEAN_RESIDUAL12}[$bi] : nan);
|
|
|
30 |
push(@out,(numberp($UPCAST{MEAN_RESIDUAL34}[$bi]) && $UPCAST{N_SAMP}[$bi]>=$minsamp) ?
|
|
|
31 |
$UPCAST{MEAN_RESIDUAL34}[$bi] : nan);
|
|
|
32 |
&antsOut(@out);
|
|
|
33 |
}
|
|
|
34 |
&antsOut('EOF'); open(STDOUT,'>&2');
|
|
|
35 |
}
|
|
|
36 |
|
|
|
37 |
1; # return true on require
|