Changeset 27277
- Timestamp:
- 09/10/22 06:13:08 (3 years ago)
- Location:
- issm/trunk-jpl/src/m/contrib/morlighem/modeldata
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpBedmachineGreenland.m
r27276 r27277 45 45 ['/Users/larour/ModelData/BedMachine/' basename '-' ncdate '.nc'],... 46 46 ['./' basename '-' ncdate '.nc'],... 47 '/media/vincent/TOSH4TB/GeorgiaTech/DataSearch/BedMachine/BedMachineGreenland-2021-04-20.nc',...48 47 }; 49 48 … … 90 89 if strcmp(string,'mask') | strcmp(string,'source'), 91 90 %Need nearest neighbor to avoid interpolation between 0 and 2 92 %output = InterpFromGrid(xdata,ydata,data,double(X),double(Y),'nearest'); 93 output = InterpFromGridToMesh(xdata,flipud(ydata),flipud(data),double(X),double(Y),NaN); %VV 91 output = InterpFromGrid(xdata,ydata,data,double(X),double(Y),'nearest'); 94 92 else 95 %output = InterpFromGrid(xdata,ydata,data,double(X),double(Y)); 96 output = InterpFromGridToMesh(xdata,flipud(ydata),flipud(data),double(X),double(Y),NaN); %VV 93 output = InterpFromGrid(xdata,ydata,data,double(X),double(Y)); 97 94 end 98 95 -
issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpSeaRISE.m
r27276 r27277 58 58 searisenc='/totten_1/ModelData/SeaRISE/Antarctica_5km_dev1.0.nc'; 59 59 end 60 case {'vincentxps159500'}61 if hemisphere==162 searisenc = '/media/vincent/TOSH4TB/GeorgiaTech/ISSM/fromSeroussi/GrIS04Feb2022/seariseData/Greenland_5km_dev1.2.nc';63 end64 60 otherwise 65 61 error('hostname not supported yet'); … … 84 80 if verbose, disp([' -- SeaRISE: interpolating ' string]); end 85 81 if strcmpi(string,'LandMask'); 86 %dataout = InterpFromGrid(xdata,ydata,data,xproj,yproj,'nearest'); 87 dataout = InterpFromGridToMesh(xdata,ydata,data,xproj,yproj,NaN); %VV 82 dataout = InterpFromGrid(xdata,ydata,data,xproj,yproj,'nearest'); 88 83 else 89 %dataout = InterpFromGrid(xdata,ydata,data,xproj,yproj); 90 dataout = InterpFromGridToMesh(xdata,ydata,data,xproj,yproj,NaN); %VV 84 dataout = InterpFromGrid(xdata,ydata,data,xproj,yproj); 91 85 end 92 86 dataout = reshape(dataout,size(X,1),size(X,2));
Note:
See TracChangeset
for help on using the changeset viewer.