Last change
on this file since 140 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.