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: Sat Feb 22 16:57:30 2014 +# dlm: Thu May 29 09:23:46 2014 # (c) 2003 A.M. Thurnherr -# uE-Info: 463 0 NIL 0 0 72 10 2 4 NIL ofnI +# uE-Info: 302 0 NIL 0 0 72 10 2 4 NIL ofnI #====================================================================== # miscellaneous RDI-specific utilities @@ -49,6 +49,7 @@ # - added set_range_lim() # Feb 22, 2014: - changed gap heuristic # - Earth coord beam-pair warning removed +# May 29, 2014: - removed unused code (disabled warning) from ref_lr_w use strict; @@ -261,8 +262,6 @@ # calculate reference-layer vertical and incident velocities -{ my($warned); # static scope - sub ref_lr_w($$$$$$$) { my($dta,$ens,$rl_b0,$rl_b1,$min_corr,$max_e,$min_pctg) = @_; @@ -296,15 +295,11 @@ $dta->{ENSEMBLE}[$ens]->{PERCENT_GOOD}[$i][2] > 0 || $dta->{ENSEMBLE}[$ens]->{PERCENT_GOOD}[$i][3] < $min_pctg); @v = @{$dta->{ENSEMBLE}[$ens]->{VELOCITY}[$i]}; - unless ($warned) { -### print(STDERR "WARNING: incident-flow & beam-pair velocities not yet implemented for earth coordinates"); - $warned = 1; - } } ### next if (!defined($v[3]) || abs($v[3]) > $max_e); # disallow 3-beam solutions next if (defined($v[3]) && abs($v[3]) > $max_e); # allow 3-beam solutions - if (defined($v[2])) { # valid velocity + if (defined($v[2])) { # valid vertical velocity $vel[2] += $v[2]; $n[2]++; # vertical velocity $vel[3] += $v[3], $n[3]++ if defined($v[3]); # error velocity push(@w,$v[2]); # save for stderr calculation @@ -407,8 +402,6 @@ } -} # static scope - sub mk_prof(...) # make profile {