Changeset 18502
- Timestamp:
- Jul 13, 2008, 3:36:28 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branch_20080706/ippScripts/scripts/detrend_reject_exp.pl
r18501 r18502 153 153 # neb:///*/gpc1/20080130 (volume not specified) 154 154 155 my $logFile;156 155 unless ($no_op) { 157 $logFile = $ipprc->file_create_append( $logName ); 158 print $logFile "Ensemble mean $ensMeanMean +/- $ensMeanStdev\n"; 159 print $logFile "Ensemble stdev $ensStdevMean +/- $ensStdevStdev\n"; 160 print $logFile "Ensemble mean rms (over imfiles) $ensMeanStdevMean +/- $ensMeanStdevStdev\n\n"; 156 print "Ensemble mean $ensMeanMean +/- $ensMeanStdev\n"; 157 print "Ensemble stdev $ensStdevMean +/- $ensStdevStdev\n"; 158 print "Ensemble mean rms (over imfiles) $ensMeanStdevMean +/- $ensMeanStdevStdev\n\n"; 161 159 } 162 160 … … 189 187 # Rejected this at an earlier stage 190 188 unless ($no_op) { 191 print $logFile"Rejecting $expID based on earlier determination.\n";189 print "Rejecting $expID based on earlier determination.\n"; 192 190 } 193 191 $reject = 1; … … 208 206 my $delta = abs($mean - $ensMeanMean); 209 207 if ($delta > ($reject_mean * $ensMeanStdev)) { 210 print $logFile"Rejecting $expID based on ensemble mean value: ";208 print "Rejecting $expID based on ensemble mean value: "; 211 209 $reject = 1; 212 210 #goto UPDATE; 213 211 } else { 214 print $logFile"$expID OK against ensemble mean: ";215 } 216 print $logFile"$mean --> $delta vs " . $reject_mean * $ensMeanStdev . "\n";212 print "$expID OK against ensemble mean: "; 213 } 214 print "$mean --> $delta vs " . $reject_mean * $ensMeanStdev . "\n"; 217 215 } else { 218 print $logFile"No rejection of $expID for ensemble mean\n";216 print "No rejection of $expID for ensemble mean\n"; 219 217 } 220 218 … … 222 220 my $delta = abs($stdev - $ensStdevMean); 223 221 if ($delta > ($reject_stdev * $ensStdevStdev)) { 224 print $logFile"Rejecting $expID based on ensemble stdev: ";222 print "Rejecting $expID based on ensemble stdev: "; 225 223 $reject = 1; 226 224 #goto UPDATE; 227 225 } else { 228 print $logFile"$expID OK against ensemble stdev: ";229 } 230 print $logFile"$stdev --> $delta sigma vs " . $reject_stdev * $ensStdevStdev . "\n";226 print "$expID OK against ensemble stdev: "; 227 } 228 print "$stdev --> $delta sigma vs " . $reject_stdev * $ensStdevStdev . "\n"; 231 229 } else { 232 print $logFile"No rejection of $expID for ensemble stdev\n";230 print "No rejection of $expID for ensemble stdev\n"; 233 231 } 234 232 … … 236 234 my $delta = abs($meanStdev - $ensMeanStdevMean); 237 235 if ($delta > ($reject_meanstdev * $ensMeanStdevStdev)) { 238 print $logFile"Rejecting $expID based on ensemble mean stdev: ";236 print "Rejecting $expID based on ensemble mean stdev: "; 239 237 $reject = 1; 240 238 #goto UPDATE; 241 239 } else { 242 print $logFile"$expID OK against ensemble mean stdev: ";243 } 244 print $logFile"$meanStdev --> $delta sigma vs " . $reject_meanstdev * $ensMeanStdevStdev. "\n";240 print "$expID OK against ensemble mean stdev: "; 241 } 242 print "$meanStdev --> $delta sigma vs " . $reject_meanstdev * $ensMeanStdevStdev. "\n"; 245 243 } else { 246 print $logFile"No rejection of $expID for ensemble mean stdev\n";244 print "No rejection of $expID for ensemble mean stdev\n"; 247 245 } 248 246 … … 256 254 if ((not $include and not $reject) or ($include and $reject)) { 257 255 unless ($no_op) { 258 print $logFile"Status of $expID has changed.\n";256 print "Status of $expID has changed.\n"; 259 257 } 260 258 $numChanges++; … … 289 287 290 288 unless ($no_op) { 291 print $logFile "Master: $master\n"; 292 print $logFile "Stop: $stop\n"; 293 close $logFile; 289 print "Master: $master\n"; 290 print "Stop: $stop\n"; 294 291 } 295 292 … … 305 302 exit($PS_EXIT_PROG_ERROR); 306 303 } 307 308 # Allow iteration to be turned off309 my $allow_iter = metadataLookupBool($ipprc->{rejection}, "ITERATION"); # Allow iteration?310 my $force_master = metadataLookupBool($ipprc->{rejection}, "MASTER"); # Force the stack to be accepted311 $stop = 1 unless $allow_iter;312 $master = 1 if $force_master;313 304 314 305 ## add the summary statistics, and request a new iteration if needed
Note:
See TracChangeset
for help on using the changeset viewer.
