| 799 |
ZP(NP)=ZM(I) |
ZP(NP)=ZM(I) |
| 800 |
ENDIF |
ENDIF |
| 801 |
ENDDO |
ENDDO |
| 802 |
|
IFLAG=0 !no debug mode |
| 803 |
CALL TRICIRCLE(NP,XP,ZP,AP,RP,CHI,XC,ZC,RADIUS,IFLAG) |
CALL TRICIRCLE(NP,XP,ZP,AP,RP,CHI,XC,ZC,RADIUS,IFLAG) |
| 804 |
print*,' circle: ',XC,ZC,RADIUS,' --- ',CHI,IFLAG |
c print*,' circle: ',XC,ZC,RADIUS,' --- ',CHI,IFLAG |
| 805 |
IF(IFLAG.NE.0)GOTO 10 !straigth fit |
IF(IFLAG.NE.0)GOTO 10 !straigth fit |
| 806 |
if(CHI.gt.100)GOTO 10 !straigth fit |
if(CHI.gt.100)GOTO 10 !straigth fit |
| 807 |
ARG = RADIUS**2-(ZINI-ZC)**2 |
ARG = RADIUS**2-(ZINI-ZC)**2 |