1 %====================================================================== |
1 %====================================================================== |
2 % L O A D S A D C P . M |
2 % L O A D S A D C P . M |
3 % doc: Sun Jun 27 23:42:04 2004 |
3 % doc: Sun Jun 27 23:42:04 2004 |
4 % dlm: Fri Mar 5 15:49:11 2010 |
4 % dlm: Tue May 22 11:01:21 2012 |
5 % (c) 2004 ladcp@ |
5 % (c) 2004 ladcp@ |
6 % uE-Info: 133 14 NIL 0 0 72 0 2 8 NIL ofnI |
6 % uE-Info: 109 4 NIL 0 0 72 0 2 8 NIL ofnI |
7 %====================================================================== |
7 %====================================================================== |
8 |
8 |
9 % CHANGES BY ANT: |
9 % CHANGES BY ANT: |
10 % Jun 30, 2004: - added Figure 9 title |
10 % Jun 30, 2004: - added Figure 9 title |
11 % Jul 9, 2004: - made lack of SADCP data in time window into real warning |
11 % Jul 9, 2004: - made lack of SADCP data in time window into real warning |
12 % Jan 7, 2009: - tightened use of exist() |
12 % Jan 7, 2009: - tightened use of exist() |
13 % Mar 4, 2010: - changed default of p.sadcp_dtok from 5min to zero |
13 % Mar 4, 2010: - changed default of p.sadcp_dtok from 5min to zero |
|
14 % May 22, 2012: - removed code that took GPS information from SADCP data, |
|
15 % because these data are unlikely to be accurate enough |
|
16 % for the ship-drift constraint; if they are, the user |
|
17 % should verify and make a GPS file during pre-processing |
14 |
18 |
15 function [di,p]=loadsadcp(f,di,p) |
19 function [di,p]=loadsadcp(f,di,p) |
16 % function [di,p]=loadsadcp(f,di,p) |
20 % function [di,p]=loadsadcp(f,di,p) |
17 % |
21 % |
18 % |
22 % |
76 xlabel('longitude') |
80 xlabel('longitude') |
77 ylabel('latitude') |
81 ylabel('latitude') |
78 streamer([p.name,' Figure 9']); |
82 streamer([p.name,' Figure 9']); |
79 pause(0.001) |
83 pause(0.001) |
80 return; |
84 return; |
81 end |
85 end % if no SADCP data in LADCP region |
82 |
86 |
83 % interpolate SADCP navigation to LADCP time series |
87 % interpolate SADCP navigation to LADCP time series |
84 di.sadcp_lon=interp1(tim_sadcp,lon_sadcp(:),di.time_jul); |
88 di.sadcp_lon=interp1(tim_sadcp,lon_sadcp(:),di.time_jul); |
85 di.sadcp_lat=interp1(tim_sadcp,lat_sadcp(:),di.time_jul); |
89 di.sadcp_lat=interp1(tim_sadcp,lat_sadcp(:),di.time_jul); |
86 |
90 |
87 % set position from SADCP nav |
91 % set position from SADCP nav |
88 if abs(p.lon)+abs(p.lat)==0 |
92 if abs(p.lon)+abs(p.lat)==0 |
|
93 error('as of version IX_9, using GPS info from SADCP data stream is no longer supported'); |
89 slat=di.sadcp_lat(1); |
94 slat=di.sadcp_lat(1); |
90 slon=di.sadcp_lon(1); |
95 slon=di.sadcp_lon(1); |
91 elat=di.sadcp_lat(end); |
96 elat=di.sadcp_lat(end); |
92 elon=di.sadcp_lon(end); |
97 elon=di.sadcp_lon(end); |
93 p.poss=[fix(slat), (slat-fix(slat))*60, fix(slon), (slon-fix(slon))*60]; |
98 p.poss=[fix(slat), (slat-fix(slat))*60, fix(slon), (slon-fix(slon))*60]; |
94 p.pose=[fix(elat), (slat-fix(elat))*60, fix(elon), (slon-fix(elon))*60]; |
99 p.pose=[fix(elat), (slat-fix(elat))*60, fix(elon), (slon-fix(elon))*60]; |
95 end |
100 end |
96 |
101 |
97 % if no other ship navigation exists, use SADCP navigation |
102 % if no other ship navigation exists, use SADCP navigation |
98 if existf(di,'slon')==0 |
103 if existf(di,'slon')==0 |
|
104 error('as of version IX_9, using GPS info from SADCP data stream is no longer supported'); |
99 di.slon=di.sadcp_lon; |
105 di.slon=di.sadcp_lon; |
100 di.slat=di.sadcp_lat; |
106 di.slat=di.sadcp_lat; |
101 else |
107 else |
102 if sum(isfinite(di.slon+di.slat))==0 |
108 if sum(isfinite(di.slon+di.slat))==0 |
|
109 error('as of version IX_9, using GPS info from SADCP data stream is no longer supported'); |
103 di.slon=di.sadcp_lon; |
110 di.slon=di.sadcp_lon; |
104 di.slat=di.sadcp_lat; |
111 di.slat=di.sadcp_lat; |
105 end |
112 end |
106 end |
113 end |
107 |
114 |