source: altifloat/matlab_toolbox/interpol_aviso.m @ 199

Last change on this file since 199 was 129, checked in by jbrlod, 10 years ago

last version of Varanth
  • Property svn:executable set to *
File size: 696 bytes
Line 
1aviso_file='../obs_aviso/lebanon/dt_med_allsat_msla_uv_20130101_20140704.nc';
2
3Lon_ref=ncread('../obs_aviso/dt_ref_med_merged_madt_uv_20080102_20080102_20100503.nc','NbLongitudes');
4Lat_ref=ncread('../obs_aviso/dt_ref_med_merged_madt_uv_20080102_20080102_20100503.nc','NbLatitudes');
5Lon_ref=Lon_ref(end-86:end);
6Lat_ref=Lat_ref(1:58);
7
8Lon=ncread(aviso_file,'lon');
9Lat=ncread(aviso_file,'lat');
10
11u=ncread(aviso_file,'u');
12v=ncread(aviso_file,'v');
13u=double(u')*0.0001; %Transpose to have the same dim as ref
14v=double(v')*0.0001;
15
16[Lon2,Lat2]=meshgrid(Lon,Lat);
17[Lon_ref2,Lat_ref2]=meshgrid(Lon_ref,Lat_ref);
18
19ui=interp2(Lon2,Lat2,u,Lon_ref2,Lat_ref2);
20vi=interp2(Lon2,Lat2,v,Lon_ref2,Lat_ref2);
Note: See TracBrowser for help on using the repository browser.