profilesMatlabProcessing/profiles_misc/aviso_get_grid.m

function grid_aviso=aviso_get_grid
% function [grid_aviso]=aviso_get_grid
%  load grid using loaddap interface to opendap website:
%  http://opendap.aviso.oceanobs.com/thredds/dodsC/dataset-duacs-dt-ref-global-merged-madt-h
%
% output: grid_aviso is a structure
%   grid_aviso.NbLatitudes : latitudes
%   grid_aviso.NbLongitudes : longitudes
%   grid_aviso.time : time (matlab julian day)
%   grid_aviso.Nlat : number of latitudes
%   grid_aviso.Nlon : number of longitudes
%   grid_aviso.Ntime : number of time samples

URL='http://opendap.aviso.oceanobs.com/thredds/dodsC/dataset-duacs-dt-ref-global-merged-madt-h?NbLatitudes[0:1:914],NbLongitudes[0:1:1079],time[0:1:911]';
grid_aviso=loaddap(URL);
if isempty(grid_aviso), error('access problem to aviso data with opendap'); end

grid_aviso.time=grid_aviso.time/24+datenum([1950 01 01 0 0 0]);
grid_aviso.Nlat=length(grid_aviso.NbLatitudes);
grid_aviso.Nlon=length(grid_aviso.NbLongitudes);
grid_aviso.Ntime=length(grid_aviso.time);
1	roquet	1.1	function grid_aviso=aviso_get_grid
2			% function [grid_aviso]=aviso_get_grid
3			% load grid using loaddap interface to opendap website:
4			% http://opendap.aviso.oceanobs.com/thredds/dodsC/dataset-duacs-dt-ref-global-merged-madt-h
5			%
6			% output: grid_aviso is a structure
7			% grid_aviso.NbLatitudes : latitudes
8			% grid_aviso.NbLongitudes : longitudes
9			% grid_aviso.time : time (matlab julian day)
10			% grid_aviso.Nlat : number of latitudes
11			% grid_aviso.Nlon : number of longitudes
12			% grid_aviso.Ntime : number of time samples
13
14			URL='http://opendap.aviso.oceanobs.com/thredds/dodsC/dataset-duacs-dt-ref-global-merged-madt-h?NbLatitudes[0:1:914],NbLongitudes[0:1:1079],time[0:1:911]';
15			grid_aviso=loaddap(URL);
16			if isempty(grid_aviso), error('access problem to aviso data with opendap'); end
17
18			grid_aviso.time=grid_aviso.time/24+datenum([1950 01 01 0 0 0]);
19			grid_aviso.Nlat=length(grid_aviso.NbLatitudes);
20			grid_aviso.Nlon=length(grid_aviso.NbLongitudes);
21			grid_aviso.Ntime=length(grid_aviso.time);