function []=MITprof_global(varargin);
%object:    take care of path and global variables (mygrid and mitprofenv), 
%           and sends global variables to caller routine workspace
%notes:     - in any call, if this has not yet been done, 
%           this routine also adds MITprof subdirectories 
%           to the matlab path, and it defines mitprofenv.
%           - this routine replaces MITprof_path

%get/define global variables:
gcmfaces_global;

%take care of path:
test0=which('MITprof_load.m');
if isempty(test0);
    test0=which('MITprof_global.m'); ii=strfind(test0,filesep);
    mydir=test0(1:ii(end));
    %
    addpath(fullfile(mydir));
    addpath(fullfile(mydir,'profiles_process_main_v2'));
    addpath(fullfile(mydir,'profiles_IO_v2'));
    addpath(fullfile(mydir,'profiles_IO_external'));
    addpath(fullfile(mydir,'profiles_misc'));
    addpath(fullfile(mydir,'profiles_stats'));
    addpath(fullfile(mydir,'ecco_v4'));
    addpath(fullfile(mydir,'profiles_devel'));
end;

%environment variables:
if ~isfield(myenv,'MITprof_dir');
    test0=which('MITprof_global.m'); ii=strfind(test0,filesep);
    myenv.MITprof_dir=test0(1:ii(end));
    %
    gridDir='';
    tmpDir=fullfile(myenv.gcmfaces_dir,'..','GRID',filesep);
    if isdir(tmpDir); gridDir=tmpDir; end;
    tmpDir=fullfile(myenv.gcmfaces_dir,'..','nctiles_grid',filesep);
    if isdir(tmpDir); gridDir=tmpDir; end;
    if isdir('GRID/'); gridDir='GRID/'; end;
    if isdir('nctiles_grid/'); gridDir='nctiles_grid/'; end;

    if isempty(gridDir); error('could not find grid'); end;
    addpath(gridDir);
    fil=which('GRID.0001.nc');
    if isempty(fil); fil=which('XC.meta'); end;
    if isempty(fil); error('could not find grid'); end;
    myenv.MITprof_griddir=[fileparts(fil) filesep];
    %
    climDir='';    
    tmpDir=fullfile(myenv.gcmfaces_dir,'sample_input','OCCAetcONv4GRID',filesep);
    if isdir(tmpDir); climDir=tmpDir; end;
    tmpDir=fullfile(myenv.MITprof_dir,'..','MITprof_climatologies',filesep);
    if isdir(tmpDir); climDir=tmpDir; end;
    tmpDir=fullfile('sample_input','OCCAetcONv4GRID',filesep);
    if isdir(tmpDir); climDir=tmpDir; end;
    tmpDir=fullfile('MITprof_climatologies',filesep);
    if isdir(tmpDir); climDir=tmpDir; end;

    addpath(climDir);
    fil=which('sigma_T_mad_feb2013.bin');
    if isempty(fil);
     fprintf('\n please indicate the climatologies directory (e.g., ''MITprof_climatologies/'') \n\n');
     fprintf('   It can be obtained as follows: \n');
     fprintf('   wget --recursive ftp://mit.ecco-group.org/gforget/OCCAetcONv4GRID .\n\n');
     fprintf('   mv mit.ecco-group.org/gforget/OCCAetcONv4GRID MITprof_climatologies \n\n');
     myenv.MITprof_climdir=input('');
    end;
 
    addpath(climDir);
    fil=which('sigma_T_mad_feb2013.bin');
    if isempty(fil); error('could not find sigma_T_mad_feb2013.bin'); end;
    myenv.MITprof_climdir=[fileparts(fil) filesep];

end;

%send to workspace:
evalin('caller','global mygrid mitprofenv');