0
|
1 |
%======================================================================
|
|
2 |
% U H _ L A D C P 2 A N T S . M
|
|
3 |
% doc: Sun Jan 22 15:19:00 2006
|
5
|
4 |
% dlm: Tue Feb 21 13:59:44 2012
|
0
|
5 |
% (c) 2006 A.M. Thurnherr
|
5
|
6 |
% uE-Info: 22 75 NIL 0 0 72 2 2 4 NIL ofnI
|
|
7 |
% uE-Info: 55 0 NIL 0 0 72 0 2 4 NIL ofnI
|
0
|
8 |
%======================================================================
|
|
9 |
%
|
|
10 |
% export LDEO LADCP output to ANTS file
|
|
11 |
%
|
3
|
12 |
% USAGE: UH_LADCP2ANTS(inFile,outBaseName)
|
0
|
13 |
%
|
|
14 |
|
|
15 |
% HISTORY:
|
3
|
16 |
% Jan 22, 2006: - created
|
|
17 |
% Nov 2, 2008: - BUG: v_var had not been exported because of typo
|
|
18 |
% Oct 12, 2009: - adapted to new struct2ANTS
|
|
19 |
% Aug 4, 2011: - added usage message on zero args
|
|
20 |
% - updated output field names
|
5
|
21 |
% Feb 21, 2012: - manually merged different versions
|
|
22 |
% - updated all output field names (dc/uc instead of dn/un)
|
0
|
23 |
|
3
|
24 |
function [] = UH_LADCP2ANTS(ifn,obn)
|
|
25 |
|
|
26 |
if nargin ~= 2
|
|
27 |
help UH_LADCP2ANTS
|
|
28 |
return
|
|
29 |
end
|
0
|
30 |
|
|
31 |
eval(sprintf('load %s',ifn));
|
|
32 |
|
3
|
33 |
prof.dayno = mean(txy_start_end(:,1));
|
0
|
34 |
|
|
35 |
prof.lon = mean(txy_start_end(find(isfinite(txy_start_end(:,2))),2));
|
|
36 |
prof.lat = mean(txy_start_end(find(isfinite(txy_start_end(:,3))),3));
|
|
37 |
|
|
38 |
good = find(sm_mn_i > 0);
|
|
39 |
prof.depth = d_samp(good);
|
|
40 |
prof.max_depth = max(prof.depth);
|
|
41 |
|
|
42 |
dn_bad = find(sm_dn_i == 0);
|
|
43 |
up_bad = find(sm_up_i == 0);
|
|
44 |
|
|
45 |
prof.u = su_mn_i(good);
|
3
|
46 |
su_dn_i(dn_bad) = NaN; prof.dc_u = su_dn_i(good);
|
|
47 |
su_up_i(up_bad) = NaN; prof.uc_u = su_up_i(good);
|
0
|
48 |
prof.v = sv_mn_i(good);
|
3
|
49 |
sv_dn_i(dn_bad) = NaN; prof.dc_v = sv_dn_i(good);
|
|
50 |
sv_up_i(up_bad) = NaN; prof.uc_v = sv_up_i(good);
|
5
|
51 |
prof.w = sw_mn_i(good);
|
|
52 |
sw_dn_i(dn_bad) = NaN; prof.dc_w = sw_dn_i(good);
|
|
53 |
sw_up_i(up_bad) = NaN; prof.uc_w = sw_up_i(good);
|
0
|
54 |
|
|
55 |
prof.u_var = su_var_mn_i(good);
|
3
|
56 |
su_var_dn_i(dn_bad) = NaN; prof.dc_u_var = su_var_dn_i(good);
|
|
57 |
su_var_up_i(up_bad) = NaN; prof.uc_u_var = su_var_up_i(good);
|
5
|
58 |
prof.v_var = sv_var_mn_i(good);
|
|
59 |
sv_var_dn_i(dn_bad) = NaN; prof.dc_v_var = sv_var_dn_i(good);
|
|
60 |
sv_var_up_i(up_bad) = NaN; prof.uc_v_var = sv_var_up_i(good);
|
|
61 |
prof.w_var = sw_var_mn_i(good);
|
|
62 |
sw_var_dn_i(dn_bad) = NaN; prof.dc_w_var = sw_var_dn_i(good);
|
|
63 |
sw_var_up_i(up_bad) = NaN; prof.uc_w_var = sw_var_up_i(good);
|
0
|
64 |
|
3
|
65 |
prof.nsamp = sn_mn_i(good);
|
|
66 |
sn_dn_i(dn_bad) = NaN; prof.dc_nsamp = sn_dn_i(good);
|
|
67 |
sn_up_i(up_bad) = NaN; prof.uc_nsamp = sn_up_i(good);
|
0
|
68 |
|
1
|
69 |
struct2ANTS(prof,ifn,sprintf('%s.prof',obn));
|