if (abs(fai-ffai)<PI/6. ! 这里,一次用这么多逻辑词,是在.or. 那里断句吗?
+ .or. (fai-ffai)>5*PI/6 .and. (fai-ffai)<PI
+ .or. (fai-ffai)>-PI .and. (fai-ffai)<-5*PI/6) then
if (ww<=w0) then
mu=4.06
else
mu=-2.34
end if
p=real(9.77*(ww/w0)**mu)
normp=2**(1-2*p)*PI*gamma(2*p+1)/gamma(p+1)**2
Direc_new=(cos((fai-ffai)/2))**(2*p)/normp
else
Direc_new=0.
end if