#======================================================================
#                    D U M P _ R E S I D U A L _ P R O F I L E S . P L 
#                    doc: Thu Mar 24 07:55:07 2016
#                    dlm: Tue May 24 16:35:30 2016
#                    (c) 2016 A.M. Thurnherr
#                    uE-Info: 46 24 NIL 0 0 72 2 2 4 NIL ofnI
#======================================================================

# HISTORY:
#	Mar 24, 2016: - created from [plot_residuals.pl]
#	Mar 29, 2016: - cleaned up
#   May 24, 2016: - calc_binDepths() -> binDepths()

sub dump_residual_profiles($)
{
	my($odir) = @_;

	unless (-d $odir) {
		warning(0,"Creating residual-profile output directory ./$odir\n");
		my(@dirs) = split('/',$odir);
		my($path) = '.';
		foreach my $d (@dirs) {
			$path .= "/$d";
			mkdir($path);
		}
    }
	
	return unless ($P{max_depth});

	@antsNewLayout = ('bin','depth','residual');

	for ($ens=$firstGoodEns; $ens<=$lastGoodEns; $ens++) {
		next unless numberp($LADCP{ENSEMBLE}[$ens]->{CTD_DEPTH});

		my($of) = sprintf('>%s/%04d.rprof',$odir,$LADCP{ENSEMBLE}[$ens]->{NUMBER});
		open(STDOUT,$of) || error("$of: $!\n");
		undef($antsActiveHeader) unless ($ANTS_TOOLS_AVAILABLE);

		&antsAddParams('ensemble',	$LADCP{ENSEMBLE}[$ens]->{NUMBER},
					   'elapsed',	$LADCP{ENSEMBLE}[$ens]->{ELAPSED},
					   'CTD_depth',	$LADCP{ENSEMBLE}[$ens]->{CTD_DEPTH},
					   'CTD_w',		$CTD{W}[$LADCP{ENSEMBLE}[$ens]->{CTD_SCAN}],
					   'CTD_accel',	$CTD{W_t}[$LADCP{ENSEMBLE}[$ens]->{CTD_SCAN}],
					   'ADCP_tilt',	$LADCP{ENSEMBLE}[$ens]->{TILT});

	  	my(@bindepth) = binDepths($ens);
		for ($bin=$LADCP_firstBin-1; $bin<=$LADCP_lastBin-1; $bin++) {
			next unless ($bin+1>=$outGrid_firstBin && $bin+1<=$outGrid_lastBin);
		  	next unless numberp($LADCP{ENSEMBLE}[$ens]->{W}[$bin]);
		  	my($bi) = $bindepth[$bin]/$opt_o;
			my($res) = ($ens < $LADCP_atbottom) ? 
						$LADCP{ENSEMBLE}[$ens]->{SSCORRECTED_OCEAN_W}[$bin] - $DNCAST{MEDIAN_W}[$bi] :
						$LADCP{ENSEMBLE}[$ens]->{SSCORRECTED_OCEAN_W}[$bin] - $UPCAST{MEDIAN_W}[$bi];
			&antsOut($bin,$bindepth[$bin],$res);
	    			 
        }
	    &antsOut('EOF'); open(STDOUT,'>&2');
    }
}

1; # return true on require