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 |