0
|
1 |
#======================================================================
|
|
2 |
# L I B W O C E . P L
|
|
3 |
# doc: Mon Aug 28 11:07:47 2000
|
|
4 |
# dlm: Sun Jul 2 00:16:46 2006
|
|
5 |
# (c) 2000 A.M. Thurnherr
|
|
6 |
# uE-Info: 53 35 NIL 0 0 72 2 2 4 NIL ofnI
|
|
7 |
#======================================================================
|
|
8 |
|
|
9 |
# HISTORY:
|
|
10 |
# Sep 04, 2000: - created from [libNODC.pl]
|
|
11 |
# Sep 19, 2000: - added &q_OK()
|
|
12 |
# Jan 2, 2002: - added optional length argument to &q_OK()
|
|
13 |
# - allowed NaN observations in &q_OK()
|
|
14 |
# Aug 14, 2002: - added &csv_q_OK()
|
|
15 |
# Jun 19, 2004: - made Y2K compatible
|
|
16 |
# Apr 4, 2006: - added flt_qual(), flt_src()
|
|
17 |
# Jul 1, 2006: - Version 3.3 [HISTORY]
|
|
18 |
|
|
19 |
require "$ANTS/libconv.pl"; # imply lat/lon conversions
|
|
20 |
|
|
21 |
#----------------------------------------------------------------------
|
|
22 |
# date/time
|
|
23 |
#----------------------------------------------------------------------
|
|
24 |
|
|
25 |
sub MMDDYY(@) # 6 digit date
|
|
26 |
{
|
|
27 |
my($MMDDYY) = &antsFunUsage(1,"6","MMDDYY",@_);
|
|
28 |
my($YY) = substr($MMDDYY,4,2);
|
|
29 |
return substr($MMDDYY,0,2) . "/" .
|
|
30 |
substr($MMDDYY,2,2) . ($YY>50 ? "/19$YY" : "/20$YY");
|
|
31 |
}
|
|
32 |
|
|
33 |
sub HHMM(@) # 4 digits
|
|
34 |
{
|
|
35 |
my($HHMM) = &antsFunUsage(1,"4","HHMM",@_);
|
|
36 |
return substr($HHMM,0,2) . ":" . substr($HHMM,2,2);
|
|
37 |
}
|
|
38 |
|
|
39 |
#----------------------------------------------------------------------
|
|
40 |
# CTD quality flags
|
|
41 |
#----------------------------------------------------------------------
|
|
42 |
|
|
43 |
sub q_OK(@) # accept only perfect measurements (table 4.10)
|
|
44 |
{
|
|
45 |
my($nflags) = (@_ == 3) ? 4 : $_[3];
|
|
46 |
my($obs,$qf,$qi) =
|
|
47 |
&antsFunUsage(-3,".${nflags}c","obs, quality flags, flag index[, nflags]",@_);
|
|
48 |
return (substr($qf,$qi,1) == 2) ? $obs : nan;
|
|
49 |
}
|
|
50 |
|
|
51 |
sub csv_q_OK(@) # exchange-format version (single flags)
|
|
52 |
{
|
|
53 |
return $_[1] == 2 ? $_[0] : nan;
|
|
54 |
}
|
|
55 |
|
|
56 |
#----------------------------------------------------------------------
|
|
57 |
# Float Quality/Source Flags [/Data/Floats/DBE/WFDAC/WDBE/quality.doc]
|
|
58 |
# - 2-digit decimal number
|
|
59 |
# - 10s: quality (0-9, with 9 being best)
|
|
60 |
# - 1s: source:
|
|
61 |
# 0 missing
|
|
62 |
# 1 interpolated (backward diff for u/v)
|
|
63 |
# 2 (forward diff for u/v)
|
|
64 |
# 3 splined
|
|
65 |
# 4 manually edited
|
|
66 |
# 5 filtered/averaged
|
|
67 |
# 9 original value
|
|
68 |
# - should accept source values >= 4
|
|
69 |
#----------------------------------------------------------------------
|
|
70 |
|
|
71 |
sub flt_qual($)
|
|
72 |
{ return int($_[0] / 10); }
|
|
73 |
|
|
74 |
sub flt_src($)
|
|
75 |
{ return $_[0] % 10; }
|
|
76 |
|
|
77 |
1;
|