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