TUGAS GEODESI FISIS

( Disusun untuk memenuhi tugas mata kuliah Geodesi Fisis)

Disusun Oleh :

Viona Yashinta 21110115120030

DEPARTEMEN TEKNIK GEODESI

FAKULTAS TEKNIK – UNIVERSITAS DIPONEGORO

Jl. Prof. Sudarto SH, Tembalang Semarang Telp. (024) 76480785, 76480788

email : jurusan@geodesi.ft.undip.ac.id

2018

1
Coding yang digunakan:

function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm);

TIDAMP=1.2;
BLATA=BLATA;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+TT*0.0
072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;
H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+BLONG)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);
C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);

2
 D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-2*H+P)+M*M*cos(2*(S-
H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(hh*100.0+mm+1200.0,2400.0);
addday=floor((hh*100+mm+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(YYYY,MM,DD);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
End

3
Tanggal 30, 10, 2018 pukul 01;00

function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,1,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

4
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(1*100.0+0+1200.0,2400.0);
addday=floor((1*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = 0.0124

5
Tanggal 30, 10, 2018 pukul 02;00

function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,2,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

6
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(2*100.0+0+1200.0,2400.0);
addday=floor((2*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = 0.0310

7
Tanggal 30, 10, 2018 pukul 03;00

function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,3,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

8
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(3*100.0+0+1200.0,2400.0);
addday=floor((3*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = 0.0331

9
Tanggal 30, 10, 2018 pukul 04;00
function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,4,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

10
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(4*100.0+0+1200.0,2400.0);
addday=floor((4*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = 0.0198

11
Tanggal 30, 10, 2018 pukul 05;00

function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,5,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);

12
 CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(5*100.0+0+1200.0,2400.0);
addday=floor((5*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = -0.0042

13
Tanggal 30, 10, 2018 pukul 06;00
function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,6,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

14
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(6*100.0+0+1200.0,2400.0);
addday=floor((6*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = -0.0319

15
Tanggal 30, 10, 2018 pukul 07;00
function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,7,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

16
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(7*100.0+0+1200.0,2400.0);
addday=floor((7*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = -0.0562

17
Tanggal 30, 10, 2018 pukul 08;00
function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,8,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

18
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(8*100.0+0+1200.0,2400.0);
addday=floor((8*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = -0.0714

19
Tanggal 30, 10, 2018 pukul 09;00
function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,9,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

20
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(9*100.0+0+1200.0,2400.0);
addday=floor((9*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = -0.0743

21
Tanggal 30, 10, 2018 pukul 10;00

function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,10,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

22
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(10*100.0+0+1200.0,2400.0);
addday=floor((10*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = -0.0649

23
Tanggal 30, 10, 2018 pukul 11;00
function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,11,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

24
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(11*100.0+0+1200.0,2400.0);
addday=floor((11*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = -0.0459

25
Tanggal 30, 10, 2018 pukul 12;00

function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,12,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

26
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(12*100.0+0+1200.0,2400.0);
addday=floor((12*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = -0.0221

27
Tanggal 30, 10, 2018 pukul 13;00

function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(6,11,2018,13,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

28
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(13*100.0+0+1200.0,2400.0);
addday=floor((13*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = 0.0005

29
Tanggal 30, 10, 2018 pukul 14;00
function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,14,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

30
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(14*100.0+0+1200.0,2400.0);
addday=floor((14*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = 0.0165

31
Tanggal 30, 10, 2018 pukul 15;00

function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,15,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);

32
 CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(15*100.0+0+1200.0,2400.0);
addday=floor((15*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = 0.021

33
Tanggal 30, 10, 2018 pukul 16;00
function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,16,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

34
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(16*100.0+0+1200.0,2400.0);
addday=floor((16*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = 0.0128

35
Tanggal 30, 10, 2018 pukul 17;00
function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,17,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

36
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(17*100.0+0+1200.0,2400.0);
addday=floor((17*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = -0.0078

37
Tanggal 30, 10, 2018 pukul 18;00
function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,18,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

38
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(18*100.0+0+1200.0,2400.0);
addday=floor((18*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = -0.0369

39
Tanggal 30, 10, 2018 pukul 19;00

function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,19,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

40
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(19*100.0+0+1200.0,2400.0);
addday=floor((19*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = -0.0683

41
Tanggal 30, 10, 2018 pukul 20;00
function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,20,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

42
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(20*100.0+0+1200.0,2400.0);
addday=floor((20*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = -0.0950

43
Tanggal 30, 10, 2018 pukul 21;00

function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,21,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

44
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(21*100.0+0+1200.0,2400.0);
addday=floor((21*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = -0.1107

45
Tanggal 30, 10, 2018 pukul 22;00
function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,22,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

46
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(22*100.0+0+1200.0,2400.0);
addday=floor((22*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = -0.1116

47
Tanggal 30, 10, 2018 pukul 23;00
function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,23,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

48
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(23*100.0+0+1200.0,2400.0);
addday=floor((23*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = -0.0974

49
Tanggal 30, 10, 2018 pukul 24;00

function TidalGravity=Tidaleffect(BLATA,BLONG,DD,MM,YYYY,hh,mm)
[TIM,JT]=TIMANDJT(30,10,2018,24,0);

TIDAMP=1.2;
BLATA=-6.9467;
ONE=1;
PI=pi;
RDN=PI/180;
RDNS=RDN/3600;
TT=JT/(36525*1440);
LAT=-BLATA*RDN;
NN=((259-mod(5*TT,ONE)*360)*3600+657.12+TT*(-
482912.63+TT*(7.58+TT*0.008)))*RDNS;

S=((270.0+mod(1336.0*TT,ONE)*360.0)*3600.0+1571.72+TT*(1108406.05+TT*(7.128+T
T*0.0072)))*RDNS;
P=((334.0+mod(11*TT,ONE)*360.0)*3600.0+1186.42+TT*(392522.51+TT*(-37.15-
TT*0.036)))*RDNS;

H=((279.0+mod(100*TT,ONE)*360.0)*3600.0+2508.05+TT*(2768.11+TT*1.08))*RDNS ;
CSN=cos(NN);
SNN=sin(NN);
I=5.145*RDN;
OMEGA=23.452*RDN;
SNI=sin(I);
CSI=cos(I);
SNO=sin(OMEGA);
CSO=cos(OMEGA);
E=0.05490;
M=0.074804;
CSII=CSO*CSI-SNO*SNI*CSN;
SNII=sqrt(ONE-CSII*CSII);
SNNU=SNI*SNN/SNII;
CSNU=sqrt(ONE-SNNU*SNNU);
SIG=S-NN+2*atan((SNO*SNN/SNII)/(ONE+CSN*CSNU+SNN*SNNU*CSO));
L=SIG+2*E*sin(S-P)+1.25*E*E*sin(2*(S-P))+3.75*M*E*sin(S-
2*H+P)+1.375*M*M*sin(2*(S-H));
E1=(16751.04+TT*(-41.80-TT*0.126))*(10^-6);
P1=((281.0*60+13)*60+15+TT*(6189.03+TT*(1.63+TT*0.012)))*RDNS;
L1=H+2*E1*sin(H-P1);
A=floor(TIM/100);
B=(TIM-A*100);
T0=A+B/60-12.0;
T=(15*(T0-12)+110.4287)*RDN;
CHI=T+H-atan(SNNU/CSNU);
CHI1=T+H;
SNLAT=sin(LAT);
CSLAT=cos(LAT);
SNL=sin(L);
CSL=cos(L);
SNL1=sin(L1);
CSL1=cos(L1);
CM=SNLAT*SNII*SNL+CSLAT*(CSL*cos(CHI)+SNL*sin(CHI)*CSII);
CS=SNLAT*SNO*SNL1+CSLAT*(CSL1*cos(CHI1)+SNL1*sin(CHI1)*CSO);

50
 C=3.84402*(10^10);
C1=1.495*(10^13);
AD=((C*(1-E*E))^-1);
AD1=((C1*(ONE-E1*E1))^-1);
D=(C^-1)+AD*(E*cos(S-P)+E*E*cos(2*(S-P))+1.875*M*E*cos(S-
2*H+P)+M*M*cos(2*(S-H)));
D1=(C1^-1)+AD1*E1*cos(H-P1);
R=6.37827*(10^8)/sqrt(ONE+0.006738*SNLAT*SNLAT);
TID=-6.670*(10^-8)*((7.3537*(10^25)*D*D*R*D*(3.0*CM*CM-
ONE+1.5*R*D*CM*(5*CM*CM-3)))+(1.993*(10^33)*D1*D1*R*D1*(3*CS*CS-1)))*(10^3);
TidalGravity=TID*TIDAMP;

function [TIM,JT]=TIMANDJT(DD,MM,YYYY,hh,mm)
TIM=mod(24*100.0+0+1200.0,2400.0);
addday=floor((24*100+0+1200)/2400);
datenum1900=datenum(1900,01,01);
datenumtoday=datenum(2018,11,6);
dayspast=datenumtoday-datenum1900+addday;
hourspast=dayspast*24;
HH1=floor(TIM/100);
MM1=TIM-HH1*100;
JT=hourspast*60.0+HH1*60.0+MM1;
end

end

ans = -0.0712

51
Dari perhitungan pasut bumi diperoleh data pasut bumi (mgal) harian dan grafik earth tide
sebagai berikut:

Tabel 1 Hasil perhitungan Earth Tide

01:00:00 0.0124
02:00:00 0.0310
03:00:00 0.0331
04:00:00 0.0198
05:00:00 -0.0042
06:00:00 -0.0319
07:00:00 -0.0562
08:00:00 -0.0714
09:00:00 -0.0743
10:00:00 -0.0649
11:00:00 -0.0459
12:00:00 -0.0221
13:00:00 0.0005
14:00:00 0.0165
15:00:00 0.0212
16:00:00 0.0128
17:00:00 -0.0078
18:00:00 -0.0369
19:00:00 -0.0683
20:00:00 -0.0950
21:00:00 -0.1107
22:00:00 -0.1116
23:00:00 -0.0974
00:00:00 -0.0712

Grafik Earth Tide (Mgal)

0.0400

0.0200

0.0000
6:00:00 AM
1:00:00 AM
2:00:00 AM
3:00:00 AM
4:00:00 AM
5:00:00 AM

7:00:00 AM
8:00:00 AM
9:00:00 AM
10:00:00 AM
11:00:00 AM

12:00:00 AM
6:00:00 PM
12:00:00 PM
1:00:00 PM
2:00:00 PM
3:00:00 PM
4:00:00 PM
5:00:00 PM

7:00:00 PM
8:00:00 PM
9:00:00 PM
10:00:00 PM
11:00:00 PM

-0.0200

-0.0400

-0.0600

-0.0800

-0.1000

-0.1200

52