1 |
|
---|
2 |
|
---|
3 |
|
---|
4 |
|
---|
5 |
|
---|
6 |
|
---|
7 |
|
---|
8 |
|
---|
9 |
|
---|
10 |
|
---|
11 |
|
---|
12 |
|
---|
13 |
|
---|
14 |
|
---|
15 |
|
---|
16 |
|
---|
17 |
|
---|
18 |
|
---|
19 |
|
---|
20 |
|
---|
21 |
|
---|
22 |
|
---|
23 |
|
---|
24 |
use lib '/opt/local/seacoos/bin/'; |
---|
25 |
use lib '/opt/local/seacoos/bin/meta_data/'; |
---|
26 |
use meta_data; |
---|
27 |
$bin_dir = '/opt/local/seacoos/bin/meta_data/'; |
---|
28 |
|
---|
29 |
|
---|
30 |
|
---|
31 |
|
---|
32 |
|
---|
33 |
|
---|
34 |
|
---|
35 |
|
---|
36 |
|
---|
37 |
|
---|
38 |
|
---|
39 |
|
---|
40 |
|
---|
41 |
|
---|
42 |
|
---|
43 |
|
---|
44 |
use NetCDF; |
---|
45 |
use Date::Manip; |
---|
46 |
use Statistics::Basic::StdDev; |
---|
47 |
use Math::BigFloat; |
---|
48 |
|
---|
49 |
|
---|
50 |
|
---|
51 |
$starttime = time; |
---|
52 |
|
---|
53 |
|
---|
54 |
|
---|
55 |
|
---|
56 |
&handle_input; |
---|
57 |
|
---|
58 |
|
---|
59 |
|
---|
60 |
|
---|
61 |
|
---|
62 |
|
---|
63 |
|
---|
64 |
|
---|
65 |
|
---|
66 |
|
---|
67 |
if ($debug) {print "\nReading data files:\n";} |
---|
68 |
$first_one = 1; |
---|
69 |
&read_data; |
---|
70 |
|
---|
71 |
|
---|
72 |
|
---|
73 |
|
---|
74 |
|
---|
75 |
&write_netcdf; |
---|
76 |
|
---|
77 |
|
---|
78 |
|
---|
79 |
$cummtime = time - $starttime; |
---|
80 |
if ($debug) { |
---|
81 |
print "Closing NetCDF files:\nALL DONE!"; |
---|
82 |
print "\nTotal time $cummtime sec\n"; |
---|
83 |
} |
---|
84 |
|
---|
85 |
if ($nemo) { |
---|
86 |
print "Attempting to scp file to nemo... (will require passwd)\n"; |
---|
87 |
system("scp $out_data_dir$out_file_name nemo.isis.unc.edu:/seacoos/". |
---|
88 |
"data/nc-coos/roof/"); |
---|
89 |
} |
---|
90 |
|
---|
91 |
|
---|
92 |
|
---|
93 |
|
---|
94 |
|
---|
95 |
|
---|
96 |
|
---|
97 |
sub read_data { |
---|
98 |
|
---|
99 |
|
---|
100 |
$time_index = -1; |
---|
101 |
$file_index = -1; |
---|
102 |
|
---|
103 |
|
---|
104 |
$fill_value = -9999; |
---|
105 |
|
---|
106 |
|
---|
107 |
|
---|
108 |
@tower_files = (); |
---|
109 |
|
---|
110 |
|
---|
111 |
|
---|
112 |
|
---|
113 |
foreach $date_string (@dates) { |
---|
114 |
$data_sub_dir = $month_dirs{$date_string}; |
---|
115 |
@t_files = glob("${d_data_dir}/$data_sub_dir/${tower_id}_${date_string}*.tx*"); |
---|
116 |
push(@tower_files,@t_files); |
---|
117 |
|
---|
118 |
} |
---|
119 |
|
---|
120 |
if ($debug) { |
---|
121 |
print "... Found ".eval($ |
---|
122 |
} |
---|
123 |
|
---|
124 |
|
---|
125 |
|
---|
126 |
$expected_vars{time} = "SAMPLE_LOOP"; |
---|
127 |
|
---|
128 |
|
---|
129 |
foreach $file (@tower_files) { |
---|
130 |
$file_time = $file; |
---|
131 |
|
---|
132 |
$file_time =~m/([0-9]{4})([0-9]{2})([0-9]{2})_([0-9]{2})([0-9]{2})[0-9]+.txt/; |
---|
133 |
$file_time_str = "$1/$2/$3 $4:$5:00"; |
---|
134 |
|
---|
135 |
|
---|
136 |
$file_u_date = &UnixDate(&ParseDate("$file_time_str UTC"),"%s"); |
---|
137 |
|
---|
138 |
|
---|
139 |
$epoch_dn = 719529; |
---|
140 |
|
---|
141 |
|
---|
142 |
$file_datenum = $file_u_date/86400 + $epoch_dn; |
---|
143 |
|
---|
144 |
$file_index++; |
---|
145 |
${$Vars{'file_time'}}[$file_index] = $file_datenum; |
---|
146 |
|
---|
147 |
open(DATA_FILE,$file); |
---|
148 |
|
---|
149 |
while ($data_line=<DATA_FILE>) { |
---|
150 |
|
---|
151 |
|
---|
152 |
|
---|
153 |
|
---|
154 |
&seperate_sensors; |
---|
155 |
|
---|
156 |
} |
---|
157 |
close(DATA_FILE); |
---|
158 |
|
---|
159 |
} |
---|
160 |
|
---|
161 |
|
---|
162 |
|
---|
163 |
|
---|
164 |
|
---|
165 |
} |
---|
166 |
|
---|
167 |
|
---|
168 |
|
---|
169 |
sub test_data { |
---|
170 |
print "Trying all keys of data array\n"; |
---|
171 |
foreach $param (@vars) { |
---|
172 |
print "...$param,"; |
---|
173 |
@temp_array = @{$Vars{$param}}; |
---|
174 |
print "Array: $#temp_array, $temp_array[0]\n"; |
---|
175 |
} |
---|
176 |
|
---|
177 |
} |
---|
178 |
|
---|
179 |
|
---|
180 |
sub write_netcdf { |
---|
181 |
|
---|
182 |
|
---|
183 |
|
---|
184 |
|
---|
185 |
|
---|
186 |
|
---|
187 |
|
---|
188 |
if ($debug) { |
---|
189 |
print "\nCreating NetCDF files:\n"; |
---|
190 |
} |
---|
191 |
$nc_name = "${out_data_dir}$out_file_name"; |
---|
192 |
$nc_sname = "${out_data_dir}$out_file_sname"; |
---|
193 |
|
---|
194 |
unless ($ |
---|
195 |
print "Exiting Now! \nFATAL ERROR:\n"; |
---|
196 |
die " There must be at least one valid record!\n". |
---|
197 |
" Please try another date range!\n"; |
---|
198 |
} |
---|
199 |
|
---|
200 |
|
---|
201 |
$ncid = NetCDF::create($nc_name,NetCDF::CLOBBER); |
---|
202 |
die "Couldn't open $nc_name\n" if $ncid<0; |
---|
203 |
|
---|
204 |
$ncid_system = NetCDF::create($nc_sname,NetCDF::CLOBBER); |
---|
205 |
die "Couldn't open $nc_sname\n" if $ncid_system<0; |
---|
206 |
|
---|
207 |
|
---|
208 |
|
---|
209 |
|
---|
210 |
|
---|
211 |
if ($debug) { |
---|
212 |
print "... Putting global variables\n"; |
---|
213 |
} |
---|
214 |
|
---|
215 |
|
---|
216 |
|
---|
217 |
foreach $global_att (@globals) { |
---|
218 |
$g_a_value = eval("$global_value{$global_att}"); |
---|
219 |
$g_a_format = $global_format{$global_att}; |
---|
220 |
|
---|
221 |
unless ($g_a_value =~/N\/A/) { |
---|
222 |
if ($debug>2) { |
---|
223 |
print "$global_att is $g_a_value\n" |
---|
224 |
} |
---|
225 |
$attid = NetCDF::attput($ncid,NetCDF::GLOBAL,"$global_att", |
---|
226 |
eval('NetCDF::'.$g_a_format), |
---|
227 |
"$g_a_value"); |
---|
228 |
$attid = NetCDF::attput($ncid_system,NetCDF::GLOBAL,"$global_att", |
---|
229 |
eval('NetCDF::'.$g_a_format), |
---|
230 |
"$g_a_value"); |
---|
231 |
$attid != -1 || die "\n Couldn't write global attribute $global_att\n"; |
---|
232 |
} |
---|
233 |
} |
---|
234 |
|
---|
235 |
|
---|
236 |
|
---|
237 |
|
---|
238 |
|
---|
239 |
@{$Vars{LONGITUDE}} = ($tower_lon{$tower_id}); |
---|
240 |
@{$Vars{LATITUDE}} = ($tower_lat{$tower_id}); |
---|
241 |
|
---|
242 |
|
---|
243 |
|
---|
244 |
%dims = ( |
---|
245 |
LONGITUDE => 1, |
---|
246 |
LATITUDE => 1, |
---|
247 |
TOWER_FILES => $ |
---|
248 |
CHARACTERS => 50, |
---|
249 |
time => $ |
---|
250 |
); |
---|
251 |
%dim_names = ( |
---|
252 |
'a' => 'LATITUDE', |
---|
253 |
'o' => 'LONGITUDE', |
---|
254 |
'f' => 'TOWER_FILES', |
---|
255 |
'c' => 'CHARACTERS', |
---|
256 |
't' => 'time' |
---|
257 |
); |
---|
258 |
|
---|
259 |
|
---|
260 |
foreach $dim (keys %dims) { |
---|
261 |
if ($dimensions{$dim}=~/[fc]/){ |
---|
262 |
$ncid{$dim} = $ncid_system; |
---|
263 |
} |
---|
264 |
elsif ($dimensions{$dim}=~/[tao]/) { |
---|
265 |
$ncid{$dim} = $ncid; |
---|
266 |
} |
---|
267 |
|
---|
268 |
unless ($dim =~/time/) { |
---|
269 |
$dimid{$dim} = NetCDF::dimdef($ncid{$dim},$dim,$dims{$dim}); |
---|
270 |
die "Couldn't define $dim dimension\n" if $dimid{$dim}<0; |
---|
271 |
} |
---|
272 |
} |
---|
273 |
|
---|
274 |
|
---|
275 |
|
---|
276 |
$dimid{time} = NetCDF::dimdef($ncid,"time",NetCDF::UNLIMITED); |
---|
277 |
die "Couldn't define record dimension\n" if $dimid{time}<0; |
---|
278 |
|
---|
279 |
|
---|
280 |
|
---|
281 |
|
---|
282 |
if ($debug) { |
---|
283 |
print "... Defining variables and putting attributes\n"; |
---|
284 |
} |
---|
285 |
|
---|
286 |
|
---|
287 |
|
---|
288 |
|
---|
289 |
@all_attributes = ('long_name', 'standard_name', 'units', |
---|
290 |
'description', 'range', 'precision' |
---|
291 |
); |
---|
292 |
|
---|
293 |
@zeros = (0,0,0); |
---|
294 |
|
---|
295 |
|
---|
296 |
foreach $param (@vars) { |
---|
297 |
if ($param =~/CLOCK/ || $param =~/BBAD/) { |
---|
298 |
|
---|
299 |
} |
---|
300 |
elsif (defined $Vars{$param}) { |
---|
301 |
&vardef_attput; |
---|
302 |
} |
---|
303 |
} |
---|
304 |
|
---|
305 |
|
---|
306 |
|
---|
307 |
$status = NetCDF::endef($ncid); |
---|
308 |
die "Couldn't end definition\n" if $status<0; |
---|
309 |
$status = NetCDF::endef($ncid_system); |
---|
310 |
die "Couldn't end definition\n" if $status<0; |
---|
311 |
|
---|
312 |
|
---|
313 |
|
---|
314 |
|
---|
315 |
if ($debug) { |
---|
316 |
print "... Putting variables:\n\n"; |
---|
317 |
} |
---|
318 |
|
---|
319 |
foreach $param (@vars) { |
---|
320 |
if ($param =~/CLOCK/ || $param =~/BBAD/) { |
---|
321 |
|
---|
322 |
} |
---|
323 |
elsif (defined $Vars{$param}) { |
---|
324 |
&varput; |
---|
325 |
} |
---|
326 |
} |
---|
327 |
|
---|
328 |
|
---|
329 |
|
---|
330 |
$status = NetCDF::close($ncid); |
---|
331 |
$status = NetCDF::close($ncid_system); |
---|
332 |
|
---|
333 |
} |
---|
334 |
|
---|
335 |
|
---|
336 |
|
---|
337 |
sub seperate_sensors { |
---|
338 |
|
---|
339 |
|
---|
340 |
|
---|
341 |
|
---|
342 |
|
---|
343 |
|
---|
344 |
|
---|
345 |
@data_line = split(', ',$data_line); |
---|
346 |
@va = (); |
---|
347 |
|
---|
348 |
|
---|
349 |
|
---|
350 |
if ($ |
---|
351 |
|
---|
352 |
if ($ |
---|
353 |
($time,$param_id,$type,$value,$units,$error_string) = @data_line; |
---|
354 |
} |
---|
355 |
|
---|
356 |
|
---|
357 |
elsif ($ |
---|
358 |
($time,$param_id,$type,$va[0],$va[1],$va[2],$va[3],$units, |
---|
359 |
$error_string) = @data_line; |
---|
360 |
$value = join(',',@va); |
---|
361 |
$value =~s/\s+//g; |
---|
362 |
} |
---|
363 |
|
---|
364 |
|
---|
365 |
elsif ($ |
---|
366 |
|
---|
367 |
|
---|
368 |
($time,$param_id,$type,$va[0],$va[1],$va[2],$va[3],$va[4], |
---|
369 |
$units,$error_string) = @data_line; |
---|
370 |
$value = join(',',@va); |
---|
371 |
$value =~s/\s+//g; |
---|
372 |
} |
---|
373 |
|
---|
374 |
|
---|
375 |
elsif ($ |
---|
376 |
|
---|
377 |
|
---|
378 |
($time,$param_id,$type,$va[0],$va[1],$va[2],$va[3],$va[4], |
---|
379 |
$va[5],$units,$error_string) = @data_line; |
---|
380 |
$value = join(',',@va); |
---|
381 |
$value =~s/\s+//g; |
---|
382 |
} |
---|
383 |
|
---|
384 |
else { |
---|
385 |
|
---|
386 |
|
---|
387 |
|
---|
388 |
($time,$param_id,$type,@junk) = @data_line; |
---|
389 |
|
---|
390 |
|
---|
391 |
|
---|
392 |
$value = "$fill_value,$fill_value,$fill_value,$fill_value"; |
---|
393 |
|
---|
394 |
} |
---|
395 |
if($param_id=~m/CTD1/) { |
---|
396 |
|
---|
397 |
|
---|
398 |
|
---|
399 |
} |
---|
400 |
|
---|
401 |
|
---|
402 |
|
---|
403 |
|
---|
404 |
|
---|
405 |
|
---|
406 |
|
---|
407 |
|
---|
408 |
|
---|
409 |
$big_type = $type; |
---|
410 |
$big_type =~s/\([0-9]\)//g; |
---|
411 |
|
---|
412 |
|
---|
413 |
|
---|
414 |
|
---|
415 |
|
---|
416 |
|
---|
417 |
@sub_vars = @{$sub_vars{$param_id}}; |
---|
418 |
|
---|
419 |
&compute_sub_vars; |
---|
420 |
|
---|
421 |
|
---|
422 |
|
---|
423 |
|
---|
424 |
|
---|
425 |
|
---|
426 |
if ($type=~/SETTING/ && $param_id=~/ACTIVE/) { |
---|
427 |
|
---|
428 |
|
---|
429 |
|
---|
430 |
|
---|
431 |
|
---|
432 |
|
---|
433 |
|
---|
434 |
|
---|
435 |
|
---|
436 |
|
---|
437 |
$value =~s/\"//g; |
---|
438 |
|
---|
439 |
@value = split("","$value"); |
---|
440 |
@this_val = @value; |
---|
441 |
|
---|
442 |
unless ($ |
---|
443 |
@value = ("-"); |
---|
444 |
} |
---|
445 |
|
---|
446 |
for ($ii=$ |
---|
447 |
@this_val[$ii] = " "; |
---|
448 |
} |
---|
449 |
|
---|
450 |
push(@{$Vars{$param_id}},@this_val); |
---|
451 |
|
---|
452 |
$expected_vars = $value; |
---|
453 |
$expected_vars=~s/\"//g; |
---|
454 |
@expected_vars = split('\s',$expected_vars); |
---|
455 |
|
---|
456 |
|
---|
457 |
my($t_type) = $param_id; |
---|
458 |
$t_type =~s/^ACTIVE\_//; |
---|
459 |
|
---|
460 |
@{$expected_vars{$t_type}} = @expected_vars; |
---|
461 |
} |
---|
462 |
elsif ($type=~/CPU/) { |
---|
463 |
|
---|
464 |
$value =~s/\"//g; |
---|
465 |
@value = split("","$value"); |
---|
466 |
@this_val = @value; |
---|
467 |
unless ($ |
---|
468 |
@value = (" "); |
---|
469 |
} |
---|
470 |
for ($ii=$ |
---|
471 |
@this_val[$ii] = " "; |
---|
472 |
} |
---|
473 |
push(@{$Vars{$param_id}},@this_val); |
---|
474 |
} |
---|
475 |
elsif ($type=~/CLOCK/) { |
---|
476 |
|
---|
477 |
$value =~s/\"//g; |
---|
478 |
@value = split("","$value"); |
---|
479 |
@this_val = @value; |
---|
480 |
unless ($ |
---|
481 |
@value = (" "); |
---|
482 |
} |
---|
483 |
for ($ii=$ |
---|
484 |
@this_val[$ii] = " "; |
---|
485 |
} |
---|
486 |
push(@{$Vars{$param_id}},@this_val); |
---|
487 |
} |
---|
488 |
|
---|
489 |
|
---|
490 |
|
---|
491 |
|
---|
492 |
|
---|
493 |
|
---|
494 |
|
---|
495 |
|
---|
496 |
elsif ($ |
---|
497 |
push(@{$loop_vars{$big_type}},$param_id); |
---|
498 |
|
---|
499 |
if ($param_id=~/LINK$/) { |
---|
500 |
if ($first_one) { |
---|
501 |
@expected_vars = ('LINK','LINK_AMPS','LINK_VOLTS','LINK_AMP_HRS', |
---|
502 |
'LINK_TEMP'); |
---|
503 |
push(@{$expected_vars{$big_type}},@expected_vars); |
---|
504 |
|
---|
505 |
@expected_vars = ('SAMPLE_LOOP','BBAD_SAMPLE_LOOP'); |
---|
506 |
push(@{$expected_vars{TIME}},@expected_vars); |
---|
507 |
|
---|
508 |
@expected_vars = ('PIR'); |
---|
509 |
push(@{$expected_vars{RLCAD}},@expected_vars); |
---|
510 |
|
---|
511 |
$first_one = 0; |
---|
512 |
} |
---|
513 |
push(@{$loop_vars{$big_type}}, 'LINK','LINK_AMPS','LINK_VOLTS', |
---|
514 |
'LINK_AMP_HRS','LINK_TEMP'); |
---|
515 |
} |
---|
516 |
} |
---|
517 |
|
---|
518 |
|
---|
519 |
|
---|
520 |
|
---|
521 |
if ($param_id=~/ATEMP/) { |
---|
522 |
$time_index++; |
---|
523 |
|
---|
524 |
|
---|
525 |
|
---|
526 |
|
---|
527 |
|
---|
528 |
|
---|
529 |
foreach $expected_type (keys %expected_vars) { |
---|
530 |
$num_expected = $ |
---|
531 |
$num_recieved = $ |
---|
532 |
|
---|
533 |
|
---|
534 |
|
---|
535 |
|
---|
536 |
if ($num_expected>$num_recieved && $num_expected && $time_index>=0) { |
---|
537 |
|
---|
538 |
|
---|
539 |
|
---|
540 |
|
---|
541 |
foreach $expected_var (@{$expected_vars{$expected_type}}) { |
---|
542 |
|
---|
543 |
|
---|
544 |
${$Vars{$expected_var}}[$time_index] = $fill_value; |
---|
545 |
} |
---|
546 |
|
---|
547 |
|
---|
548 |
|
---|
549 |
|
---|
550 |
|
---|
551 |
|
---|
552 |
|
---|
553 |
unless(${$Vars{"CLOCK"}}[$time_index]) { |
---|
554 |
${$Vars{"CLOCK"}}[$time_index] = $fill_value; |
---|
555 |
} |
---|
556 |
} |
---|
557 |
} |
---|
558 |
|
---|
559 |
%loop_vars = (); |
---|
560 |
|
---|
561 |
|
---|
562 |
|
---|
563 |
|
---|
564 |
|
---|
565 |
|
---|
566 |
$platform_time[$time_index] = $time; |
---|
567 |
$u_date = &UnixDate(&ParseDate("$time"),"%Y,%m,%d,%H,%M,%S"); |
---|
568 |
($year,$month,$day,$hour,$minute,$second) = split(',',$u_date); |
---|
569 |
|
---|
570 |
$u_date = &UnixDate(&ParseDate("$time UTC"),"%s"); |
---|
571 |
$um1_date = &UnixDate(&ParseDate("$platform_time[$time_index -1] UTC"), |
---|
572 |
"%s"); |
---|
573 |
$sample_interval = $u_date - $um1_date; |
---|
574 |
$hrs_per_sample = $sample_interval/3600; |
---|
575 |
|
---|
576 |
unless ($time_index) { |
---|
577 |
$hrs_per_sample = 1/30; |
---|
578 |
} |
---|
579 |
|
---|
580 |
$epoch_dn = 719529; |
---|
581 |
|
---|
582 |
|
---|
583 |
$datenum = $u_date/86400 + $epoch_dn; |
---|
584 |
|
---|
585 |
${$Vars{'year'}}[$time_index] = $year; |
---|
586 |
${$Vars{'month'}}[$time_index] = $month; |
---|
587 |
${$Vars{'day'}}[$time_index] = $day; |
---|
588 |
${$Vars{'hour'}}[$time_index] = $hour; |
---|
589 |
${$Vars{'minute'}}[$time_index] = $minute; |
---|
590 |
${$Vars{'second'}}[$time_index] = $second; |
---|
591 |
${$Vars{'time'}}[$time_index] = $datenum; |
---|
592 |
|
---|
593 |
} |
---|
594 |
|
---|
595 |
|
---|
596 |
&parse_vars; |
---|
597 |
|
---|
598 |
if ($param_id=~/GPS1/) { |
---|
599 |
$param_id = 'GPSLON'; |
---|
600 |
&convert_vars; |
---|
601 |
$param_id = 'GPSLAT'; |
---|
602 |
&convert_vars; |
---|
603 |
$param_id = 'GPS1'; |
---|
604 |
} |
---|
605 |
|
---|
606 |
|
---|
607 |
|
---|
608 |
|
---|
609 |
&convert_vars; |
---|
610 |
|
---|
611 |
|
---|
612 |
|
---|
613 |
|
---|
614 |
|
---|
615 |
|
---|
616 |
if ($param_id =~/PIRC/ | $param_id =~/PIRA/) { |
---|
617 |
|
---|
618 |
if ($Vars{PIRA}[-1]!=~/$fill_value/ && |
---|
619 |
$Vars{PIRC}[-1]!=~/$fill_value/){ |
---|
620 |
$temporary_param = $param_id; |
---|
621 |
$param_id = 'PIR'; |
---|
622 |
|
---|
623 |
$pir_ct = $Vars{PIRA}[-1]; |
---|
624 |
$pir_th = $Vars{PIRC}[-1]; |
---|
625 |
|
---|
626 |
&convert_vars; |
---|
627 |
$param_id = $temporary_param; |
---|
628 |
} |
---|
629 |
else { |
---|
630 |
$Vars{PIR} = $fill_value; |
---|
631 |
} |
---|
632 |
} |
---|
633 |
$value = undef; |
---|
634 |
|
---|
635 |
} |
---|
636 |
|
---|
637 |
|
---|
638 |
else { |
---|
639 |
|
---|
640 |
|
---|
641 |
} |
---|
642 |
} |
---|
643 |
|
---|
644 |
|
---|
645 |
|
---|
646 |
sub parse_vars { |
---|
647 |
|
---|
648 |
|
---|
649 |
|
---|
650 |
|
---|
651 |
|
---|
652 |
|
---|
653 |
$output = $parse_output{$param_id}; |
---|
654 |
@output = split(',',$output); |
---|
655 |
|
---|
656 |
|
---|
657 |
if ($units{$param_id}=~/string/) { |
---|
658 |
|
---|
659 |
eval("$output[0] = '$value';"); |
---|
660 |
} |
---|
661 |
else { |
---|
662 |
|
---|
663 |
$old_v = $value; |
---|
664 |
|
---|
665 |
|
---|
666 |
|
---|
667 |
|
---|
668 |
if ($parse_info{$param_id}=~/\@avg/) { |
---|
669 |
|
---|
670 |
$old_v =~s/\ \-9999//g; |
---|
671 |
@value = split(' ',$old_v); |
---|
672 |
eval("$output[0]".' = @value;'); |
---|
673 |
} |
---|
674 |
elsif ($parse_info{$param_id}=~/\@gva/) { |
---|
675 |
|
---|
676 |
@comp = @pitch = @roll = (); |
---|
677 |
|
---|
678 |
@value = split('\$',$old_v); |
---|
679 |
foreach $vals (@value) { |
---|
680 |
$vals =~m/$C([0-9\.\-]*)P([0-9\.\-]*)R([0-9\.\-]*)T.*/; |
---|
681 |
push(@comp,$1); |
---|
682 |
push(@pitch,$2); |
---|
683 |
push(@roll,$3); |
---|
684 |
|
---|
685 |
} |
---|
686 |
|
---|
687 |
} |
---|
688 |
|
---|
689 |
|
---|
690 |
elsif ($ |
---|
691 |
|
---|
692 |
|
---|
693 |
|
---|
694 |
|
---|
695 |
|
---|
696 |
|
---|
697 |
$value =~m/$parse_info{$param_id}/; |
---|
698 |
@value = $1; |
---|
699 |
|
---|
700 |
|
---|
701 |
if (defined $value[0]) { |
---|
702 |
eval("$output[0]".' = $value[0];'); |
---|
703 |
} |
---|
704 |
else { |
---|
705 |
eval("$output[0] = $fill_value;"); |
---|
706 |
} |
---|
707 |
|
---|
708 |
} |
---|
709 |
|
---|
710 |
|
---|
711 |
elsif ($ |
---|
712 |
if ($param_id =~m/CTD/) { |
---|
713 |
|
---|
714 |
if ($value=~m/^\"\ |
---|
715 |
$pi = $parse_info{$param_id}; |
---|
716 |
} |
---|
717 |
else { |
---|
718 |
$pi = 'NOBAD!!'; |
---|
719 |
} |
---|
720 |
} |
---|
721 |
else { |
---|
722 |
|
---|
723 |
$pi = $parse_info{$param_id}; |
---|
724 |
} |
---|
725 |
if ($param_id=~/LINK/) { |
---|
726 |
|
---|
727 |
} |
---|
728 |
elsif ($param_id=~/GPS1/) { |
---|
729 |
|
---|
730 |
@v = split('',$value); |
---|
731 |
$mygps_lat = join('',@v[11..17]); |
---|
732 |
$mygps_lon = join('',@v[19..25]); |
---|
733 |
|
---|
734 |
} |
---|
735 |
$value =~m/$pi/; |
---|
736 |
|
---|
737 |
|
---|
738 |
|
---|
739 |
|
---|
740 |
|
---|
741 |
|
---|
742 |
|
---|
743 |
|
---|
744 |
for ($out_num=1; $out_num<=$ |
---|
745 |
eval("$output[$out_num-1] =\$$out_num;"); |
---|
746 |
} |
---|
747 |
if ($param_id=~/LINK/) { |
---|
748 |
|
---|
749 |
} |
---|
750 |
elsif ($param_id=~/GPS1/) { |
---|
751 |
if ($mygps_lat+0) {$gps_lat = $mygps_lat;} |
---|
752 |
else {$gps_lat = $fill_value;} |
---|
753 |
if ($mygps_lon+0) {$gps_lon = $mygps_lon;} |
---|
754 |
else {$gps_lon = $fill_value;} |
---|
755 |
} |
---|
756 |
} |
---|
757 |
} |
---|
758 |
} |
---|
759 |
|
---|
760 |
|
---|
761 |
|
---|
762 |
sub convert_vars { |
---|
763 |
our($fill_value); |
---|
764 |
|
---|
765 |
|
---|
766 |
|
---|
767 |
|
---|
768 |
|
---|
769 |
|
---|
770 |
|
---|
771 |
|
---|
772 |
|
---|
773 |
eval ('$c_coeff = $'.$tower_id.'_conversion_coeff{$param_id};'); |
---|
774 |
eval ('$c_formula = $'.$tower_id.'_conversion_formula{$param_id};'); |
---|
775 |
|
---|
776 |
|
---|
777 |
@v = split(',',$c_coeff); |
---|
778 |
|
---|
779 |
|
---|
780 |
|
---|
781 |
if ($ |
---|
782 |
$temp_var = ''; |
---|
783 |
if ($c_formula=~/^\&/) { |
---|
784 |
|
---|
785 |
eval("$c_formula"); |
---|
786 |
} |
---|
787 |
else { |
---|
788 |
|
---|
789 |
my ($out) = eval("$output[0]"); |
---|
790 |
unless($out=~/$fill_value/) { |
---|
791 |
$temp_var = eval("$c_formula"); |
---|
792 |
} |
---|
793 |
|
---|
794 |
|
---|
795 |
|
---|
796 |
else { |
---|
797 |
$temp_var = $fill_value; |
---|
798 |
} |
---|
799 |
|
---|
800 |
|
---|
801 |
&set_index; |
---|
802 |
|
---|
803 |
|
---|
804 |
if ($param_id=~/RAIN/) { |
---|
805 |
|
---|
806 |
$temp_var = $temp_var/$hrs_per_sample; |
---|
807 |
} |
---|
808 |
|
---|
809 |
|
---|
810 |
|
---|
811 |
unless ($index +1) { |
---|
812 |
print "Warning. $param_id has index $index. Adjusting...\n"; |
---|
813 |
${$Vars{$param_id}}[$index+1] = $temp_var; |
---|
814 |
} |
---|
815 |
else { |
---|
816 |
${$Vars{$param_id}}[$index] = $temp_var; |
---|
817 |
} |
---|
818 |
} |
---|
819 |
} |
---|
820 |
|
---|
821 |
else { |
---|
822 |
|
---|
823 |
} |
---|
824 |
} |
---|
825 |
|
---|
826 |
|
---|
827 |
|
---|
828 |
|
---|
829 |
sub handle_input { |
---|
830 |
|
---|
831 |
|
---|
832 |
|
---|
833 |
open(DEFAULTS,"${bin_dir}tower_parse.defaults"); |
---|
834 |
while (<DEFAULTS>) { |
---|
835 |
|
---|
836 |
chop; |
---|
837 |
|
---|
838 |
$_=~s/ |
---|
839 |
|
---|
840 |
if ($_) { |
---|
841 |
my($var_name, $default_value) = split(','); |
---|
842 |
eval("$var_name = $default_value;"); |
---|
843 |
} |
---|
844 |
} |
---|
845 |
|
---|
846 |
($date_string,%varargin) = @ARGV; |
---|
847 |
|
---|
848 |
|
---|
849 |
%varargout = ( |
---|
850 |
'-out' => '$out_file_name', |
---|
851 |
'-d' => '$debug', |
---|
852 |
'-t' => '$tower_id', |
---|
853 |
'-p' => '$package', |
---|
854 |
); |
---|
855 |
|
---|
856 |
|
---|
857 |
|
---|
858 |
|
---|
859 |
foreach $argin (keys %varargin) { |
---|
860 |
eval("$varargout{$argin} = '$varargin{$argin}';"); |
---|
861 |
} |
---|
862 |
|
---|
863 |
$d_data_dir = "$d_data_dir$tower_id/"; |
---|
864 |
$out_data_dir = "$out_data_dir$tower_id/"; |
---|
865 |
|
---|
866 |
|
---|
867 |
|
---|
868 |
if ($date_string) { |
---|
869 |
&handle_date; |
---|
870 |
} |
---|
871 |
|
---|
872 |
|
---|
873 |
else { |
---|
874 |
die "Must at least input a date string!\n". |
---|
875 |
"./tower_parse.perl '200406' -out 2004_06.nc -d 1\n"; |
---|
876 |
} |
---|
877 |
|
---|
878 |
|
---|
879 |
unless ($out_file_name=~/.nc$/) { |
---|
880 |
if ($out_file_name) { |
---|
881 |
$out_file_sname = "${out_file_name}_sys.nc"; |
---|
882 |
$out_file_name = "${out_file_name}.nc"; |
---|
883 |
} |
---|
884 |
else { |
---|
885 |
$out_file_name = $d_out_file_name; |
---|
886 |
@do = split('',$d_out_file_name); |
---|
887 |
$out_file_sname = join('',@do[0..$ |
---|
888 |
} |
---|
889 |
} |
---|
890 |
else { |
---|
891 |
unless ($out_file_sname) { |
---|
892 |
@do = split('',$out_file_name); |
---|
893 |
$out_file_sname = join('',@do[0..$ |
---|
894 |
} |
---|
895 |
} |
---|
896 |
|
---|
897 |
|
---|
898 |
|
---|
899 |
if ($debug) { |
---|
900 |
print "\n\nData dir : $d_data_dir\nOutput dir : $out_data_dir\n". |
---|
901 |
"Outfile name : $out_file_name and $out_file_sname\n"; |
---|
902 |
if ($ |
---|
903 |
print "\nLooking for data between $dates[0] and $dates[-1]\n"; |
---|
904 |
} |
---|
905 |
} |
---|
906 |
|
---|
907 |
|
---|
908 |
} |
---|
909 |
|
---|
910 |
|
---|
911 |
|
---|
912 |
|
---|
913 |
sub handle_date { |
---|
914 |
|
---|
915 |
if ($date_string=~/\ /) { |
---|
916 |
($st_date,$en_date) = split(' ',$date_string); |
---|
917 |
@st_date = split('',$st_date); |
---|
918 |
@en_date = split('',$en_date); |
---|
919 |
|
---|
920 |
|
---|
921 |
$st_date = join('',(@st_date[0..3],' ',@st_date[4..5],' ',@st_date[6..7])); |
---|
922 |
$en_date = join('',(@en_date[0..3],' ',@en_date[4..5],' ',@en_date[6..7])); |
---|
923 |
|
---|
924 |
|
---|
925 |
$this_date = $st_date; |
---|
926 |
$this_parsed = &ParseDate($st_date); |
---|
927 |
$en_parsed = &ParseDate($en_date); |
---|
928 |
|
---|
929 |
@dates = (); |
---|
930 |
%month_dirs = (); |
---|
931 |
|
---|
932 |
while ($this_parsed <= $en_parsed) { |
---|
933 |
$ymd = &UnixDate($this_parsed,"%Y%m%d"); |
---|
934 |
push(@dates,$ymd); |
---|
935 |
$month_dirs{$ymd} = &UnixDate($this_parsed,"%Y_%m"); |
---|
936 |
$this_date = &UnixDate(&DateCalc($this_date,"+ 1 day"), "%Y %m %d"); |
---|
937 |
$this_parsed = &ParseDate($this_date); |
---|
938 |
} |
---|
939 |
|
---|
940 |
} |
---|
941 |
elsif ($date_string=~/latest/) { |
---|
942 |
$today = &ParseDate("today"); |
---|
943 |
$date1 = &UnixDate(&DateCalc("$today","- 2 day"),"%Y%m%d"); |
---|
944 |
$date2 = &UnixDate(&DateCalc("$today","- 1 day"),"%Y%m%d"); |
---|
945 |
$date3 = &UnixDate("$today","%Y%m%d"); |
---|
946 |
$mond1 = &UnixDate(&DateCalc("$today","- 2 day"),"%Y_%m"); |
---|
947 |
$mond2 = &UnixDate(&DateCalc("$today","- 1 day"),"%Y_%m"); |
---|
948 |
$mond3 = &UnixDate("$today","%Y_%m"); |
---|
949 |
@dates = ($date1,$date2,$date3); |
---|
950 |
%month_dirs = ($date1,$mond1,$date2,$mond2,$date3,$mond3); |
---|
951 |
|
---|
952 |
} |
---|
953 |
elsif ($date_string =~/help/) { |
---|
954 |
|
---|
955 |
|
---|
956 |
print "\nUsage: ./tower_parse.perl date_string\n\n". |
---|
957 |
"Optional argument pairs:\n". |
---|
958 |
" -out output file name\n". |
---|
959 |
" -d debug\n". |
---|
960 |
|
---|
961 |
" -t tower id\n". |
---|
962 |
" -nemo scp data file to nemo\n\n"; |
---|
963 |
print "Date string can have any of the following formats:\n". |
---|
964 |
" yyyymm Finds all data for the month\n". |
---|
965 |
" yyyymmdd Finds data just for one day\n". |
---|
966 |
" 'yyyymmdd yyyymmdd' Finds data within the date range\n". |
---|
967 |
" latest Finds today, yesterday and the previous day\n"; |
---|
968 |
die "\n"; |
---|
969 |
} |
---|
970 |
else { |
---|
971 |
@dates = ($date_string); |
---|
972 |
@date_str = split('',$date_string); |
---|
973 |
$fn = join('',(@date_str[0..3],'_',@date_str[4..5])); |
---|
974 |
|
---|
975 |
unless($out_file_name) { |
---|
976 |
|
---|
977 |
$out_file_name = "${tower_id}_${fn}.nc"; |
---|
978 |
$out_file_sname = "${tower_id}_${fn}_sys.nc"; |
---|
979 |
} |
---|
980 |
|
---|
981 |
|
---|
982 |
$month_dirs{$date_string} = $fn; |
---|
983 |
} |
---|
984 |
@month_dirs = values (%month_dirs); |
---|
985 |
} |
---|
986 |
|
---|
987 |
|
---|
988 |
|
---|
989 |
|
---|
990 |
sub parse_adcp { |
---|
991 |
|
---|
992 |
|
---|
993 |
|
---|
994 |
|
---|
995 |
} |
---|
996 |
|
---|
997 |
|
---|
998 |
|
---|
999 |
|
---|
1000 |
sub vardef_attput { |
---|
1001 |
|
---|
1002 |
|
---|
1003 |
|
---|
1004 |
$format = 'NetCDF::'.$format{$param}; |
---|
1005 |
if ($debug>2) { |
---|
1006 |
print "$param is of format $format\n"; |
---|
1007 |
} |
---|
1008 |
|
---|
1009 |
@dimensions = split('\s',$dimensions{$param}); |
---|
1010 |
@dim_names = @dim_names{@dimensions}; |
---|
1011 |
@dim_ids = @dimid{@dim_names}; |
---|
1012 |
|
---|
1013 |
$dimensions = join('',@dimensions); |
---|
1014 |
|
---|
1015 |
|
---|
1016 |
|
---|
1017 |
if ($dimensions=~/f/) { |
---|
1018 |
$ncid{$param} = $ncid_system; |
---|
1019 |
} |
---|
1020 |
elsif ($dimensions=~/t/) { |
---|
1021 |
$ncid{$param} = $ncid; |
---|
1022 |
} |
---|
1023 |
|
---|
1024 |
|
---|
1025 |
$varid{$param} = NetCDF::vardef($ncid{$param},$param, |
---|
1026 |
eval("$format"), |
---|
1027 |
[@dim_ids]); |
---|
1028 |
|
---|
1029 |
die "Couldn't define variable $param\n" if $varid{$param}<0; |
---|
1030 |
|
---|
1031 |
|
---|
1032 |
|
---|
1033 |
foreach $attribute (@all_attributes) { |
---|
1034 |
$att_val = eval('\$'.$attribute.'{$param}'); |
---|
1035 |
if (defined $att_val && $att_val!=~/^N\/A/) { |
---|
1036 |
$attid = NetCDF::attput($ncid{$param},$varid{$param},"$attribute", |
---|
1037 |
NetCDF::CHAR, |
---|
1038 |
eval('\$'.$attribute.'{$param}')); |
---|
1039 |
} |
---|
1040 |
} |
---|
1041 |
} |
---|
1042 |
|
---|
1043 |
|
---|
1044 |
|
---|
1045 |
sub varput { |
---|
1046 |
$format = 'NetCDF::'.$format{$param}; |
---|
1047 |
|
---|
1048 |
@dimensions = split('\s',$dimensions{$param}); |
---|
1049 |
@dim_names = @dim_names{@dimensions}; |
---|
1050 |
|
---|
1051 |
@dim_counts = @dims{@dim_names}; |
---|
1052 |
@dim_starts = @zeros[0..$ |
---|
1053 |
|
---|
1054 |
unless ($param =~/CLOCK/ || $param =~/BBAD/ || |
---|
1055 |
$param =~/GPS$/ || $param =~/LINK$/) { |
---|
1056 |
|
---|
1057 |
|
---|
1058 |
|
---|
1059 |
|
---|
1060 |
@temp_array = @{$Vars{$param}}; |
---|
1061 |
for ($ii=0; $ii<=$ |
---|
1062 |
unless (defined($temp_array[$ii])) { |
---|
1063 |
${$Vars{$param}}[$ii] = $fill_value; |
---|
1064 |
} |
---|
1065 |
} |
---|
1066 |
|
---|
1067 |
if ($debug>1) { |
---|
1068 |
print "$param length ".eval($ |
---|
1069 |
", id $varid{$param}\n"; |
---|
1070 |
} |
---|
1071 |
|
---|
1072 |
unless ($format=~/CHAR/) { |
---|
1073 |
|
---|
1074 |
$status = NetCDF::varput($ncid{$param},$varid{$param}, |
---|
1075 |
\@dim_starts,\@dim_counts, |
---|
1076 |
\@{$Vars{$param}}); |
---|
1077 |
} |
---|
1078 |
else { |
---|
1079 |
|
---|
1080 |
if ($ |
---|
1081 |
$status = NetCDF::varput($ncid{$param},$varid{$param}, |
---|
1082 |
\@dim_starts,\@dim_counts, |
---|
1083 |
\@{$Vars{$param}}); |
---|
1084 |
} |
---|
1085 |
elsif ($debug>1) { |
---|
1086 |
print "Didn't have enough values for $param\n"; |
---|
1087 |
} |
---|
1088 |
} |
---|
1089 |
} |
---|
1090 |
else { |
---|
1091 |
print "$param being ignored\n"; |
---|
1092 |
} |
---|
1093 |
} |
---|
1094 |
|
---|
1095 |
|
---|
1096 |
|
---|
1097 |
sub compute_sub_vars { |
---|
1098 |
|
---|
1099 |
if ($ |
---|
1100 |
$orig_param = $param_id; |
---|
1101 |
|
---|
1102 |
foreach $sub_param (@sub_vars) { |
---|
1103 |
$param_id = $sub_param; |
---|
1104 |
&convert_vars; |
---|
1105 |
} |
---|
1106 |
$param_id = $orig_param; |
---|
1107 |
} |
---|
1108 |
} |
---|
1109 |
|
---|
1110 |
|
---|
1111 |
|
---|
1112 |
sub set_index { |
---|
1113 |
|
---|
1114 |
|
---|
1115 |
if ($dimensions{$param_id}=~/t/) { |
---|
1116 |
$index = $time_index; |
---|
1117 |
} |
---|
1118 |
|
---|
1119 |
elsif ($dimensions{$param_id}=~/f/) { |
---|
1120 |
$index = $file_index; |
---|
1121 |
} |
---|
1122 |
} |
---|
1123 |
|
---|
1124 |
|
---|
1125 |
|
---|
1126 |
|
---|
1127 |
sub avg_wind { |
---|
1128 |
|
---|
1129 |
|
---|
1130 |
|
---|
1131 |
|
---|
1132 |
|
---|
1133 |
|
---|
1134 |
|
---|
1135 |
|
---|
1136 |
|
---|
1137 |
|
---|
1138 |
|
---|
1139 |
|
---|
1140 |
|
---|
1141 |
|
---|
1142 |
|
---|
1143 |
|
---|
1144 |
|
---|
1145 |
|
---|
1146 |
my(@junk) = (); |
---|
1147 |
my(@input) = (@_); |
---|
1148 |
%names = ('spda'=>'wi_spd_a', |
---|
1149 |
'dira'=>'wi_dir_a', |
---|
1150 |
'spdb'=>'wi_spd_b', |
---|
1151 |
'dirb'=>'wi_dir_b', |
---|
1152 |
'spdn'=>'junk', |
---|
1153 |
'dirn'=>'junk', |
---|
1154 |
'cmp1'=>'comp_1', |
---|
1155 |
'cmp2'=>'comp_2'); |
---|
1156 |
|
---|
1157 |
|
---|
1158 |
|
---|
1159 |
if ($input[0]=~/spda/) { |
---|
1160 |
|
---|
1161 |
|
---|
1162 |
|
---|
1163 |
if ($dira[0]) { |
---|
1164 |
$real_last_var = $last_var; |
---|
1165 |
|
---|
1166 |
if ($ |
---|
1167 |
|
---|
1168 |
@cmp1 = (0); |
---|
1169 |
} |
---|
1170 |
} |
---|
1171 |
else { |
---|
1172 |
|
---|
1173 |
$real_last_var = 'NOPE'; |
---|
1174 |
} |
---|
1175 |
|
---|
1176 |
|
---|
1177 |
@spda = @dira = @spdb = @dirb = @cmp1 = @cmp2 = (); |
---|
1178 |
} |
---|
1179 |
elsif ($input[0]=~/dir/) { |
---|
1180 |
if ($ |
---|
1181 |
|
---|
1182 |
|
---|
1183 |
$offset = $input[1]; |
---|
1184 |
} |
---|
1185 |
else { |
---|
1186 |
$offset = 0; |
---|
1187 |
} |
---|
1188 |
} |
---|
1189 |
|
---|
1190 |
eval("\@$input[0] = \@$names{$input[0]};"); |
---|
1191 |
|
---|
1192 |
if ($input[0]=~/$real_last_var/) { |
---|
1193 |
|
---|
1194 |
&handle_wind; |
---|
1195 |
&write_wind; |
---|
1196 |
} |
---|
1197 |
|
---|
1198 |
$last_var = $input[0]; |
---|
1199 |
} |
---|
1200 |
|
---|
1201 |
|
---|
1202 |
|
---|
1203 |
|
---|
1204 |
sub handle_wind { |
---|
1205 |
|
---|
1206 |
|
---|
1207 |
|
---|
1208 |
my($time); |
---|
1209 |
our($u_mean_a,$u_max_a,$u_std_a,$u_mean_b,$u_max_b,$u_std_b, |
---|
1210 |
$v_mean_a,$v_max_a,$v_std_a,$v_mean_b,$v_max_b,$v_std_b, |
---|
1211 |
$mean_spd_a,$mean_dir_a,$mean_spd_b,$mean_dir_b,$spd_std_a, |
---|
1212 |
$dir_std_a,$spd_std_b,$dir_std_b,$gust_a,$gust_b) = $fill_value; |
---|
1213 |
|
---|
1214 |
my(@u_a,@v_a,@u_b,@v_b,@temp_cmp,@comp_off) = (); |
---|
1215 |
|
---|
1216 |
|
---|
1217 |
|
---|
1218 |
for ($time=0; $time<=$ |
---|
1219 |
|
---|
1220 |
|
---|
1221 |
if ($ |
---|
1222 |
|
---|
1223 |
if ($cmp1[0]=~/\ /) { |
---|
1224 |
@temp_cmp = split(' ',$cmp1[0]); |
---|
1225 |
@cmp1 = @temp_cmp; |
---|
1226 |
$comp_off[$time] = &volts2comp($cmp1[$time]) + $offset; |
---|
1227 |
} |
---|
1228 |
|
---|
1229 |
else { |
---|
1230 |
$comp_off[$time] = &volts2comp($cmp1[0]) + $offset; |
---|
1231 |
} |
---|
1232 |
} |
---|
1233 |
else{ |
---|
1234 |
|
---|
1235 |
$comp_off[$time] = &volts2comp($cmp1[$time]) + $offset; |
---|
1236 |
} |
---|
1237 |
|
---|
1238 |
|
---|
1239 |
|
---|
1240 |
($spda[$time],$dira[$time]) = &volts2spddir($spda[$time],$dira[$time]); |
---|
1241 |
($spdb[$time],$dirb[$time]) = &volts2spddir($spdb[$time],$dirb[$time]); |
---|
1242 |
|
---|
1243 |
|
---|
1244 |
|
---|
1245 |
|
---|
1246 |
$bada = ($spda[$time]>100 | $spda[$time]==$fill_value | |
---|
1247 |
$dira[$time]==$fill_value); |
---|
1248 |
if ($bada) { |
---|
1249 |
$spda[$time] = $dira[$time] = $u_a[$time] = $v_a[$time] = $fill_value; |
---|
1250 |
} |
---|
1251 |
else { |
---|
1252 |
|
---|
1253 |
$dira[$time] = &modulo($dira[$time] + $comp_off[$time], 360); |
---|
1254 |
|
---|
1255 |
($u_a[$time],$v_a[$time]) = &spddir2uv($spda[$time],$dira[$time]); |
---|
1256 |
} |
---|
1257 |
|
---|
1258 |
|
---|
1259 |
unless ($ |
---|
1260 |
$badb = ($spdb[$time]>100 | $spdb[$time]==$fill_value | |
---|
1261 |
$dirb[$time]==$fill_value); |
---|
1262 |
if ($badb) { |
---|
1263 |
$spdb[$time] = $dirb[$time] = $u_b[$time] = |
---|
1264 |
$v_b[$time] = $fill_value; |
---|
1265 |
} |
---|
1266 |
else { |
---|
1267 |
$dirb[$time] = &modulo($dirb[$time] + $comp_off[$time], 360); |
---|
1268 |
($u_b[$time],$v_b[$time]) = &spddir2uv($spdb[$time],$dirb[$time]); |
---|
1269 |
} |
---|
1270 |
} |
---|
1271 |
} |
---|
1272 |
|
---|
1273 |
|
---|
1274 |
|
---|
1275 |
($u_mean_a,$u_std_a,$u_max_a) = &fill_stats(@u_a); |
---|
1276 |
($v_mean_a,$v_std_a,$v_max_a) = &fill_stats(@v_a); |
---|
1277 |
|
---|
1278 |
(undef,$spd_std_a,$gust_a) = &fill_stats(@spda); |
---|
1279 |
|
---|
1280 |
(undef,$dir_std_a,undef) = &fill_stats(@dira); |
---|
1281 |
|
---|
1282 |
($mean_spd_a,$mean_dir_a) = &uv2spddir($u_mean_a,$v_mean_a); |
---|
1283 |
|
---|
1284 |
($comp_mean,$comp_std,undef) = &fill_stats(@comp_off); |
---|
1285 |
|
---|
1286 |
|
---|
1287 |
unless ($ |
---|
1288 |
($u_mean_b,$u_std_b,$u_max_b) = &fill_stats(@u_b); |
---|
1289 |
($v_mean_b,$v_std_b,$v_max_b) = &fill_stats(@v_b); |
---|
1290 |
|
---|
1291 |
(undef,$spd_std_b,$gust_b) = &fill_stats(@spdb); |
---|
1292 |
|
---|
1293 |
(undef,$dir_std_b,undef) = &fill_stats(@dirb); |
---|
1294 |
|
---|
1295 |
($mean_spd_b,$mean_dir_b) = &uv2spddir($u_mean_b,$v_mean_b); |
---|
1296 |
} |
---|
1297 |
} |
---|
1298 |
|
---|
1299 |
|
---|
1300 |
|
---|
1301 |
sub write_wind { |
---|
1302 |
|
---|
1303 |
unless ($mean_dir_a == $fill_value || |
---|
1304 |
$mean_spd_a == $fill_value) { |
---|
1305 |
|
---|
1306 |
${$Vars{WSPDA}}[$index] = $mean_spd_a; |
---|
1307 |
${$Vars{WDIRA}}[$index] = $mean_dir_a; |
---|
1308 |
${$Vars{WEASTA}}[$index] = $u_mean_a; |
---|
1309 |
${$Vars{WNORTHA}}[$index] = $v_mean_a; |
---|
1310 |
${$Vars{WSPDSTDA}}[$index] = $spd_std_a; |
---|
1311 |
${$Vars{WDIRSTDA}}[$index] = $dir_std_a; |
---|
1312 |
${$Vars{WGSTA}}[$index] = $gust_a; |
---|
1313 |
} |
---|
1314 |
else { |
---|
1315 |
${$Vars{WSPDA}}[$index] = $fill_value; |
---|
1316 |
${$Vars{WDIRA}}[$index] = $fill_value; |
---|
1317 |
${$Vars{WEASTA}}[$index] = $fill_value; |
---|
1318 |
${$Vars{WNORTHA}}[$index] = $fill_value; |
---|
1319 |
${$Vars{WSPDSTDA}}[$index] = $fill_value; |
---|
1320 |
${$Vars{WDIRSTDA}}[$index] = $fill_value; |
---|
1321 |
${$Vars{WGSTA}}[$index] = $fill_value; |
---|
1322 |
} |
---|
1323 |
unless ($mean_dir_b == $fill_value || |
---|
1324 |
$mean_spd_b == $fill_value) { |
---|
1325 |
|
---|
1326 |
${$Vars{WSPDB}}[$index] = $mean_spd_b; |
---|
1327 |
${$Vars{WDIRB}}[$index] = $mean_dir_b; |
---|
1328 |
${$Vars{WEASTB}}[$index] = $u_mean_b; |
---|
1329 |
${$Vars{WNORTHB}}[$index] = $v_mean_b; |
---|
1330 |
${$Vars{WSPDSTDB}}[$index] = $spd_std_b; |
---|
1331 |
${$Vars{WDIRSTDB}}[$index] = $dir_std_b; |
---|
1332 |
${$Vars{WGSTB}}[$index] = $gust_b; |
---|
1333 |
} |
---|
1334 |
else { |
---|
1335 |
${$Vars{WSPDB}}[$index] = $fill_value; |
---|
1336 |
${$Vars{WDIRB}}[$index] = $fill_value; |
---|
1337 |
${$Vars{WEASTB}}[$index] = $fill_value; |
---|
1338 |
${$Vars{WNORTHB}}[$index] = $fill_value; |
---|
1339 |
${$Vars{WSPDSTDB}}[$index] = $fill_value; |
---|
1340 |
${$Vars{WDIRSTDB}}[$index] = $fill_value; |
---|
1341 |
${$Vars{WGSTB}}[$index] = $fill_value; |
---|
1342 |
} |
---|
1343 |
|
---|
1344 |
unless ($comp_std == 0 || |
---|
1345 |
$comp_mean == -9999) { |
---|
1346 |
${$Vars{COMP1}}[$index] = $comp_mean; |
---|
1347 |
} |
---|
1348 |
else { |
---|
1349 |
${$Vars{COMP1}}[$index] = $fill_value; |
---|
1350 |
} |
---|
1351 |
} |
---|
1352 |
|
---|
1353 |
|
---|
1354 |
|
---|
1355 |
sub handle_comp { |
---|
1356 |
|
---|
1357 |
|
---|
1358 |
|
---|
1359 |
|
---|
1360 |
|
---|
1361 |
|
---|
1362 |
|
---|
1363 |
|
---|
1364 |
our($compd_mean,$compd_std,$compd_max,$pitch_mean,$pitch_std, |
---|
1365 |
$pitch_max,$roll_mean,$roll_std,$roll_max) = $fill_value; |
---|
1366 |
|
---|
1367 |
my(@u_c,@v_c,@u_p,@v_p,@u_r,@v_r) = (); |
---|
1368 |
|
---|
1369 |
|
---|
1370 |
if ($ |
---|
1371 |
|
---|
1372 |
|
---|
1373 |
|
---|
1374 |
for ($time=0; $time<=$ |
---|
1375 |
|
---|
1376 |
|
---|
1377 |
|
---|
1378 |
|
---|
1379 |
|
---|
1380 |
|
---|
1381 |
$comp[$time] = &modulo($comp[$time], 360); |
---|
1382 |
$pitch[$time] = &modulo($pitch[$time], 360); |
---|
1383 |
$roll[$time] = &modulo($roll[$time], 360); |
---|
1384 |
|
---|
1385 |
|
---|
1386 |
|
---|
1387 |
|
---|
1388 |
($u_c[$time],$v_c[$time]) = &spddir2uv(1,$comp[$time]); |
---|
1389 |
($u_p[$time],$v_p[$time]) = &spddir2uv(1,$pitch[$time]); |
---|
1390 |
($u_r[$time],$v_r[$time]) = &spddir2uv(1,$roll[$time]); |
---|
1391 |
|
---|
1392 |
|
---|
1393 |
|
---|
1394 |
} |
---|
1395 |
|
---|
1396 |
|
---|
1397 |
($u_c_mean,$u_c_std,$u_c_max) = &fill_stats(@u_c); |
---|
1398 |
($u_p_mean,$u_p_std,$u_p_max) = &fill_stats(@u_p); |
---|
1399 |
($u_r_mean,$u_r_std,$u_r_max) = &fill_stats(@u_r); |
---|
1400 |
|
---|
1401 |
($v_c_mean,$v_c_std,$v_c_max) = &fill_stats(@v_c); |
---|
1402 |
($v_p_mean,$v_p_std,$v_p_max) = &fill_stats(@v_p); |
---|
1403 |
($v_r_mean,$v_r_std,$v_r_max) = &fill_stats(@v_r); |
---|
1404 |
|
---|
1405 |
|
---|
1406 |
(undef,$compd_mean) = &uv2spddir($u_c_mean,$v_c_mean); |
---|
1407 |
(undef,$compd_std) = ($u_c_std**2+$v_c_std**2)**(1/2); |
---|
1408 |
(undef,$compd_max) = &uv2spddir($u_c_max,$v_c_max); |
---|
1409 |
|
---|
1410 |
(undef,$pitch_mean) = &uv2spddir($u_p_mean,$v_p_mean); |
---|
1411 |
(undef,$pitch_std) = ($u_p_std**2+$v_p_std**2)**(1/2); |
---|
1412 |
(undef,$pitch_max) = &uv2spddir($u_p_max,$v_p_max); |
---|
1413 |
|
---|
1414 |
(undef,$roll_mean) = &uv2spddir($u_r_mean,$v_r_mean); |
---|
1415 |
(undef,$roll_std) = ($u_r_std**2+$v_r_std**2)**(1/2); |
---|
1416 |
(undef,$roll_max) = &uv2spddir($u_r_max,$v_r_max); |
---|
1417 |
|
---|
1418 |
|
---|
1419 |
|
---|
1420 |
|
---|
1421 |
|
---|
1422 |
|
---|
1423 |
${$Vars{COMPMEAN}}[$file_index] = $compd_mean; |
---|
1424 |
${$Vars{COMPSTD}}[$file_index] = $compd_std; |
---|
1425 |
${$Vars{COMPMAX}}[$file_index] = $compd_max; |
---|
1426 |
${$Vars{PITCHMEAN}}[$file_index] = $pitch_mean; |
---|
1427 |
${$Vars{PITCHSTD}}[$file_index] = $pitch_std; |
---|
1428 |
${$Vars{PITCHMAX}}[$file_index] = $pitch_max; |
---|
1429 |
${$Vars{ROLLMEAN}}[$file_index] = $roll_mean; |
---|
1430 |
${$Vars{ROLLSTD}}[$file_index] = $roll_std; |
---|
1431 |
${$Vars{ROLLMAX}}[$file_index] = $roll_max; |
---|
1432 |
} |
---|
1433 |
else { |
---|
1434 |
|
---|
1435 |
${$Vars{COMPMEAN}}[$file_index] = $fill_value; |
---|
1436 |
${$Vars{COMPSTD}}[$file_index] = $fill_value; |
---|
1437 |
${$Vars{COMPMAX}}[$file_index] = $fill_value; |
---|
1438 |
${$Vars{PITCHMEAN}}[$file_index] = $fill_value; |
---|
1439 |
${$Vars{PITCHSTD}}[$file_index] = $fill_value; |
---|
1440 |
${$Vars{PITCHMAX}}[$file_index] = $fill_value; |
---|
1441 |
${$Vars{ROLLMEAN}}[$file_index] = $fill_value; |
---|
1442 |
${$Vars{ROLLSTD}}[$file_index] = $fill_value; |
---|
1443 |
${$Vars{ROLLMAX}}[$file_index] = $fill_value; |
---|
1444 |
} |
---|
1445 |
} |
---|
1446 |
|
---|
1447 |
|
---|
1448 |
|
---|
1449 |
|
---|
1450 |
|
---|
1451 |
sub volts2spddir { |
---|
1452 |
|
---|
1453 |
my($volts_spd,$volts_dir) = @_; |
---|
1454 |
my($wind_spd,$wind_dir) = undef; |
---|
1455 |
|
---|
1456 |
unless($volts_spd=~/$fill_value/) { |
---|
1457 |
|
---|
1458 |
$wind_spd = $volts_spd*60; |
---|
1459 |
} |
---|
1460 |
else { $wind_spd = $fill_value; } |
---|
1461 |
|
---|
1462 |
unless($volts_dir=~/$fill_value/) { |
---|
1463 |
|
---|
1464 |
$wind_dir = $volts_dir*360; |
---|
1465 |
} |
---|
1466 |
else { $wind_dir = $fill_value; } |
---|
1467 |
|
---|
1468 |
return ($wind_spd,$wind_dir); |
---|
1469 |
} |
---|
1470 |
|
---|
1471 |
|
---|
1472 |
|
---|
1473 |
sub volts2comp { |
---|
1474 |
|
---|
1475 |
my($volts_comp) = @_; |
---|
1476 |
|
---|
1477 |
unless($volts_comp=~/$fill_value/) { |
---|
1478 |
|
---|
1479 |
return $volts_comp*360; |
---|
1480 |
} |
---|
1481 |
else { return $fill_value; } |
---|
1482 |
} |
---|
1483 |
|
---|
1484 |
|
---|
1485 |
|
---|
1486 |
|
---|
1487 |
|
---|
1488 |
sub spddir2uv { |
---|
1489 |
|
---|
1490 |
my($wind_spd,$wind_dir) = @_; |
---|
1491 |
my($deg2rad) = .0174533; |
---|
1492 |
my($wind_u,$wind_v) = undef; |
---|
1493 |
|
---|
1494 |
|
---|
1495 |
unless($wind_spd=~/$fill_value/ || $wind_dir=~/$fill_value/) { |
---|
1496 |
|
---|
1497 |
$wind_v = ($wind_spd)*cos($wind_dir*$deg2rad); |
---|
1498 |
$wind_u = ($wind_spd)*sin($wind_dir*$deg2rad); |
---|
1499 |
return ($wind_u,$wind_v); |
---|
1500 |
} |
---|
1501 |
else { |
---|
1502 |
|
---|
1503 |
return ($fill_value,$fill_value); |
---|
1504 |
} |
---|
1505 |
} |
---|
1506 |
|
---|
1507 |
|
---|
1508 |
|
---|
1509 |
sub uv2spddir { |
---|
1510 |
|
---|
1511 |
my($wind_u,$wind_v) = @_; |
---|
1512 |
my($deg2rad) = .0174533; |
---|
1513 |
my($wind_spd,$wind_dir) = undef; |
---|
1514 |
|
---|
1515 |
unless($wind_u=~/$fill_value/ || $wind_v=~/$fill_value/) { |
---|
1516 |
$wind_spd = ($wind_u**2 + $wind_v**2)**(1/2); |
---|
1517 |
$wind_dir = &modulo(atan2($wind_u, $wind_v)/$deg2rad,360); |
---|
1518 |
return ($wind_spd,$wind_dir); |
---|
1519 |
} |
---|
1520 |
else { |
---|
1521 |
return ($fill_value,$fill_value); |
---|
1522 |
} |
---|
1523 |
} |
---|
1524 |
|
---|
1525 |
|
---|
1526 |
|
---|
1527 |
sub avg_rain { |
---|
1528 |
|
---|
1529 |
|
---|
1530 |
|
---|
1531 |
|
---|
1532 |
|
---|
1533 |
|
---|
1534 |
|
---|
1535 |
|
---|
1536 |
|
---|
1537 |
|
---|
1538 |
|
---|
1539 |
|
---|
1540 |
|
---|
1541 |
|
---|
1542 |
|
---|
1543 |
|
---|
1544 |
|
---|
1545 |
|
---|
1546 |
|
---|
1547 |
|
---|
1548 |
|
---|
1549 |
|
---|
1550 |
|
---|
1551 |
|
---|
1552 |
|
---|
1553 |
my($maxrain); |
---|
1554 |
$maxrain = &max(@cumrain); |
---|
1555 |
|
---|
1556 |
|
---|
1557 |
unless ($maxrain=~/$fill_value/) { |
---|
1558 |
$maxrain = $maxrain*10; |
---|
1559 |
} |
---|
1560 |
else { |
---|
1561 |
$maxrain = $fill_value; |
---|
1562 |
} |
---|
1563 |
${$Vars{RAIN}}[$index] = $maxrain; |
---|
1564 |
} |
---|
1565 |
|
---|
1566 |
|
---|
1567 |
|
---|
1568 |
|
---|
1569 |
|
---|
1570 |
|
---|
1571 |
sub fill_stats(@array) { |
---|
1572 |
|
---|
1573 |
my(@array,$mean_value); |
---|
1574 |
@array = @_; |
---|
1575 |
|
---|
1576 |
$array = join(',',@array); |
---|
1577 |
$array =~s/$fill_value//g; |
---|
1578 |
$array =~s/\,\,/\,/g; |
---|
1579 |
@array = split(',',$array); |
---|
1580 |
|
---|
1581 |
unless ($ |
---|
1582 |
($mean) = &mean(@array); |
---|
1583 |
($std) = &std_dev(@array); |
---|
1584 |
($max) = &max(@array); |
---|
1585 |
return ($mean,$std,$max); |
---|
1586 |
} |
---|
1587 |
else { |
---|
1588 |
return ($fill_value,$fill_value,$fill_value); |
---|
1589 |
} |
---|
1590 |
} |
---|
1591 |
|
---|
1592 |
|
---|
1593 |
|
---|
1594 |
sub mean(@array) { |
---|
1595 |
my(@array,$mean_value); |
---|
1596 |
@array = @_; |
---|
1597 |
|
---|
1598 |
$sum_array = join('+',@array); |
---|
1599 |
|
---|
1600 |
if ($sum_array =~m/\+\+/) { |
---|
1601 |
return $fill_value; |
---|
1602 |
} |
---|
1603 |
|
---|
1604 |
|
---|
1605 |
eval('$mean_value = ('."$sum_array".')/($#array+1);'); |
---|
1606 |
return $mean_value; |
---|
1607 |
} |
---|
1608 |
|
---|
1609 |
|
---|
1610 |
|
---|
1611 |
sub std_dev(@array) { |
---|
1612 |
my(@array,$std_dev); |
---|
1613 |
@array = @_; |
---|
1614 |
|
---|
1615 |
my $temp_std = new Statistics::Basic::StdDev(\@array); |
---|
1616 |
$std_dev = $temp_std -> query; |
---|
1617 |
|
---|
1618 |
return $std_dev; |
---|
1619 |
} |
---|
1620 |
|
---|
1621 |
|
---|
1622 |
|
---|
1623 |
sub max { |
---|
1624 |
|
---|
1625 |
my $max = shift(@_); |
---|
1626 |
foreach $foo (@_) { |
---|
1627 |
$max = $foo if $max < $foo; |
---|
1628 |
} |
---|
1629 |
return $max; |
---|
1630 |
} |
---|
1631 |
|
---|
1632 |
|
---|
1633 |
|
---|
1634 |
sub modulo($element,$modulus) { |
---|
1635 |
|
---|
1636 |
|
---|
1637 |
my($element,$modulus,$el,$floor) = @_; |
---|
1638 |
|
---|
1639 |
if ($element < 0) { |
---|
1640 |
return $element + $modulus; |
---|
1641 |
} |
---|
1642 |
elsif ($element <= $modulus) { |
---|
1643 |
return $element; |
---|
1644 |
} |
---|
1645 |
elsif ($element > 3*$modulus) { |
---|
1646 |
|
---|
1647 |
return $fill_value; |
---|
1648 |
} |
---|
1649 |
elsif ($element > 2*$modulus) { |
---|
1650 |
return $element - 2*$modulus; |
---|
1651 |
} |
---|
1652 |
elsif ($element > $modulus) { |
---|
1653 |
return $element - $modulus; |
---|
1654 |
} |
---|
1655 |
} |
---|
1656 |
|
---|
1657 |
|
---|
1658 |
|
---|
1659 |
|
---|