Télécharger crepo2.eso

Retour à la liste

Numérotation des lignes :

crepo2
  1. C CREPO2 SOURCE PV 20/03/24 21:16:41 10554
  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.  
  39. -INC PPARAM
  40. -INC CCOPTIO
  41. -INC SMCOORD
  42.  
  43. IF ((IDIM.LT.1).OR.(IDIM.GT.3)) THEN
  44. CALL ERREUR (219)
  45. RETURN
  46. END IF
  47.  
  48. idimp1=IDIM+1
  49. segact mcoord*mod
  50. NBPTS=nbpts+1
  51. SEGADJ MCOORD
  52. IPoin=(NBPTS-1)*idimp1
  53. XCOOR(IPoin+1)=R1
  54. IF (IDIM.GE.2) XCOOR(IPoin+2)=R2
  55. IF (IDIM.GE.3) XCOOR(IPoin+3)=R3
  56. XCOOR(IPoin+idimp1)=DENS
  57. IPOIN=NBPTS
  58.  
  59. RETURN
  60. END
  61.  
  62.  
  63.  
  64.  
  65.  

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