%
% A simple Matlab script to read the wind data netcdf file.
%
% Author: Nathan Bindoff
% Date: 5 June 2000
% modified by GBR 13/06/2000
% The script uses Jim Mansbridge's netcdf routines for matlab.
%
% Documetation and description of these routines can be found at 
% http://www.marine.csiro.au/sw/matlab-netcdf.html.
%
%
%
% Examine one wind data file
%
% To show the extracted parameter, you should use the folowing commands :
%    pcolor ( parameter / Scale) ; 
%    example : pcolor (Zonal_wind / attnc(file,'Zonal_wind','Scale'));
%    shading flat                ;
%    colormap ( jet )            ;
%    colorbar                    ;
%    caxis([0 3000])             ;
% and see the result on the Matlab's figure.
%

inqnc(file)
attnc(file)
%
%
Latitude=getnc(file,'Latitude')*attnc(file,'Latitude','Scale');
Longitude=getnc(file,'Longitude')*attnc(file,'Longitude','Scale');
Num_Sample=getnc(file,'Num_Sample')*attnc(file,'Num_Sample','Scale');
Wind_speed=getnc(file,'Wind_speed')*attnc(file,'Wind_speed','Scale');
Zonal_wind=getnc(file,'Zonal_wind')*attnc(file,'Zonal_wind','Scale');
Meridional_wind=getnc(file,'Meridional_wind')*attnc(file,'Meridional_wind','Scale');
Wind_error=getnc(file,'Wind_error')*attnc(file,'Wind_error','Scale');
Zonal_wind_error=getnc(file,'Zonal_wind_error')*attnc(file,'Zonal_wind_error','Scale');
Meridional_wind_error=getnc(file,'Meridional_wind_error')*attnc(file,'Meridional_wind_error','Scale');
Wind_stress=getnc(file,'Wind_stress')*attnc(file,'Wind_stress','Scale');
Zonal_stress=getnc(file,'Zonal_stress')*attnc(file,'Zonal_stress','Scale');
Meridional_stress=getnc(file,'Meridional_stress')*attnc(file,'Meridional_stress','Scale');
Wind_stress_error=getnc(file,'Wind_stress_error')*attnc(file,'Wind_stress_error','Scale');
Zonal_stress_error=getnc(file,'Zonal_stress_error')*attnc(file,'Zonal_stress_error','Scale');
Meridional_stress_error=getnc(file,'Meridional_stress_error')*attnc(file,'Meridional_stress_error','Scale');
Divergence=getnc(file,'Divergence')*attnc(file,'Divergence','Scale');
Wind_curl=getnc(file,'Wind_curl')*attnc(file,'Wind_curl','Scale');

%The 2 dimensions are  1) Latitude = 160  2) Longitude = 360.
%None of the dimensions is unlimited
% 
%----- Get further information about the following variables -----
% 
%  -1) None of them (no further information)
%   0) All of the variables
%   1) Latitude               2) Longitude              3) Num_Sample
%   4) Wind_speed             5) Zonal_wind             6) Meridional_wind
%   7) Wind_error             8) Zonal_wind_error       9) Meridional_wind_error
%  10) Wind_stress           11) Zonal_stress          12) Meridional_stress
%  13) Wind_stress_error     14) Zonal_stress_error    15) Meridional_stress_error
%  16) Divergence            17) Wind_curl      
%

