Changeset 17590 for issm/trunk-jpl/test/Par/ISMIPE.par
- Timestamp:
- 03/28/14 09:15:19 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/Par/ISMIPE.par
r16167 r17590 4 4 data=transpose(ncread('../Data/ISMIPE.nc','data')); 5 5 md.geometry.surface=zeros(md.mesh.numberofvertices,1); 6 md.geometry.b ed=zeros(md.mesh.numberofvertices,1);6 md.geometry.base=zeros(md.mesh.numberofvertices,1); 7 7 for i=1:md.mesh.numberofvertices 8 8 y=md.mesh.y(i); … … 10 10 point2=min(point1+1,51); 11 11 coeff=(y-(point1-1.)*100.)/100.; 12 md.geometry.b ed(i)=(1.-coeff)*data(point1,2)+coeff*data(point2,2);12 md.geometry.base(i)=(1.-coeff)*data(point1,2)+coeff*data(point2,2); 13 13 md.geometry.surface(i)=(1.-coeff)*data(point1,3)+coeff*data(point2,3); 14 14 end 15 md.geometry.thickness=md.geometry.surface-md.geometry.b ed;15 md.geometry.thickness=md.geometry.surface-md.geometry.base; 16 16 md.geometry.thickness(find(~md.geometry.thickness))=0.01; 17 md.geometry.b ed=md.geometry.surface-md.geometry.thickness;17 md.geometry.base=md.geometry.surface-md.geometry.thickness; 18 18 19 19 disp(' creating drag');
Note:
See TracChangeset
for help on using the changeset viewer.