function [a]=swan_read(fn,nx,ny)
% Read an output array from SWAN run
% [a]=swan_read(fn,nx,ny)
% 
disp(['Reading ', fn])
fid = fopen(fn,'r');
atr = fscanf(fid, '%f', [nx, ny]);
fclose(fid);
a = atr';
