Changeset 27966


Ignore:
Timestamp:
10/23/23 08:10:03 (17 months ago)
Author:
Cheng Gong
Message:

remove author names and dates in the script

Location:
issm/trunk-jpl/src/m/contrib/chenggong
Files:
13 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/analyzeCalvingFront.m

    r26975 r27966  
    55%                   [Xmain, positionx, positiony, time], icemask, calving rate,
    66%                   melting rate and velocity magnitude, sigma, thickness
    7 % Author: Cheng Gong
    8 % Last modified: 2020-09-25
    97
    108transient = isfield(md.results,'TransientSolution');

  • issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/averageOverTime.m

    r26975 r27966  
    1313% if within the range, do linear interpolation.
    1414%
    15 %   Author: Cheng Gong
    16 %   Last modified: 2020-09-09
    1715if nargin < 4
    1816    endP = startP;

  • issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/fillInNan.m

    r26975 r27966  
    33%               data should have the same size as md.mesh.x
    44%
    5 % Author: Cheng Gong
    6 % Last modified: 2021-12-08
    75nanvflag = find(isnan(data));
    86NNanv = length(nanvflag);

  • issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/interpZeroPos.m

    r26975 r27966  
    1414%                x0 = -------------
    1515%                        y1-y2
    16 % Author: Cheng Gong
    17 % Last modified: 2020-10-20
    1816[m, n] = size(X);
    1917[my, k] = size(Y);

  • issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/npsd.m

    r26975 r27966  
    77%   freq: frequency(non-negative half)
    88%
    9 %   Author: Cheng Gong
    10 %   Date: 2021-08-17
    119
    1210% length of data

  • issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/projectToFlowlines.m

    r26975 r27966  
    66%       fy                      -       y coordinates of the flowline
    77%
    8 % Author: Cheng Gong
    9 % Last modified: 2021-01-27
    108
    119function valueC = projectToFlowlines(md, pValue, fx, fy)

  • issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/psd.m

    r26975 r27966  
    77%   freq: frequency(non-negative half)
    88%
    9 %   Author: Cheng Gong
    10 %   Date: 2021-08-17
    119
    1210% length of data

  • issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/wfDistance.m

    r26975 r27966  
    99%   w2:  Wasserstein-Fourier Distance, or W2 distance of sx and sy
    1010%
    11 %   Author: Cheng Gong
    12 %   Date: 2021-08-18
    1311N = length(x);
    1412Sinv=linspace(0,1,N);

  • issm/trunk-jpl/src/m/contrib/chenggong/modeldata/interpFromITSLIVE.m

    r27901 r27966  
    1 function [vx_out, vy_out, vx_err, vy_err] = interpFromITSLIVE(X,Y,Tstart,Tend,varargin)
     1function [vx_out, vy_out] = interpFromITSLIVE(X,Y,Tstart,Tend,varargin)
    22        %interpFromITSLIVE:
    33        %       Interpolate ITS_LIVE velocity data to the given mesh
    … …  
    2828                dataTime(i) = date2decyear(datenum(tempConv{end}, 'yyyy'));
    2929        end
    30         disp(['  Found ', num2str(Ndata), ' records in ', foldername]);
    31         disp(['    from ', datestr(decyear2date(min(dataTime)),'yyyy-mm-dd'), ' to ', datestr(decyear2date(max(dataTime)),'yyyy-mm-dd') ]);
    32 
    33         % 0000 is a ref velocity
    34 
    3530        % find all the data files with Tstart<=t<=Tend
    3631        dataInd = (dataTime>=Tstart) & (dataTime<=Tend);
    … …  
    6358        vx_out = zeros(numel(X)+1, numel(timeToload));
    6459        vy_out = zeros(numel(X)+1, numel(timeToload));
    65         vx_err = zeros(numel(X)+1, numel(timeToload));
    66         vy_err = zeros(numel(X)+1, numel(timeToload));
    6760        for i = 1:length(dataToLoad)
    6861
    … …  
    7063                vx = (ncread(filename,'vx',[id1xh id1yh],[id2xh-id1xh+1 id2yh-id1yh+1],[1 1]));
    7164                vy = (ncread(filename,'vy',[id1xh id1yh],[id2xh-id1xh+1 id2yh-id1yh+1],[1 1]));
    72                 vxerr = (ncread(filename,'vx_err',[id1xh id1yh],[id2xh-id1xh+1 id2yh-id1yh+1],[1 1]));
    73                 vyerr = (ncread(filename,'vy_err',[id1xh id1yh],[id2xh-id1xh+1 id2yh-id1yh+1],[1 1]));
    7465
    7566                vx(vx<-32760) = nan;
    … …  
    7970                vy_out(1:end-1,i) = InterpFromGrid(xh, yh, double(vy'), X, Y);
    8071                vy_out(end, i) = timeToload(i);
    81                 vx_err(1:end-1,i) = InterpFromGrid(xh, yh, double(vxerr'), X, Y);
    82                 vx_err(end, i) = timeToload(i);
    83                 vy_err(1:end-1,i) = InterpFromGrid(xh, yh, double(vyerr'), X, Y);
    84                 vy_err(end, i) = timeToload(i);
    8572        end
    8673end

  • issm/trunk-jpl/src/m/contrib/chenggong/modeldata/interpISMIP6Temp.m

    r26975 r27966  
    44%       X and Y are the coordinates of the mesh
    55%
    6 %  Author: Cheng Gong
    7 %  Last modified: 2021-12-06
    8 
    96filename = '/totten_1/ModelData/Greenland/ISMIP6/GreenlandISMIP6-Morlighem-2020-10-01.nc';
    107

  • issm/trunk-jpl/src/m/contrib/chenggong/modeldata/interpMonthlyIceMaskGreene.m

    r27694 r27966  
    1313%   - optional 5th input argument: path to the data file, by default it is the path on totten
    1414%
    15 %   Author: Cheng Gong
    16 %   Last modified: 2023-04-19
    1715
    1816% set icemask=-1 for the region with rocks

  • issm/trunk-jpl/src/m/contrib/chenggong/visualization/imageNonUni.m

    r26975 r27966  
    22%imageNonUni - draw an image with non uniform grid
    33%
    4 % Author: Cheng Gong
    5 % Last modified: 2020-08-24
    64x = varargin{1};
    75y = varargin{2};

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

    r27789 r27966  
    1515%   yl:                 ylim value
    1616%
    17 %   Author: Cheng Gong
    18 %   Date: 2021-12-06
    1917
    2018N = length(velList);
Note: See TracChangeset for help on using the changeset viewer.