Changeset 222 for altifloat/exp_nemed/plot_article.m
- Timestamp:
- 03/17/16 11:52:56 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
altifloat/exp_nemed/plot_article.m
r191 r222 13 13 robs=load([ indir 'obs.dat']); 14 14 rsim=load([ indir 'rfloat_nodiv.dat']); 15 %rsim=load([ indir 'rfloat_div2.dat']); 15 16 uvb=load([ indir 'uv_back_aviso.dat']); 16 17 17 18 uvr=load([indir 'uv_tot_nodiv.dat']); 18 19 %uvr=load([indir 'uv_tot_div2.dat']); 19 20 meshg=load([indir 'meshgrid_aviso.dat']); 20 21 mask = load([indir 'mask_aviso.dat']); … … 125 126 %Ajout des trajectoires réelles 126 127 axis(LIM); 127 mapshow(S,'Color',[0 0 0]); 128 mapshow(S,'Color',[0 0 0]); 129 xlabel('Degrees East') 130 ylabel('Degrees North') 128 131 129 132 legend([Hb Hr Hobs],'Averaged over 9 days : background','Corrected','Drifter trajectory','Location','SouthWest'); … … 175 178 owr = okubo_weiss(Uxr,Uyr,Vxr,Vyr); 176 179 title('Okubo-Weiss parameter of the Background'); 180 xlabel('Degrees East') 181 ylabel('Degrees North') 182 177 183 if osave 178 184 print('-dpng','-r300',[outdir 'okubo_weiss_aviso.png']); … … 192 198 colormap(cmap); 193 199 title('Okubo-Weiss parameter of the corrected field'); 200 xlabel('Degrees East') 201 ylabel('Degrees North') 194 202 195 203 if osave
Note: See TracChangeset
for help on using the changeset viewer.