Télécharger dyne25.eso

Retour à la liste

Numérotation des lignes :

  1. C DYNE25 SOURCE CHAT 05/01/12 23:16:03 5004
  2. SUBROUTINE DYNE25(IBAS,KCPR,IK)
  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. * Recherche des points de r{f{rence d{finis par les points *
  11. * rep}res de la base modale. *
  12. * *
  13. * Param}tres: *
  14. * *
  15. * e IBAS Table repr{sentant une base modale *
  16. * s KCPR Segment des points *
  17. * es IK Compteur *
  18. * *
  19. * Auteur, date de cr{ation: *
  20. * *
  21. * Lionel VIVAN, le 24 octobre 1989. *
  22. * *
  23. *--------------------------------------------------------------------*
  24. -INC CCOPTIO
  25. -INC SMCOORD
  26. *
  27. SEGMENT,ICPR(XCOOR(/1)/(IDIM+1))
  28. *
  29. LOGICAL L0,L1
  30. CHARACTER*4 MO2
  31. CHARACTER*8 TYPRET,CHARRE
  32. *
  33. ICPR = KCPR
  34. *
  35. IMODE = 0
  36. 10 CONTINUE
  37. ITBMOD = 0
  38. IMODE = IMODE + 1
  39. TYPRET = ' '
  40. CALL ACCTAB(IBAS,'ENTIER',IMODE,X0,' ',L0,IP0,
  41. & TYPRET,I1,X1,CHARRE,L1,ITBMOD)
  42. IF (IERR.NE.0) RETURN
  43. IF (ITBMOD.NE.0) THEN
  44. CALL ACCTAB(ITBMOD,'MOT',I0,X0,'SOUSTYPE',L0,IP0,
  45. & 'MOT',I1,X1,MO2,L1,IP1)
  46. IF (IERR.NE.0) RETURN
  47. IF (MO2.EQ.'MODE') THEN
  48. *
  49. * Le point rep}re nous {choit
  50. *
  51. CALL ACCTAB(ITBMOD,'MOT',I0,X0,'POINT_REPERE',L0,IP0,
  52. & 'POINT',I1,X1,' ',L1,KNOE)
  53. IF (IERR.NE.0) RETURN
  54. IF (KNOE.NE.0) THEN
  55. IF (ICPR(KNOE).EQ.0) THEN
  56. IK = IK + 1
  57. ICPR(KNOE) = IK
  58. IF (IIMPI.EQ.333) THEN
  59. WRITE(IOIMP,*)'DYNE25 : basemo. ICPR(',KNOE,')=',ICPR(KNOE)
  60. ENDIF
  61. ENDIF
  62. ENDIF
  63. ELSE
  64. CALL ERREUR(314)
  65. RETURN
  66. ENDIF
  67. GOTO 10
  68. ENDIF
  69. *
  70. END
  71.  
  72.  
  73.  

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