Télécharger la plaquette Cast3M


fr en

Télécharger mocur.eso

Retour à la liste

Numérotation des lignes : 

mocur

  1. C MOCUR     SOURCE    CB215821  24/04/12    21:16:45     11897          
  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.  
  17. -INC PPARAM
  18. -INC CCOPTIO
  19. -INC SMLREEL
  20. -INC CCGEOME
  21. -INC CCREEL
  22. -INC SMMODEL
  23. -INC SMCHAML
  24.       POINTEUR IPCURY.MLREEL,IPCURZ.MLREEL,IPAXIA.MLREEL
  25.       POINTEUR IPMOMY.MLREEL,IPMOMZ.MLREEL,IPEPAX.MLREEL
  26.       CHARACTER*5 MOVERI(1)
  27.       DATA MOVERI/'VERIF'/
  28. C
  29.       NMAXIT=50
  30.       NSTRS2 = 6
  31.       zveri=.false.
  32. C
  33.       ICOUL1=IDCOUL
  34. C
  35. C
  36. C------- READ THE INPUT --
  37.       IF ((IFOUR.EQ.-2).OR.(IFOUR.EQ.-1)) THEN
  38. C
  39.         CALL LIROBJ ('LISTREEL',IPCURZ,1,IRET)
  40.         IF(IRET.EQ.0) RETURN
  41. C
  42.         CALL LIROBJ ('LISTREEL',IPAXIA,1,IRET0)
  43.         IF(IRET.EQ.0) RETURN
  44. C
  45.         CALL LIROBJ ('MMODEL  ',IPMODL,1,IRET)
  46.         IF(IRET.EQ.0) RETURN
  47. C
  48.         CALL LIROBJ ('MCHAML  ',IPCAR,1,IRET)
  49.         IF(IRET.EQ.0) RETURN
  50. C
  51.         CALL LIRREE (TOL,1,IRET)
  52.         IF(IRET.EQ.0) RETURN
  53. C
  54.         CALL LIRMOT(MOVERI,1,iret,0)
  55.         if (iret.eq.1) zveri=.true.
  56. C
  57.         CALL MOCU2D(IPCURZ,IPAXIA,IPMODL,IPCAR,IPEPAX,
  58.      $       IPMOMZ,TOL,mtable,zveri)
  59. C
  60.         if (zveri) Call ecrobj('TABLE', mtable)
  61.         CALL ECROBJ('LISTREEL',IPEPAX)
  62.         CALL ECROBJ('LISTREEL',IPMOMZ)
  63.  
  64.       ELSE
  65. C
  66.         CALL LIROBJ ('LISTREEL',IPCURY,1,IRET)
  67.         IF(IRET.EQ.0) RETURN
  68. C
  69.         CALL LIROBJ ('LISTREEL',IPCURZ,1,IRET)
  70.         IF(IRET.EQ.0) RETURN
  71. C
  72.         CALL LIROBJ ('LISTREEL',IPAXIA,1,IRET0)
  73.         IF(IRET0.EQ.0) RETURN
  74. C
  75.         CALL LIROBJ ('MMODEL  ',IPMODL,1,IRET)
  76.         IF(IRET.EQ.0) RETURN
  77. C
  78.         CALL LIROBJ ('MCHAML  ',IPCAR,1,IRET)
  79.         IF(IRET.EQ.0) RETURN
  80. C
  81.         CALL LIRREE (TOL,1,IRET)
  82.         IF(IRET.EQ.0) RETURN
  83. C
  84.         CALL LIRMOT(MOVERI,1,iret,0)
  85.         if (iret.eq.1) zveri=.true.
  86.  
  87. C
  88.          CALL MOCU3D(IPCURY,IPCURZ,IPAXIA,IPMODL,IPCAR,IPEPAX,
  89.      $ IPMOMY,IPMOMZ,TOL,mtable,zveri)
  90. C
  91.         if (zveri) Call ecrobj('TABLE', mtable)
  92.         CALL ECROBJ('LISTREEL',IPEPAX)
  93.         CALL ECROBJ('LISTREEL',IPMOMZ)
  94.         CALL ECROBJ('LISTREEL',IPMOMY)
  95.  
  96. C
  97.       ENDIF
  98. C
  99.       RETURN
  100.       END
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  

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