Télécharger dyne21.eso

Retour à la liste

Numérotation des lignes :

  1. C DYNE21 SOURCE BP208322 15/07/22 21:15:33 8586
  2. SUBROUTINE DYNE21(ITLA,PDT,MTRA,NLIAA,NXPALA,NPLAA,NIPALA)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. *--------------------------------------------------------------------*
  6. * *
  7. * Op{rateur DYNE : algorithme de Fu - de Vogelaere *
  8. * ________________________________________________ *
  9. * *
  10. * D{termination des param}tres de liaison pour la base A. *
  11. * *
  12. * Param}tres: *
  13. * *
  14. * e ITLA Table rassemblant la description des liaisons *
  15. * e PDT Valeur du pas de temps *
  16. * es MTRA Segment de travail temporaire (indique si liaison POL) *
  17. * s NLIAA Nombre total de liaisons. *
  18. * s NXPALA Maxi du nombre de param}tres d{finissant une liaison. *
  19. * s NPLAA Maxi du nombre de points intervenant dans une liaison. *
  20. * s NIPALA Maxi du nombre de param}tres d{finissant une liaison. *
  21. * *
  22. * Auteur, date de cr{ation: *
  23. * *
  24. * Lionel VIVAN, le 24 aout 1989. *
  25. * *
  26. *--------------------------------------------------------------------*
  27. -INC CCOPTIO
  28. -INC SMTABLE
  29. *
  30. SEGMENT,MTRA
  31. INTEGER IPLA(NTRA)
  32. ENDSEGMENT
  33. *
  34. LOGICAL L0,L1
  35. CHARACTER*8 TYPRET,CHARRE
  36. CHARACTER*40 CMOT
  37. *
  38. * Boucle sur le nombre de liaisons
  39. *
  40. NXPALA = 0
  41. NIPALA = 20
  42. NPLAA = 0
  43. NLIAA = 0
  44. IL = 0
  45. 10 CONTINUE
  46. IL = IL + 1
  47. TYPRET = ' '
  48. CALL ACCTAB(ITLA,'ENTIER',IL,X0,' ',L0,IP0,
  49. & TYPRET,I0,X0,CHARRE,L1,ITLIAI)
  50. IF (IERR.NE.0) RETURN
  51. IF (ITLIAI.NE.0) THEN
  52. NLIAA = NLIAA + 1
  53. CALL ACCTAB(ITLIAI,'MOT',I0,X0,'TYPE_LIAISON',L0,IP0,
  54. & 'MOT',I1,X0,CMOT,L1,IP1)
  55. IF (IERR.NE.0) RETURN
  56. *
  57. * ------ liaison {l{mentaire POINT_PLAN_FLUIDE
  58. *
  59. IF (CMOT(1:17).EQ.'POINT_PLAN_FLUIDE') THEN
  60. IPLA(IL) = 0
  61. KIPALA = 3
  62. KXPALA = 6
  63. KPLAA = 1
  64. NXPALA = MAX(NXPALA,KXPALA)
  65. NIPALA = MAX(NIPALA,KIPALA)
  66. NPLAA = MAX(NPLAA,KPLAA)
  67. *
  68. * ------ liaison {l{mentaire POINT_PLAN
  69. *
  70. ELSE IF (CMOT(1:10).EQ.'POINT_PLAN') THEN
  71. IPLA(IL) = 0
  72. TYPRET = ' '
  73. CALL ACCTAB(ITLIAI,'MOT',I0,X0,'AMORTISSEMENT',L0,IP0,
  74. & TYPRET,I1,XAMO,CHARRE,L1,IP1)
  75. IF (IERR.NE.0) RETURN
  76. IF (TYPRET.EQ.'FLOTTANT') THEN
  77. KXPALA = 3
  78. ELSE
  79. KXPALA = 2
  80. ENDIF
  81. KIPALA = 3
  82. KPLAA = 1
  83. NXPALA = MAX(NXPALA,KXPALA)
  84. NIPALA = MAX(NIPALA,KIPALA)
  85. NPLAA = MAX(NPLAA,KPLAA)
  86. *
  87. * ------ liaison {l{mentaire de couplage en vitesse
  88. *
  89. ELSE IF (CMOT(1:16).EQ.'COUPLAGE_VITESSE') THEN
  90. IPLA(IL) = 0
  91. KXPALA = 1
  92. KPLAA = 2
  93. KIPALA = 3
  94. NXPALA = MAX(NXPALA,KXPALA)
  95. NIPALA = MAX(NIPALA,KIPALA)
  96. NPLAA = MAX(NPLAA,KPLAA)
  97. *
  98. * ------ liaison {l{mentaire de couplage en d{placement
  99. *
  100. ELSE IF (CMOT(1:20).EQ.'COUPLAGE_DEPLACEMENT') THEN
  101. IPLA(IL) = 0
  102. cbp KXPALA = 1
  103. KXPALA = 2
  104. KPLAA = 2
  105. KIPALA = 3
  106. NXPALA = MAX(NXPALA,KXPALA)
  107. NIPALA = MAX(NIPALA,KIPALA)
  108. NPLAA = MAX(NPLAA,KPLAA)
  109. *+*
  110. * ------ liaison elementaire POLYNOMIALE
  111. *
  112. ELSE IF (CMOT(1:11).EQ.'POLYNOMIALE') THEN
  113. IPLA(IL) = 1
  114. CALL DYPOL0(ITLIAI,PDT,NX,NI)
  115. IF (IERR.NE.0) RETURN
  116. KXPALA = NX
  117. KPLAA = NI
  118. KIPALA = 3
  119. NXPALA = MAX(NXPALA,KXPALA)
  120. NIPALA = MAX(NIPALA,KIPALA)
  121. NPLAA = MAX(NPLAA,KPLAA)
  122. *
  123. * ------ liaison {l{mentaire ...
  124. *
  125. * ELSE IF (CMOT(1: ).EQ.' ') THEN
  126. * IPLA(IL) = 0
  127. * KXPALA = ...
  128. * KIPALA = ...
  129. * KPLAA = ...
  130. * NXPALA = MAX(NXPALA,KXPALA)
  131. * NIPALA = MAX(NIPALA,KIPALA)
  132. * NPLAA = MAX(NPLAA,KPLAA)
  133. ELSE
  134. CALL ERREUR(490)
  135. RETURN
  136. ENDIF
  137. GOTO 10
  138. ENDIF
  139. *
  140. END
  141.  
  142.  
  143.  
  144.  

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