%------------------------------------------------------------------
% Script for plotting geostrophic currents and statistics from reference
% levels. 
% 
% PARAMETERS
%
%   year_base
%
%   file_name
%      Root name of data file. Functions called within this script will add
%      automatically the suffixes to get name of grid and statistics files.
%
%   x_var
%      Index to x variable. This is a number corresponding to the x column in
%      the data files. Possible values are:
%
%         1: time
%         2: longitude
%         3: latitude
%
%      If negative, the axis direction will be reversed.
%
%   n_filter
%      Number (>=2) used to filter data. Set it to one to disable filtering.
%
%   depth
%      Vector containing list of depth from which grd data should be plotted.
%      If empty, all available depths are plotted.
%
%   plot_title
%
%   institute
%
%   info
%      String for additionnal information on plot
%
%   plot_grd
%   plot_sta
%      Flags (1 or 0) for enabling plotting of grd and sta data, respectively.
%
% OUTPUT
%
%   <file_name>_sta.ps
%   <file_name>_grd.ps
%      Graphical output files.
%
%------------------------------------------------------------------

% Edit following parameters
%------------------------------------------------------------------

year_base  = 1997;
file_name  = 'vgeo';
x_var      = 3
depth      = [];
n_filter   = 10;
plot_title = 'AGMASCO-AZORES-SSR/ADCP Geostrophic Current (m/s)';
institute  = 'P.Jaccard, Geophysical Institute, University of Bergen';
info       = 'Positive westward';
plot_grd   = 1;
plot_sta   = 1;

%------------------------------------------------------------------
% No more editing!!!

if plot_grd,
	vgeo_grd(year_base, file_name, x_var, depth, n_filter, ...
			plot_title, institute, info);
end

if plot_sta,
	vgeo_sta(year_base, file_name, x_var, plot_title, institute, info);
end

