Ignore:
Timestamp:
03/28/14 09:15:19 (11 years ago)
Author:
Mathieu Morlighem
Message:

CHG: WARNING: md.geometry.bed is now md.geometry.base and bathymetry is bed

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/test/Par/ISMIPE.par

    r16167 r17590  
    44data=transpose(ncread('../Data/ISMIPE.nc','data'));
    55md.geometry.surface=zeros(md.mesh.numberofvertices,1);
    6 md.geometry.bed=zeros(md.mesh.numberofvertices,1);
     6md.geometry.base=zeros(md.mesh.numberofvertices,1);
    77for i=1:md.mesh.numberofvertices
    88        y=md.mesh.y(i);
     
    1010        point2=min(point1+1,51);
    1111        coeff=(y-(point1-1.)*100.)/100.;
    12         md.geometry.bed(i)=(1.-coeff)*data(point1,2)+coeff*data(point2,2);
     12        md.geometry.base(i)=(1.-coeff)*data(point1,2)+coeff*data(point2,2);
    1313        md.geometry.surface(i)=(1.-coeff)*data(point1,3)+coeff*data(point2,3);
    1414end
    15 md.geometry.thickness=md.geometry.surface-md.geometry.bed;
     15md.geometry.thickness=md.geometry.surface-md.geometry.base;
    1616md.geometry.thickness(find(~md.geometry.thickness))=0.01;
    17 md.geometry.bed=md.geometry.surface-md.geometry.thickness;
     17md.geometry.base=md.geometry.surface-md.geometry.thickness;
    1818
    1919disp('      creating drag');
Note: See TracChangeset for help on using the changeset viewer.