0001 function wb=basevert(md)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011 alpha=zeros(md.numberofelements,3);
0012 beta=zeros(md.numberofelements,3);
0013 gamma=zeros(md.numberofelements,3);
0014
0015 for n=1:md.numberofelements
0016 X=inv([md.x(md.elements(n,:)) md.y(md.elements(n,:)) ones(3,1)]);
0017 alpha(n,:)=X(1,:);
0018 beta(n,:)=X(2,:);
0019 gamma(n,:)=X(3,:);
0020 end
0021
0022 hu=md.thickness.*md.vx;
0023 hv=md.thickness.*md.vy;
0024
0025 summation=[1;1;1];
0026 hux=(hu(md.elements).*alpha)*summation;
0027 hvy=(hv(md.elements).*beta)*summation;
0028
0029 uelem=md.vx(md.elements)*summation/3;
0030 velem=md.vy(md.elements)*summation/3;
0031
0032 dbdx=(md.bed(md.elements).*alpha)*summation;
0033 dbdy=(md.bed(md.elements).*beta)*summation;
0034
0035 wb=-md.rho_ice/md.rho_water*(hux+hvy)+uelem.*dbdx+velem.*dbdy;