Changeset 86 for altifloat/obs_external

Timestamp:

04/01/13 19:31:21 (11 years ago)

Author:

leila_ocean

Message:

Friday's afternoon modifications

Location:

altifloat/obs_external

Files:

• read_output.m (modified) (1 diff)
• read_traj.m (modified) (1 diff)

altifloat/obs_external/read_output.m

@@ -1 +1 @@
 close all; clear all;
-load('../obs_float/uzero.dat');
-Ut=reshape(uzero(:,4),58,87);
-Vt=reshape(uzero(:,4),58,87);
+%load('../obs_float/uzero.dat');
+%Ut=reshape(uzero(:,4),58,87);
+%Vt=reshape(uzero(:,4),58,87);
 
 
- load('../obs_float/uoptimtwin.dat');
- Urec=reshape(uoptimtwin(:,5),58,87);
- Vrec=reshape(uoptimtwin(:,5),58,87);
+ load('../obs_float/uoptim7.dat');
+ Urec=reshape(uoptim7(:,5),58,87);
+ %Vrec=reshape(uoptimtwin(:,5),58,87);
  figure;
- imagesc(Urec-Ut); axis xy; colorbar;
+ imagesc(Urec); axis xy; colorbar;
 %imagesc(Ut); axis xy; colorbar;
 

altifloat/obs_external/read_traj.m

@@ -1 +1 @@
 clear all;
-RR=load('../obs_float/r_cout_d_optim_1.dat');
+RR=load('../obs_float/r_float_optim7.dat');
 R1=RR(3:2:end,end);
 R2=RR(4:2:end,end);
+hold on
+%figure;
+plot(R1+1,R2+1,'k-');
+hold on
+plot(R1(1)+1,R2(1)+1,'*r');
 
-figure;
-plot(R1,R2);
-hold
-plot(R1(1),R2(1),'*r');
+obs=load('../obs_float/obs.dat')
+plot(obs(1,end)+1,obs(2,end)+1,'*g');