Télécharger crepo1.eso

Retour à la liste

Numérotation des lignes :

  1. C CREPO1 SOURCE BP208322 16/11/18 21:16:09 9177
  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 --------------
  15. C CALL CREPO1 (R1,R2,R3, IPOIN)
  16. C
  17. C Parametres : (E)=Entree (S)=Sortie
  18. C ------------
  19. C R1 REEL DP (E) 1ere coordonnee du point
  20. C R2 REEL DP (E) 2e coordonnee du point (si IDIM=2 ou 3)
  21. C R3 REEL DP (E) 3e coordonnee du point (si IDIM=3)
  22. C IPOIN ENTIER (S) numero du point cree (dans SMCOORD)
  23. C
  24. C Auteur, date de creation :
  25. C --------------------------
  26. C PASCAL MANIGOT 19 AVRIL 1985
  27. C 10/2003 : Modifications pour traiter le cas IDIM=1
  28. C
  29. C Langage : FORTRAN77 + ESOPE
  30. C-----------------------------------------------------------------------
  31.  
  32. SUBROUTINE CREPO1 (R1,R2,R3, IPOIN)
  33.  
  34. IMPLICIT INTEGER(I-N)
  35. IMPLICIT REAL*8 (A-H,O-Z)
  36.  
  37. -INC CCGEOME
  38. -INC CCOPTIO
  39. -INC SMCOORD
  40.  
  41. IF ((IDIM.LT.1).OR.(IDIM.GT.3)) THEN
  42. CALL ERREUR (219)
  43. RETURN
  44. ENDIF
  45.  
  46. idimp1=IDIM+1
  47. NBPTS=XCOOR(/1)/idimp1+1
  48. SEGADJ MCOORD
  49. IPoin=(NBPTS-1)*idimp1
  50. XCOOR(IPoin+1)=R1
  51. IF (IDIM.GE.2) XCOOR(IPoin+2)=R2
  52. IF (IDIM.GE.3) XCOOR(IPoin+3)=R3
  53. XCOOR(IPoin+idimp1)=DENSIT
  54. IPOIN=NBPTS
  55.  
  56. RETURN
  57. END
  58.  
  59.  
  60.  
  61.  
  62.  
  63.  
  64.  
  65.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales