vpl_contour.m

来自「ＷＡＡＳ ｍａｔｌａｂ仿真程序 ｍａａｓｔ」· M 代码 · 共 49 行

function vpl_contour(lats, lons, vpl,percent)
%*************************************************************************
%*     Copyright c 2001 The board of trustees of the Leland Stanford     *
%*                      Junior University. All rights reserved.          *
%*     This script file may be distributed and used freely, provided     *
%*     this copyright notice is always kept with it.                     *
%*                                                                       *
%*     Questions and comments should be directed to Todd Walter at:      *
%*     [e-mail]                                              *
%*************************************************************************
%
%VPL_CONTOUR contour plot of the vertical protection level
%   VPL_CONTOUR(LATS, LONS, VPL)
%   Given the user latitudes in LATS (n_lats,1) and longitudes in LONS 
%   (n_lons,1) both in degrees between -180 and 180, and a vector of vpls in
%   VPL(n_lats*n_lons,1) in meters, this function will plot the contour lines
%   according to intervals specified in INIT_GRAPH using the color scheme
%   specified there as well.  The colors can be changed afterward with calls
%   to COLORMAP.
%
%   SEE ALSO SVM_CONTOUR INIT_GRAPH

%  Created by Todd Walter 02 Apr 2001

global GRAPH_VPL_CONTOURS GRAPH_VPL_COLORS

n_lats=length(lats);
n_lons=length(lons);
n_levels=size(GRAPH_VPL_CONTOURS,2);

temp=num2str(GRAPH_VPL_CONTOURS');
ticklabels(n_levels,:)=['> ' temp(n_levels,:)];
for idx=1:n_levels-1
  ticklabels(idx,:)=['< ' temp(idx+1,:)];
end


clf
if nargin > 3
  bartext = ['VPL (m) - ' num2str(percent*100,2) '%'];
else
  bartext = 'VPL (m)';
end
svm_contour(lons,lats,reshape(vpl,n_lons,n_lats)', GRAPH_VPL_CONTOURS, ...
            ticklabels, GRAPH_VPL_COLORS, bartext);


title('VPL as a function of user location', 'FontSize', 14);