Changeset 26975

Timestamp:

04/29/22 04:22:16 (3 years ago)

Author:

Cheng Gong

Message:

add post-processing codes

Location:

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

Files:

• dataprocessing (added)
• dataprocessing/README.md (added)
• dataprocessing/analyzeCalvingFront.m (added)
• dataprocessing/averageOverTime.m (added)
• dataprocessing/computeGrad.m (added)
• dataprocessing/extractTransientSolutions.m (added)
• dataprocessing/fillInNan.m (added)
• dataprocessing/integrateOverDomain.m (added)
• dataprocessing/integrateOverTime.m (added)
• dataprocessing/interpZeroPos.m (added)
• dataprocessing/inverseHistograms.m (added)
• dataprocessing/movingAverage.m (added)
• dataprocessing/npsd.m (added)
• dataprocessing/projectToFlowlines.m (added)
• dataprocessing/psd.m (added)
• dataprocessing/rescalegradientNan.m (added)
• dataprocessing/wfDistance.m (added)
• expxy2shpll.m (added)
• interpFromMEaSUREsGeotiff.m (modified) (1 diff)
• modeldata (added)
• modeldata/interpISMIP6Temp.m (added)
• visualization (added)
• visualization/imageNonUni.m (added)
• visualization/plotCompareTransientFlowline.m (added)

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

@@ -1 @@
-% This function calls src/m/contrib/morlighem/modeldata/interpFromGeotiff.m for multiple times to load all avaliable 
-% tif data in  /totten_1/ModelData/Greenland/VelMEaSUREs/Jakobshavn_2008_2021/ within the given time period (in decimal years)
-% For some reason, each .tif file in this folder contains two sets of data, only the first dataset is useful
+function dataout = interpFromMEaSUREsGeotiff(X,Y,Tstart,Tend,varargin)
+%interpFromMEaSUREsGeotiff: 
+%       This function calls src/m/contrib/morlighem/modeldata/interpFromGeotiff.m for multiple times to load all avaliable 
+%       tif data in  /totten_1/ModelData/Greenland/VelMEaSUREs/Jakobshavn_2008_2021/ within the given time period (in decimal years)
+%       For some reason, each .tif file in this folder contains two sets of data, only the first dataset is useful
+%
+%   Usage:
+%                dataout = interpFromMEaSUREsGeotiff(X,Y,Tstart,Tend, varargin)
+%
+%       X, Y are the coordinates of the mesh 
+%       Tstart and Tend decimal year of the start and end time
+%
+%   Example:
+%                       obsData = interpFromMEaSUREsGeotiff(md.mesh.x,md.mesh.y, tstart, tend);
+%
+%   Options:
+%      - 'glacier':  which glacier to look for
+options    = pairoptions(varargin{:});
+glacier    = getfieldvalue(options,'glacier','Jakobshavn');
 
-function dataout = interpFromMEaSUREsGeotiff(X,Y,Tstart,Tend)
-
-foldername = '/totten_1/ModelData/Greenland/VelMEaSUREs/Jakobshavn_2008_2021/';
+if strcmp(glacier, 'Jakobshavn')
+        foldername = '/totten_1/ModelData/Greenland/VelMEaSUREs/Jakobshavn_2008_2021/';
+elseif strcmp(glacier, 'Kangerlussuaq')
+        foldername = '/totten_1/ModelData/Greenland/VelMEaSUREs/Kangerlussuaq_2006_2021/';
+elseif strcmp(glacier, 'Store')
+        foldername = '/totten_1/ModelData/Greenland/VelMEaSUREs/Store_2008_2021/';
+elseif strcmp(glacier, 'Rink')
+        foldername = '/totten_1/ModelData/Greenland/VelMEaSUREs/Rink_2008_2022/';
+else
+        error(['The velocity data for ', glacier, ' is not available, please download from NSIDC first.']);
+end
 
 % get the time info from file names