Changeset 27517

Timestamp:

01/17/23 06:13:52 (3 years ago)

Author:

Mathieu Morlighem

Message:

CHG: upating doc

Location:

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

Files:

• interpAdusumilliIceShelfMelt.m (modified) (1 diff)
• interpRignotIceShelfMelt.m (modified) (1 diff)

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

@@ -1 +1 @@
 function output = interpAdusumilliIceShelfMelt(X,Y)
-%INTERPADUSUMILLIICESHELFMELT imports basal melt rates from (Adusumilli et al., 2020).
-%       About the data: "Average basal melt rates for Antarctic ice shelves for the 2010–2018 period at 
-%       high spatial resolution, estimated using CryoSat-2 data. This data file was last updated on 2020-06-11."
+%INTERPADUSUMILLIICESHELFMELT - imports basal melt rates from (Adusumilli et al., 2020).
+%   About the data: "Average basal melt rates for Antarctic ice shelves for the 2010–2018 period at 
+%   high spatial resolution, estimated using CryoSat-2 data. This data file was last updated on 2020-06-11."
 %
-%       Citation: Adusumilli, Susheel; Fricker, Helen A.; Medley, Brooke C.; Padman, Laurie; Siegfried, Matthew R. (2020). 
-%               Data from: Interannual variations in meltwater input to the Southern Ocean from Antarctic ice shelves. 
-%               UC San Diego Library Digital Collections. https://doi.org/10.6075/J04Q7SHT
+%   Citation: Adusumilli, Susheel; Fricker, Helen A.; Medley, Brooke C.; Padman, Laurie; Siegfried, Matthew R. (2020). 
+%   Data from: Interannual variations in meltwater input to the Southern Ocean from Antarctic ice shelves. 
+%   UC San Diego Library Digital Collections. https://doi.org/10.6075/J04Q7SHT
 %
-%       Written by Benjamin Getraer Jan 2023
-%       Adapted from interpRignotIceShelfMelt
+%   Usage:
+%      output = interpAdusumilliIceShelfMelt(X,Y)
 
-   % define path and filename for this machine
-   switch (oshostname()),
-      case {'totten'}
-         filename ='/totten_1/ModelData/Antarctica/Adusumilli2020IceShelfMelt/ANT_iceshelf_melt_rates_CS2_2010-2018_v0.h5';
-      otherwise
-         error('hostname not supported yet');
-   end
+% define path and filename for this machine
+switch (oshostname()),
+        case {'totten'}
+                filename ='/totten_1/ModelData/Antarctica/Adusumilli2020IceShelfMelt/ANT_iceshelf_melt_rates_CS2_2010-2018_v0.h5';
+        otherwise
+                error('hostname not supported yet');
+end
 
-        disp(['   -- Adusumilli Ice Shelf Melt: loading melt data']);
-   % read in coordinates:
-        %       coordinates are in Polar Stereographic projection 'PS-71'
-   xdata = double(h5read(filename,'/x'));
-   ydata = double(h5read(filename,'/y'));
+disp(['   -- Adusumilli Ice Shelf Melt: loading melt data']);
+% read in coordinates:
+%       coordinates are in Polar Stereographic projection 'PS-71'
+xdata = double(h5read(filename,'/x'));
+ydata = double(h5read(filename,'/y'));
 
-        % read in data:
-   % 'Basal melt rate (2010–2018), in meters of ice equivalent per year, positive is melting'
-        % 'For ice shelf areas where CryoSat-2 data were not available, w_b_interp provides the 
-        %  mean melt rate measured at the same ice draft as the grid cell elsewhere on the ice shelf. 
-        %  Ice draft was estimated using BedMachine data.'
-   data = double(h5read(filename,'/w_b'));
-        data_interp = double(h5read(filename,'/w_b_interp'));
-        data = data';
-        disp(['   -- Adusumilli Ice Shelf Melt: interpolating melt data']);
-        data(isnan(data)) = data_interp(isnan(data));
-        output = InterpFromGrid(xdata,ydata,data,X(:),Y(:));
-        output = reshape(output,size(X,1),size(X,2));
+% read in data:
+% 'Basal melt rate (2010–2018), in meters of ice equivalent per year, positive is melting'
+% 'For ice shelf areas where CryoSat-2 data were not available, w_b_interp provides the 
+%  mean melt rate measured at the same ice draft as the grid cell elsewhere on the ice shelf. 
+%  Ice draft was estimated using BedMachine data.'
+data = double(h5read(filename,'/w_b'));
+data_interp = double(h5read(filename,'/w_b_interp'));
+data = data';
+disp(['   -- Adusumilli Ice Shelf Melt: interpolating melt data']);
+data(isnan(data)) = data_interp(isnan(data));
+output = InterpFromGrid(xdata,ydata,data,X(:),Y(:));
+output = reshape(output,size(X,1),size(X,2));

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

@@ -1 @@
-function output = interpRignotIceShelfMelt(X,Y,string),
+function output = interpRignotIceShelfMelt(X,Y,string)
+%INTERPRIGNOTICESHELFMELT - interp melt rates from Rignot et al. 2013
+%
+%   Usage:
+%      output = interpRignotIceShelfMelt(X,Y)
 
-switch (oshostname()),
+switch (oshostname())
         case {'ronne'}
                 rignotmelt='/home/ModelData/Antarctica/RignotMeltingrate/Ant_MeltingRate.nc';