%% dstplot.m plots the dst values for the requested months in the given year
%% of interest.
%%
%% Author: Phillip A. Webb
%% Last updated: January 9 2004                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           5 2003 
%%

clear all
close all

yeari=input('Year of interest? ');
monthv=input('Input months of interest (max of six) (vector)? ');
nop=length(monthv);
if nop<=2
 ny=1;
else
 ny=2;
end
nx=ceil(nop/ny);

eval(['load dst',int2str(yeari)])
eval(['dstm=dst',int2str(yeari),';'])
yearv=dstm(:,1);
monthv=dstm(:,2);
dayv=dstm(:,3);

figure(1)
for i=1:nop
 subplot(ny,nx,i)
 
end 
 
