Context Navigation

script_plot.m @ 160

Last change on this file since 160 was 140, checked in by jbrlod, 9 years ago
modify matlab function to prepare exp model
File size: 548 bytes

Line
1	if ~exist('uvb')
2	uvb=load(['../exp_forw_dan/uv_bck_dan.dat']);
3	end
4	filelonlat='../../data/data/MODEL_CYCO/uveta_cyco1_2013090112.nc';
5	lon=ncread(filelonlat,'lon');
6	lat=ncread(filelonlat,'lat');
7	[Lon,Lat]=meshgrid(lon,lat);
8	it=1:2:25
9	nlon=length(unique(uvb(:,2)));
10	nlat=length(unique(uvb(:,3)));
11	ntime=length(unique(uvb(:,1)));
12	Ub=reshape(uvb(:,end-1),nlat,nlon,ntime);
13	Vb=reshape(uvb(:,end),nlat,nlon,ntime);
14
15	for i=it
16	Hb=quiver(Lon,Lat,squeeze(Ub(:,:,i)),squeeze(Vb(:,:,i)),'b');
17	axis([35.5 36 34.5 35.5])
18	title(['it=' int2str(i)]);
19	pause(0.2)
20	end

Note: See TracBrowser for help on using the repository browser.