div_curl_2d Subroutine

public subroutine div_curl_2d(r, t, ur, vt, divr, curl, undef)

Calculation of rotation and divergence from radial and tangential winds
divr = $\dfrac{\partial ru_r}{r\partial r} + \dfrac{\partial v_t}{r\partial \theta}$
curl = $\dfrac{\partial rv_t}{r\partial r} - \dfrac{\partial u_r}{r\partial \theta}$

Arguments

Type	Intent	Optional	Attributes		Name
double precision,	intent(in)			::	r(:)	radius [m]
double precision,	intent(in)			::	t(:)	angle [rad]
double precision,	intent(in)			::	ur(size(r),size(t))	Ur [m/s]
double precision,	intent(in)			::	vt(size(r),size(t))	Vt [m/s]
double precision,	intent(out)			::	divr(size(r),size(t))	divergence [1/s]
double precision,	intent(out)			::	curl(size(r),size(t))	rotation [1/s]
double precision,	intent(in),	optional		::	undef	Undefined value