Télécharger d2vint.eso

Retour à la liste

Numérotation des lignes :

d2vint
  1. C D2VINT SOURCE BP208322 22/09/20 21:15:02 11462
  2. C DEVINT SOURCE PLAF 90/05/03 21:11:36
  3.  
  4. *--------------------------------------------------------------------*
  5. * *
  6. * Operateur DYNE : algorithme de Fu - de Vogelaere *
  7. * ________________________________________________ *
  8. * *
  9. * Interpolation d'un liste REEL. *
  10. * Remplissage du segment de travail KTRAV. *
  11. * *
  12. * Param}tres: *
  13. * *
  14. * e ICHATE pointeur sur un LISTREEL, abscisse *
  15. * e ICHAFO pointeur sur un LISTREEL, ordonnee *
  16. * e KTNUM segment contenant les param}tres numeriques *
  17. * s KTRAV segment de travail *
  18. * e ICHAR numero du chargement *
  19. * *
  20. * Auteur, date de creation: *
  21. * Lionel VIVAN, le 17 avril 1990 *
  22. * *
  23. * Parallélisation : BP, 2022-09-19 *
  24. * *
  25. *--------------------------------------------------------------------*
  26.  
  27. SUBROUTINE D2VINT(ICHATE,ICHAFO,KTNUM,KTRAV,ICHAR)
  28. IMPLICIT INTEGER(I-N)
  29. IMPLICIT REAL*8(A-H,O-Z)
  30.  
  31. -INC PPARAM
  32. -INC CCOPTIO
  33. -INC SMLREEL
  34.  
  35. SEGMENT,MTNUM
  36. REAL*8 XDT(NPC1),XTEMPS(NPC1)
  37. ENDSEGMENT
  38. SEGMENT,MTRAV
  39. REAL*8 FTCHG(NCHAR,NPC1)
  40. REAL*8 XFORCA(NPREF,NCHAR)
  41. ENDSEGMENT
  42.  
  43. *--------------------------------------------------------------------*
  44. * lecture arguments d'entree/sortie
  45. *--------------------------------------------------------------------*
  46.  
  47. * discretisation temporelle de DYNE
  48. MTNUM = KTNUM
  49. NPC1 = XDT(/1)
  50. PASTEM = XDT(1)
  51.  
  52. * objets du chargement fourni
  53. MLREE1 = ICHATE
  54. SEGACT MLREE1
  55. NF = MLREE1.PROG(/1)
  56. MLREE2 = ICHAFO
  57. SEGACT MLREE2
  58.  
  59. cbp,2022 SEGINI MTRAV
  60. MTRAV = KTRAV
  61. *
  62. *--------------------------------------------------------------------*
  63. * interpolation du chargement
  64. *--------------------------------------------------------------------*
  65. *
  66. N1 = 1
  67. N2 = 2
  68.  
  69. DO 10 IT = 1,NPC1
  70. TEMP = XTEMPS(IT)
  71. CALL INTLIN(TEMP,ICHATE,ICHAFO,NF,N1,N2, FTEM,IRET)
  72. IF (IRET.EQ.0) THEN
  73. * temps en dehors de la liste
  74. CALL ERREUR(208)
  75. RETURN
  76. ENDIF
  77. FTCHG(ICHAR,IT) = FTEM
  78. 10 CONTINUE
  79.  
  80. SEGDES,MLREE1,MLREE2
  81.  
  82. END
  83.  
  84.  
  85.  
  86.  

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