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 |