% This function provides a Gaussian vorticity distrubution
% which will make a vortex with the specified vmax and rmw


function zeta=gaussianvort(r,vmax,rmw);

zeta=zeros(size(r));

%a = (3.51296*vmax)/rmw;
%a = (3.518*vmax)/rmw;
%b = rmw/1.1209;
%b = rmw/1.12;

%clear;clc;clf;VfromVort;diff = vmax - (max(vt).');test = abs(diff);max(test)


%Nolan and Montgomery 2002 parameters
a = 1.5*10.^(-3);
b = 45*1000;

zeta = a*exp( -(r./b).^2 );
