Télécharger la plaquette Cast3M


fr en

Télécharger dyne24.eso

Retour à la liste

Numérotation des lignes : 

dyne24

  1. C DYNE24    SOURCE    CHAT      05/01/12    23:15:54     5004
  2.       SUBROUTINE DYNE24(ITABL,ITLIA,KTRAV,II,NLS,NVAR)
  3.       IMPLICIT INTEGER(I-N)
  4.       IMPLICIT REAL*8(A-H,O-Z)
  5. *--------------------------------------------------------------------*
  6. *                                                                    *
  7. *     Operateur DYNE : algorithme de Fu - de Vogelaere               *
  8. *     ________________________________________________               *
  9. *                                                                    *
  10. *     Dimensionne les variables NLS et NVAR, gestion des demandes de *
  11. *     sorties pour les liaisons.                                     *
  12. *                                                                    *
  13. *     Parametres:                                                    *
  14. *                                                                    *
  15. * e   ITABL   Table definissant les liaisons en sortie.              *
  16. * e   ITLIA   Table definissant toutes les liaisons sur une base.    *
  17. * e   KTRAV   Segment de travail qui va remplir MTRES dans DYNE15    *
  18. * es  II      Compteur                                               *
  19. * s   NLS     Nombre de demandes de liaisons en sortie               *
  20. * s   NVAR    Nombre maxi de variables internes de liaisons          *
  21. *                                                                    *
  22. *                                                                    *
  23. *     Auteur, date de cr{ation:                                      *
  24. *                                                                    *
  25. *     Lionel VIVAN, le 2 octobre 1989.                               *
  26. *                                                                    *
  27. *--------------------------------------------------------------------*
  28.  
  29. -INC PPARAM
  30. -INC CCOPTIO
  31. *
  32.       SEGMENT,MTRAV
  33.          INTEGER ICHRE2(NBLS),MPOLA(NBLS),MNULA(NBLS)
  34.          INTEGER MPOLB(NBLS),MNULB(NBLS)
  35.          INTEGER IVLIAA(NBLS,NTVAR),INLIAA(NBLS,NTVAR)
  36.          INTEGER IVLIAB(NBLS,NTVAR),INLIAB(NBLS,NTVAR)
  37.       ENDSEGMENT
  38. *
  39.       LOGICAL L0,L1,LVAR
  40.       CHARACTER*8  MONMOT,TYPRET,CHARRE
  41.       CHARACTER*40 CMOT
  42. *
  43.       MTRAV = KTRAV
  44.       NLS  = 0
  45.       NVAR = 0
  46. *
  47. *     recherche des liaisons en sortie
  48. *
  49.       INDICE = 0
  50.       IVAR   = 0
  51.  100  CONTINUE
  52.       INDICE = INDICE + 1
  53.       TYPRET = ' '
  54.       CALL ACCTAB(ITLIA,'ENTIER',INDICE,X0,' ',L0,IP0,
  55.      &                  TYPRET,I1,X1,CHARRE,L1,ITTL)
  56.       IF (TYPRET.EQ.'TABLE   ' .AND. ITTL.NE.0) THEN
  57.          TYPRET = ' '
  58.          CALL ACCTAB(ITABL,'TABLE',I0,X0,' ',L0,ITTL,
  59.      &                     TYPRET,I1,X1,CHARRE,LVAR,ITVAR)
  60.          IF (TYPRET.EQ.'LOGIQUE ' .AND. LVAR) THEN
  61.             NLS = NLS + 1
  62.             MNULA(NLS) = INDICE
  63.             MPOLA(NLS) = ITTL
  64.             CALL ACCTAB(ITTL,'MOT',I0,X0,'TYPE_LIAISON',L0,IP0,
  65.      &                       'MOT',I1,X1,CMOT,L1,IP1)
  66.             IF (CMOT(1:17).EQ.'POINT_PLAN_FLUIDE') THEN
  67.                IVAR = 6
  68.                NVAR = MAX(NVAR,IVAR)
  69.                II = II + 1
  70.                ICHRE2(II) = IVAR
  71.                IVLIAA(NLS,1) = 1
  72.                INLIAA(NLS,1) = 7
  73.                IVLIAA(NLS,2) = 1
  74.                INLIAA(NLS,2) = 8
  75.                IVLIAA(NLS,3) = 1
  76.                INLIAA(NLS,3) = 69
  77.                IVLIAA(NLS,4) = 1
  78.                INLIAA(NLS,4) = 70
  79.                IVLIAA(NLS,5) = 1
  80.                INLIAA(NLS,5) = 71
  81.                IVLIAA(NLS,6) = 1
  82.                INLIAA(NLS,6) = 72
  83. *
  84.             ELSE IF (CMOT(1:10).EQ.'POINT_PLAN') THEN
  85.                MONMOT = ' '
  86.                CALL ACCTAB(ITTL,'MOT',I0,X0,'AMORTISSEMENT',L0,IP0,
  87.      &                          MONMOT,I1,XAMO,CHARRE,L1,IP1)
  88. *
  89.                IF (MONMOT.EQ.'FLOTTANT') THEN
  90.                   IVAR = 3
  91.                   IVLIAA(NLS,3) = 1
  92.                   INLIAA(NLS,3) = 39
  93.                ELSE
  94.                   IVAR = 2
  95.                ENDIF
  96.                NVAR = MAX(NVAR,IVAR)
  97.                II = II + 1
  98.                ICHRE2(II) = IVAR
  99.                IVLIAA(NLS,1) = 1
  100.                INLIAA(NLS,1) = 15
  101.                IVLIAA(NLS,4) = 1
  102.                INLIAA(NLS,4) = 7
  103. *
  104.             ELSE IF (CMOT(1:16).EQ.'COUPLAGE_VITESSE') THEN
  105.                IVAR = 3
  106.                NVAR = MAX(NVAR,IVAR)
  107.                II = II + 1
  108.                ICHRE2(II) = IVAR
  109.                IVLIAA(NLS,1) = 1
  110.                INLIAA(NLS,1) = 73
  111.                IVLIAA(NLS,3) = 1
  112.                INLIAA(NLS,3) = 8
  113.                IVLIAA(NLS,4) = 1
  114.                INLIAA(NLS,4) = 7
  115. *
  116.             ELSE IF (CMOT(1:20).EQ.'COUPLAGE_DEPLACEMENT') THEN
  117.                IVAR = 2
  118.                NVAR = MAX(NVAR,IVAR)
  119.                II = II + 1
  120.                ICHRE2(II) = IVAR
  121.                IVLIAA(NLS,1) = 1
  122.                INLIAA(NLS,1) = 74
  123.                IVLIAA(NLS,4) = 1
  124.                INLIAA(NLS,4) = 7
  125. *
  126.             ELSE IF (CMOT(1:11).EQ.'POLYNOMIALE') THEN
  127.                IVAR = 1
  128.                NVAR = MAX(NVAR,IVAR)
  129.                II = II + 1
  130.                ICHRE2(II) = IVAR
  131.                IVLIAA(NLS,1) = 1
  132.                INLIAA(NLS,1) = 77
  133. *
  134.             ELSE
  135.                CALL ERREUR(490)
  136.                RETURN
  137.             ENDIF
  138. *
  139.          ELSE IF (TYPRET.EQ.'TABLE   ' .AND. ITVAR.NE.0) THEN
  140.             NLS = NLS + 1
  141.             MNULA(NLS) = INDICE
  142.             MPOLA(NLS) = ITTL
  143.             CALL ACCTAB(ITTL,'MOT',I0,X0,'TYPE_LIAISON',L0,IP0,
  144.      &                       'MOT',I1,X1,CMOT,L1,IT1)
  145.             IVAR = 0
  146. *
  147.             IF (CMOT(1:17).EQ.'POINT_PLAN_FLUIDE') THEN
  148.                TYPRET = ' '
  149.                CALL ACCTAB(ITVAR,'MOT',I0,X0,'DEPLACEMENT',L0,IP0,
  150.      &                           TYPRET,I1,X1,CHARRE,LVAR,IT1)
  151.                IF (TYPRET.EQ.'LOGIQUE ' .AND. LVAR ) THEN
  152.                   IVAR = IVAR + 1
  153.                   IVLIAA(NLS,1) = 1
  154.                   INLIAA(NLS,1) = 7
  155.                ENDIF
  156.                TYPRET = ' '
  157.                CALL ACCTAB(ITVAR,'MOT',I0,X0,'VITESSE',L0,IP0,
  158.      &                           TYPRET,I1,X1,CHARRE,LVAR,IT1)
  159.                IF (TYPRET.EQ.'LOGIQUE ' .AND. LVAR ) THEN
  160.                   IVAR = IVAR + 1
  161.                   IVLIAA(NLS,2) = 1
  162.                   INLIAA(NLS,2) = 8
  163.                ENDIF
  164.                TYPRET = ' '
  165.                CALL ACCTAB(ITVAR,'MOT',I0,X0,'MASSE_AJOUTEE',L0,IP0,
  166.      &                           TYPRET,I1,X1,CHARRE,LVAR,IT1)
  167.                IF (TYPRET.EQ.'LOGIQUE ' .AND. LVAR ) THEN
  168.                   IVAR = IVAR + 1
  169.                   IVLIAA(NLS,3) = 1
  170.                   INLIAA(NLS,3) = 69
  171.                ENDIF
  172.                TYPRET = ' '
  173.                CALL ACCTAB(ITVAR,'MOT',I0,X0,'FORCE_CONVECTION',L0,IP0,
  174.      &                           TYPRET,I1,X1,CHARRE,LVAR,IT1)
  175.                IF (TYPRET.EQ.'LOGIQUE ' .AND. LVAR ) THEN
  176.                   IVAR = IVAR + 1
  177.                   IVLIAA(NLS,4) = 1
  178.                   INLIAA(NLS,4) = 70
  179.                ENDIF
  180.                TYPRET = ' '
  181.                CALL ACCTAB(ITVAR,'MOT',I0,X0,'FORCE_VISCOSITE',L0,IP0,
  182.      &                           TYPRET,I1,X1,CHARRE,LVAR,IT1)
  183.                IF (TYPRET.EQ.'LOGIQUE ' .AND. LVAR ) THEN
  184.                   IVAR = IVAR + 1
  185.                   IVLIAA(NLS,5) = 1
  186.                   INLIAA(NLS,5) = 71
  187.                ENDIF
  188.                TYPRET = ' '
  189.                CALL ACCTAB(ITVAR,'MOT',I0,X0,'FORCE_PERTE_DE_CHARGE',
  190.      &                           L0,IP0,TYPRET,I1,X1,CHARRE,LVAR,IT1)
  191.                IF (TYPRET.EQ.'LOGIQUE ' .AND. LVAR ) THEN
  192.                   IVAR = IVAR + 1
  193.                   IVLIAA(NLS,6) = 1
  194.                   INLIAA(NLS,6) = 72
  195.                ENDIF
  196.                IF (IVAR.EQ.0) THEN
  197.                   CALL ERREUR(496)
  198.                   RETURN
  199.                ENDIF
  200. *
  201.             ELSE IF (CMOT(1:10).EQ.'POINT_PLAN') THEN
  202.                TYPRET = ' '
  203.                CALL ACCTAB(ITVAR,'MOT',I0,X0,'DEPLACEMENT',L0,IP0,
  204.      &                           TYPRET,I1,X1,CHARRE,LVAR,IT1)
  205.                IF (TYPRET.EQ.'LOGIQUE ' .AND. LVAR ) THEN
  206.                   IVAR = IVAR + 1
  207.                   IVLIAA(NLS,4) = 1
  208.                   INLIAA(NLS,4) = 7
  209.                ENDIF
  210.                TYPRET = ' '
  211.                CALL ACCTAB(ITVAR,'MOT',I0,X0,'VITESSE_NORMALE',L0,IP0,
  212.      &                           TYPRET,I1,X1,CHARRE,LVAR,IT1)
  213.                IF (TYPRET.EQ.'LOGIQUE ' .AND. LVAR ) THEN
  214.                   IVAR = IVAR + 1
  215.                   IVLIAA(NLS,3) = 1
  216.                   INLIAA(NLS,3) = 39
  217.                ENDIF
  218.                TYPRET = ' '
  219.                CALL ACCTAB(ITVAR,'MOT',I0,X0,'FORCE_DE_CHOC',L0,IP0,
  220.      &                           TYPRET,I1,X1,CHARRE,LVAR,IT1)
  221.                IF (TYPRET.EQ.'LOGIQUE ' .AND. LVAR ) THEN
  222.                   IVAR = IVAR + 1
  223.                   IVLIAA(NLS,1) = 1
  224.                   INLIAA(NLS,1) = 15
  225.                ENDIF
  226.                IF (IVAR.EQ.0) THEN
  227.                   CALL ERREUR(496)
  228.                   RETURN
  229.                ENDIF
  230. *
  231.             ELSE IF (CMOT(1:16).EQ.'COUPLAGE_VITESSE') THEN
  232.                TYPRET = ' '
  233.                CALL ACCTAB(ITVAR,'MOT',I0,X0,'DEPLACEMENT',L0,IP0,
  234.      &                           TYPRET,I1,X1,CHARRE,LVAR,IT1)
  235.                IF (TYPRET.EQ.'LOGIQUE ' .AND. LVAR ) THEN
  236.                   IVAR = IVAR + 1
  237.                   IVLIAA(NLS,4) = 1
  238.                   INLIAA(NLS,4) = 7
  239.                ENDIF
  240.                TYPRET = ' '
  241.                CALL ACCTAB(ITVAR,'MOT',I0,X0,'VITESSE',L0,IP0,
  242.      &                           TYPRET,I1,X1,CHARRE,LVAR,IT1)
  243.                IF (TYPRET.EQ.'LOGIQUE ' .AND. LVAR ) THEN
  244.                   IVAR = IVAR + 1
  245.                   IVLIAA(NLS,3) = 1
  246.                   INLIAA(NLS,3) = 8
  247.                ENDIF
  248.                TYPRET = ' '
  249.                CALL ACCTAB(ITVAR,'MOT',I0,X0,
  250.      &                     'FORCE_DE_COUPLAGE_VITESSE',L0,IP0,
  251.      &                     TYPRET,I1,X1,CHARRE,LVAR,IT1)
  252.                IF (TYPRET.EQ.'LOGIQUE ' .AND. LVAR ) THEN
  253.                   IVAR = IVAR + 1
  254.                   IVLIAA(NLS,1) = 1
  255.                   INLIAA(NLS,1) = 73
  256.                ENDIF
  257.                IF (IVAR.EQ.0) THEN
  258.                   CALL ERREUR(496)
  259.                   RETURN
  260.                ENDIF
  261. *
  262.             ELSE IF (CMOT(1:20).EQ.'COUPLAGE_DEPLACEMENT') THEN
  263.                TYPRET = ' '
  264.                CALL ACCTAB(ITVAR,'MOT',I0,X0,'DEPLACEMENT',L0,IP0,
  265.      &                           TYPRET,I1,X1,CHARRE,LVAR,IT1)
  266.                IF (TYPRET.EQ.'LOGIQUE ' .AND. LVAR ) THEN
  267.                   IVAR = IVAR + 1
  268.                   IVLIAA(NLS,4) = 1
  269.                   INLIAA(NLS,4) = 7
  270.                ENDIF
  271.                TYPRET = ' '
  272.                CALL ACCTAB(ITVAR,'MOT',I0,X0,
  273.      &                     'FORCE_DE_COUPLAGE_DEPLACEMENT',L0,IP0,
  274.      &                     TYPRET,I1,X1,CHARRE,LVAR,IT1)
  275.                IF (TYPRET.EQ.'LOGIQUE ' .AND. LVAR ) THEN
  276.                   IVAR = IVAR + 1
  277.                   IVLIAA(NLS,1) = 1
  278.                   INLIAA(NLS,1) = 74
  279.                ENDIF
  280.                IF (IVAR.EQ.0) THEN
  281.                   CALL ERREUR(496)
  282.                   RETURN
  283.                ENDIF
  284. *
  285.             ELSE IF (CMOT(1:11).EQ.'POLYNOMIALE') THEN
  286.                TYPRET = ' '
  287.                CALL ACCTAB(ITVAR,'MOT',I0,X0,'FORCE_POLYNOMIALE',
  288.      &                     L0,IP0,TYPRET,I1,X1,CHARRE,LVAR,IT1)
  289.                IF (TYPRET.EQ.'LOGIQUE '.AND.LVAR) THEN
  290.                   IVAR = IVAR + 1
  291.                   IVLIAA(NLS,1) = 1
  292.                   INLIAA(NLS,1) = 77
  293.                ENDIF
  294. *
  295.             ELSE
  296.                CALL ERREUR(490)
  297.                RETURN
  298.             ENDIF
  299.             NVAR = MAX(NVAR,IVAR)
  300.             II = II + 1
  301.             ICHRE2(II) = IVAR
  302.             IF (IIMPI.EQ.333) THEN
  303.                WRITE(IOIMP,*)'DYNE24 : ICHRE2(',II,')=',IVAR
  304.             ENDIF
  305.          ENDIF
  306.          GOTO 100
  307.       ENDIF
  308.       IF (IIMPI.EQ.333) THEN
  309.       WRITE(IOIMP,*) 'DYNE24 : NVAR = ',NVAR
  310.       WRITE(IOIMP,*) 'DYNE24 : NLS  = ',NLS
  311.       WRITE(IOIMP,*) 'DYNE24 : IVAR = ',IVAR
  312.       ENDIF
  313. *
  314.       END
  315.  
  316.  
  317.  
  318.  

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