% This script loads a NetCDF file named wind_stress_zonal_monthly_maps.nc and plots the 
% zonal monthly wind stress field for January.
%
% wind_stress_zonal_monthly_maps.nc contains global ocean zonal monthly wind stress fields 
% for January through December. These fields form part of the SCOW wind atlas, which is avaiable at 
% http://cioss.coas.oregonstate.edu/scow/.
% 
% This script was written by Craig Risien on 31 May 2007 and tested using Matlab 5.3.1.29215a (R11.1) 
% and the NetCDF Toolbox For Matlab-5 Version of 30-Apr-2003 11:16:19.


clear all
close all

dummy = netcdf('wind_stress_zonal_monthly_maps.nc','nowrite');

ncdump(dummy)	% the ncdump command will give you a listing of the NetCDF file headers.

% extract NetCDF variables of interest from wind_stress_zonal_monthly_maps.nc.

temp = dummy{'january'};	% to plot the zonal wind stress for other months replace 'january' with 'february', 'march', ..., or 'december'.
january = squeeze(temp(:,:));
january(find(january==-9999))=nan;	% missing data are flagged as -9999.

temp = dummy{'latitude'};
latitude = squeeze(temp(:,:));

temp = dummy{'longitude'};
longitude = squeeze(temp(:,:));

% plot the zonal monthly wind stress field for January.

pcolor(longitude,latitude,january)
shading flat 
caxis([-.4 .4])
colorbar
title('SCOW January Zonal Wind Stress (N/m^2)')
xlabel('longitude')
ylabel('latitude')
orient landscape
print -dtiff -r150 January_Zonal_Wind_Stress