/[MITgcm]/MITgcm/verification/lab_sea/matlab/myquiver.m
ViewVC logotype

Contents of /MITgcm/verification/lab_sea/matlab/myquiver.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.4 - (show annotations) (download)
Fri May 23 20:19:16 2003 UTC (20 years, 10 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint62v, checkpoint57m_post, checkpoint52l_pre, checkpoint62u, checkpoint57g_pre, checkpoint62t, checkpoint57s_post, checkpoint58b_post, checkpoint57b_post, checkpoint52d_pre, checkpoint57g_post, checkpoint56b_post, checkpoint57y_post, checkpoint52j_pre, checkpoint51o_pre, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint54d_post, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint54e_post, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint62c, checkpoint51l_post, checkpoint57r_post, checkpoint57d_post, checkpoint57i_post, checkpoint52k_post, checkpoint59, checkpoint58, checkpoint55, checkpoint54, checkpoint57, checkpoint56, checkpoint51, checkpoint53, checkpoint52, checkpoint58f_post, checkpoint52f_post, checkpoint57n_post, checkpoint58d_post, checkpoint62s, checkpoint58a_post, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint57z_post, checkpoint54f_post, checkpoint51f_post, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62w, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint58y_post, checkpoint51d_post, checkpoint51t_post, checkpoint58t_post, checkpoint51n_post, checkpoint55i_post, checkpoint58m_post, checkpoint57l_post, checkpoint52i_pre, checkpoint51s_post, checkpoint57t_post, checkpoint55c_post, checkpoint51j_post, checkpoint63g, checkpoint52e_pre, checkpoint57v_post, checkpoint57f_post, checkpoint52e_post, checkpoint51n_pre, checkpoint53d_post, checkpoint64, checkpoint65, checkpoint60, checkpoint61, checkpoint62, checkpoint63, checkpoint57a_post, checkpoint57h_pre, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint52b_pre, checkpoint54b_post, checkpoint58w_post, checkpoint57h_post, checkpoint51l_pre, checkpoint52m_post, checkpoint57y_pre, checkpoint55g_post, checkpoint51q_post, checkpoint51b_pre, checkpoint52b_post, checkpoint52c_post, checkpoint51h_pre, checkpoint58o_post, checkpoint57c_post, checkpoint58p_post, checkpoint58q_post, checkpoint52f_pre, checkpoint55d_post, checkpoint58e_post, checkpoint54a_pre, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint53c_post, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint55d_pre, checkpoint57c_pre, checkpoint58r_post, checkpoint55j_post, branchpoint-genmake2, checkpoint54a_post, checkpoint55h_post, checkpoint58n_post, checkpoint51r_post, checkpoint51i_post, checkpoint57e_post, checkpoint55b_post, checkpoint51b_post, checkpoint51c_post, checkpoint53a_post, checkpoint65o, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint55f_post, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint52d_post, checkpoint53g_post, checkpoint57p_post, checkpint57u_post, checkpoint50g_post, checkpoint57q_post, eckpoint57e_pre, checkpoint58k_post, checkpoint52a_pre, checkpoint62b, checkpoint58v_post, checkpoint50h_post, checkpoint52i_post, checkpoint50i_post, checkpoint51i_pre, checkpoint52h_pre, checkpoint56a_post, checkpoint64y, checkpoint64x, checkpoint58l_post, checkpoint64z, checkpoint53f_post, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint57h_done, checkpoint52j_post, checkpoint57j_post, checkpoint57f_pre, checkpoint61f, checkpoint58g_post, branch-netcdf, checkpoint52l_post, checkpoint58x_post, checkpoint61n, checkpoint52n_post, checkpoint53b_pre, checkpoint59j, checkpoint58h_post, checkpoint56c_post, checkpoint58j_post, checkpoint51e_post, checkpoint57a_pre, checkpoint55a_post, checkpoint57o_post, checkpoint51o_post, checkpoint61q, checkpoint57k_post, checkpoint51f_pre, checkpoint53b_post, checkpoint52a_post, checkpoint57w_post, checkpoint61e, checkpoint58i_post, checkpoint51g_post, ecco_c52_e35, checkpoint57x_post, checkpoint58c_post, checkpoint58u_post, checkpoint51m_post, checkpoint53d_pre, checkpoint58s_post, checkpoint55e_post, checkpoint61g, checkpoint61d, checkpoint54c_post, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint51a_post, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint51p_post, checkpoint61z, checkpoint61x, checkpoint61y, checkpoint51u_post, HEAD
Branch point for: branch-genmake2, branch-nonh, tg2-branch, netcdf-sm0, checkpoint51n_branch
Changes since 1.3: +0 -0 lines
checkpoint50g_post
o merged with release1_p17 (pkg/seaice and verification/lab_sea)
  - added SEAICE_MULTILEVEL for 8-category sea-ice thermodynamics
  - LSR sea-ice dynamic solver moved to SouthWest B-grid location and
    made the default because of faster convergence than ADI

1 function hh = myquiver(arg1,arg2,arg3,arg4,arg5,arg6)
2 %MYQUIVER Quiver (or velocity) plot.
3 % MYQUIVER(X,Y,U,V) plots the velocity vectors with components
4 % (u,v) at the points (x,y). The matrices X,Y,U,V must all be
5 % the same size and contain the cooresponding position and
6 % velocity components (X and Y can also be vectors to specify a
7 % uniform grid). MYQUIVER automatically scales the velocity
8 % vectors to fit within the grid.
9 %
10 % MYQUIVER(U,V) plots the velocity vectors at equally spaced
11 % points in the x-y plane.
12 %
13 % MYQUIVER(X,Y,S) or MYQUIVER(X,Y,U,V,S,...) automatically scales
14 % the velocity vectors to fit within the grid and then multiplies
15 % them by S. Use S=0 to plot the velocity vectors without the
16 % automatic scaling.
17 %
18 % MYQUIVER(...,STYLE) uses the plot linestyle specified by the string
19 % STYLE for the velocity vectors. See PLOT for other linestyles.
20 %
21 % H = MYQUIVER(...) returns a vector of line handles.
22 %
23 % Example:
24 % [x,y] = meshgrid(-2:.2:2,-1:.15:1);
25 % z = x .* exp(-x.^2 - y.^2); [px,py] = gradient(z,.2,.15);
26 % contour(x,y,z), hold on
27 % myquiver(x,y,px,py), hold off, axis image
28 %
29 % See also: FEATHER, PLOT, QUIVER
30
31 % Clay M. Thompson 3-3-94
32 % Copyright (c) 1994 by The MathWorks, Inc.
33 % $Revision: 1.1.2.1 $
34
35 % modified D Menemenlis 28 mar 95
36 % modified from QUIVER for more consistent scaling
37 % search for "modified" in text for details
38
39 error(nargchk(2,6,nargin));
40
41 % Arrow head parameters
42 alpha = 0.33; % Size of arrow head relative to the length of the vector
43 beta = 0.33; % Width of the base of the arrow head relative to the length
44 autoscale = 1; % Autoscale if ~= 0 then scale by this.
45
46 % Check numeric input arguments
47 if nargin<4
48 [msg,x,y,u,v] = xyzchk(arg1,arg2);
49 elseif nargin==4
50 if isstr(arg4)
51 [msg,x,y,u,v] = xyzchk(arg1,arg2);
52 else
53 [msg,x,y,u,v] = xyzchk(arg1,arg2,arg3,arg4);
54 end
55 else
56 [msg,x,y,u,v] = xyzchk(arg1,arg2,arg3,arg4);
57 end
58 if ~isempty(msg), error(msg); end
59
60 if nargin==2, % myquiver(u,v)
61 lo = get(gca,'LineStyleOrder'); sym = lo(1,:);
62 elseif nargin==3, % myquiver(u,v,s) or myquiver(u,v,'style')
63 if isstr(arg3),
64 sym = arg3;
65 else
66 autoscale = arg3;
67 lo = get(gca,'LineStyleOrder'); sym = lo(1,:);
68 end
69 elseif nargin==4, % myquiver(x,y,u,v) or myquiver(x,y,s,'style')
70 if isstr(arg4),
71 autoscale = arg3;
72 sym = arg4;
73 else
74 lo = get(gca,'LineStyleOrder'); sym = lo(1,:);
75 end
76 elseif nargin==5, % myquiver(x,y,u,v,s) or myquiver(x,y,u,v,'style')
77 if isstr(arg5),
78 sym = arg5;
79 else
80 autoscale = arg5;
81 lo = get(gca,'LineStyleOrder'); sym = lo(1,:);
82 end
83 elseif nargin==6, % myquiver(x,y,u,v,s,style)
84 autoscale = arg5;
85 sym = arg6;
86 end
87
88 if autoscale,
89 % Base autoscale value on average spacing in the x and y
90 % directions. Estimate number of points in each direction as
91 % either the size of the input arrays or the effective square
92 % spacing if x and y are vectors.
93 if min(size(x))==1, n=sqrt(prod(size(x))); m=n; else [m,n]=size(x); end
94
95 % modified D Menemenlis 28 mar 95
96 % replaced "n" and "m" by "(n-1)" and "(m-1)" for more consistent scaling as
97 % size of domain is increased
98 delx = diff([min(x(find(~isnan(x)))) max(x(find(~isnan(x))))])/(n-1);
99 dely = diff([min(y(find(~isnan(y)))) max(y(find(~isnan(y))))])/(m-1);
100
101 len = sqrt((u/delx).^2 + (v/dely).^2);
102 autoscale = autoscale*0.9 / max(len(find(~isnan(len))));
103 u = u*autoscale; v = v*autoscale;
104 end
105
106 ax = newplot;
107 next = lower(get(ax,'NextPlot'));
108 hold_state = ishold;
109
110 % Make velocity vectors
111 x = x(:).'; y = y(:).';
112 u = u(:).'; v = v(:).';
113
114 % modified D Menemenlis 28 mar 95
115 % added "+0*u" and "+0*v" to remove dots on plot for nans
116
117 uu = [x+0*u;x+u;NaN*ones(size(u))];
118 vv = [y+0*v;y+v;NaN*ones(size(u))];
119
120 h = plot(uu(:),vv(:),sym);
121
122 % Make arrow heads and plot them
123 hu = [x+u-alpha*(u+beta*(v+eps));x+u; ...
124 x+u-alpha*(u-beta*(v+eps));NaN*ones(size(u))];
125 hv = [y+v-alpha*(v-beta*(u+eps));y+v; ...
126 y+v-alpha*(v+beta*(u+eps));NaN*ones(size(v))];
127 hold on
128 h = [h;plot(hu(:),hv(:),sym)];
129
130 if ~hold_state, hold off, view(2); set(ax,'NextPlot',next); end
131
132 if nargout>0, hh = h; end

  ViewVC Help
Powered by ViewVC 1.1.22