load 'vlbarain2003' load 'vlbarain2004' load 'vlbarain2005' load 'vlbarain2006' load 'vlbarain2007' load 'vlbarain2008' load 'vlbarain2009' load 'vlbarain2010' vlbarain2003jan = extractjan (vlbarain2003); vlbarain2004jan = extractjan (vlbarain2004); vlbarain2005jan = extractjan (vlbarain2005); vlbarain2006jan = extractjan (vlbarain2006); vlbarain2007jan = extractjan (vlbarain2007); vlbarain2008jan = extractjan (vlbarain2008); vlbarain2009jan = extractjan (vlbarain2009); vlbarain2010jan = extractjan (vlbarain2010); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2003jan','vlbarain2003jan'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2004jan','vlbarain2004jan'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2005jan','vlbarain2005jan'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2006jan','vlbarain2006jan'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2007jan','vlbarain2007jan'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2008jan','vlbarain2008jan'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2009jan','vlbarain2009jan'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2010jan','vlbarain2010jan'); vlbarainjan = [vlbarain2003jan;vlbarain2004jan;vlbarain2005jan;vlbarain2006jan;vlbarain2007jan;vlbarain2008jan;vlbarain2009jan;vlbarain2010jan]; save ('/Users/sarahstuart/Desktop/precipitation report/vlbarainjan','vlbarainjan'); vlbarain2003feb = extractfeb (vlbarain2003); vlbarain2004feb = extractfeb (vlbarain2004); vlbarain2005feb = extractfeb (vlbarain2005); vlbarain2006feb = extractfeb (vlbarain2006); vlbarain2007feb = extractfeb (vlbarain2007); vlbarain2008feb = extractfeb (vlbarain2008); vlbarain2009feb = extractfeb (vlbarain2009); vlbarain2010feb = extractfeb (vlbarain2010); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2003feb','vlbarain2003feb'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2004feb','vlbarain2004feb'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2005feb','vlbarain2005feb'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2006feb','vlbarain2006feb'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2007feb','vlbarain2007feb'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2008feb','vlbarain2008feb'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2009feb','vlbarain2009feb'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2010feb','vlbarain2010feb'); vlbarainfeb = [vlbarain2003feb;vlbarain2004feb;vlbarain2005feb;vlbarain2006feb;vlbarain2007feb;vlbarain2008feb;vlbarain2009feb;vlbarain2010feb]; save ('/Users/sarahstuart/Desktop/precipitation report/vlbarainfeb','vlbarainfeb'); vlbarain2003mar = extractmar (vlbarain2003); vlbarain2004mar = extractmar (vlbarain2004); vlbarain2005mar = extractmar (vlbarain2005); vlbarain2006mar = extractmar (vlbarain2006); vlbarain2007mar = extractmar (vlbarain2007); vlbarain2008mar = extractmar (vlbarain2008); vlbarain2009mar = extractmar (vlbarain2009); vlbarain2010mar = extractmar (vlbarain2010); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2003mar','vlbarain2003mar'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2004mar','vlbarain2004mar'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2005mar','vlbarain2005mar'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2006mar','vlbarain2006mar'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2007mar','vlbarain2007mar'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2008mar','vlbarain2008mar'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2009mar','vlbarain2009mar'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2010mar','vlbarain2010mar'); vlbarainmar = [vlbarain2003mar;vlbarain2004mar;vlbarain2005mar;vlbarain2006mar;vlbarain2007mar;vlbarain2008mar;vlbarain2009mar;vlbarain2010mar]; save ('/Users/sarahstuart/Desktop/precipitation report/vlbarainmar','vlbarainmar'); vlbarain2003apr = extractapr (vlbarain2003); vlbarain2004apr = extractapr (vlbarain2004); vlbarain2005apr = extractapr (vlbarain2005); vlbarain2006apr = extractapr (vlbarain2006); vlbarain2007apr = extractapr (vlbarain2007); vlbarain2008apr = extractapr (vlbarain2008); vlbarain2009apr = extractapr (vlbarain2009); vlbarain2010apr = extractapr (vlbarain2010); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2003apr','vlbarain2003apr'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2004apr','vlbarain2004apr'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2005apr','vlbarain2005apr'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2006apr','vlbarain2006apr'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2007apr','vlbarain2007apr'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2008apr','vlbarain2008apr'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2009apr','vlbarain2009apr'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2010apr','vlbarain2010apr'); vlbarainapr = [vlbarain2003apr;vlbarain2004apr;vlbarain2005apr;vlbarain2006apr;vlbarain2007apr;vlbarain2008apr;vlbarain2009apr;vlbarain2010apr]; save ('/Users/sarahstuart/Desktop/precipitation report/vlbarainapr','vlbarainapr'); vlbarain2003may = extractmay (vlbarain2003); vlbarain2004may = extractmay (vlbarain2004); vlbarain2005may = extractmay (vlbarain2005); vlbarain2006may = extractmay (vlbarain2006); vlbarain2007may = extractmay (vlbarain2007); vlbarain2008may = extractmay (vlbarain2008); vlbarain2009may = extractmay (vlbarain2009); vlbarain2010may = extractmay (vlbarain2010); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2003may','vlbarain2003may'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2004may','vlbarain2004may'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2005may','vlbarain2005may'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2006may','vlbarain2006may'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2007may','vlbarain2007may'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2008may','vlbarain2008may'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2009may','vlbarain2009may'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2010may','vlbarain2010may'); vlbarainmay = [vlbarain2003may;vlbarain2004may;vlbarain2005may;vlbarain2006may;vlbarain2007may;vlbarain2008may;vlbarain2009may;vlbarain2010may]; save ('/Users/sarahstuart/Desktop/precipitation report/vlbarainmay','vlbarainmay'); vlbarain2003jun = extractjun (vlbarain2003); vlbarain2004jun = extractjun (vlbarain2004); vlbarain2005jun = extractjun (vlbarain2005); vlbarain2006jun = extractjun (vlbarain2006); vlbarain2007jun = extractjun (vlbarain2007); vlbarain2008jun = extractjun (vlbarain2008); vlbarain2009jun = extractjun (vlbarain2009); vlbarain2010jun = extractjun (vlbarain2010); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2003jun','vlbarain2003jun'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2004jun','vlbarain2004jun'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2005jun','vlbarain2005jun'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2006jun','vlbarain2006jun'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2007jun','vlbarain2007jun'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2008jun','vlbarain2008jun'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2009jun','vlbarain2009jun'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2010jun','vlbarain2010jun'); vlbarainjun = [vlbarain2003jun;vlbarain2004jun;vlbarain2005jun;vlbarain2006jun;vlbarain2007jun;vlbarain2008jun;vlbarain2009jun;vlbarain2010jun]; save ('/Users/sarahstuart/Desktop/precipitation report/vlbarainjun','vlbarainjun'); vlbarain2003jul = extractjul (vlbarain2003); vlbarain2004jul = extractjul (vlbarain2004); vlbarain2005jul = extractjul (vlbarain2005); vlbarain2006jul = extractjul (vlbarain2006); vlbarain2007jul = extractjul (vlbarain2007); vlbarain2008jul = extractjul (vlbarain2008); vlbarain2009jul = extractjul (vlbarain2009); vlbarain2010jul = extractjul (vlbarain2010); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2003jul','vlbarain2003jul'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2004jul','vlbarain2004jul'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2005jul','vlbarain2005jul'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2006jul','vlbarain2006jul'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2007jul','vlbarain2007jul'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2008jul','vlbarain2008jul'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2009jul','vlbarain2009jul'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2010jul','vlbarain2010jul'); vlbarainjul = [vlbarain2003jul;vlbarain2004jul;vlbarain2005jul;vlbarain2006jul;vlbarain2007jul;vlbarain2008jul;vlbarain2009jul;vlbarain2010jul]; save ('/Users/sarahstuart/Desktop/precipitation report/vlbarainjul','vlbarainjul'); vlbarain2003aug = extractaug (vlbarain2003); vlbarain2004aug = extractaug (vlbarain2004); vlbarain2005aug = extractaug (vlbarain2005); vlbarain2006aug = extractaug (vlbarain2006); vlbarain2007aug = extractaug (vlbarain2007); vlbarain2008aug = extractaug (vlbarain2008); vlbarain2009aug = extractaug (vlbarain2009); vlbarain2010aug = extractaug (vlbarain2010); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2003aug','vlbarain2003aug'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2004aug','vlbarain2004aug'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2005aug','vlbarain2005aug'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2006aug','vlbarain2006aug'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2007aug','vlbarain2007aug'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2008aug','vlbarain2008aug'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2009aug','vlbarain2009aug'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2010aug','vlbarain2010aug'); vlbarainaug = [vlbarain2003aug;vlbarain2004aug;vlbarain2005aug;vlbarain2006aug;vlbarain2007aug;vlbarain2008aug;vlbarain2009aug;vlbarain2010aug]; save ('/Users/sarahstuart/Desktop/precipitation report/vlbarainaug','vlbarainaug'); vlbarain2003sep = extractsep (vlbarain2003); vlbarain2004sep = extractsep (vlbarain2004); vlbarain2005sep = extractsep (vlbarain2005); vlbarain2006sep = extractsep (vlbarain2006); vlbarain2007sep = extractsep (vlbarain2007); vlbarain2008sep = extractsep (vlbarain2008); vlbarain2009sep = extractsep (vlbarain2009); vlbarain2010sep = extractsep (vlbarain2010); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2003sep','vlbarain2003sep'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2004sep','vlbarain2004sep'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2005sep','vlbarain2005sep'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2006sep','vlbarain2006sep'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2007sep','vlbarain2007sep'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2008sep','vlbarain2008sep'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2009sep','vlbarain2009sep'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2010sep','vlbarain2010sep'); vlbarainsep = [vlbarain2003sep;vlbarain2004sep;vlbarain2005sep;vlbarain2006sep;vlbarain2007sep;vlbarain2008sep;vlbarain2009sep;vlbarain2010sep]; save ('/Users/sarahstuart/Desktop/precipitation report/vlbarainsep','vlbarainsep'); vlbarain2003oct = extractoct (vlbarain2003); vlbarain2004oct = extractoct (vlbarain2004); vlbarain2005oct = extractoct (vlbarain2005); vlbarain2006oct = extractoct (vlbarain2006); vlbarain2007oct = extractoct (vlbarain2007); vlbarain2008oct = extractoct (vlbarain2008); vlbarain2009oct = extractoct (vlbarain2009); vlbarain2010oct = extractoct (vlbarain2010); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2003oct','vlbarain2003oct'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2004oct','vlbarain2004oct'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2005oct','vlbarain2005oct'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2006oct','vlbarain2006oct'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2007oct','vlbarain2007oct'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2008oct','vlbarain2008oct'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2009oct','vlbarain2009oct'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2010oct','vlbarain2010oct'); vlbarainoct = [vlbarain2003oct;vlbarain2004oct;vlbarain2005oct;vlbarain2006oct;vlbarain2007oct;vlbarain2008oct;vlbarain2009oct;vlbarain2010oct]; save ('/Users/sarahstuart/Desktop/precipitation report/vlbarainoct','vlbarainoct'); vlbarain2003nov = extractnov (vlbarain2003); vlbarain2004nov = extractnov (vlbarain2004); vlbarain2005nov = extractnov (vlbarain2005); vlbarain2006nov = extractnov (vlbarain2006); vlbarain2007nov = extractnov (vlbarain2007); vlbarain2008nov = extractnov (vlbarain2008); vlbarain2009nov = extractnov (vlbarain2009); vlbarain2010nov = extractnov (vlbarain2010); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2003nov','vlbarain2003nov'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2004nov','vlbarain2004nov'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2005nov','vlbarain2005nov'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2006nov','vlbarain2006nov'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2007nov','vlbarain2007nov'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2008nov','vlbarain2008nov'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2009nov','vlbarain2009nov'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2010nov','vlbarain2010nov'); vlbarainnov = [vlbarain2003nov;vlbarain2004nov;vlbarain2005nov;vlbarain2006nov;vlbarain2007nov;vlbarain2008nov;vlbarain2009nov;vlbarain2010nov]; save ('/Users/sarahstuart/Desktop/precipitation report/vlbarainnov','vlbarainnov'); vlbarain2003dec = extractdec (vlbarain2003); vlbarain2004dec = extractdec (vlbarain2004); vlbarain2005dec = extractdec (vlbarain2005); vlbarain2006dec = extractdec (vlbarain2006); vlbarain2007dec = extractdec (vlbarain2007); vlbarain2008dec = extractdec (vlbarain2008); vlbarain2009dec = extractdec (vlbarain2009); vlbarain2010dec = extractdec (vlbarain2010); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2003dec','vlbarain2003dec'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2004dec','vlbarain2004dec'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2005dec','vlbarain2005dec'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2006dec','vlbarain2006dec'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2007dec','vlbarain2007dec'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2008dec','vlbarain2008dec'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2009dec','vlbarain2009dec'); save ('/Users/sarahstuart/Desktop/precipitation report/vlbarain2010dec','vlbarain2010dec'); vlbaraindec = [vlbarain2003dec;vlbarain2004dec;vlbarain2005dec;vlbarain2006dec;vlbarain2007dec;vlbarain2008dec;vlbarain2009dec;vlbarain2010dec]; save ('/Users/sarahstuart/Desktop/precipitation report/vlbaraindec','vlbaraindec'); vlbarainjan = [vlbarain2003jan;vlbarain2004jan;vlbarain2005jan;vlbarain2006jan;vlbarain2007jan;vlbarain2008jan;vlbarain2009jan;vlbarain2010jan]; vlbarainfeb = [vlbarain2003feb;vlbarain2004feb;vlbarain2005feb;vlbarain2006feb;vlbarain2007feb;vlbarain2008feb;vlbarain2009feb;vlbarain2010feb]; vlbarainmar = [vlbarain2003mar;vlbarain2004mar;vlbarain2005mar;vlbarain2006mar;vlbarain2007mar;vlbarain2008mar;vlbarain2009mar;vlbarain2010mar]; vlbarainapr = [vlbarain2003apr;vlbarain2004apr;vlbarain2005apr;vlbarain2006apr;vlbarain2007apr;vlbarain2008apr;vlbarain2009apr;vlbarain2010apr]; vlbarainmay = [vlbarain2003may;vlbarain2004may;vlbarain2005may;vlbarain2006may;vlbarain2007may;vlbarain2008may;vlbarain2009may;vlbarain2010may]; vlbarainjun = [vlbarain2003jun;vlbarain2004jun;vlbarain2005jun;vlbarain2006jun;vlbarain2007jun;vlbarain2008jun;vlbarain2009jun;vlbarain2010jun]; vlbarainjul = [vlbarain2003jul;vlbarain2004jul;vlbarain2005jul;vlbarain2006jul;vlbarain2007jul;vlbarain2008jul;vlbarain2009jul;vlbarain2010jul]; vlbarainaug = [vlbarain2003aug;vlbarain2004aug;vlbarain2005aug;vlbarain2006aug;vlbarain2007aug;vlbarain2008aug;vlbarain2009aug;vlbarain2010aug]; vlbarainsep = [vlbarain2003sep;vlbarain2004sep;vlbarain2005sep;vlbarain2006sep;vlbarain2007sep;vlbarain2008sep;vlbarain2009sep;vlbarain2010sep]; vlbarainoct = [vlbarain2003oct;vlbarain2004oct;vlbarain2005oct;vlbarain2006oct;vlbarain2007oct;vlbarain2008oct;vlbarain2009oct;vlbarain2010oct]; vlbarainnov = [vlbarain2003nov;vlbarain2004nov;vlbarain2005nov;vlbarain2006nov;vlbarain2007nov;vlbarain2008nov;vlbarain2009nov;vlbarain2010nov]; vlbaraindec = [vlbarain2003dec;vlbarain2004dec;vlbarain2005dec;vlbarain2006dec;vlbarain2007dec;vlbarain2008dec;vlbarain2009dec;vlbarain2010dec]; vlbarainjanm = nanmean (vlbarainjan); vlbarainfebm = nanmean (vlbarainfeb); vlbarainmarm = nanmean (vlbarainmar); vlbarainaprm = nanmean (vlbarainapr); vlbarainmaym = nanmean (vlbarainmay); vlbarainjunm = nanmean (vlbarainjun); vlbarainjulm = nanmean (vlbarainjul); vlbarainaugm = nanmean (vlbarainaug); vlbarainsepm = nanmean (vlbarainsep); vlbarainoctm = nanmean (vlbarainoct); vlbarainnovm = nanmean (vlbarainnov); vlbaraindecm = nanmean (vlbaraindec); vlbarainannual = [vlbarainjanm;vlbarainfebm;vlbarainmarm;vlbarainaprm;vlbarainmaym;vlbarainjunm;vlbarainjulm;vlbarainaugm;vlbarainsepm;vlbarainoctm;vlbarainnovm;vlbaraindecm]; save ('/Users/sarahstuart/Desktop/precipitation report/vlbarainannual','vlbarainannual'); m = vlbarainannual (:,2); rain = vlbarainannual (:,5); vlbarainannualfig = figure; hold on vlbarainannualfig = bar (m,rain,'blue'); hold on axis ([0 13 0 1.5]); hold on xlabel ('Months','Fontsize',12); hold on ylabel ('Precipitation (mm)','Fontsize',12); hold on grid on hold on set(gca,'XTick',1:1:12) set(gca,'XTickLabel',{'Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep','Oct','Nov','Dec'}) saveas (vlbarainannualfig,'/Users/sarahstuart/Desktop/precipitation report/vlbarainannualfig','png'); load 'vlbarainannual.mat' M = vlbarainannual; wsws = M (:,5); wsmed = mean (wsws); s = size (M); n = s(1); Sa = 0; Sb = 0; WS = []; for i = 1 : n; ws = M (i,5); a = ws * cos ((2*pi*i)/12); b = ws * sin ((2*pi*i)/12); Sa = Sa + a; Sb = Sb + b; end A = 1/6 * Sa; B = 1/6 * Sb; amp = sqrt (A.^2+B.^2) phase = atan (B/A) for i = 1 : n m = M (i,2); ws = M (i,5); wsf = wsmed + amp * cos ((2*pi*i)/12 - phase); WSS = [m,ws,wsf]; WS = [WS;WSS]; end m = WS(:,1); ws = WS (:,2); wsf = WS (:,3); amp a = figure; a = bar (m,ws,'blue'); hold on a = plot (m,wsf,'r','Linewidth',2); hold on grid on xlabel ('Months','Fontsize',12); ylabel ('Precipitation (mm)'); hold on set(gca,'XTick',1:1:12) set(gca,'XTickLabel',{'Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep','Oct','Nov','Dec'}) saveas (a,'/users/sarahstuart/desktop/vlbarainharmana.png');