--- MITgcm_contrib/gael/profilesMatlabProcessing/MITprof_global.m	2014/07/28 19:34:16	1.1
+++ MITgcm_contrib/gael/profilesMatlabProcessing/MITprof_global.m	2016/01/31 15:52:18	1.8
@@ -1,13 +1,11 @@
 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
+% 
+% 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:
-global mitprofenv;
+gcmfaces_global;
 
 %take care of path:
 test0=which('MITprof_load.m');
@@ -26,16 +24,52 @@
 end;
 
 %environment variables:
-if isempty(mitprofenv);
+if ~isfield(myenv,'MITprof_dir');
     test0=which('MITprof_global.m'); ii=strfind(test0,filesep);
-    mitprofenv.MITprof_dir=test0(1:ii(end));
-    mitprofenv.verbose=0;
-    mitprofenv.lessplot=0;
-    mitprofenv.lesstest=0;
-    mitprofenv.useNativeMatlabNetcdf = ~isempty(which('netcdf.open'));
-    %... check for gcmfaces
+    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 mitprofenv');
+evalin('caller','global mygrid myenv');