Changeset 2296
- Timestamp:
- 09/23/09 15:21:29 (15 years ago)
- Location:
- issm/trunk/src/m
- Files:
-
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/classes/public/SectionValues.m
r1181 r2296 72 72 73 73 %Interpolation of data on specified points 74 data_interp=InterpFromMesh 2d(md.elements,md.x,md.y,data,X,Y,NaN);74 data_interp=InterpFromMeshToMesh2d(md.elements,md.x,md.y,data,X,Y,NaN); 75 75 76 76 %Compute index … … 83 83 %Get bed and surface for each 2d point, offset to make sure that it is inside the glacier system 84 84 offset=10^-10; 85 bed=InterpFromMesh 2d(md.elements2d,md.x2d,md.y2d,project2d(md,md.bed,1),X,Y,NaN)+offset;86 surface=InterpFromMesh 2d(md.elements2d,md.x2d,md.y2d,project2d(md,md.surface,1),X,Y,NaN)-offset;85 bed=InterpFromMeshToMesh2d(md.elements2d,md.x2d,md.y2d,project2d(md,md.bed,1),X,Y,NaN)+offset; 86 surface=InterpFromMeshToMesh2d(md.elements2d,md.x2d,md.y2d,project2d(md,md.surface,1),X,Y,NaN)-offset; 87 87 88 88 %Some useful parameters … … 109 109 110 110 %Interpolation of data on specified points 111 data_interp=InterpFromMesh 3d(md.elements,md.x,md.y,md.z,data,X3,Y3,Z3,NaN);111 data_interp=InterpFromMeshToMesh3d(md.elements,md.x,md.y,md.z,data,X3,Y3,Z3,NaN); 112 112 113 113 %build outputs -
issm/trunk/src/m/classes/public/mesh/meshyams.m
r2173 r2296 54 54 %interpolate velocities onto mesh 55 55 disp(' interpolating velocities...'); 56 vx_obs=InterpFromGrid (x_m,y_m,vx,md.x,md.y,0);57 vy_obs=InterpFromGrid (x_m,y_m,vy,md.x,md.y,0);56 vx_obs=InterpFromGridToMesh(x_m,y_m,vx,md.x,md.y,0); 57 vy_obs=InterpFromGridToMesh(x_m,y_m,vy,md.x,md.y,0); 58 58 field=sqrt(vx_obs.^2+vy_obs.^2); 59 59 … … 96 96 md.gridonwater=zeros(md.numberofgrids,1); 97 97 end 98 md.vx_obs=InterpFromGrid (x_m,y_m,vx,md.x,md.y,0);99 md.vy_obs=InterpFromGrid (x_m,y_m,vy,md.x,md.y,0);98 md.vx_obs=InterpFromGridToMesh(x_m,y_m,vx,md.x,md.y,0); 99 md.vy_obs=InterpFromGridToMesh(x_m,y_m,vy,md.x,md.y,0); 100 100 md.vel_obs=sqrt(md.vx_obs.^2+md.vy_obs.^2); -
issm/trunk/src/m/classes/public/modeldefault/defaultparams.m
r2216 r2296 101 101 disp(' reading geothermal flux'); 102 102 load(heatfluxpath); 103 md.geothermalflux=InterpFromGrid (x_m,y_m,heatflux_Antarctica,md.x,md.y,80);103 md.geothermalflux=InterpFromGridToMesh(x_m,y_m,heatflux_Antarctica,md.x,md.y,80); 104 104 pos=find(md.geothermalflux==0);md.geothermalflux(pos)=80; 105 105 md.geothermalflux=md.geothermalflux/1000; %map is given in mW/m^2, we need it in W/m^2 -
issm/trunk/src/m/classes/public/plugvelocities.m
r1411 r2296 94 94 95 95 %interpolate 96 md.vx_obs=InterpFromGrid (x,y,vx,md.x,md.y,default_value);97 md.vy_obs=InterpFromGrid (x,y,vy,md.x,md.y,default_value);96 md.vx_obs=InterpFromGridToMesh(x,y,vx,md.x,md.y,default_value); 97 md.vy_obs=InterpFromGridToMesh(x,y,vy,md.x,md.y,default_value); 98 98 md.vx=md.vx_obs; 99 99 md.vy=md.vy_obs; -
issm/trunk/src/m/classes/public/plugvelocitiesraw.m
r1318 r2296 94 94 95 95 %interpolate 96 md.vx_obs_raw=InterpFromGrid (x,y,vx,md.x,md.y,default_value);97 md.vy_obs_raw=InterpFromGrid (x,y,vy,md.x,md.y,default_value);96 md.vx_obs_raw=InterpFromGridToMesh(x,y,vx,md.x,md.y,default_value); 97 md.vy_obs_raw=InterpFromGridToMesh(x,y,vy,md.x,md.y,default_value); 98 98 md.vel_obs_raw=sqrt(md.vx_obs_raw.^2+md.vy_obs_raw.^2); -
issm/trunk/src/m/utils/Interp/InterpFromFile.m
r1183 r2296 198 198 199 199 %interpolate 200 data_out=InterpFromMesh 2d(index_data,x_data(:),y_data(:),data(:),x(:),y(:),default_value);200 data_out=InterpFromMeshToMesh2d(index_data,x_data(:),y_data(:),data(:),x(:),y(:),default_value); 201 201 202 202 else … … 216 216 217 217 %interpolate 218 data_out=InterpFromGrid (x_data(:),y_data(:),data,x(:),y(:),default_value);219 end 218 data_out=InterpFromGridToMesh(x_data(:),y_data(:),data,x(:),y(:),default_value); 219 end -
issm/trunk/src/m/utils/LatLong/stereomap.m
r2271 r2296 23 23 %ok, we have the transformed coordinates in the current projection, go pick up the corresponding 24 24 %values. 25 new_data_line=InterpFromGrid (x_m,y_m,data,new_x,new_y,NaN);25 new_data_line=InterpFromGridToMesh(x_m,y_m,data,new_x,new_y,NaN); 26 26 27 27 new_data=zeros(M,N);
Note:
See TracChangeset
for help on using the changeset viewer.