|
Revision 495
(checked in by cbc, 12 years ago)
|
Initial import of Stark code.
|
| Line | |
|---|
| 1 |
% |
|---|
| 2 |
% DIVG - Compute the divergence of a FEM 2-D vector field |
|---|
| 3 |
% |
|---|
| 4 |
% DIVG(e,x,y,u,v) computes the divergence of a 2-D |
|---|
| 5 |
% vector field (u,v) over the FEM domain specified by the |
|---|
| 6 |
% element list e and the corresponding horizontal node |
|---|
| 7 |
% coordinates (x,y). The result is a 1-D scalar field |
|---|
| 8 |
% returned to the workspace. |
|---|
| 9 |
% |
|---|
| 10 |
% Call as: dv=divg(e,x,y,u,v); |
|---|
| 11 |
% |
|---|
| 12 |
% All arguments are REQUIRED. |
|---|
| 13 |
% |
|---|
| 14 |
% WRitten by:: |
|---|
| 15 |
% Brian O. Blanton |
|---|
| 16 |
% Curr. in Marine Sciences |
|---|
| 17 |
% 15-1A Venable Hall |
|---|
| 18 |
% CB# 3300 |
|---|
| 19 |
% Uni. of North Carolina |
|---|
| 20 |
% Chapel Hill, NC |
|---|
| 21 |
% 27599-3300 |
|---|
| 22 |
% |
|---|
| 23 |
% 919-962-4466 |
|---|
| 24 |
% blanton@marine.unc.edu |
|---|
| 25 |
% |
|---|
| 26 |
function dv=divg(e,x,y,u,v) |
|---|
| 27 |
|
|---|
| 28 |
if nargout~=1 |
|---|
| 29 |
error('DIVG must have 1 (and only 1) output argument.') |
|---|
| 30 |
else |
|---|
| 31 |
dv=divgmex5(x,y,e,u,v); |
|---|
| 32 |
end |
|---|
