Changeset 27694

Timestamp:

04/19/23 12:01:03 (2 years ago)

Author:

Cheng Gong

Message:

ADD: a function to interpolate ice mask from Greene monthly reconstruction data

Location:

issm/trunk-jpl/src/m/contrib/chenggong

Files:

• dataprocessing/computeGrad.m (modified) (1 diff)
• dataprocessing/extractTransientSolutions.m (modified) (1 diff)
• dataprocessing/integrateOverDomain.m (modified) (2 diffs)
• interpFromMEaSUREsGeotiff.m (modified) (1 diff)
• modeldata/interpMonthlyIceMaskGreene.m (added)
• visualization/plotCompareTransientFlowline.m (modified) (2 diffs)

issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/computeGrad.m

@@ -1 +1 @@
 function [gradx, grady]=computeGrad(index,x,y,field)
-%COMPUTEHESSIAN - compute the gradient from a field
+%COMPUTEGRAD - compute the gradient from a field
 
 %some variables

issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/extractTransientSolutions.m

@@ -9 +9 @@
 transientSolutions.thickness = cell2mat({md.results.TransientSolution(:).Thickness});
 transientSolutions.SigmaVM = cell2mat({md.results.TransientSolution(:).SigmaVM});
-transientSolutions.smb = cell2mat({md.results.TransientSolution(:).SmbMassBalance});
+if (isfield(md.results.TransientSolution, 'SmbMassBalance'))
+        transientSolutions.smb = cell2mat({md.results.TransientSolution(:).SmbMassBalance});
+end
 transientSolutions.ice_levelset = cell2mat({md.results.TransientSolution(:).MaskIceLevelset});
 transientSolutions.calvingRate = cell2mat({md.results.TransientSolution(:).CalvingCalvingrate});

issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/integrateOverDomain.m

@@ -8 +8 @@
         weights = ones(size(data));
         if nargin<3
-                masked = [];
+                masked = logical(zeros(size(data)));
         end
 end
 
+masked = masked | isnan(data) | isnan(weights);
 % Set the area with masked=1 to nan
 data(masked) = nan;
 weights(masked) =nan;
+
 
 % get the mesh
@@ -28 +30 @@
 eleAreas = 1/3*eleAreas.*(weights(elements(:,1),:)+weights(elements(:,2),:)+weights(elements(:,3),:));
 
-intData = sum(eleData(:),'omitnan');
-areas = sum(eleAreas(:),'omitnan');
-meanData = intData / areas;
+intData = sum(eleData, 1, 'omitnan');
+areas = sum(eleAreas, 1, 'omitnan');
+meanData = intData ./ areas;

issm/trunk-jpl/src/m/contrib/chenggong/interpFromMEaSUREsGeotiff.m

@@ -27 +27 @@
 elseif strcmp(glacier, 'Rink')
         foldername = '/totten_1/ModelData/Greenland/VelMEaSUREs/Rink_2008_2022/';
+elseif strcmp(glacier, 'Upernavik')
+        foldername = '/totten_1/ModelData/Greenland/VelMEaSUREs/Upernavik_2008_2022/';
+elseif strcmp(glacier, 'Helheim')
+        foldername = '/totten_1/ModelData/Greenland/VelMEaSUREs/Helheim_2008_2023/';
 else
         error(['The velocity data for ', glacier, ' is not available, please download from NSIDC first.']);

issm/trunk-jpl/src/m/contrib/chenggong/visualization/plotCompareTransientFlowline.m

@@ -16 +16 @@
 %
 %   Author: Cheng Gong
-%   Date: 2021-12-06
+%   Date: 2023-04-10
 
 N = length(velList);
@@ -46 +46 @@
     h=colorbar;
     title(h,'m/a')
+         % add a vertical line to indicate initial ice front position
+         id = max(find(~isnan(cumsum(vel_flowline{p}(:,1)))));
+         icefront = flowline.Xmain(id);
+         plot([icefront, icefront], [2007,2020], '-r')
 end