Visibility of Luna and Lunar Crescent

)
{
ip nt(iprr4,""/od:%d{'/o8.3lo\t",in(thou),in((lbou'in(lhou))'60),nma8-len)i
'
{
falt=malt+.1 ; fam=hem+.1:
kndisl=9o.o-fali:
Glomonlaco(sin(matl'PUl 80 0)rsin(fal1'Ptl80-0)+co(malt'Pv1 80 0)'co(i
ah.P/180.0.)r@s((n'a-faa)rPvl80.)))r l80.o/lli
relonsui=(acos(si nGan' P V 1 80.0)* sin(lalt'P V l E 0.0)+cos(salr' PI/ 1 80 0)'cos(ral(
iPl/180.0)'co(G@-ran)'Pvl80.)))1 180.o/PI;
sacon=1.0/(cos(zendist'P/180.0)10.0286rq(-10-51@s(andist1Pvl80.0)));
aacom-1.o/(cos(rndisl.Ptl80.0)+0.0l2lrexpc24.5xcos(4ndistiPI/180.0)));
oacon=pow(1.0-posGin(Endisl'PVl 80.0)/(1.0+20.0/6378.0),2.0),-0.5)i
fo(i=ori<5;i++)
{
kFo.1066'ex(-palr8200 0)rpow(Mchlil/0.55.4.0);
lc=o l'pus(*6chlil/0.55.-l.r)'dptpalL/1<00 0\
ka=kaipow(1.0-0.12llog(phunVl 00.0),1.i3)
'(1.0+0.13'si.Galpha'Pl/180 0).(pla/aab(plat)))i
ko{2$htir(3.or0.4'(plarr(PV180.0)1@s(elph.tPUl 80.0)-
@s(3.orplaCPvl 80.0))Y3.0;
kw:Brehlil*0.94.(ph$!i 100.0)*cxp(pteDp/15.0)iexpcparr/E200.0)j
li$h[i]=kd.ka+ko+kw:
dnshlil=krSeom+k .dcon+ko.olcom+kw.gacom;
)
npo$i=1.o/Go(EndistxPVl80.0)+0.025rexp(-t
ii(mah<=o.o) mhpos€f=4o.o;
l.0icos(zbdist*Pv180.0))j
nnposeF I -0lcos(90.Gtu1r). PVI 80.0)+0.025rexpt I L0tco((90.0-
mlrxPyl80.0)));
ir(elt<=o.o) snposef=4o.0i
snposcf=1.0(cos((90.onatrPr/l 80.0)r{.025.ex(-t 1.01@(90.0-
ex)1PVl 80.0)));
lb(i=oj<s:i+)
{
.idrb=boschtil.(1.or0.3'co<6.283.(ye&j992.0y1 1.0));
nigb|b=nisntb*(0.4+o.6/pow(1.G0.961po{sir(andistrpvl8o.o),2.0),-o.s));
nishrbrightb.pow(l0.0,c0_41t$htiirslpost)i
oMas=-12.?3+0.0261rabs(l 80.o-etonsp)+4.0.pow(180.0-
elonsp),4.0).(!ow(t 0.0,-9.0))i
mmag=nnmas+.nshlil;
co'@=pow(t 0.0,c0.4.k$blil.@pceo)
229