equal
deleted
inserted
replaced
1 #!/usr/bin/perl |
1 #!/usr/bin/perl |
2 #====================================================================== |
2 #====================================================================== |
3 # L I S T W |
3 # L I S T W |
4 # doc: Wed Mar 24 06:45:09 2004 |
4 # doc: Wed Mar 24 06:45:09 2004 |
5 # dlm: Thu Jul 30 17:42:33 2009 |
5 # dlm: Mon Nov 25 18:30:33 2013 |
6 # (c) 2004 A.M. Thurnherr |
6 # (c) 2004 A.M. Thurnherr |
7 # uE-Info: 205 53 NIL 0 0 72 2 2 4 NIL ofnI |
7 # uE-Info: 23 0 NIL 0 0 72 2 2 4 NIL ofnI |
8 #====================================================================== |
8 #====================================================================== |
9 |
9 |
10 # dump vertical velocities |
10 # dump vertical velocities |
11 |
11 |
12 # NB: currently broken |
12 # NB: currently broken |
17 # - floatized time |
17 # - floatized time |
18 # Apr 3, 2004: - cosmetics |
18 # Apr 3, 2004: - cosmetics |
19 # Nov 8, 2005: - UNIXTIME => UNIX_TIME |
19 # Nov 8, 2005: - UNIXTIME => UNIX_TIME |
20 # Sep 19, 2007: - adapted to new [RDI_BB_Read.pl] (not tested) |
20 # Sep 19, 2007: - adapted to new [RDI_BB_Read.pl] (not tested) |
21 # Jul 30, 2009: - NaN => nan |
21 # Jul 30, 2009: - NaN => nan |
|
22 # Nov 25, 2013: - checkEnsemble() expunged |
22 |
23 |
23 $0 =~ m{(.*)/[^/]+}; |
24 $0 =~ m{(.*)/[^/]+}; |
24 require "$1/WorkhorseBinRead.pl"; |
25 require "$1/WorkhorseBinRead.pl"; |
25 require "$1/WorkhorseCoords.pl"; |
26 require "$1/WorkhorseCoords.pl"; |
26 require "$1/WorkhorseUtils.pl"; |
27 require "$1/WorkhorseUtils.pl"; |
121 #---------------------------------------------------------------------- |
122 #---------------------------------------------------------------------- |
122 |
123 |
123 print(STDERR "Generating profile by integrating w..."); |
124 print(STDERR "Generating profile by integrating w..."); |
124 |
125 |
125 for ($e=0; $e<=$#{$dta{ENSEMBLE}}; $e++) { |
126 for ($e=0; $e<=$#{$dta{ENSEMBLE}}; $e++) { |
126 checkEnsemble(\%dta,$e); # sanity checks |
|
127 filterEnsemble(\%dta,$e) # filter ensemble |
127 filterEnsemble(\%dta,$e) # filter ensemble |
128 if (defined($opt_F) && |
128 if (defined($opt_F) && |
129 $dta{ENSEMBLE}[$e]->{PERCENT_GOOD}[0][0] > 0); |
129 $dta{ENSEMBLE}[$e]->{PERCENT_GOOD}[0][0] > 0); |
130 |
130 |
131 $dta{ENSEMBLE}[$e]->{REFW} = ref_lr_w($e); |
131 $dta{ENSEMBLE}[$e]->{REFW} = ref_lr_w($e); |