| 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 |
|