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); |
40 |
fprintf('\n please indicate the ECCO v4 grid directory (e.g., ''nctiles_grid/'') \n\n'); |
41 |
fprintf(' It can be obtained as follows: \n'); |
42 |
fprintf(' wget --recursive ftp://mit.ecco-group.org/ecco_for_las/version_4/release2/nctiles_grid/ .\n'); |
43 |
fprintf(' mv mit.ecco-group.org/ecco_for_las/version_4/release2/nctiles_grid/ . \n\n'); |
44 |
gridDir=input(''); |
45 |
end; |
46 |
|
47 |
addpath(gridDir); |
48 |
fil=which('GRID.0001.nc'); |
49 |
if isempty(fil); fil=which('XC.meta'); end; |
50 |
if isempty(fil); error('could not find grid'); end; |
51 |
myenv.MITprof_griddir=[fileparts(fil) filesep]; |
52 |
% |
53 |
climDir=''; |
54 |
tmpDir=fullfile(myenv.gcmfaces_dir,'sample_input','OCCAetcONv4GRID',filesep); |
55 |
if isdir(tmpDir); climDir=tmpDir; end; |
56 |
tmpDir=fullfile(myenv.MITprof_dir,'..','gcmfaces_climatologies',filesep); |
57 |
if isdir(tmpDir); climDir=tmpDir; end; |
58 |
tmpDir=fullfile('sample_input','OCCAetcONv4GRID',filesep); |
59 |
if isdir(tmpDir); climDir=tmpDir; end; |
60 |
tmpDir=fullfile('gcmfaces_climatologies',filesep); |
61 |
if isdir(tmpDir); climDir=tmpDir; end; |
62 |
|
63 |
addpath(climDir); |
64 |
fil=which('sigma_T_mad_feb2013.bin'); |
65 |
if isempty(fil); |
66 |
fprintf('\n please indicate the climatologies directory (e.g., ''gcmfaces_climatologies/'') \n\n'); |
67 |
fprintf(' It can be obtained as follows: \n'); |
68 |
fprintf(' wget --recursive ftp://mit.ecco-group.org/gforget/OCCAetcONv4GRID .\n'); |
69 |
fprintf(' mv mit.ecco-group.org/gforget/OCCAetcONv4GRID gcmfaces_climatologies \n\n'); |
70 |
climDir=input(''); |
71 |
end; |
72 |
|
73 |
addpath(climDir); |
74 |
fil=which('sigma_T_mad_feb2013.bin'); |
75 |
if isempty(fil); error('could not find sigma_T_mad_feb2013.bin'); end; |
76 |
myenv.MITprof_climdir=[fileparts(fil) filesep]; |
77 |
|
78 |
end; |
79 |
|
80 |
%send to workspace: |
81 |
evalin('caller','global mygrid myenv'); |
82 |
|