def terbit(y,m,d,B,Li,Z,H):
waktuter = waktu_sholat.transit(y,m,d,B,Li,Z,H)
delta = waktu_sholat.delta(jdconverter.masehikejd(y,m,d))
L = Li*np.pi/180
ht = (-0.8333 - 0.0347*np.sqrt(H))*np.pi/180
hat = np.arccos((np.sin(ht)-np.sin(L)*np.sin(delta))/(np.cos(L)*np.cos(delta))) * 180/np.pi
waktu_terbit = waktu_sholat.wkt(waktuter - hat/15)
return waktu_terbit
def subuh(y,m,d,B,Li,Z,H):
waktusu = waktu_sholat.transit(y,m,d,B,Li,Z,H)
delta = waktu_sholat.delta(jdconverter.masehikejd(y,m,d))
L = Li*np.pi/180
hs = -20*np.pi/180
has = np.arccos((np.sin(hs)-np.sin(L)*np.sin(delta))/(np.cos(L)*np.cos(delta))) * 180/np.pi
waktu_subuh = waktu_sholat.wkt(waktusu - has/15)
return waktu_subuh
def isya(y,m,d,B,Li,Z,H):
waktuis = waktu_sholat.transit(y,m,d,B,Li,Z,H)
delta = waktu_sholat.delta(jdconverter.masehikejd(y,m,d))
L = Li*np.pi/180
hi = -18*np.pi/180
hai = np.arccos((np.sin(hi)-np.sin(L)*np.sin(delta))/(np.cos(L)*np.cos(delta))) * 180/np.pi
waktu_isya = waktu_sholat.wkt(waktuis + hai/15)
return waktu_isya
def ashar(y,m,d,B,Li,Z,H):
waktuas = waktu_sholat.transit(y,m,d,B,Li,Z,H)
delta = waktu_sholat.delta(jdconverter.masehikejd(y,m,d))
L = Li*np.pi/180
ha = np.arctan(1/(1+np.tan(abs(delta-L))))
haa = np.arccos((np.sin(ha)-np.sin(L)*np.sin(delta))/(np.cos(L)*np.cos(delta))) * 180/np.pi
waktu_ashar = waktu_sholat.wkt(waktuas + haa/15)
return waktu_ashar
def maghrib(y,m,d,B,Li,Z,H):
waktuma = waktu_sholat.transit(y,m,d,B,Li,Z,H)
akarhm = np.sqrt(H)
hm = (-0.0347*akarhm - 0.8333)*np.pi/180
delta = waktu_sholat.delta(jdconverter.masehikejd(y,m,d))
L = Li*np.pi/180
ham = np.arccos((np.sin(hm)-np.sin(L)*np.sin(delta))/(np.cos(L)*np.cos(delta))) * 180/np.pi
waktu_maghrib = waktu_sholat.wkt(waktuma + ham/15)
return waktu_maghrib
def transit(y,m,d,B,Li,Z,H):
jd1 = jdconverter.masehikejd(y,m,d)
T = 2*np.pi*(jd1-2451545)/365.25
delta = (0.37877 + 23.264*np.sin(57.297*T*np.pi/180-79.547*np.pi/180)
+ 0.3812*np.sin(2*57.297*T*np.pi/180-82.682*np.pi/180)
+ 0.17132*np.sin(3*57.297*T*np.pi/180-59.722*np.pi/180))*np.pi/180
U = (jd1 - 2451545)/36525
L0 = (280.46607 + 36000.7698*U)*mt.pi/180
ET = (-(1789 + 237*U)*mt.sin(L0)-(7146 - 62*U)*mt.cos(L0) + (9934-14*U)*mt.sin(2*L0) - (29+ 5*U)
*mt.cos(2*L0) + (74 + 10*U)*mt.sin(3*L0) + (320 - 4*U)*mt.cos(3*L0) - 212*mt.sin(4*L0))/1000
transit1 = 12 + Z - B/15 - ET/60
L = Li*np.pi/180
#waktu dzuhur
#dzuhur = transit1
#waktu_dzuhur = waktu_sholat.wkt(dzuhur)
#waktu ashar
#ha = np.arctan(1/(1+np.tan(abs(delta-L))))
#haa = np.arccos((np.sin(ha)-np.sin(L)*np.sin(delta))/(np.cos(L)*np.cos(delta))) * 180/np.pi
#waktu_ashar = waktu_sholat.wkt(transit1 + haa/15)
#waktu maghrib
#akarhm = np.sqrt(H)
#hm = (-0.0347*akarhm - 0.8333)*np.pi/180
#ham = np.arccos((np.sin(hm)-np.sin(L)*np.sin(delta))/(np.cos(L)*np.cos(delta))) * 180/np.pi
#waktu_maghrib = waktu_sholat.wkt(transit1 + ham/15)
#waktu isya
#hi = -18*np.pi/180
#hai = np.arccos((np.sin(hi)-np.sin(L)*np.sin(delta))/(np.cos(L)*np.cos(delta))) * 180/np.pi
#waktu_isya = waktu_sholat.wkt(transit1 + hai/15)
#waktu subuh
#hs = -20*np.pi/180
#has = np.arccos((np.sin(hs)-np.sin(L)*np.sin(delta))/(np.cos(L)*np.cos(delta))) * 180/np.pi
#waktu_subuh = waktu_sholat.wkt(transit1 - has/15)
#terbit matahari
#ht = (-0.8333 - 0.0347*np.sqrt(H))*np.pi/180
#hat = np.arccos((np.sin(ht)-np.sin(L)*np.sin(delta))/(np.cos(L)*np.cos(delta))) * 180/np.pi
#waktu_terbit = waktu_sholat.wkt(transit1 - hat/15)
#waktu_array = ['subuh', 'terbit', 'dzuhur', 'ashar', 'maghrib', 'isya']
#waktu_array2=[waktu_subuh, waktu_terbit, waktu_dzuhur, waktu_ashar, waktu_maghrib, waktu_isya]
#gunakan dictionary
return transit1