diff --git a/RDI_Utils.pl b/RDI_Utils.pl --- a/RDI_Utils.pl +++ b/RDI_Utils.pl @@ -1,9 +1,9 @@ #====================================================================== # R D I _ U T I L S . P L # doc: Wed Feb 12 10:21:32 2003 -# dlm: Tue Aug 8 16:25:05 2017 +# dlm: Mon Nov 27 10:32:50 2017 # (c) 2003 A.M. Thurnherr -# uE-Info: 58 55 NIL 0 0 72 2 2 4 NIL ofnI +# uE-Info: 59 75 NIL 0 0 72 2 2 4 NIL ofnI #====================================================================== # miscellaneous RDI-specific utilities @@ -56,6 +56,7 @@ # May 19, 2016: - adapted to new velBeamToInstrument() usage # Aug 7, 2017: - added abmiguity velocity # Aug 8, 2017: - changed transducer frequency to kHz +# Nov 27, 2017: - BUG: profile-restart heuristic did not work with P6#001 use strict; @@ -467,10 +468,11 @@ $dta->{ENSEMBLE}[$lastgood]->{UNIX_TIME}; # ... last good ens if ($dt > $max_gap) { -# if ($dta->{ENSEMBLE}[$lastgood]->{UNIX_TIME} - # heuristic changed Feb 22, 2014 -# $dta->{ENSEMBLE}[$firstgood]->{UNIX_TIME} > 15*60) { - if (@{$dta->{ENSEMBLE}}-$e < @{$dta->{ENSEMBLE}}/2) { + # 2nd heuristic test added Nov 2017 for P06 profile #001 + if ((@{$dta->{ENSEMBLE}}-$e < @{$dta->{ENSEMBLE}}/2) && + ($maxz > 25 && $z < $maxz/2)) { printf(STDERR "WARNING: %.1f-s gap in 2nd half of profile is too long; profile ended at ensemble $lastgood\n",$dt); +# printf(STDERR "\t[#ens = %d, end-of-gap = $e]\n",scalar(@{$dta->{ENSEMBLE}})); last; } printf(STDERR "WARNING: %.1f-s gap in first half of profile is too long; profile restarted at ensemble $e\n",$dt);