Télécharger la plaquette Cast3M


fr en

Télécharger active.eso

Retour à la liste

Numérotation des lignes : 

  1. C ACTIVE    SOURCE    FANDEUR   10/12/14    21:15:23     6812
  2.       SUBROUTINE ACTIVE
  3. C--------------------------------------------------------------------
  4. C              OPERATEUR ACCELERATION DE CONVERGENCE
  5. C
  6. C    SYNTAXE : VOIR NOTICE
  7. C--------------------------------------------------------------------
  8.       IMPLICIT INTEGER(I-N)
  9.       IMPLICIT REAL*8(A-H,O-Z)
  10. C
  11. -INC CCOPTIO
  12. C
  13.       CHARACTER*4 MOCOMP
  14.       CHARACTER*4 LISMO(2)
  15.       DATA LISMO/'GEOM','SECA'/
  16. C
  17.       CALL LIRMOT(LISMO,2,IDET,1)
  18.       IF(IERR.NE.0) RETURN
  19. C
  20.       IF(IDET.EQ.1) THEN
  21.       FLOT=0.98D0
  22.       CALL LIRREE(FLOT,0,IRETOU)
  23.       ENDIF
  24. C
  25. C  CAS DES CHPOINTS   - ACCELERATIONS GEOMETRIQUE OU SECANTE
  26. C
  27.       CALL LIROBJ('CHPOINT',MCHPO1,0,IRETOU)
  28.       IF(IERR.NE.0) RETURN
  29.       IF(IRETOU.EQ.1) THEN
  30.         CALL LIROBJ('CHPOINT',MCHPO2,1,IRETOU)
  31.         IF(IERR.NE.0) RETURN
  32.         CALL LIROBJ('CHPOINT',MCHPO3,1,IRETOU)
  33.         IF(IERR.NE.0) RETURN
  34. C
  35.         IF(IDET.EQ.2) THEN
  36.           CALL LIROBJ('CHPOINT',MCHPO4,1,IRETOU)
  37.           IF(IERR.NE.0) RETURN
  38.         ENDIF
  39. C
  40.         CALL ACTIPO(FLOT,IDET,MCHPO1,MCHPO2,MCHPO3,MCHPO4)
  41. C
  42.         RETURN
  43.       ENDIF
  44. C
  45. C  CAS DES MCHAMLS   - ACCELERATION GEOMETRIQUE UNIQUEMENT
  46. C
  47.       IF(IDET.NE.1) THEN
  48.         CALL ERREUR(21)
  49.         RETURN
  50.       ENDIF
  51. C
  52.       CALL LIROBJ('MCHAML',IPCHE1,1,IRETOU)
  53.       IF(IERR.NE.0)RETURN
  54.       CALL LIROBJ('MCHAML',IPCHE2,1,IRETOU)
  55.       IF(IERR.NE.0) RETURN
  56.       CALL LIROBJ('MCHAML',IPCHE3,1,IRETOU)
  57.       IF(IERR.NE.0)RETURN
  58.       MOCOMP=' '
  59.       CALL LIRCHA(MOCOMP,0,IRETOU)
  60.       IF(IERR.NE.0) RETURN
  61.       CALL ACTICH(FLOT,IPCHE1,IPCHE2,IPCHE3,MOCOMP,IPCHE4)
  62. C
  63.       IF (IERR.NE.0) CALL ECROBJ('MCHAML',IPCHE4)
  64. C
  65.       RETURN
  66.       END
  67.  
  68.  
  69.  

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