1 : C CREPO1 SOURCE CHAT 07/10/29 21:17:35 5940 2 : 3 : C----------------------------------------------------------------------- 4 : C C R E P O 1 5 : C ----------- 6 : C 7 : C Fonction : 8 : C ---------- 9 : C Creation d'un POINT (equivalent au sous-programme CREPOI sans 10 : C les lectures et ecritures). La densite du point cree est la 11 : C valeur courante (stockee dans CCOPTION). 12 : C 13 : C Mode d'appel : 14 : C -------------- 14 : C CALL CREPO1 (R1,R2,R3, IPOIN) 15 : C 16 : C Parametres : (E)=Entree (S)=Sortie 17 : C ------------ 18 : C R1 REEL DP (E) 1ere coordonnee du point 19 : C R2 REEL DP (E) 2e coordonnee du point (si IDIM=2 ou 3) 20 : C R3 REEL DP (E) 3e coordonnee du point (si IDIM=3) 21 : C IPOIN ENTIER (S) numero du point cree (dans SMCOORD) 22 : C 23 : C Auteur, date de creation : 24 : C -------------------------- 25 : C PASCAL MANIGOT 19 AVRIL 1985 26 : C 10/2003 : Modifications pour traiter le cas IDIM=1 27 : C 28 : C Langage : FORTRAN77 + ESOPE 29 : C----------------------------------------------------------------------- 30 : 31 : SUBROUTINE CREPO1 (R1,R2,R3, IPOIN) 32 : 33 : IMPLICIT INTEGER(I-N) 34 : IMPLICIT REAL*8 (A-H,O-Z) 35 : 36 : -INC CCGEOME 37 : -INC CCOPTIO 38 : -INC SMCOORD 39 : 40 : IF ((IDIM.LT.1).OR.(IDIM.GT.3)) THEN 40 : CALL ERREUR (219) 41 : RETURN 42 : ENDIF 43 : 44 : idimp1=IDIM+1 45 : NBPTS=XCOOR(/1)/idimp1+1 46 : SEGADJ MCOORD 47 : IPoin=(NBPTS-1)*idimp1 48 : XCOOR(IPoin+1)=R1 49 : IF (IDIM.GE.2) XCOOR(IPoin+2)=R2 50 : IF (IDIM.GE.3) XCOOR(IPoin+3)=R3 51 : XCOOR(IPoin+idimp1)=DENSIT 52 : IPOIN=NBPTS 53 : 54 : RETURN 55 : END 56 : 57 : 58 : 59 : 60 :
© Cast3M 2003 - All rights reserved.
Disclaimer