function transitioncal
% n w a phe ro r1 rk
h1=findobj(0,'tag','now');
m=findobj(h1,'tag','EditText2');
phe=get(m,'string');
phe=str2num(phe)
m=findobj(h1,'tag','EditText3');
AK=get(m,'string');
AK=str2num(AK)
h2=findobj(0,'tag','tra');
tkk=findobj(h2,'tag','EditText1hour');
wk=get(tkk,'string');
w=str2num(wk)
t=findobj(h2,'tag','EditText2day');
nk=get(t,'string');
n=str2num(nk)
del=23.45*sin((360*(n+79+284)/365)*3.14/180);
delrad=del*pi/180;
wdegree=(w-12)*15;
wrad=wdegree*(pi/180);
zenith=(180/pi)*acos(sin(delrad)*sin(phe*pi/180)+cos(delrad)*cos(phe*pi/180)*cos(wrad));
a0star=0.4237-0.00821*((6-AK)^2);
a1star=0.5055+0.00595*((6.5-AK)^2);
kstar=0.2711+0.01858*((2.5-AK)^2);
load a:\data1
pop1=findobj(0,'tag','PopupMenutransition');
so=get(pop1,'value');
st=get(pop1,'string');
%st{so}
eval(['kk.' st{so}])
eval(['ro=kk.' st{so} '.ro'])
eval(['r1=kk.' st{so} '.r1'])
eval(['rk=kk.' st{so} '.rk'])
a0=a0star*ro
a1=a1star*r1
k=kstar*rk
costetaz=cos(zenith*pi/180);
tb=a0+a1*(exp(-k/costetaz))
msgbox(num2str(tb),'(tb)atmospher beam transition')
td=0.271-0.294*tb;
msgbox(num2str(td),'(td)atmospher beam transition')
