Télécharger mocur.eso

Retour à la liste

Numérotation des lignes :

  1. C MOCUR SOURCE BP208322 16/11/18 21:19:19 9177
  2. SUBROUTINE MOCUR
  3. **********************************************************************
  4. *
  5. * OPERATEUR MOCU (MOment/CUrvature)
  6. *
  7. **********************************************************************
  8. * Pierre Pegon (ISPRA) Juillet/Aout 1993
  9. **********************************************************************
  10. IMPLICIT INTEGER(I-N)
  11. IMPLICIT REAL*8(A-H,O-Z)
  12. C
  13. DIMENSION DEPSI(6),SIGMA(6),CRIGI(12),CMASS(12)
  14. logical zveri
  15. C
  16. -INC CCOPTIO
  17. -INC SMLREEL
  18. -INC CCGEOME
  19. -INC CCREEL
  20. -INC SMMODEL
  21. -INC SMCHAML
  22. POINTEUR IPCURY.MLREEL,IPCURZ.MLREEL,IPAXIA.MLREEL
  23. POINTEUR IPMOMY.MLREEL,IPMOMZ.MLREEL,IPEPAX.MLREEL
  24. C
  25. NMAXIT=50
  26. NSTRS2 = 6
  27. zveri=.false.
  28. C
  29. ICOUL1=IDCOUL
  30. C
  31. C
  32. C------- READ THE INPUT --
  33. IF ((IFOUR.EQ.-2).OR.(IFOUR.EQ.-1)) THEN
  34. C
  35. CALL LIROBJ ('LISTREEL',IPCURZ,1,IRET)
  36. IF(IRET.EQ.0) RETURN
  37. C
  38. CALL LIROBJ ('LISTREEL',IPAXIA,1,IRET0)
  39. IF(IRET.EQ.0) RETURN
  40. C
  41. CALL LIROBJ ('MMODEL ',IPMODL,1,IRET)
  42. IF(IRET.EQ.0) RETURN
  43. C
  44. CALL LIROBJ ('MCHAML ',IPCAR,1,IRET)
  45. IF(IRET.EQ.0) RETURN
  46. C
  47. CALL LIRREE (TOL,1,IRET)
  48. IF(IRET.EQ.0) RETURN
  49. C
  50. CALL LIRMOT('VERIF',5,iret,0)
  51. if (iret.eq.1) zveri=.true.
  52. C
  53. CALL MOCU2D(IPCURZ,IPAXIA,IPMODL,IPCAR,IPEPAX,
  54. $ IPMOMZ,TOL,mtable,zveri)
  55. C
  56. if (zveri) Call ecrobj('TABLE', mtable)
  57. CALL ECROBJ('LISTREEL',IPEPAX)
  58. CALL ECROBJ('LISTREEL',IPMOMZ)
  59.  
  60. ELSE
  61. C
  62. CALL LIROBJ ('LISTREEL',IPCURY,1,IRET)
  63. IF(IRET.EQ.0) RETURN
  64. C
  65. CALL LIROBJ ('LISTREEL',IPCURZ,1,IRET)
  66. IF(IRET.EQ.0) RETURN
  67. C
  68. CALL LIROBJ ('LISTREEL',IPAXIA,1,IRET0)
  69. IF(IRET0.EQ.0) RETURN
  70. C
  71. CALL LIROBJ ('MMODEL ',IPMODL,1,IRET)
  72. IF(IRET.EQ.0) RETURN
  73. C
  74. CALL LIROBJ ('MCHAML ',IPCAR,1,IRET)
  75. IF(IRET.EQ.0) RETURN
  76. C
  77. CALL LIRREE (TOL,1,IRET)
  78. IF(IRET.EQ.0) RETURN
  79. C
  80. CALL LIRMOT('VERIF',5,iret,0)
  81. if (iret.eq.1) zveri=.true.
  82.  
  83. C
  84. CALL MOCU3D(IPCURY,IPCURZ,IPAXIA,IPMODL,IPCAR,IPEPAX,
  85. $ IPMOMY,IPMOMZ,TOL,mtable,zveri)
  86. C
  87. if (zveri) Call ecrobj('TABLE', mtable)
  88. CALL ECROBJ('LISTREEL',IPEPAX)
  89. CALL ECROBJ('LISTREEL',IPMOMZ)
  90. CALL ECROBJ('LISTREEL',IPMOMY)
  91.  
  92. C
  93. ENDIF
  94. C
  95. RETURN
  96. END
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  

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