Changeset 25735

Timestamp:

11/03/20 12:32:32 (4 years ago)

Author:

Mathieu Morlighem

Message:

CHG: merging with Dan's code

File:

• issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpBedmachineGreenland.m (modified) (2 diffs)

issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpBedmachineGreenland.m

@@ -3 +3 @@
 if nargin<3, string = 'bed'; end
 if nargin<4,
-        %ncdate='2013-05-21';
-        %ncdate='2013-06-27';
-        %ncdate='2013-07-18';
-        %ncdate='2013-11-15';
-        %ncdate='2013-12-03';
-        %ncdate='2014-02-26';
-        %ncdate='2014-03-24';
-        %ncdate='2014-07-31';
-        %ncdate='2014-11-14';
-        %ncdate='2015-03-03';
-        %ncdate='2015-03-10';
-        %ncdate='2015-03-26';
-        %ncdate='2015-03-30';
         %ncdate='2015-04-27'; %BedMachine v2
-        %ncdate='2015-07-30';
-        %ncdate='2015-10-02';
-        %ncdate='2016-03-21';
-        %ncdate='2016-05-12';
-        ncdate='2016-07-06';
-        ncdate='2016-08-04';
-        ncdate='2016-10-26';
-        ncdate='2016-11-23';
-        ncdate='2016-12-21';
-        ncdate='2017-01-19';
-        ncdate='2017-03-28';
-        ncdate='2017-05-10';
-        ncdate='2017-07-21';
         ncdate='2017-09-25'; %BedMachine v3
-        ncdate='2018-06-01';
-        ncdate='2019-04-18';
         ncdate='2020-04-14';
 end
 
-if exist('datetime','file') 
-        date1 = sscanf(ncdate,'%d-%d-%d');
-        date2 = datetime(date1(1),date1(2),date1(3));
-        if date2<datetime(2016,10,24),
-                basename = 'MCdataset'; 
-        else
-                basename = 'BedMachineGreenland';
+basename = 'BedMachineGreenland';
+
+%List of common paths to try
+paths = {...
+        ['/u/astrid-r1b/ModelData/ModelData/MCdataset-' ncdate '.nc'],...
+        ['/home/ModelData/Greenland/BedMachine/' basename '-' ncdate '.nc'],...
+        ['/Users/larour/ModelData/BedMachine/' basename '-' ncdate '.nc'],...
+        ['./' basename '-' ncdate '.nc'],...
+        };
+
+found = 0;
+for i=1:numel(paths)+1
+        if exist(paths{i},'file')
+                found = 1;
+                ncfile = parths{i}
         end
-else
-  basename = 'BedMachineGreenland';
 end
 
-switch oshostname(),
-        case {'murdo','thwaites','astrid'}
-                morlighem2013nc=['/u/astrid-r1b/ModelData/ModelData/MCdataset-' ncdate '.nc']';
-        case {'ronne'}
-                morlighem2013nc=['/home/ModelData/Greenland/BedMachine/' basename '-' ncdate '.nc'];
-        otherwise
-                error('machine not supported yet');
+if ~found
+        error(['Could not find ' basename '-' ncdate '.nc, you can add the path to the list');
 end
 
+
 disp(['   -- BedMachine Greenland version: ' ncdate]);
-xdata = double(ncread(morlighem2013nc,'x'));
-ydata = double(ncread(morlighem2013nc,'y'));
+xdata = double(ncread(ncfile,'x'));
+ydata = double(ncread(ncfile,'y'));
 
 offset=2;
@@ -76 +50 @@
 
 disp(['   -- BedMachine Greenland: loading ' string]);
-data  = double(ncread(morlighem2013nc,string,[id1x id1y],[id2x-id1x+1 id2y-id1y+1],[1 1]))';
+data  = double(ncread(ncfile,string,[id1x id1y],[id2x-id1x+1 id2y-id1y+1],[1 1]))';
 xdata=xdata(id1x:id2x);
 ydata=ydata(id1y:id2y);