LAMPIRAN
21
Lampiran 1 Time Height Section Angin Zonal di Kototabang (H = 2 – 20 km) Periode 1 Desember 2007 – 31 Desember 2008.
Lampiran 2 Time Height Section Angin Meridional di Kototabang (H = 2 – 20 km) Periode 1 Desember 2007 – 31 Desember 2008.
22
Lampiran 3 Time Height Section Angin Zonal di Kototabang (H = 2 – 20 km) Periode Desember 2007 – Februari 2008.
Lampiran 4 Time Height Section Angin Meridional di Kototabang (H = 2–20 km) Periode Desember 2007 – Februari 2008.
23
Lampiran 5 Time Height Section Angin Meridional di Kototabang (H = 14–20 km) Periode Desember 2007 – Februari 2008.
Lampiran 6 Time Height Section Angin Zonal di Kototabang (H = 2 – 20 km) Periode Juni – Agustus 2008.
24
Lampiran 7 Time Height Section Angin Meridional di Kototabang (H = 2–20 km) Periode Juni – Agustus 2008.
Lampiran 8 Time Height Section Angin Meridional di Kototabang (H = 14–20 km) Periode Juni – Agustus 2008.
25
Lampiran 9 Struktur vertikal gelombang Kelvin di Atas Kototabang Periode 1 Desember 2007 – 31 Desember 2008.
Lampiran 10 Struktur vertikal gelombang Kelvin .di Atas Kototabang Periode Desember 2007 – Februari 2008
26
Lampiran 11 Struktur vertikal gelombang Kelvin .di Atas Kototabang Periode Juni – Agustus 2008
Lampiran 12 Time Height Section Vertikal Echo Intensity di Kototabang (H = 2– 20 km) Periode 1 Desember 2007 – 31 Desember 2008.
27
Lampiran 13 Time Height Section Vertikal Echo Intensity di Kototabang (H = 2 – 20 km) Periode Desember 2007 – Februari 2008.
Lampiran 14 Time Height Section Vertikal Echo Intensity di Kototabang (H = 2 – 20 km) Periode 1 Juni - Agustus 2008.
28
Lampiran 15 Power Spektral Density angin zonal pada ketinggian 2 - 20 km di Kototabang periode 1 Desember 2007 - 31 Desember 2008
Lampiran 16 Power Spektral Density angin meridional pada ketinggian 2 - 20 km di Kototabang periode 1 Desember 2007 - 31 Desember 2008
29
Lampiran 17 Power Spektral Density angin zonal pada ketinggian 2 - 20 km di Kototabang periode 1 Desember 2007 - 31 Desember 2008
Lampiran 18 Power Spektral Density angin meridional pada ketinggian 2 - 20 km di Kototabang periode 1 Desember 2007 - 31 Desember 2008
30
Lampiran 19 Power Spektral Density angin zonal pada ketinggian 2 - 20 km di Kototabang periode 1 Desember 2007 - 31 Desember 2008
Lampiran 20 Power Spektral Density angin meridional pada ketinggian 2 - 20 km di Kototabang periode 1 Desember 2007 - 31 Desember 2008
Lampiran 21 Daftar Nilai Power Spektral Density Angin Zonal Di Atas Kototabang (H = 15.05 - 18.00 km) pada beberapa periode puncak selama 1 Desember 2007 – 1 Desember 2008 H (km)
Periode (hari)
15.05
7 7,884,901.36
11 43,306,129.48
16 45,162,149.89
18 21,025,513.89
25 21,609,989.29
31 1,919,177.87
66 101,227,218.32
199 146,550,052.06
15.19
13,031,424.59
50,087,586.50
53,222,616.19
19,188,509.26
47,027,994.02
6,513,070.86
129,531,948.41
151,227,445.17
15.34
14,469,398.46
48,197,692.77
48,608,473.27
12,353,872.76
73,560,580.15
20,875,505.30
168,597,052.20
157,528,617.83
15.49
34,305,069.27
56,692,061.46
44,642,176.48
14,872,248.73
88,161,961.86
30,674,538.60
168,576,518.20
142,732,546.94
15.64
56,939,566.81
61,276,539.36
59,283,945.94
8,590,159.96
89,478,581.35
37,349,638.32
189,145,252.53
131,897,939.52
15.79 15.93
54,229,422.16 56,430,556.23
49,641,882.14 63,514,055.67
61,471,027.06 23,798,116.85
4,289,356.88 7,609,921.35
76,364,622.77 47,731,129.02
57,073,414.18 73,398,974.55
171,911,284.28 178,079,955.66
129,065,543.34 103,094,544.13
16.08
28,584,980.00
28,079,925.32
12,475,873.04
3,638,224.77
29,215,806.23
84,449,068.36
237,903,746.63
123,128,876.16
16.23
5,664,996.93
17,801,130.33
10,431,179.82
7,653,341.60
23,426,536.58
104,779,676.57
251,130,778.69
140,479,108.33
16.38
1,486,843.43
16,057,709.78
18,607,140.08
20,282,030.33
22,438,101.89
97,890,949.36
243,936,636.74
161,073,085.08
16.52
1,930,861.90
16,589,066.99
42,687,502.06
42,254,455.58
22,841,170.18
117,091,422.82
244,473,107.12
173,650,420.63
16.67
4,949,321.53
30,266,909.98
82,876,684.79
44,305,238.03
24,965,869.17
123,466,363.37
245,765,482.30
134,637,756.88
16.82 16.97
1,911,001.58 2,981,197.87
27,867,966.66 10,068,092.56
62,576,067.87 52,075,983.36
49,546,536.73 58,149,863.50
36,928,810.41 36,958,792.49
71,256,957.15 20,939,463.67
254,532,066.83 252,792,098.13
87,538,815.87 74,350,478.21
17.11
11,458,095.62
4,002,762.82
47,657,564.60
69,465,953.56
30,869,734.38
3,988,700.86
224,712,873.05
72,025,116.55
17.26 17.41
10,614,522.59 4,654,539.35
171,059.62 3,162,601.95
34,744,351.70 33,205,363.84
94,779,334.61 106,706,456.39
21,195,606.36 13,610,036.46
11,563,941.81 28,101,967.80
197,467,008.02 134,737,329.75
75,977,595.08 83,870,331.59
17.56
2,277,657.90
5,583,949.80
25,202,852.01
87,684,869.92
11,570,246.84
23,149,996.73
59,568,304.76
87,883,143.55
17.71
3,523,905.60
2,659,871.75
14,190,187.15
58,021,066.10
11,265,481.92
9,845,730.44
26,384,794.97
87,641,923.76
17.85
4,955,020.43
7,671,016.66
17,702,644.14
50,974,916.67
10,835,467.67
3,285,951.16
14,369,115.90
65,716,989.89
18.00
2,573,751.65
8,345,697.73
22,969,799.15
44,687,630.39
11,662,318.63
1,539,852.43
5,055,793.42
35,542,440.50
31
Lampiran 22 Daftar Nilai Power Spektral Density Angin Zonal Di Atas Kototabang (H = 15.05 - 18.00 km) pada beberapa periode puncak selama Desember 2007 – Februari 2008. Periode (Hari)
H (km)
9
11
15
46
15.05 15.19 15.34 15.49 15.64 15.79 15.93 16.08 16.23 16.38 16.52 16.67 16.82 16.97 17.11 17.26 17.41 17.56 17.71 17.85 18.00
1,113,830.12 1,466,164.70 2,386,214.52 2,423,947.83 5,166,843.39 7,400,192.11 7,189,551.99 11,742,870.58 27,655,091.37 24,141,424.32 17,940,141.19 13,415,350.23 3,770,537.48 682,908.51 866,515.49 241,821.69 608,280.70 1,489,931.22 3,733,039.02 4,534,750.24 1,992,377.36
9,177,071.74 5,649,362.35 4,255,208.18 2,597,919.97 3,889,757.56 1,924,625.04 1,686,489.04 2,059,316.54 11,193,502.46 18,965,257.46 27,581,586.59 37,269,398.82 40,915,478.92 25,642,167.49 11,559,179.85 299,758.28 3,786,567.58 8,137,368.67 6,363,862.84 8,242,604.08 5,116,479.04
12,972,372.84 12,002,398.18 9,448,017.76 5,745,399.02 5,517,533.47 3,395,572.83 1,931,609.71 2,267,525.64 2,588,797.51 1,520,081.40 4,004,357.99 19,002,772.15 25,150,716.34 36,461,085.01 35,712,442.52 27,256,499.28 34,446,988.20 34,313,520.61 21,821,308.74 15,515,801.96 15,922,083.46
37,924,957.33 53,823,711.23 51,522,759.13 63,707,784.90 70,679,619.37 57,373,097.27 38,857,843.04 49,644,668.14 53,040,720.23 44,302,884.76 22,329,553.97 12,049,192.55 5,568,259.08 9,188,802.85 13,357,119.83 25,283,255.40 29,949,447.88 23,669,688.99 18,084,828.42 14,536,673.83 9,579,458.47
Lampiran 23 Daftar Nilai Power Spektral Density Angin Zonal Di Atas Kotabang (H = 15.05 - 18.00 km) pada beberapa periode puncak selama Juni – Agustus 2008 H (km)
Periode (Hari) 9
18
15.05
717,711.11
2,211,697.15
15.19
1,975,035.59
1,687,106.40
15.34
2,165,804.61
1,341,553.71
15.49
2,423,337.34
651,322.84
15.64
1,580,638.67
16,564.04
15.79
1,910,753.22
373,641.12
15.93
1,525,769.58
1,936,997.07
16.08
1,234,246.20
4,630,867.62
16.23
2,554,082.53
7,789,722.50
16.38
7,859,213.88
13,995,218.24
16.52
11,584,532.75
19,786,282.76
16.67
12,990,470.95
17,454,639.12
16.82
11,341,467.42
15,589,512.15
16.97
7,455,169.08
15,125,432.65
17.11
4,991,090.73
12,233,953.80
17.26
2,172,658.97
12,533,607.62
17.41
66,210.05
13,082,202.98
17.56
570,063.74
11,336,978.76
17.71
1,077,611.54
9,007,474.95
17.85
842,942.58
6,413,961.72
18
692,793.96
4,030,736.52
32
33
Lampiran 24
Script untuk pengolahan data dengan menggunakan software Matlab 7.2.0.232 (R2006a)
%========================================================= %Program membaca data EAR periode Desember 2007- Desember 2008 %untuk melihat profil kontur %Modified by : Widya Ningrum % Departemen Geofisika dan Meteorologi % Institut Pertanian Bogor %========================================================= %LOAD DATA.xls %========================================================= clc;clear;close all % LOAD DATA .xls data=xlsread('LONG_TERM',2); [m,n]=size(data); h=data(1,4:n); data=data(2:m,4:n); [m,n]=size(data); [a,b]=find(data==999); for j=1:length(a) data(a(j),b(j))=nan; end % RATA-RATA================================================== for i=1:n temp=data(:,i); index=find(isnan(temp)==1); temp(index)=[]; avg(i)=mean(temp); warning off end % PLOT ===================================================== x=data(1,:); xmin=min(x); subplot(2,3,[1 2 4 5]) H1=plot(x,h); hold off for i=2:m x=(i-1)*1+data(i,:); H1=plot(x,h); set(H1,'color',[rand(1) rand(1) rand(1)]) end xmax=max(x); xmin=min(data'); h1=pcolor(data');shading interp; contourcmap([2:0.5:1],'jet','colorbar','on','location','vertical') xmax=max(data'); title('Angin Meridional di Kototabang Periode 1 Desember 2007 - 31 Desember2008','fontweight','bold','fontsize',20)
34
set(gca,'xtick',[1000 2000 3000 4000 5000 6000 7000 8000 9000]); set(gca,'xticklabel',{'11 Jan' '22 Feb' '3 Apr' '15 Mei' '26 Jun' '6 Ags' '17 Sep' '29 Okt' '9 Des'}) xlabel('Tanggal','fontweight','bold','fontsize',18) set(gca,'ytick',[1 14 28 41 55 68 82 95 109 120]); set(gca,'yticklabel',{'2' '4' '6' '8' '10' '12' '14' '16' '18' '20'}) ylabel('Ketinggian (km)','fontweight','bold','fontsize',18); subplot(2,3,[3 6]) plot(avg,h) xmin=min(avg); xmax=max(avg); set(gca,'xlim',[xmin xmax]) title('Rata-rata','fontweight','bold','fontsize',19) xlabel('Kecepatan (m/det)','fontweight','bold','fontsize',18) ylabel('Ketinggian (km)','fontweight','bold','fontsize',18)
angin
%========================================================= %Program membaca data EAR periode Desember 2007- Desember 2008 %untuk melihat profil vertikal %Modified by : Widya Ningrum % Departemen Geofisika dan Meteorologi % Institut Pertanian Bogor %========================================================= %LOAD DATA.xls %========================================================= clc;clear;close all % LOAD DATA .xls data=xlsread('LONG',2); [m,n]=size(data); h=data(1,4:n); data=data(2:m,4:n); [m,n]=size(data); [a,b]=find(data==999); for j=1:length(a) data(a(j),b(j))=nan; end % PLOT .................... x=data(1,:); xmin=min(x); subplot(2,3,[1 2 4 5]); H1=plot(x,h,'-k'); hold on for i=2:m x=(i-1)*1+data(i,:); H1=plot(x,h); end xmax=max(x);
35
set(gca,'xlim',[xmin xmax]) title('Profil Vertikal Angin Zonal di Kototabang Periode 1 Januari 2008 - 31 Desember 2008','fontweight','bold','fontsize',20) set(gca,'xtick',[1000 2000 3000 4000 5000 6000 7000 8000]); set(gca,'xticklabel',{'11 Feb' '24 Mar' '4 Mei' '15 Jun' '27 Jul' '6 Sep' '18 Okt' '29 Nov'}) xlabel('Tanggal','fontweight','bold','fontsize',18) ylabel('Ketinggian (km)','fontweight','bold','fontsize',18); subplot(2,3,[3 6]) for i=1:n temp=data(:,i); index=find(isnan(temp)==1); temp(index)=[]; avg(i)=mean(temp); warning off end plot(avg,h) xmin=min(avg); xmax=max(avg); set(gca,'xlim',[xmin xmax]) title('Rata-rata','fontweight','bold','fontsize',18) xlabel('Kecepatan (m/det)','fontweight','bold','fontsize',18) ylabel('Ketinggian (km)','fontweight','bold','fontsize',18) set(gca,'ytick',[2 4 6 8 10 12 14 16 18 20])
angin
%========================================================= %Program membaca data EAR periode Desember 2007-Desember 2008 %untuk melihat anomali angin %Modified by : Widya Ningrum % Departemen Geofisika dan Meteorologi % Institut Pertanian Bogor %========================================================= %LOAD DATA.xls %========================================================= clear all; i=1; data=csvread('LONGII.csv'); [m,n]=size(data); [a,b]=find(data==0); for j=1:length (a); data (a(j),b(j))=nan; end tinggi=data(1:m,1); [mh nh]=size(tinggi); z1=csvread('LONGII.csv'); [m n]=size(z1); for j=1:m for i=1:n if(z1(j,i)==999.0) z1(j,i)=NaN;
36
end end
end j=1; sc=2; for i=1:m plot(sc*i+z1(i,:),'black');hold on; mx(i,:)=max(sc*i+z1(i,:)); mn(i,:)=min(sc*i+z1(i,:)); set(gca, 'Tickdir', 'in'); set(gca, 'Ticklength',[0.02,0.02]); end m0=(mn+mx)/2; m1=min(m0); m2=max(m0); ml1=min(mn); ml2=max(mx); m3=(m2-m1)/5; ms=m1:m3:m2; h1=tinggi(2,:); h2=tinggi(mh-1,:); h3=(h2-h1)/5; hs=h1:h3:h2; ylim([ml1 ml2]); set(gca,'YTick',ms); d=1:406; xlim([min(d) max(d)]); set(gca,'yTickLabel',{hs}); set(gca, 'Fontweight', 'bold'); set(gca, 'Fontsize',14); xlabel('Day Number','fontweight','bold','fontsize',20); ylabel('Height (km)','fontweight','bold','fontsize',20); title('Zonal Wind Anomaly di Kototabang 100,32 BT, 0,23 LS Periode 1 Desember 2007 - 31 Desember 2008')
%========================================================= %Program membaca data angin zonal periode Desember 2007- Desember 2008 %untuk melihat power spectral density dengan metode FFT %Modified by : Widya Ningrum % Departemen Geofisika dan Meteorologi % Institut Pertanian Bogor %========================================================= %LOAD DATA.xls %========================================================= %load data dari excel data=xlsread('LONG_TERM',3,'A2:U9529'); [m,n]=size(data); t=1:length(data); y=data; [spec,f]= fftrl(y,t); spec=real(spec).^2+imag(spec).^2; %rms frekuensi f=1./f;
37
figure;semilogx(f,spec);grid on set(h,'Interpreter','none') set (gca,'xtick',[0 10^0 275 425 750 1440]) xlabel('Periode (jam)','fontweight','bold','fontsize',20) ylabel('Energi Spektral (J/Hz)','fontweight','bold','fontsize',20) title('Power Spectral Density (PSD) Angin Zonal Harian di Kototabang Periode 1 Januari 2008 31 Desember 2008','fontweight','bold','fontsize',22)
%========================================================= %Program membaca data angin zonal periode Desember 2007- Desember 2008 %untuk melihat power spectral density dengan metode wavelet %Modified by : Widya Ningrum % Departemen Geofisika dan Meteorologi % Institut Pertanian Bogor %========================================================= %LOAD DATA.xls %========================================================= % % % % % % % % % %
WAVETEST Example Matlab script for WAVELET, using NINO3 SST dataset See "http://paos.colorado.edu/research/wavelets/" Written January 1998 by C. Torrence Modified Oct 1999, changed Global Wavelet Spectrum (GWS) to be sideways, changed all "log" to "log2", changed logarithmic axis on GWS to a normal axis. normalize by standard deviation (not necessary, but makes it easier to compare with plot on Interactive Wavelet page, at "http://paos.colorado.edu/research/wavelets/plot/" Modified by Widya Ningrum on March 2009
% ------------------------ loading data -------------------------% load 'enambelaskmsembln.txt' ; % input zonal wind series % ---------------------------------------------------------------%------------------------- computation --------------------------% madden_julian = enambelaskmsembln(:,:); variance = std(madden_julian)^2; madden_julian =(madden_julian-mean(madden_julian))/sqrt(variance); n dt time xlim pad dj s0 j1 lag1 mother
= = = = = = = = = =
length(madden_julian); 1 ; [0:length(madden_julian)-1]*dt + 1.0 ; % construct time array [1,397]; % plotting range 1; % pad the time series with zeroes (recommended) 0.25; % this will do 4 sub-octaves per octave 2*dt; % this says start at a scale of 6 months 7/dj; % this says do 7 powers-of-two with dj sub-octaves each % lag-1 autocorrelation for red noise background 0.72; 'Morlet';
% Wavelet transform: [wave,period,scale,coi]= wavelet(madden_julian,dt,pad,dj,s0,j1,mother; power = (abs(wave)).^2 ; % compute wavelet power spectrum
38
% Significance levels: (variance=1 for the normalized SST) [signif,fft_theor] = wave_signif(1.0,dt,scale,0,lag1,-1,-1,mother); sig95 = (signif')*(ones(1,n)); % expand signif --> (J+1)x(N) array
sig95 = power ./ sig95 ;
%
where ratio > 1, power is significant
% Global wavelet spectrum & significance levels: global_ws = variance*(sum(power')/n); % time-average over all times dof = n - scale; % the -scale corrects for padding at edges global_signif = wave_signif(variance,dt,scale,1,lag1,-1,dof,mother);
% Scale-average between Madden_Julian periods of 30--60days avg = find((scale >= 10) & (scale < 20)); Cdelta = 0.776; % this is for the MORLET wavelet scale_avg = (scale')*(ones(1,n));% expand scale --> (J+1)x(N) array scale_avg = power ./ scale_avg; % [Eqn(24)] scale_avg
= variance*dj*dt/Cdelta*sum(scale_avg(avg,:));
% [Eqn(24)]
scaleavg_signif = wave_signif(variance,dt,scale,2,lag1,-1,[2,7.9],mother);
whos %-------------------------------------Plotting %--- Plot time series subplot('position',[0.08 0.75 0.56 0.18]) plot(time,madden_julian) grid set(gca,'XLim',xlim(:)) xlabel('Tanggal','fontweight','bold','fontsize',14); set(gca,'xtick',[50 100 150 200 250 300 350]); set(gca,'xticklabel',{'19 Jan' '9 Mar' '28 Apr' '17 Jun' '6 Aug' '25 Sep' '14 Nov'}) ylabel('Kecepatan (m/det)','fontweight','bold','fontsize',14) title('a) Time Series Kecepatan Angin Zonal Harian pada Ketinggian 17.67 km di Atas Kototabang Periode 1 Desember 2007 - 31 Desember 2008','fontweight','bold','fontsize',16) hold off %--- Contour plot wavelet power spectrum subplot('position',[0.08 0.38 0.66 0.25]) levels = [0.0625,0.125,0.25,0.5,1,2,4,8,16] ; Yticks = 2.^(fix(log2(min(period))):fix(log2(max(period)))); % contour(time,log2(period),log2(power),log2(levels)); use 'contourfill'
%*** or
[C,h]=contourf(time,log2(period),log2(power),log2(levels)); %*** or use 'contourfill' colormap colorbar %imagesc(time,log2(period),log2(power));
%*** uncomment for 'image' plot
xlabel('Tanggal','fontweight','bold','fontsize',14) set(gca,'xtick',[50 100 150 200 250 300 350]); set(gca,'xticklabel',{'19 Jan' '9 Mar' '28 Apr' '17 Jun' '6 Aug' '25 Sep' '14 Nov'}) ylabel('Periode (hari)','fontweight','bold','fontsize',14) title('b) Wavelet Power Spektrum','fontweight','bold','fontsize',16)
set(gca,'XLim',xlim(:)) set(gca,'YLim',log2([min(period),max(period)]), ... 'YDir','reverse', ...
39
'YTick',log2(Yticks(:)), ... 'YTickLabel',Yticks)
% 95% significance contour, levels at -99 (fake) and 1 (95% signif)
hold on contour(time,log2(period),sig95,[-99,1],'k'); hold on % cone-of-influence, anything "below" is dubious plot(time,log2(coi),'k') hold off
%--- Plot global wavelet spectrum subplot('position',[0.78 0.37 0.2 0.25]) plot(global_ws,log2(period)) grid hold on plot(global_signif,log2(period),'--') grid hold off xlabel('Power (m/det)^2','fontweight','bold','fontsize',14) Yticks = 2.^(fix(log2(min(period))):fix(log2(max(period)))); title('c) Global Wavelet Spektrum','fontweight','bold','fontsize',16)
set(gca,'YLim',log2([min(period),max(period)]), ... 'YDir','reverse', ... 'YTick',log2(Yticks(:)), ... 'YTickLabel',Yticks) set(gca,'XLim',[0,1.25*max(global_ws)])
%--- Plot 30--60 days scale-average time series subplot('position',[0.08 0.07 0.56 0.18]) plot(time,scale_avg) grid set(gca,'XLim',xlim(:)) xlabel('Tanggal','fontweight','bold','fontsize',14) set(gca,'xtick',[50 100 150 200 250 300 350]); set(gca,'xticklabel',{'19 Jan' '9 Mar' '28 Apr' '17 Jun' '6 Aug' '25 Sep' '14 Nov'}) ylabel('Rata-rata varians (m/det)^2','fontweight','bold','fontsize',14) title('d) Rata-rata Time Series 10 - 20 Harian','fontweight','bold','fontsize',16)
hold on plot(xlim,scaleavg_signif+[0,0],'--') hold off
40
Lampiran 25 Hasil spektral silang (cross spectrum) antara angin zonal (H=16.97 km) dengan Kototabang Tropopause Layer periode Desember 2007 – Februari 2008 Period 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45
90.00 45.00 30.00 22.50 18.00 15.00 12.86 11.25 10.00 9.00 8.18 7.50 6.92 6.43 6.00 5.63 5.29 5.00 4.74 4.50 4.29 4.09 3.91 3.75 3.60 3.46 3.33 3.21 3.10 3.00 2.90 2.81 2.73 2.65 2.57 2.50 2.43 2.37 2.31 2.25 2.20 2.14 2.09 2.05 2.00
Cross Quad
Cross Amplit.
Squared Coherncy
0.000 26.299 33.730 35.467 18.587 20.397 35.558 6.901 -12.868 -3.860 -1.206 -0.931 -3.545 -0.642 8.360 7.235 -3.869 -0.164 10.772 6.745 6.143 9.010 -2.037 -3.149 7.363 0.744 -5.990 -3.344 -5.214 1.753 5.577 -8.478 -11.943 -2.064 -0.854 -4.559 -6.241 -6.715 -0.424 5.207 3.481 4.348 4.972 10.679 15.325 0.000
106.936 137.064 127.368 75.463 50.240 35.988 40.089 9.877 13.247 8.921 2.574 8.810 20.330 7.631 11.631 11.592 15.843 1.029 11.592 8.039 9.399 9.493 9.399 25.801 15.387 2.502 9.863 7.666 7.550 4.209 7.938 12.058 11.943 4.600 9.278 6.462 6.464 7.645 5.381 5.247 3.494 7.448 9.692 13.918 17.436 6.581
0.962 0.966 0.805 0.566 0.382 0.138 0.200 0.030 0.110 0.067 0.008 0.104 0.373 0.035 0.059 0.073 0.232 0.002 0.227 0.140 0.208 0.108 0.056 0.282 0.156 0.010 0.182 0.139 0.147 0.034 0.070 0.139 0.177 0.047 0.277 0.192 0.230 0.280 0.119 0.100 0.047 0.247 0.524 0.631 0.713 0.126
Gain val X over Y 0.134 0.134 0.124 0.116 0.094 0.031 0.021 0.007 0.017 0.016 0.004 0.010 0.021 0.009 0.019 0.028 0.048 0.004 0.047 0.022 0.016 0.015 0.017 0.045 0.028 0.005 0.019 0.018 0.019 0.009 0.019 0.033 0.032 0.012 0.021 0.016 0.017 0.017 0.011 0.012 0.008 0.020 0.026 0.030 0.038 0.017
Gain val Y over X 7.156 7.224 6.490 4.895 4.081 4.466 9.726 4.181 6.521 4.247 1.936 10.932 18.122 3.785 3.072 2.623 4.807 0.387 4.833 6.447 12.872 7.259 3.274 6.256 5.502 1.865 9.724 7.876 7.949 3.750 3.750 4.223 5.485 4.059 12.979 12.323 13.392 16.060 10.850 8.431 5.781 12.647 19.959 21.384 18.748 7.544
Phase Spectrum
Hamming Weights
0.000 2.949 2.874 2.652 2.763 2.539 2.051 2.368 -1.811 -2.694 -2.654 -0.106 -0.175 -0.084 2.340 0.674 -0.247 -0.160 1.949 2.146 2.429 1.891 -2.923 -3.019 2.643 2.840 -2.489 -2.690 -2.379 2.712 0.779 -0.780 -1.561 -0.465 -0.092 -0.783 -1.307 -2.069 -3.063 1.694 1.655 2.518 2.603 2.267 2.068 0.000
0.036 0.241 0.446 0.241 0.036
41
Lampiran 26 Hasil spektral silang antara angin zonal (H=16.52 km) dengan Kototabang Tropopause Layer periode Juni - Agustus 2008 Period 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46
92.00 46.00 30.67 23.00 18.40 15.33 13.14 11.50 10.22 9.20 8.36 7.67 7.08 6.57 6.13 5.75 5.41 5.11 4.84 4.60 4.38 4.18 4.00 3.83 3.68 3.54 3.41 3.29 3.17 3.07 2.97 2.88 2.79 2.71 2.63 2.56 2.49 2.42 2.36 2.30 2.24 2.19 2.14 2.09 2.04 2.00
Cross Quad
Cross Amplit.
Squared Coherncy
0.000 -1.987 -0.172 -8.656 -17.939 -16.415 -19.538 -20.407 -10.215 -1.163 11.496 13.238 5.027 0.929 -2.324 -6.968 -5.240 -4.533 -3.479 -1.851 -2.209 -1.682 -0.849 -0.024 0.632 0.380 0.859 2.367 3.201 4.152 0.938 -4.249 -4.094 -2.297 -0.916 -0.387 -0.392 0.055 0.611 -0.216 -1.435 -0.119 1.453 1.757 0.864 0.153 0.000
20.172 28.932 40.398 51.514 77.039 92.129 51.143 20.411 12.685 6.609 12.789 13.595 5.131 3.725 6.692 7.987 5.366 4.721 4.057 1.857 3.334 1.854 1.838 0.299 1.256 0.492 2.275 4.647 4.720 5.406 3.042 4.596 4.271 2.327 1.262 0.773 0.505 0.085 1.164 3.233 4.126 0.537 1.719 1.780 1.069 0.262 0.037
0.932 0.857 0.676 0.653 0.801 0.871 0.506 0.212 0.208 0.072 0.176 0.244 0.074 0.049 0.080 0.100 0.085 0.163 0.360 0.090 0.208 0.096 0.209 0.012 0.250 0.024 0.161 0.424 0.395 0.308 0.084 0.307 0.715 0.415 0.138 0.067 0.025 0.002 0.236 0.566 0.422 0.011 0.249 0.315 0.304 0.249 0.183
Gain val X over Y 0.126 0.122 0.133 0.157 0.145 0.138 0.105 0.069 0.062 0.037 0.077 0.130 0.057 0.035 0.050 0.047 0.047 0.056 0.051 0.033 0.093 0.112 0.222 0.043 0.126 0.035 0.193 0.463 0.364 0.219 0.074 0.126 0.228 0.235 0.193 0.168 0.049 0.005 0.080 0.163 0.151 0.030 0.162 0.124 0.070 0.032 0.008
Gain val Y over X 7.392 7.024 5.083 4.164 5.516 6.318 4.808 3.080 3.377 1.963 2.282 1.876 1.294 1.404 1.588 2.145 1.823 2.899 7.021 2.754 2.242 0.859 0.944 0.269 1.987 0.693 0.833 0.916 1.086 1.403 1.124 2.435 3.140 1.770 0.714 0.401 0.521 0.306 2.957 3.479 2.795 0.376 1.535 2.540 4.363 7.800 21.653
Phase Spectrum
Hamming Weights
0.000 -3.073 -3.137 -2.973 -2.907 -2.962 -2.750 -1.550 -0.936 -0.177 1.117 1.341 1.370 0.252 -0.355 -1.060 -1.354 -1.854 -2.111 -1.654 -0.724 -1.137 -2.661 -3.061 0.527 0.881 2.755 2.607 2.396 2.266 2.828 -1.962 -1.860 -1.731 -0.812 -0.525 -0.889 0.707 2.589 -3.075 -2.786 -2.919 1.007 1.733 2.200 2.521 0.000
0.036 0.241 0.446 0.241 0.036
