lat2=rad2deg(0.33181)
long2= rad2deg(-1.4144)
az=0:10:350;
rge=100:25:400
u=-20*cos(deg2rad(az));
u=ones(13,1)*u;
v=20*sin(deg2rad(az));
v=ones(13,1)*v;
mod_az=az;
w=find(mod_az >=270 & mod_az <360);
a=mod_az;
a(w)=360-a(w)+90;
w=find(mod_az>=0 & mod_az<270);
a(w)=90-a(w);
mod_az=a;
[th,r]=meshgrid(deg2rad(mod_az),rge);
[x,y]=pol2cart(th,r);
quiver(x,y,u,v)
for i=1:13
for j=1:36
[lat_b(i,j),long_b(i,j)]=reckon(lat2,long2,distdim(rge(i),'km','deg'),az(j));
end
quiver(long_b,lat_b,u,v)
end
for i=1:13
for j=1:36
[lat_b(i,j),long_b(i,j)]=reckon(lat2,long2,distdim(rge(i),'km','deg'),az(j));
end
quiver(long_b,lat_b,u,v)
axis equal
hold on
plot(long2,lat2,'+r')
lat3=reshape(lat_b,468,1);
long3=reshape(long_b,468,1);
vr=data(:,5);
pcolor(lat3,long3,vr)
vr=reshape(data(:,5),13,36);
figure;pcolor(lat_b,long_b,vr)
contourf(long_b,lat_b,vr,30)
axis equal
contourf(long_b,lat_b,vr,30)
shading flat
hold on
lat1 = rad2deg(0.32)
long1 = rad2deg(-1.44)
plot(lat1,long1,'or')
quiver(long_b,lat_b,u,v)
axis equal
hold on
plot(long1,lat1,'or')
plot(long2,lat2,'+r')
contourf(long_b,lat_b,vr,30)
quiver(long_b,lat_b,u,v)
axis equal
shading flat
plot(long1,lat1,'or')