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    PV090527  24/09/03    21:15:02     11998          
  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.         CALL ACTOBJ('MMODEL',IPMODL,1)
  48. C
  49.         CALL LIROBJ ('MCHAML  ',IPCAR,1,IRET)
  50.         IF(IRET.EQ.0) RETURN
  51.         CALL ACTOBJ('MCHAML',IPCAR ,1)
  52. C
  53.         CALL LIRREE (TOL,1,IRET)
  54.         IF(IRET.EQ.0) RETURN
  55. C
  56.         CALL LIRMOT(MOVERI,1,iret,0)
  57.         if (iret.eq.1) zveri=.true.
  58. C
  59.         CALL MOCU2D(IPCURZ,IPAXIA,IPMODL,IPCAR,IPEPAX,
  60.      $       IPMOMZ,TOL,mtable,zveri)
  61. C
  62.         if (zveri) Call ecrobj('TABLE', mtable)
  63.         CALL ECROBJ('LISTREEL',IPEPAX)
  64.         CALL ECROBJ('LISTREEL',IPMOMZ)
  65.  
  66.       ELSE
  67. C
  68.         CALL LIROBJ ('LISTREEL',IPCURY,1,IRET)
  69.         IF(IRET.EQ.0) RETURN
  70. C
  71.         CALL LIROBJ ('LISTREEL',IPCURZ,1,IRET)
  72.         IF(IRET.EQ.0) RETURN
  73. C
  74.         CALL LIROBJ ('LISTREEL',IPAXIA,1,IRET0)
  75.         IF(IRET0.EQ.0) RETURN
  76. C
  77.         CALL LIROBJ ('MMODEL  ',IPMODL,1,IRET)
  78.         IF(IRET.EQ.0) RETURN
  79.         CALL ACTOBJ('MMODEL',IPMODL,1)
  80. C
  81.         CALL LIROBJ ('MCHAML  ',IPCAR,1,IRET)
  82.         IF(IRET.EQ.0) RETURN
  83.         CALL ACTOBJ('MCHAML',IPCAR ,1)
  84. C
  85.         CALL LIRREE (TOL,1,IRET)
  86.         IF(IRET.EQ.0) RETURN
  87. C
  88.         CALL LIRMOT(MOVERI,1,iret,0)
  89.         if (iret.eq.1) zveri=.true.
  90.  
  91. C
  92.          CALL MOCU3D(IPCURY,IPCURZ,IPAXIA,IPMODL,IPCAR,IPEPAX,
  93.      $ IPMOMY,IPMOMZ,TOL,mtable,zveri)
  94. C
  95.         if (zveri) Call ecrobj('TABLE', mtable)
  96.         CALL ECROBJ('LISTREEL',IPEPAX)
  97.         CALL ECROBJ('LISTREEL',IPMOMZ)
  98.         CALL ECROBJ('LISTREEL',IPMOMY)
  99.  
  100. C
  101.       ENDIF
  102. C
  103.       RETURN
  104.       END
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  

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