Télécharger crepo2.eso

Retour à la liste

Numérotation des lignes :

  1. C CREPO2 SOURCE CHAT 05/01/12 22:29:54 5004
  2.  
  3. C-----------------------------------------------------------------------
  4. C C R E P O 2
  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) en donnant egalement la densite
  11. C associee au point a creer.
  12. C
  13. C Mode d'appel :
  14. C --------------
  15. C CALL CREPO2 (R1,R2,R3,DENS, 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 DENS REEL DP (E) densite du point a creer
  23. C IPOIN ENTIER (S) numero du point cree (dans SMCOORD)
  24. C
  25. C Auteur, date de creation :
  26. C --------------------------
  27. C PASCAL MANIGOT 19 AVRIL 1985
  28. C 10/2003 : Modifications pour traiter le cas IDIM=1
  29. C
  30. C Langage : FORTRAN77 + ESOPE
  31. C-----------------------------------------------------------------------
  32.  
  33. SUBROUTINE CREPO2 (R1,R2,R3,DENS, IPOIN)
  34.  
  35. IMPLICIT INTEGER(I-N)
  36. IMPLICIT REAL*8 (A-H,O-Z)
  37.  
  38. -INC CCOPTIO
  39. -INC SMCOORD
  40.  
  41. IF ((IDIM.LT.1).OR.(IDIM.GT.3)) THEN
  42. CALL ERREUR (219)
  43. RETURN
  44. END IF
  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)=DENS
  54. IPOIN=NBPTS
  55.  
  56. RETURN
  57. END
  58.  
  59.  
  60.  
  61.  

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