| 460 |
|
|
| 461 |
endif |
endif |
| 462 |
|
|
| 463 |
100 return |
c 100 return |
| 464 |
|
return |
| 465 |
end |
end |
| 466 |
|
|
| 467 |
*** * * * *** * * * *** * * * *** * * * *** * * * *** * * * *** |
*** * * * *** * * * *** * * * *** * * * *** * * * *** * * * *** |
| 508 |
|
|
| 509 |
endif |
endif |
| 510 |
|
|
| 511 |
100 return |
c 100 return |
| 512 |
|
return |
| 513 |
end |
end |
| 514 |
*** * * * *** * * * *** * * * *** * * * *** * * * *** * * * *** |
*** * * * *** * * * *** * * * *** * * * *** * * * *** * * * *** |
| 515 |
c real function riseta(ic,angle) |
c real function riseta(ic,angle) |
| 559 |
endif |
endif |
| 560 |
|
|
| 561 |
|
|
| 562 |
100 return |
c 100 return |
| 563 |
|
return |
| 564 |
end |
end |
| 565 |
|
|
| 566 |
*** * * * *** * * * *** * * * *** * * * *** * * * *** * * * *** |
*** * * * *** * * * *** * * * *** * * * *** * * * *** * * * *** |
| 721 |
$ ,cog2-iadd,' -->',pfaeta2 |
$ ,cog2-iadd,' -->',pfaeta2 |
| 722 |
|
|
| 723 |
|
|
| 724 |
100 return |
c 100 return |
| 725 |
|
return |
| 726 |
end |
end |
| 727 |
|
|
| 728 |
*** * * * *** * * * *** * * * *** * * * *** * * * *** * * * *** |
*** * * * *** * * * *** * * * *** * * * *** * * * *** * * * *** |
| 825 |
if(DEBUG.EQ.1)print*,'ETA3 (ic ',ic,' ang',angle,')' |
if(DEBUG.EQ.1)print*,'ETA3 (ic ',ic,' ang',angle,')' |
| 826 |
$ ,cog3-iadd,' -->',pfaeta3 |
$ ,cog3-iadd,' -->',pfaeta3 |
| 827 |
|
|
| 828 |
100 return |
c 100 return |
| 829 |
|
return |
| 830 |
end |
end |
| 831 |
|
|
| 832 |
*** * * * *** * * * *** * * * *** * * * *** * * * *** * * * *** |
*** * * * *** * * * *** * * * *** * * * *** * * * *** * * * *** |
| 935 |
if(DEBUG.EQ.1)print*,'ETA4 (ic ',ic,' ang',angle,')' |
if(DEBUG.EQ.1)print*,'ETA4 (ic ',ic,' ang',angle,')' |
| 936 |
$ ,cog4-iadd,' -->',pfaeta4 |
$ ,cog4-iadd,' -->',pfaeta4 |
| 937 |
|
|
| 938 |
100 return |
c 100 return |
| 939 |
|
return |
| 940 |
end |
end |
| 941 |
|
|
| 942 |
|
|
| 1014 |
if(sl1.gt.sr1)then |
if(sl1.gt.sr1)then |
| 1015 |
if((sl1+sc).ne.0)COG = -sl1/(sl1+sc) |
if((sl1+sc).ne.0)COG = -sl1/(sl1+sc) |
| 1016 |
elseif(sl1.lt.sr1)then |
elseif(sl1.lt.sr1)then |
| 1017 |
if((sc+sr1).ne.0)COG = sr1/(sc+sr1) |
if((sc+sr1).ne.0)COG = sr1/(sc+sr1) |
| 1018 |
elseif( sl1.eq.sr1.and.sl1.ne.-9999.)then |
elseif( sl1.eq.sr1.and.sl1.ne.-9999.)then |
| 1019 |
if( clsigma(indmax(ic)-1).lt.clsigma(indmax(ic)+1) |
if( clsigma(indmax(ic)-1).lt.clsigma(indmax(ic)+1) |
| 1020 |
$ .and.(sl1+sc).ne.0 )cog = -sl1/(sl1+sc) |
$ .and.(sl1+sc).ne.0 )cog = -sl1/(sl1+sc) |
| 2217 |
if(DEBUG.eq.1)print*,'LANDI (ic ',ic,' ang',angle,') -->',pfacorr |
if(DEBUG.eq.1)print*,'LANDI (ic ',ic,' ang',angle,') -->',pfacorr |
| 2218 |
|
|
| 2219 |
|
|
| 2220 |
100 return |
c 100 return |
| 2221 |
|
return |
| 2222 |
end |
end |