Changeset 25733

Timestamp:

11/03/20 12:17:53 (4 years ago)

Author:

dlcheng

Message:

CHG: dlcheng contrib.

Location:

issm/trunk-jpl/src/m/contrib

Files:

• dlcheng (added)
• dlcheng/MeshMergedOutlineToShp.m (added)
• dlcheng/MeshOutlineToShp.m (added)
• dlcheng/VxVyToVv.m (added)
• dlcheng/interpBedmachineGreenland.m (added)
• morlighem/modeldata/interpBedmachineGreenland.m (modified) (2 diffs)

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

@@ -54 +54 @@
                 morlighem2013nc=['/home/ModelData/Greenland/BedMachine/' basename '-' ncdate '.nc'];
         otherwise
-                error('machine not supported yet');
+                morlighem2013nc=['/Users/larour/ModelData/BedMachine/' basename '-' ncdate '.nc'];
+                disp('machine not supported yet, using default path /Users/larour/ModelData/BedMachine/');
 end
 
@@ -81 +82 @@
 data(find(data==-9999))=NaN;
 
+ydata=flipud(ydata);
+data=flipud(data);
+    
 disp(['   -- BedMachine Greenland: interpolating ' string]);
 if strcmp(string,'mask') | strcmp(string,'source'),
         %Need nearest neighbor to avoid interpolation between 0 and 2
-        output = InterpFromGrid(xdata,ydata,data,double(X),double(Y),'nearest');
+        output = InterpFromGridToMesh(xdata,ydata,data,double(X),double(Y),'nearest');
 else
-        output = InterpFromGrid(xdata,ydata,data,double(X),double(Y));
+    
+        output = InterpFromGridToMesh(xdata,ydata,data,double(X),double(Y),NaN);
 end