
hola='/home/hrd/annane/AL092012_0828_1630'
% Process HWind file
[hwind_xyres,hwind_slon,hwind_slat,hwind_dx,hwind_dy,hwind_dlon,hwind_dlat,hwind_u10,hwind_v10] = read_gridded_hwind_data(hola);

% Get 10-m wind speed and its maximum
hwind_wspd = sqrt(hwind_u10.^2+hwind_v10.^2);
hwind_wspdmax = nanmax(nanmax(hwind_wspd));

% Convert to polar coordinates
[hwind_dxp,hwind_dyp,hwind_azm,hwind_rge,hwind_wspdpol] = do_polar_conversion2(hwind_dx,hwind_dy,hwind_wspd,0:360,0:500);

% Compute wavenumber components
[hwind_wspdpol_wvnum,hwind_varexp,hwind_ampval,hwind_ampratio] = get_wavenumber(hwind_wspdpol,5);

