function []=plot_std_field();
%
%PLOT_STD_FIELD computes and displays a standard deviation field
%
%stand-alone call: addpath gcmfaces/sample_analysis/; plot_std_field;
%
%needed input files:
%  mkdir release1
%  wget --recursive ftp://mit.ecco-group.org/ecco_for_las/version_4/release1/nctiles_climatology/ETAN
%  mv mit.ecco-group.org/ecco_for_las/version_4/release1/nctiles_climatology release1/.

p = genpath('gcmfaces/'); addpath(p);
p = genpath('MITprof/'); addpath(p);

gcmfaces_global;

input_list_check('plot_std_field',nargin);

%expected location:
myenv.nctilesdir=fullfile(pwd,'/release1/nctiles_climatology/ETAN/');
%if ETAN is not found then try old location:
if ~isdir(myenv.nctilesdir);
    %if not found then try old location:
    tmpdir=fullfile(pwd,'/gcmfaces/sample_input/nctiles_climatology/ETAN/');
    if isdir(tmpdir); myenv.nctilesdir=tmpdir; end;
end;
%if ETAN is still not found then issue warning and skip plot_std_field
if ~isdir(myenv.nctilesdir);
    help plot_std_field;
    warning('plot_std_field requires release1/nctiles_climatology/ETAN that was not found ---> abort!');
    return;
end;

%%%%%%%%%%%%%%%%%
%load grid:
%%%%%%%%%%%%%%%%%

mygrid=[];
if isempty(mygrid);
   grid_load;
end;

%%%%%%%%%%%
%get field:
%%%%%%%%%%%
fld=read_nctiles([myenv.nctilesdir 'ETAN'],'ETAN');
fld=std(fld,[],3); fld(find(fld==0))=NaN;
cc=[0:0.1:1]*0.10;

%%%%%%%%%%%%
%plot field:
%%%%%%%%%%%%
if ~myenv.lessplot;
    figureL; gcmfaces_sphere(fld,cc,[],[],1);
end;