quiverwcolorbar adds functionality to the original quiver function by assigning color to the arrows according to the vector magnitude. The corresponding colorbar spans the minimum and maximum values of the dataset or can be specified by the user. Any changes to the colorbar boundaries change the vector colors accordingly.
Syntax:
quiverwcolorbar(x,y,u,v,scale,'bounds',[colormin colormax])

Inputs:
x, y: Vector locations
u, v: Vector directions (East-West, North-South)
scale: Scalar value to set vector lengths
'bounds': Optional setting to specify colorbar axis minimum and maximum values

Example:
x = rand(1,50).*100;
y = rand(1,50).*100;
u = rand(1,50) .* 10;
v = rand(1,50) .* 10;
scale = 0;
figure; quiverwcolorbar(x',y',u',v',scale); %compare to:
figure; quiverwcolorbar(x',y',u',v',scale,'bounds',[0 10]);