


GETB_VERT - compute the B matrix for vertical velocity of diagnostic
Compute B matrix. B=[dh1/dz dh2/dz dh3/dz dh4/dz dh5/dz dh6/dz];
where hi is the interpolation function for grid i.
Usage:
B=GetB_vert(pentaelem,xyz_list,gauss_coord)
See also GETBPRIME_VERT, CREATEKMATRIX

0001 function B=GetB_vert(pentaelem,xyz_list,gauss_coord) 0002 %GETB_VERT - compute the B matrix for vertical velocity of diagnostic 0003 % 0004 % Compute B matrix. B=[dh1/dz dh2/dz dh3/dz dh4/dz dh5/dz dh6/dz]; 0005 % where hi is the interpolation function for grid i. 0006 % 0007 % Usage: 0008 % B=GetB_vert(pentaelem,xyz_list,gauss_coord) 0009 % 0010 % See also GETBPRIME_VERT, CREATEKMATRIX 0011 0012 dh1dh6_basic=GetNodalFunctionsDerivativesBasic(pentaelem,xyz_list,gauss_coord); 0013 %Build B: 0014 B=zeros(1,6); 0015 for i=1:6, 0016 B(1,i)=dh1dh6_basic(3,i); 0017 end