real(kind=8) :: the,al,om,detApi,wart,y2,wykl,radmax
real(kind=8) :: tant,zz1,W1i,radius,zk,fac,dV,sum,sum1
real(kind=8) :: which_lobe
lcheck=.false.
lprint=.false.
fail=.false.
real(kind=8) :: the,al,om,detApi,wart,y2,wykl,radmax
real(kind=8) :: tant,zz1,W1i,radius,zk,fac,dV,sum,sum1
real(kind=8) :: which_lobe
lcheck=.false.
lprint=.false.
fail=.false.