--- MITgcm_contrib/gael/profilesMatlabProcessing/MITprof_global.m 2014/12/03 01:17:50 1.4 +++ MITgcm_contrib/gael/profilesMatlabProcessing/MITprof_global.m 2016/01/17 19:28:27 1.7 @@ -29,8 +29,47 @@ if ~isfield(myenv,'MITprof_dir'); test0=which('MITprof_global.m'); ii=strfind(test0,filesep); myenv.MITprof_dir=test0(1:ii(end)); - myenv.MITprof_griddir=fullfile(myenv.gcmfaces_dir,'..','GRID',filesep); - myenv.MITprof_climdir=fullfile(myenv.gcmfaces_dir,'sample_input','OCCAetcONv4GRID',filesep); + % + 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: