NCAR

LAMPIRAN

Lampiran 1.

Curah hujan bulanan (mm/hari) sekitar Indonesia (90°BT-140°BT dan 12°LU12°LS) periode Maret 2007-Februari 2008 dengan menggunakan data NCEP/NCAR Reanalysis (curah hujan, angin zonal dan angin meridional).

Lampiran 2. Propagasi SCC (Super Cloud Cluster) dari OLR Periode Maret 2007-Februari 2008 di 90°BT-140°BT

Lampiran 3.

Curah hujan bulanan (mm/hari) sekitar Indonesia (90°BT-140°BT dan 12°LU12°LS) periode Maret 2007-Februari 2008 dengan menggunakan data GPCP

Lampiran 4. Script untuk pengolahan dengan menggunakan sofware Matlab 7.1 %==================================================================== % Program Membaca Data Curah Hujan Bulanan Daerah Kototabang,Pontianak, dan % % Biak % Periode Januari 2007-Februari 2008 % By: Ining Sunarsih % Departemen Geofisika dan Meteorologi % Institut Pertanian Bogor %==================================================================== % Load Data %==================================================================== clear all; y=xlsread('Curah hujan',1,'B2:D13'); x=1:length(y); bar(x,y); xlabel('Bulan','fontweight','bold','fontsize',16); set(gca,'xticklabel',{'Mar' 'Apr' 'Mei' 'Jun' 'Jul' 'Aug' 'Sept' 'Okt' 'Nov' 'Des' 'Jan' 'Feb'}); ylabel('Curah Hujan (mm)','fontweight','bold','fontsize',16); set(gca,'ytick',[50 100 150 200 250 300 350 400 450 500 550]); title('Curah Hujan Bulanan Daerah Kototabang, Pontianak, dan Biak Periode Maret 2007 Februari 2008','fontweight','bold','fontsize',16)

%==================================================================== % Program Membaca Data Curah Hujan Bulanan Daerah Kototabang,Pontianak, dan % Biak % Periode Desember 2007 % By: Ining Sunarsih % Departemen Geofisika dan Meteorologi % Institut Pertanian Bogor %==================================================================== % Load Data %==================================================================== clear all; y=xlsread('Curah hujan',1,'B21:D51'); x=1:length(y); bar(x,y); xlabel('Tanggal','fontweight','bold','fontsize',20); set(gca,'xtick',[3 6 9 12 15 18 21 24 27 30]); ylabel('Curah Hujan (mm)','fontweight','bold','fontsize',20); set(gca,'ytick',[15 30 45 60 75 90 105 120 135]); title('Curah Hujan pada Bulan Desember Daerah Kototabang, Pontianak, dan Biak','fontweight','bold','fontsize',20)

%==================================================================== % Program Membaca Data Curah Hujan Bulanan Daerah Kototabang,Pontianak, dan % Biak % Periode Juni 2007 % By: Ining Sunarsih % Departemen Geofisika dan Meteorologi % Institut Pertanian Bogor %==================================================================== % Load Data %==================================================================== clear all; y=xlsread('Curah hujan',1,'B56:D85'); x=1:length(y);

bar(x,y); xlabel('Tanggal','fontweight','bold','fontsize',20); set(gca,'xtick',[3 6 9 12 15 18 21 24 27 30]); ylabel('Curah Hujan (mm)','fontweight','bold','fontsize',20); set(gca,'ytick',[15 30 45 60 75 90 105 120 135]); title('Curah Hujan pada Bulan Juni Daerah Kototabang, Pontianak, dan Biak','fontweight','bold','fontsize',20)

%==================================================================== % Program Profil Angin Zonal di Kototabang Periode 2 Maret 2007-29 Februari 2008 %==================================================================== clc;clear;close all % LOAD DATA .xls data=xlsread('Koto_rev'); [m,n]=size(data); h=data(1,2:n); data=data(2:m,2:n); [m,n]=size(data); [a,b]=find(data==999); for j=1:length(a) data(a(j),b(j))=0; 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); set(gca,'xlim',[xmin xmax]) xmin=min(data'); h1=pcolor(data');shading interp;colorbar xmax=max(data'); %H1=plot(x,h); %hold on %set(gca,'xlim',[xmin xmax]) title('Angin Zonal di Kototabang Periode 2 Maret 2007-29 Februari 2008','fontweight','bold','fontsize',20) xlabel('Tanggal','fontweight','bold','fontsize',20) set(gca,'xticklabel',{'26 Apr' '15 Jun' '21 Sep' '10 Nov' '30 Des' '18 Feb'}); set(gca,'ytick',[1 7 14 21 28 35 41]); ylabel('Ketinggian (km)','fontweight','bold','fontsize',20); set(gca,'yticklabel',{'2' '3' '4' '5' '6' '7' '8'})

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',20) xlabel('Kecepatan angin (m/det)','fontweight','bold','fontsize',20) ylabel('Ketinggian (km)','fontweight','bold','fontsize',20) set(gca,'ytick',[2 3 4 5 6 7 8])

%==================================================================== % Program Profil Angin Zonal di Pontianak Periode 1 Maret 2007-29 Februari 2008 %==================================================================== clc;clear;close all % LOAD DATA .xls data=xlsread('Ponti_rev'); [m,n]=size(data); h=data(1,2:n); data=data(2:m,2:n); [m,n]=size(data); [a,b]=find(data==999); for j=1:length(a) data(a(j),b(j))=0; 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); set(gca,'xlim',[xmin xmax]) xmin=min(data'); h1=pcolor(data');shading interp;colorbar xmax=max(data'); %H1=plot(x,h); %hold on %set(gca,'xlim',[xmin xmax]) title('Angin Zonal di Pontianak Periode 1 Maret 2007 - 29 Februari 2008','fontweight','bold','fontsize',20) xlabel('Tanggal','fontweight','bold','fontsize',20) set(gca,'xticklabel',{'26 Apr' '15 Jun' '21 Sep' '10 Nov' '30 Des' '18 Feb'}); ylabel('Ketinggian (km)','fontweight','bold','fontsize',20); set(gca,'ytick',[9 20 30 40 51 60 70 80]);

set(gca,'yticklabel',{'1' '2' '3' '4' '5' '6' '7' '8'}) 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',20) xlabel('Kecepatan angin (m/det)','fontweight','bold','fontsize',20) ylabel('Ketinggian (km)','fontweight','bold','fontsize',20) set(gca,'ytick',[1 2 3 4 5 6 7 8])

%==================================================================== % Program Profil Angin Zonal di Biak Periode 11 Maret 2007-13 Februari 2008 %==================================================================== clc;clear;close all % LOAD DATA .xls data=xlsread('Biak_rev'); [m,n]=size(data); h=data(1,2:n); data=data(2:m,2:n); [m,n]=size(data); [a,b]=find(data==999); for j=1:length(a) data(a(j),b(j))=0; 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); set(gca,'xlim',[xmin xmax]) xmin=min(data'); h1=pcolor(data');shading interp;colorbar xmax=max(data'); %H1=plot(x,h); %hold on %set(gca,'xlim',[xmin xmax]) title('Angin Zonal di Biak Periode 11 Maret 2007-13 Februari 2008','fontweight','bold','fontsize',16) xlabel('Tanggal','fontweight','bold','fontsize',16)

set(gca,'xticklabel',{'29 Apr' '25 Jun' '19 Agu' '10 Okt' '1 Des' '24 Jan'}); ylabel('Ketinggian (km)','fontweight','bold','fontsize',16); set(gca,'ytick',[8 19 29 39 50 60 70 79]); set(gca,'yticklabel',{'1' '2' '3' '4' '5' '6' '7' '8'}) 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',16) xlabel('Kecepatan angin (m/det)','fontweight','bold','fontsize',16) ylabel('Ketinggian (km)','fontweight','bold','fontsize',16) set(gca,'ytick',[1 2 3 4 5 6 7 8])

%==================================================================== % Program Profil Angin Zonal di Kototabang pada Tanggal 5-13 Desember 2007 %==================================================================== clc;clear;close all % LOAD DATA .xls data=xlsread('Des_Koto'); % ganti data sesuai lokasi [m,n]=size(data); h=data(1,2:n); data=data(2:m,2:n); [m,n]=size(data); [a,b]=find(data==999); for j=1:length(a) data(a(j),b(j))=0; 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); set(gca,'xlim',[xmin xmax]) xmin=min(data'); h1=pcolor(data');shading interp;colorbar xmax=max(data'); %H1=plot(x,h); %hold on %set(gca,'xlim',[xmin xmax])

title('Angin Zonal di Kototabang pada Tanggal 5-13 Desember 2007','fontweight','bold','fontsize',16); & ganti judul sesuai lokasi xlabel('Tanggal','fontweight','bold','fontsize',16) set(gca,'xtick',[1 2 3 4 5 6 7 8 9]); set(gca,'xticklabel',{'5' '6' '7' '8' '9' '10' '11' '12' '13'}); set(gca,'ytick',[1 7 14 21 28 35 41]); ylabel('Ketinggian (km)','fontweight','bold','fontsize',16); set(gca,'yticklabel',{'2' '3' '4' '5' '6' '7' '8'}) 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',16) xlabel('Kecepatan angin (m/det)','fontweight','bold','fontsize',16) ylabel('Ketinggian (km)','fontweight','bold','fontsize',16) set(gca,'ytick',[2 3 4 5 6 7 8])

%==================================================================== % Program Profil Angin Zonal di Kototabang pada Tanggal 14-17 Juni 2007 %==================================================================== clc;clear;close all % LOAD DATA .xls data=xlsread('Jun_Koto'); % ganti data sesuai lokasi [m,n]=size(data); h=data(1,2:n); data=data(2:m,2:n); [m,n]=size(data); [a,b]=find(data==999); for j=1:length(a) data(a(j),b(j))=0; 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); set(gca,'xlim',[xmin xmax]) xmin=min(data'); h1=pcolor(data');shading interp;colorbar

xmax=max(data'); %H1=plot(x,h); %hold on %set(gca,'xlim',[xmin xmax]) title('Angin Zonal di Kototabang Tanggal 14-17 Juni 2007','fontweight','bold','fontsize',20); %ganti judul sesuai lokasi xlabel('Tanggal','fontweight','bold','fontsize',20) set(gca,'xtick',[1 2 3 4]); set(gca,'xticklabel',{'14' '15' '16' '17'}); set(gca,'ytick',[1 7 14 21 28 35 41]); ylabel('Ketinggian (km)','fontweight','bold','fontsize',20); set(gca,'yticklabel',{'2' '3' '4' '5' '6' '7' '8'}) 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',20) xlabel('Kecepatan angin (m/det)','fontweight','bold','fontsize',20) ylabel('Ketinggian (km)','fontweight','bold','fontsize',20) set(gca,'ytick',[2 3 4 5 6 7 8])

%==================================================================== % Program membaca data kecepatan angin pada reversal wind %==================================================================== data=xlsread('BarKot',1,'B191:B282'); % ganti data dan range data ining=flipud(data); [m,n]=size(ining); x=1:length(ining); y=ining; barh(x,y);grid on ylim([1 92]); % ganti jumlah data xlabel('Kecepatan angin (m/det)','fontweight','bold','fontsize',16); ylabel('Tanggal','fontweight','bold','fontsize',16) title('Kecepatan Angin Zonal Harian pada Ketinggian 5.1 km di Atas Kototabang 1 November 2007 - 31 Januari 2008','fontweight','bold','fontsize',16); % ganti judul sesuai lokasi, ketinggian dan periode data set(gca,'ytick',[23 46 69 92]); % ganti urutan data set(gca,'xtick',[-10 -8 -6 -4 -2 0 2 4 6 8 10 12]); set(gca,'yticklabel',{'9 Jan' '17 Des' '24 Nov' '1 Nov'}); %ganti label sesuai data

%==================================================================== %Program membaca data angin periode 2 Maret 2007 – 29 Februari 2008 di Kototabang dengan menggunakan teknik FFT % ==================================================================== %load data dari excel %==================================================================== data=xlsread('reversal_EAR',1,'B1:C311'); % ganti data dan range data (Data Pontianak dan Biak [m,n]=size(data); t=1:length(data); y=data;

figure;plot(t,y);grid on xlabel('Tanggal','fontweight','bold','fontsize',16); set(gca,'XTickLabel',{'','26 Apr','15 Jun','3 Sep','23 Okt','13 Des','2 Feb'}) ylabel('Kecepatan angin (m/det)','fontweight','bold','fontsize',16); title('Time Series Angin Zonal Harian pada Ketinggian 2.9 dan 5 km di Kototabang 2 Maret 2007 29 Februari 2008','fontweight','bold','fontsize',16); % ganti judul sesuai lokasi [spec,f]= fftrl(y,t); spec=real(spec).^2+imag(spec).^2; %rms frekuensi f=1./f; figure;semilogx(f,spec);grid on set (gca,'xtick',[0 10^0 5 10^1 15 20 30 40 50 60 10^2 500]) xlabel('Periode (hari)','fontweight','bold','fontsize',16) ylabel('Energi Spektral (W/Hz)','fontweight','bold','fontsize',16) title('Power Spectral Density (PSD) Angin Zonal Harian pada Ketinggian 2.9 dan 5 km di Kototabang 2 Maret 2007 - 29 Februari 2008','fontweight','bold','fontsize',16) ; % ganti judul sesuai lokasi

%==================================================================== %Program membaca data angin periode 2 Maret 2007 – 29 Februari 2008 dengan menggunakan teknik wavelet %==================================================================== % 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 Ining Sunarsih on Juni 19, 2008 % ------------------------ loading data -------------------------------- % load 'reversal_koto.txt' ; % input zonal wind series % ---------------------------------------------------------------------%------------------------- computation ----------------------------- % madden_julian = reversal_koto(:,:); variance = std(madden_julian)^2; madden_julian = (madden_julian-mean(madden_julian))/sqrt(variance) ; n = length(madden_julian); dt = 1 ; time = [0:length(madden_julian)-1]*dt + 1.0 ; % construct time array xlim = [1,311]; % plotting range pad = 1; % pad the time series with zeroes (recommended) dj = 0.25; % this will do 4 sub-octaves per octave s0 = 2*dt; % this says start at a scale of 6 months j1 = 7/dj; % this says do 7 powers-of-two with dj sub-octaves each lag1 = 0.72; % lag-1 autocorrelation for red noise background mother = 'Morlet'; % Wavelet transform: [wave,period,scale,coi] = wavelet(madden_julian,dt,pad,dj,s0,j1,mother);

power

= (abs(wave)).^2 ; % compute wavelet power spectrum

% 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 >= 30) & (scale < 60)); 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(:)) %set(gca,'XTickLabel',a) xlabel('Tanggal','fontweight','bold','fontsize',14); set(gca,'xticklabel',{'26 Apr' '15 Jun' '21 Sep' '10 Nov' '30 Des' '18 Feb'}); ylabel('Kecepatan (m/det)','fontweight','bold','fontsize',14) title('a) Time Series Kecepatan Angin Zonal Harian pada Ketinggian 5.1 km di Kototabang Periode 2 Maret 2007 - 29 Februari 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)); %*** or use 'contourfill' [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,'xticklabel',{'26 Apr' '15 Jun' '21 Sep' '10 Nov' '30 Des' '18 Feb'}); 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', ... 'YTick',log2(Yticks(:)), ...o '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,'xticklabel',{'26 Apr' '15 Jun' '21 Sep' '10 Nov' '30 Des' '18 Feb'}); ylabel('Rata-rata varians (m/det)^2','fontweight','bold','fontsize',14) title('d) Rata-rata Time Series 30-60 Harian','fontweight','bold','fontsize',16) hold on plot(xlim,scaleavg_signif+[0,0],'--') hold off

Lampiran 5. Script untuk pengolahan data dengan menggunakan GrADS %==================================================================== %Program membaca data angin dan data curah hujan tahun 2007-2008 %Data dapat didapatkan secara online dari www.cdc.noaa.gov %Oleh : Ining Sunarsih % Departemen Geofisika dan Meteorologi % Institut Pertanian Bogor (IPB) %==================================================================== 'sdfopen d:\prate.sfc.gauss.2007.nc' %ganti data dengan data CH tahun lainnya 'sdfopen d:\uwnd.sig995.2007.nc' %ganti data dengan data CH tahun lainnya 'sdfopen d:\vwnd.sig995.2007.nc' %ganti data dengan data CH tahun lainnya 'set lat -12 12' 'set lon 90 140' 'a=ave(prate,t=60,t=90)' %nilai t ganti untuk bulan Maret 2007-Februari 2008 'b=a*86400' 'c=ave(uwnd.2,t=60,t=90)' %nilai t ganti untuk bulan Maret 2007-Februari 2008 'd=ave(vwnd.3,t=60,t=90)' %nilai t ganti untuk bulan Maret 2007-Februari 2008 'set rgb 16 170 255 255' 'set rgb 17 50 220 255' 'set rgb 18 0 160 255' 'set rgb 19 0 100 255' 'set rgb 20 0 20 255' 'set rgb 21 255 189 189' 'set rgb 22 255 140 140' 'set rgb 23 0 20 200' 'set rgb 24 0 0 170' 'set rgb 25 0 0 110' 'set rgb 26 0 0 0' 'set rgb 27 0 0 0' 'set csmooth on' 'set grid off' 'set grads off' 'set display color white' 'clear' 'set gxout shaded' 'set clevs 10 13 16 19 21 24 27 30 33' 'set ccols 0 16 17 18 19 20 23 24 25 26' 'd b' 'd c;d' 'cbarn' 'set font 0' 'set string 1 1' 'set strsiz 0.14' 'draw string 9.7 0.7 mm' 'set font 5' 'set string 1 1' 'set strsiz 0.14' 'draw string 9.4 1.0 m/det' 'draw title Pola Angin dan Curah Hujan pada Bulan Maret 2007' %ganti judul untuk setiap bulan 'draw xlab Bujur' 'printim d:\0307.gif white' %ganti nama gambar untuk setiap bulan 'reinit'

%==================================================================== %Program membaca data OLR %Data OLR dari www.cdc.noaa.gov % Oleh : Ining Sunarsih % Departemen Geofisika dan Meteorologi % Institut Pertanian Bogor (IPB) %==================================================================== 'sdfopen d:\olr.day.mean.nc' 'set lat 0.0' 'set lon 90 135' 'set t 11962 12327' 'set rgb 16 150 0 220' 'set rgb 17 100 0 200' 'set rgb 18 0 0 255' 'set rgb 19 0 150 220' 'set rgb 20 0 200 220' 'set rgb 21 0 225 150' 'set rgb 23 150 255 0' 'set rgb 24 230 230 0' 'set rgb 25 230 150 0' 'set rgb 26 230 100 0' 'set rgb 27 255 0 0' 'set grads off' 'set display color white' 'clear' 'd olr' 'set yflip on' 'c' 'set yflip on' 'set gxout shaded' 'set clevs 100 120 140 160 180 200 220 240 260 280' 'set ccols 16 17 18 19 20 21 23 24 25 26 27' 'd olr' 'draw title Diagram Hovmoller Periode Maret 2007-Februari 2008' 'draw ylab Tanggal' 'set font 0' 'set string 1 1' 'set strsiz 0.15' 'draw string 10.4 0.5 watt/m2' 'set string 1 1' 'set strsiz 0.16' 'draw string 6 0.75 Bujur' 'cbarn' 'printim d:\hovmoller.gif white' 'reinit'

Lampiran 6. Hasil korelasi silang kecepatan angin zonal pada Ketinggian 5.1 km dengan curah hujan Periode 1 November 2007 – 29 Februari 2008 a.

Kototabang Model Description Model Name Series Name

MOD_1 Angin_zOnal

1 2

Hujan

Transformation Non-Seasonal Differencing

None 0

Seasonal Differencing Length of Seasonal Period Range of Lags

0 No periodicity

From To

-7 7

Display and Plot

All lags

Case Processing Summary Series Length Number of Excluded User-Missing Value Cases Due to System-Missing Value Number of Valid Cases Number of Computable Zero-Order Correlations After Differencing

Cross Correlations Series Pair: Angin_zOnal with Hujan Lag -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7

Cross Correlation .006 -.030 .014 .151 .175 .049 -.005 .116 .234 .255 .240 .240 .143 .059 .068

Std.Error( a) .094 .093 .093 .092 .092 .092 .091 .091 .091 .092 .092 .092 .093 .093 .094

121 0 0 121 121

b. Pontianak Model Description Model Name Series Name

MOD_2 Angin_zOnal

1 2

Hujan

Transformation Non-Seasonal Differencing

None 0

Seasonal Differencing Length of Seasonal Period Range of Lags

0 No periodicity

From To

-7 7

Display and Plot

All lags Case Processing Summary

Series Length Number of Excluded Cases Due to

User-Missing Value System-Missing Value

Number of Valid Cases Number of Computable Zero-Order Correlations After Differencing

Cross Correlations Series Pair: Angin_zOnal with Hujan Lag -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7

Cross Correlation .147 .182 .105 -.066 .051 .015 .035 -.050 .047 .097 .037 -.064 -.039 -.109 .018

Std.Error( a) .094 .093 .093 .092 .092 .092 .091 .091 .091 .092 .092 .092 .093 .093 .094

121 0 0 121 121

c.

Biak Model Description Model Name Series Name

MOD_6 Angin_Zonal

1 2

Hujan

Transformation Non-Seasonal Differencing

None 0

Seasonal Differencing Length of Seasonal Period Range of Lags

0 No periodicity

From To

-7 7

Display and Plot

All lags

Case Processing Summary Series Length Number of Excluded User-Missing Value Cases Due to System-Missing Value Number of Valid Cases Number of Computable Zero-Order Correlations After Differencing

Cross Correlations Series Pair: Tanggal with Angin_Zonal Lag -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7

Cross Correlation -.023 -.037 -.031 -.040 -.124 -.072 .199 .103 -.078 -.073 -.036 -.078 -.058 .076 .084

Std.Error( a) .094 .093 .093 .092 .092 .092 .091 .091 .091 .092 .092 .092 .093 .093 .094

121 0 0 121 121

Lampiran 7. Data Curah Hujan Harian Daerah Kototabang, Pontianak, dan Biak Periode 1 Maret 2007-29 Februari 2008 a. Data Curah Hujan Daerah Kototabang Tanggal Mar Apr May 1 3.72 0.03 2 3.91 0.50 3 33.19 3.17 4 1.92 2.51 5 0 3.39 0 6 0 0.07 0 7 0 18.78 0.04 8 8.05 0.01 23.93 9 0.02 4.12 17.04 10 0.02 18.15 0.76 11 0 3.90 0.02 12 0.28 8.89 0.49 13 21.60 19.49 0.96 14 2.30 13.27 3.41 15 21.47 20.43 1.99 16 10.36 0.17 0 17 0.05 0.77 4.33 18 0.81 0.05 0 19 66.52 17.56 0 20 32.23 1.14 0.03 21 11.60 0 24.15 22 0 0 23 0.04 24 0 25 0.45 26 0 0.29

Jun 5.90 18.64 30.74 14.56 1.16 0.06 0 17.99 0.00 14.02 4.34 0 0.01 41.50 26.11 7.96 17.93 0.26 0.86 0 0 0.72 0.11 0 0 0

Jul 0 0 0 0 0 1.86 0 0.41 1.57 0 9.95 0.93 3.92 0 41.22 0.37 8.97 12.55 0 9.40 8.11 22.50 42.27 1.27 0.01 40.45

Aug 0.25 0.01 0 0 0.73 0 5.13 5.839 0 0.03 1.11 0.02 9.96 0.56 2.41 0.19 0 0 9.75 4.54 1.36 1.79 2.03 9.45 0.01 9.21

Sep 0 7.51 50.24 1.19 0.06 0 49.71 4.917 0.163 0 4.01 0.01 0 9.18 5.81 3.13 0 3.87 0.92 0.11 2.50 1.41 11.57 0.23 0 7.36

Oct 23.07 0.05 0.00 0 0.10 5.73 34.05 18.82 17.48 41.35 3.24 36.11 32.31 2.57 -

Nov 0 5.52 21.09 21.34 0.53 1.231 0 2.10 7.38 0.01 0.68 0.00 0.04 0.04 0.11 8.42

Dec 0 0 0 0.12 31.51 47.32 17.86 12.42 23.61 21.06 0 46.80 49.88 8.44 0.71 3.94 11.4 7.40 0 1.98 18.30 0.26 4.66 14.90 10.05 0.20

Jan 0.02 0 0 0 1.48 0 0 10.10 7.59 1.41 0.02 0.00 0 11.44 0 0 0 9.36 15.33 0.84 4.20 0 0 0 1.65 2.97

Feb 16.74 0.33 0.51 0 0 0 0 0.50 0 2.24 5.76 0 0 0.02 0 0 0 0.12 6.61 1.57 20.50 0.85 0.88 9.65 6.66 25.04

Lanjutan Lampiran 7a 27 0.59 28 0.91 29 3.55 30 1.10 31 2.36

0

74.18 21.38 8.99 0.47 3.63

0 0 0 0

0.58 0.03 0.02 0 0

4.89 0.66 0.95 0.25 0

b. Data Curah Hujan Pontianak Periode 1 Maret 2007 – 29 Februari 2008 Tanggal Mar Apr May Jun Jul Aug 1 6 0 8 3 0 2 30 19 0 0 3 0.5 0 10 32 0 0 4 0 3 7 0 0 0 5 0 0 16 0 3 0 6 0 7 3 5 0 7 0 0 0 2 0 0 8 0 3 0 3 0 9 0 24 103 9 18 0 10 0 0 6 1 11 0 11 0 29 22 0.2 0 12 59 1 3 3 0 31 13 6 92 0 0 0 0 14 0 0.6 6 0 0 15 0 77 47 125 0 15 16 19 0 66 34 8 0 17 14 0 13 11 0 18 6 0 1 49 40 0 19 0 15 4 12 5 0 20 0 2 32 14 1

0 0 0.11 0.51

-

8.56 18.96 0.23 0.03

1.07 3.48 18.12 0 0

0.42 4.28 18.08 8.89 55.64

10.92 0.56 1.98

Sep 18 26 73 24 26 0 0 0 2 0 0 0 11 0 0 2 16 0 0

Oct 0 0 0 0 0 0 4 105 7 34 0 28 66 0 28 67 4 7 24 15

Nov 5 16 20 1 2 44 11 19 13 3 16 0.4 0 0 0 1 14 16 1 0

Dec 0 33 0 38 8 10 46 15 40 20 0.4 31 0.8 4 1 7 69 3 15

Jan 0 20 0 1 1 0 3 0 0 0 0 0 0 0 0 0 0 0 0

Feb 10 27 12 0 15 0 0 2 4 0 0 0 0 0 0 0 0

Lanjutan Lampiran 7b 21 0.8 22 0 23 36 24 10 25 0 26 0 27 0 28 29 30 2 31 21

1 5 3 1 8 11 0 0 0 8

0.7 39 38 0 37 7 2 0 7 6 32

48 0 0 0 0 0 5 0 2

31 21 0.5 4 3 45 0 2 0 16 81

c. Data Curah Hujan Biak Periode 1 Maret 2007 – 29 Februari 2008 Tanggal Mar Apr May Jun Jul 1 7 0.2 32 41 2 0 6 4 12 3 26 4 0 0 43 4 26 11 0 2 5 0 0 1 6 16 16 2 0 7 22 0 7 0 8 12 0 0 0 0 9 0 9 6 10 0.2 0 1 0 11 73 0 11 0 0 12 0 15 2 10 17 13 0 1 0 38 14 2 0 0 24 30

31 19 0 2 7 3 4 1 0 3

0 0 0 0 0 0 37 1 0 0

1 81 6 6 13 28 12 -

20 33 0 0 0 1 9 1 2 4

2 4 0 0 1 12 5 0 0 0 2

0 0 0 0 4 45 3 2 0.5 1 34

3 0 0 0 0 0 0 6

Aug 0 8 4 22 4 0.4 4 26 2 11 0.5 -

Sep 3 3 3 27 9 0 6 6 0 3 69

Oct 0 0 0 0 33 0 0 1 32 9 2 0 6 -

Nov 12 6 0 23 20 0 0.4 20 0 0.4 0.2 3

Dec 10 0 0 48 13 0 15 2 0 0 1

Jan 43 44 19 0 33 19 16 29 10 0.3 5 0 0

Feb 38 3 38 0 20 43 1 2 23 0 22 0.2

Lanjutan Lampiran 7c 15 16 0.5 17 0.5 18 0.8 19 7 20 0.4 21 60 22 0 23 0 24 11 25 5 26 0.5 27 28 29 3 30 31 0 Keterangan:

30 1 0 0 9 15 0 0 5 0 0

-

4 0 0 19 25 4 2 6 25 60 10 4 -

: Tidak ada hujan

24 14 0.2 0.3 1 9 15 0 6 4 28 11

1 9 14 41 54 7 15 35 28 6 10 0.2 5 51

5 1 18 0 0 3 46 0.9 0 3 13 2 18 18 12

29 22 10 7 4 2 9 1 0 0 13 0 0 0

0.2 22 0 1 0.2 5 0.8 5 3 9 0 -

1 25 9 0.1 0 0 3 18 0 0 0 0 0.7 -

21 12 0 3 17 12 12 26 0.3 1 16 20 0.3 -

0 20 3 8 1 2 0 65 0 1 58 0 1 6

31 8 0 35 5 26 3 1 0.2 3 0