procon
C PROCON SOURCE PV090527 23/07/13 21:15:04 11708 IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) DIMENSION XX(*),YY(*) data ISAV /0/ save ISAV data sqxgr /0/ save sqxgr -INC CCREEL PROCO = 0.D0 IF(ISAV .EQ. 0)THEN sqxgr=sqrt(xgrand*xzprec) isav=1 ENDIF IF (NSTRS .GT. 3) GOTO 1 DO I=1,NSTRS XI = min(sqxgr,max(XX(i),-sqxgr)) YI = min(sqxgr,max(YY(i),-sqxgr)) PROCO = PROCO + XI*YI ENDDO PROCON = PROCO RETURN 1 DO I=4,NSTRS XI = min(sqxgr,max(XX(i),-sqxgr)) YI = min(sqxgr,max(YY(i),-sqxgr)) PROCO = PROCO + XI*YI ENDDO PROCO = PROCO * 2.D0 DO I=1,3 XI = min(sqxgr,max(XX(i),-sqxgr)) YI = min(sqxgr,max(YY(i),-sqxgr)) PROCO = PROCO + XI*YI ENDDO PROCON = PROCO END
© Cast3M 2003 - Tous droits réservés.
Mentions légales