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	function []=MITprof_global(varargin);
2	%
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
7	%get/define global variables:
8	gcmfaces_global;
9
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	if ~isfield(myenv,'MITprof_dir');
28	test0=which('MITprof_global.m'); ii=strfind(test0,filesep);
29	myenv.MITprof_dir=test0(1:ii(end));
30	%
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	tmpDir=fullfile(myenv.MITprof_dir,'..','gcmfaces_climatologies',filesep);
50	if isdir(tmpDir); climDir=tmpDir; end;
51	tmpDir=fullfile('sample_input','OCCAetcONv4GRID',filesep);
52	if isdir(tmpDir); climDir=tmpDir; end;
53	tmpDir=fullfile('gcmfaces_climatologies',filesep);
54	if isdir(tmpDir); climDir=tmpDir; end;
55
56	addpath(climDir);
57	fil=which('sigma_T_mad_feb2013.bin');
58	if isempty(fil);
59	fprintf('\n please indicate the climatologies directory (e.g., ''gcmfaces_climatologies/'') \n\n');
60	fprintf(' It can be obtained as follows: \n');
61	fprintf(' wget --recursive ftp://mit.ecco-group.org/gforget/OCCAetcONv4GRID .\n\n');
62	fprintf(' mv mit.ecco-group.org/gforget/OCCAetcONv4GRID gcmfaces_climatologies \n\n');
63	climDir=input('');
64	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	end;
72
73	%send to workspace:
74	evalin('caller','global mygrid myenv');
75