gael/profilesMatlabProcessing/MITprof_global.m

function []=MITprof_global(varargin);
% 
% MITPROF_GLOBAL calls gcmfaces_global, adds MITprof paths, adds
%    defines MITprof_climdir and MITprof_griddir in myenv, and 
%    adds myenv (global variable) to caller routine workspace

%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,'..','gcmfaces_climatologies',filesep);
    if isdir(tmpDir); climDir=tmpDir; end;
    tmpDir=fullfile('sample_input','OCCAetcONv4GRID',filesep);
    if isdir(tmpDir); climDir=tmpDir; end;
    tmpDir=fullfile('gcmfaces_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., ''gcmfaces_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 gcmfaces_climatologies \n\n');
     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 myenv');

1	gforget	1.1	function []=MITprof_global(varargin);
2	gforget	1.8	%
3			% MITPROF_GLOBAL calls gcmfaces_global, adds MITprof paths, adds
4			% defines MITprof_climdir and MITprof_griddir in myenv, and
5			% adds myenv (global variable) to caller routine workspace
6	gforget	1.1
7			%get/define global variables:
8	gforget	1.2	gcmfaces_global;
9	gforget	1.1
10			%take care of path:
11			test0=which('MITprof_load.m');
12			if isempty(test0);
13			test0=which('MITprof_global.m'); ii=strfind(test0,filesep);
14			mydir=test0(1:ii(end));
15			%
16			addpath(fullfile(mydir));
17			addpath(fullfile(mydir,'profiles_process_main_v2'));
18			addpath(fullfile(mydir,'profiles_IO_v2'));
19			addpath(fullfile(mydir,'profiles_IO_external'));
20			addpath(fullfile(mydir,'profiles_misc'));
21			addpath(fullfile(mydir,'profiles_stats'));
22			addpath(fullfile(mydir,'ecco_v4'));
23			addpath(fullfile(mydir,'profiles_devel'));
24			end;
25
26			%environment variables:
27	gforget	1.2	if ~isfield(myenv,'MITprof_dir');
28	gforget	1.1	test0=which('MITprof_global.m'); ii=strfind(test0,filesep);
29	gforget	1.2	myenv.MITprof_dir=test0(1:ii(end));
30	gforget	1.5	%
31			gridDir='';
32			tmpDir=fullfile(myenv.gcmfaces_dir,'..','GRID',filesep);
33			if isdir(tmpDir); gridDir=tmpDir; end;
34			tmpDir=fullfile(myenv.gcmfaces_dir,'..','nctiles_grid',filesep);
35			if isdir(tmpDir); gridDir=tmpDir; end;
36			if isdir('GRID/'); gridDir='GRID/'; end;
37			if isdir('nctiles_grid/'); gridDir='nctiles_grid/'; end;
38
39			if isempty(gridDir); error('could not find grid'); end;
40			addpath(gridDir);
41			fil=which('GRID.0001.nc');
42			if isempty(fil); fil=which('XC.meta'); end;
43			if isempty(fil); error('could not find grid'); end;
44			myenv.MITprof_griddir=[fileparts(fil) filesep];
45			%
46			climDir='';
47			tmpDir=fullfile(myenv.gcmfaces_dir,'sample_input','OCCAetcONv4GRID',filesep);
48			if isdir(tmpDir); climDir=tmpDir; end;
49	gforget	1.7	tmpDir=fullfile(myenv.MITprof_dir,'..','gcmfaces_climatologies',filesep);
50	gforget	1.5	if isdir(tmpDir); climDir=tmpDir; end;
51			tmpDir=fullfile('sample_input','OCCAetcONv4GRID',filesep);
52			if isdir(tmpDir); climDir=tmpDir; end;
53	gforget	1.7	tmpDir=fullfile('gcmfaces_climatologies',filesep);
54	gforget	1.5	if isdir(tmpDir); climDir=tmpDir; end;
55
56			addpath(climDir);
57			fil=which('sigma_T_mad_feb2013.bin');
58			if isempty(fil);
59	gforget	1.7	fprintf('\n please indicate the climatologies directory (e.g., ''gcmfaces_climatologies/'') \n\n');
60	gforget	1.5	fprintf(' It can be obtained as follows: \n');
61			fprintf(' wget --recursive ftp://mit.ecco-group.org/gforget/OCCAetcONv4GRID .\n\n');
62	gforget	1.7	fprintf(' mv mit.ecco-group.org/gforget/OCCAetcONv4GRID gcmfaces_climatologies \n\n');
63	gforget	1.6	climDir=input('');
64	gforget	1.5	end;
65
66			addpath(climDir);
67			fil=which('sigma_T_mad_feb2013.bin');
68			if isempty(fil); error('could not find sigma_T_mad_feb2013.bin'); end;
69			myenv.MITprof_climdir=[fileparts(fil) filesep];
70
71	gforget	1.1	end;
72
73			%send to workspace:
74	gforget	1.8	evalin('caller','global mygrid myenv');
75	gforget	1.1