cycling/quick_mesh.m

   1 function [ ] = quick_mesh(varargin)
   2 %msg = quick_mesh(matrix)
   3 %   function meshes only a subset of a matrix, for quick handling
   4
   5 if nargin == 1
   6     matrix = varargin{1};
   7     figure;
   8     mesh(matrix(1:5:end,1:5:end));
   9     return
  10 end
  11
  12
  13 if nargin == 2
  14     matrix = varargin{1};
  15     fig_num = varargin{2};
  16     figure(fig_num);
  17     mesh(matrix(1:10:end,1:10:end));
  18     return
  19 end
  20
  21
  22 if nargin == 3
  23     skip = 7;
  24     figure
  25     x =  varargin{1}; % x,y = meshgrid
  26     y =  varargin{2};
  27     z =  varargin{3}; % same size as x and y
  28     lons = x(1,1:skip:end);
  29     lats = y(1:skip:end,1);
  30     small_z = z(1:skip:end,1:skip:end);
  31     [m,n] = size(small_z);
  32     %c = lons.*lats;
  33     CO = zeros(m,n,3);
  34     CO(:,:,1) = zeros(m,n); % red
  35     CO(:,:,2) = small_z/max(small_z(:));
  36     %CO(:,:,2) = ones(m,1).*linspace(0.5,0.6,m); % green
  37     %CO(:,:,3) = ones(m,1).*linspace(0,1,m); % blue
  38     CO(:,:,3) = small_z/max(small_z(:))/2;
  39     c = small_z;
  40     %mesh(lons,lats,small_z,CO)
  41     mesh(lons,lats,small_z,'FaceAlpha',0.5,'EdgeColor','interp','FaceColor','interp','LineStyle',':')
  42     xlabel('Lon (degrees)')
  43     ylabel('Lat (degrees)')
  44 %    figure
  45     hold on
  46     z_level = [1.7279e+05 1.7279e+05];
  47     contour3(lons,lats,small_z,z_level,'k')
  48     hold off
  49 %     xlabel('Lon (degrees)')
  50 %     ylabel('Lat (degrees)')
  51
  52
  53 end
  54
  55
  56 %msg = 'mesh(matrix(1:10:end,1:10:end))'
  57
  58 end
  59