1 | function [ output_args ] = plotongraph(nboffloater, jtlagt, totaltime) |
---|
2 | %%%%%%%%%% |
---|
3 | |
---|
4 | |
---|
5 | |
---|
6 | close all |
---|
7 | |
---|
8 | figure(1) |
---|
9 | |
---|
10 | fid = fopen('../obs_float/obs_float_test.dat', 'rb'); |
---|
11 | fseek(fid, 0, 'eof'); |
---|
12 | fileSize = ftell(fid); |
---|
13 | frewind(fid); |
---|
14 | data = fread(fid, fileSize, 'uint8'); |
---|
15 | numLines = sum(data == 10) + 1; |
---|
16 | fclose(fid); |
---|
17 | m=numLines-1; |
---|
18 | n=m/nboffloater; |
---|
19 | |
---|
20 | |
---|
21 | load ../obs_float/obs_float_test.dat; |
---|
22 | counterobs=0; |
---|
23 | for g=1:nboffloater |
---|
24 | |
---|
25 | x=obs_float_test (counterobs+1:2:counterobs+n,4); |
---|
26 | y=obs_float_test (counterobs+2:2:counterobs+n,4); |
---|
27 | plot(x(1),y(1),'+') |
---|
28 | hold on; |
---|
29 | plot(x(n/2),y(n/2),'>') |
---|
30 | hold on; |
---|
31 | plot(x(1:n/2),y(1:n/2),'+'); |
---|
32 | hold on; |
---|
33 | plot(x(1:n/2),y(1:n/2),'-r') |
---|
34 | hold on; |
---|
35 | counterobs=counterobs+n; |
---|
36 | end |
---|
37 | |
---|
38 | |
---|
39 | fid = fopen('../obs_float/rfloat_total.dat', 'rb'); |
---|
40 | fseek(fid, 0, 'eof'); |
---|
41 | fileSize = ftell(fid); |
---|
42 | frewind(fid); |
---|
43 | data = fread(fid, fileSize, 'uint8'); |
---|
44 | numLines = sum(data == 10) + 1; |
---|
45 | fclose(fid); |
---|
46 | p=numLines-1; |
---|
47 | u=p/nboffloater; |
---|
48 | |
---|
49 | countercorr=0; |
---|
50 | |
---|
51 | |
---|
52 | load ../obs_float/rfloat_total.dat |
---|
53 | ligneobs=0; |
---|
54 | for d=1:nboffloater |
---|
55 | x2=rfloat_total(countercorr+1:countercorr+u,4); |
---|
56 | y2=rfloat_total(countercorr+1:countercorr+u,3); |
---|
57 | ligne=ligneobs+1; |
---|
58 | for h=1:jtlagt:u |
---|
59 | x2(h)=obs_float_test(ligne,4); |
---|
60 | y2(h)=obs_float_test(ligne+1,4); |
---|
61 | ligne=ligne+2; |
---|
62 | end |
---|
63 | [k,h]=size(x2); |
---|
64 | plot(x2(1),y2(1),'+') |
---|
65 | hold on; |
---|
66 | plot(x2(k),y2(k),'>') |
---|
67 | hold on; |
---|
68 | plot(x2(1:k),y2(1:k),'+') |
---|
69 | hold on; |
---|
70 | plot(x2(1:k),y2(1:k),'-g') |
---|
71 | hold on; |
---|
72 | ligneobs=ligneobs+n; |
---|
73 | countercorr=countercorr+u; |
---|
74 | end |
---|
75 | |
---|
76 | fid = fopen('../obs_float/third_observation.dat', 'rb'); |
---|
77 | fseek(fid, 0, 'eof'); |
---|
78 | fileSize = ftell(fid); |
---|
79 | frewind(fid); |
---|
80 | data = fread(fid, fileSize, 'uint8'); |
---|
81 | numLines = sum(data == 10) + 1; |
---|
82 | fclose(fid); |
---|
83 | m5=numLines-1; |
---|
84 | n5=m5/2; |
---|
85 | load ../obs_float/third_observation.dat; |
---|
86 | |
---|
87 | x5=third_observation (1:2:m5,4); |
---|
88 | y5=third_observation (2:2:m5,4); |
---|
89 | x=x+1; |
---|
90 | y=y+1; |
---|
91 | plot(x5(1),y5(1),'+') |
---|
92 | hold on; |
---|
93 | plot(x5(n5),y5(n5),'>') |
---|
94 | hold on; |
---|
95 | plot(x5(1:n5),y5(1:n5),'+'); |
---|
96 | hold on; |
---|
97 | plot(x5(1:n5),y5(1:n5),'-m') |
---|
98 | hold on; |
---|
99 | |
---|
100 | %%% To show the same path but forwarded with the bck that is beeing used |
---|
101 | fid = fopen('../obs_float/generated_fwd.dat', 'rb'); |
---|
102 | fseek(fid, 0, 'eof'); |
---|
103 | fileSize = ftell(fid); |
---|
104 | frewind(fid); |
---|
105 | data = fread(fid, fileSize, 'uint8'); |
---|
106 | numLines = sum(data == 10) + 1; |
---|
107 | fclose(fid); |
---|
108 | p3=numLines-1; |
---|
109 | u3=p3/3; |
---|
110 | |
---|
111 | load ../obs_float/generated_fwd.dat |
---|
112 | x11=generated_fwd (1:2:u3,4); |
---|
113 | y11=generated_fwd (2:2:u3,4); |
---|
114 | |
---|
115 | plot(x11(1),y11(1),'+') |
---|
116 | hold on; |
---|
117 | plot(x11(u3/2),y11(u3/2),'>') |
---|
118 | hold on; |
---|
119 | plot(x11(1:u3/2),y11(1:u3/2),'+'); |
---|
120 | hold on; |
---|
121 | plot(x11(1:u3/2),y11(1:u3/2),'-b') |
---|
122 | hold on; |
---|
123 | |
---|
124 | x13=generated_fwd (u3+1:2:2*u3,4); |
---|
125 | y13=generated_fwd (u3+2:2:2*u3,4); |
---|
126 | |
---|
127 | plot(x13(1),y13(1),'+') |
---|
128 | hold on; |
---|
129 | plot(x13(u3/2),y13(u3/2),'>') |
---|
130 | hold on; |
---|
131 | plot(x13(1:u3/2),y13(1:u3/2),'+'); |
---|
132 | hold on; |
---|
133 | plot(x13(1:u3/2),y13(1:u3/2),'-b') |
---|
134 | hold on; |
---|
135 | |
---|
136 | x15=generated_fwd (2*u3+1:2:p3,4); |
---|
137 | y15=generated_fwd (2*u3+2:2:p3,4); |
---|
138 | |
---|
139 | plot(x15(1),y15(1),'+') |
---|
140 | hold on; |
---|
141 | plot(x15(u3/2),y15(u3/2),'>') |
---|
142 | hold on; |
---|
143 | plot(x15(1:u3/2),y15(1:u3/2),'+'); |
---|
144 | hold on; |
---|
145 | plot(x15(1:u3/2),y15(1:u3/2),'-b') |
---|
146 | hold on; |
---|
147 | |
---|
148 | %%%% To show the same path but forwarded with the mean that is beeing used |
---|
149 | fid = fopen('../obs_float/forwardedbymean.dat', 'rb'); |
---|
150 | fseek(fid, 0, 'eof'); |
---|
151 | fileSize = ftell(fid); |
---|
152 | frewind(fid); |
---|
153 | data = fread(fid, fileSize, 'uint8'); |
---|
154 | numLines = sum(data == 10) + 1; |
---|
155 | fclose(fid); |
---|
156 | p4=numLines-1; |
---|
157 | u4=p4/3; |
---|
158 | % |
---|
159 | load ../obs_float/forwardedbymean.dat |
---|
160 | x11=forwardedbymean (1:2:u4,4); |
---|
161 | y11=forwardedbymean (2:2:u4,4); |
---|
162 | % x=x+1; |
---|
163 | % y=y+1; |
---|
164 | plot(x11(1),y11(1),'+') |
---|
165 | hold on; |
---|
166 | plot(x11(u4/2),y11(u4/2),'>') |
---|
167 | hold on; |
---|
168 | plot(x11(1:u4/2),y11(1:u4/2),'+'); |
---|
169 | hold on; |
---|
170 | plot(x11(1:u4/2),y11(1:u4/2),'-r') |
---|
171 | hold on; |
---|
172 | |
---|
173 | x13=forwardedbymean (u4+1:2:2*u4,4); |
---|
174 | y13=forwardedbymean (u4+2:2:2*u4,4); |
---|
175 | % x3=x3+1; |
---|
176 | % y3=y3+1; |
---|
177 | plot(x13(1),y13(1),'+') |
---|
178 | hold on; |
---|
179 | plot(x13(u4/2),y13(u4/2),'>') |
---|
180 | hold on; |
---|
181 | plot(x13(1:u4/2),y13(1:u4/2),'+'); |
---|
182 | hold on; |
---|
183 | plot(x13(1:u4/2),y13(1:u4/2),'-r') |
---|
184 | hold on; |
---|
185 | |
---|
186 | x15=forwardedbymean (2*u4+1:2:p4,4); |
---|
187 | y15=forwardedbymean (2*u4+2:2:p4,4); |
---|
188 | % x3=x3+1; |
---|
189 | % y3=y3+1; |
---|
190 | plot(x15(1),y15(1),'+') |
---|
191 | hold on; |
---|
192 | plot(x15(u4/2),y15(u4/2),'>') |
---|
193 | hold on; |
---|
194 | plot(x15(1:u4/2),y15(1:u4/2),'+'); |
---|
195 | hold on; |
---|
196 | plot(x15(1:u4/2),y15(1:u4/2),'-r') |
---|
197 | hold on; |
---|
198 | |
---|
199 | |
---|
200 | |
---|
201 | load ../obs_float/uv_total.dat; |
---|
202 | whos; |
---|
203 | X1=uv_total(1:5046,end-3); |
---|
204 | Y1=uv_total(1:5046,end-2); |
---|
205 | X1=X1+2; |
---|
206 | Y1=Y1+2; |
---|
207 | %%%%%%change si pas de temps change |
---|
208 | U1=zeros(5046,1); |
---|
209 | V1=zeros(5046,1); |
---|
210 | e=1; |
---|
211 | for r=1:1:totaltime |
---|
212 | |
---|
213 | Unew=zeros(1,5046); |
---|
214 | Vnew=zeros(1,5046); |
---|
215 | for count=1:5046 |
---|
216 | Unew(count)=uv_total(e,end-1); |
---|
217 | Vnew(count)=uv_total(e,end); |
---|
218 | e=e+1; |
---|
219 | end |
---|
220 | U1=U1+Unew'; |
---|
221 | V1=V1+Vnew'; |
---|
222 | |
---|
223 | |
---|
224 | end |
---|
225 | |
---|
226 | U1=U1./totaltime; |
---|
227 | V1=V1./totaltime; |
---|
228 | %%% Because if the jtlag is 3 we have a set of zeros each 3 time steps in |
---|
229 | %%% the uv_total. |
---|
230 | U1=U1; |
---|
231 | V1=V1; |
---|
232 | % U1=roundn(U1,-6); |
---|
233 | % V1=roundn(V1,-6); |
---|
234 | |
---|
235 | quiver(X1,Y1,U1,V1,2,'r','Autoscale','off'); |
---|
236 | |
---|
237 | |
---|
238 | load ../obs_float/vel_bck_all.dat; |
---|
239 | X2=vel_bck_all(1:5046,end-3); |
---|
240 | Y2=vel_bck_all(1:5046,end-2); |
---|
241 | X2=X2+1; |
---|
242 | Y2=Y2+1; |
---|
243 | %change si pas de temps change |
---|
244 | U2=zeros(5046,1); |
---|
245 | V2=zeros(5046,1); |
---|
246 | ebcknew=1; |
---|
247 | for r=1:1:28 |
---|
248 | |
---|
249 | Ubcknew=zeros(1,5046); |
---|
250 | Vnew=zeros(1,5046); |
---|
251 | for count=1:5046 |
---|
252 | Ubcknew(count)=vel_bck_all(ebcknew,end-1); |
---|
253 | Vbcknew(count)=vel_bck_all(ebcknew,end); |
---|
254 | ebcknew=ebcknew+1; |
---|
255 | end |
---|
256 | U2=U2+Ubcknew'; |
---|
257 | V2=V2+Vbcknew'; |
---|
258 | |
---|
259 | |
---|
260 | end |
---|
261 | |
---|
262 | U2=U2./28; |
---|
263 | V2=V2./28; |
---|
264 | |
---|
265 | quiver(X2,Y2,U2,V2,2,'b','Autoscale','off'); |
---|
266 | |
---|
267 | |
---|
268 | end |
---|
269 | |
---|