#======================================================================
# W R I T E _ R E S I D U A L _ P R O F S . P L
# doc: Fri May 15 20:22:54 2020
# dlm: Fri May 15 20:56:06 2020
# (c) 2020 A.M. Thurnherr
# uE-Info: 21 53 NIL 0 0 72 2 2 4 NIL ofnI
#======================================================================
# HISTORY:
# May 25, 2020: - created from [plot_residual_profs.pl]
sub write_residual_profs($)
{
my($ofn) = @_;
@antsNewLayout = ('depth','dc_residual12','dc_residual34','uc_residual12','uc_residual34');
open(STDOUT,">$ofn") || error("$ofn: $!\n");
undef($antsActiveHeader) unless ($ANTS_TOOLS_AVAILABLE);
for (my($bi)=0; $bi<=$#{$DNCAST{MEAN_RESIDUAL12}}; $bi++) {
my(@out);
push(@out,($bi+0.5) * $opt_o); # depth
push(@out,(numberp($DNCAST{MEAN_RESIDUAL12}[$bi]) && $DNCAST{N_SAMP}[$bi]>=$minsamp) ?
$DNCAST{MEAN_RESIDUAL12}[$bi] : nan);
push(@out,(numberp($DNCAST{MEAN_RESIDUAL34}[$bi]) && $DNCAST{N_SAMP}[$bi]>=$minsamp) ?
$DNCAST{MEAN_RESIDUAL34}[$bi] : nan);
push(@out,(numberp($UPCAST{MEAN_RESIDUAL12}[$bi]) && $UPCAST{N_SAMP}[$bi]>=$minsamp) ?
$UPCAST{MEAN_RESIDUAL12}[$bi] : nan);
push(@out,(numberp($UPCAST{MEAN_RESIDUAL34}[$bi]) && $UPCAST{N_SAMP}[$bi]>=$minsamp) ?
$UPCAST{MEAN_RESIDUAL34}[$bi] : nan);
&antsOut(@out);
}
&antsOut('EOF'); open(STDOUT,'>&2');
}
1; # return true on require