Esope source files

   1 : C CFL       SOURCE    CHAT      06/08/24    21:19:19     5529           
   2 :       SUBROUTINE CFL
   3 : *-----------------------------------------------------------------------
   4 : *
   5 : *   chapeau de l'opérateur cfl
   6 : *
   7 : *   appelle la routine clf1.eso qui est aussi appelé par taille et cson
   8 : *
   9 : *-----------------------------------------------------------------------
  10 : *
  11 :       IMPLICIT INTEGER(I-N)
  12 :       IMPLICIT REAL*8(A-H,O-Z)
  13 :       CHARACTER*4 MOTCLE(2)
  14 :       DATA MOTCLE /'CSON','TAIL'/
  15 : *
  16 :       IRET = 0
  17 :       IRET1 = 0
  18 :       IRET2 = 0
  19 :       IRET3 = 0
  20 :       IPCHA1 = 0
  21 :       IPCHA2 = 0
  22 :       IPCHA3 = 0
  23 : *
  23 :       CALL LIROBJ('MMODEL',IPMODL,1,IRET)
  24 :       IF (IRET.EQ.0) RETURN
  25 : c
  26 : *
  27 : *  determination des trois option
  28 : *
  28 :       CALL LIRMOT(MOTCLE,2,IVAL,0)
  29 :       IF (IVAL .EQ. 0 ) THEN
  30 : *      cas du calcul global
  31 :           ICAS = 1
  31 :           CALL LIROBJ('MCHAML',IPCHA1,1,IRET1)
  32 :           IF (IRET1.EQ.0) RETURN
  33 : *
  34 :       ELSE IF (IVAL .EQ. 1 ) THEN
  35 : *       cas ou la vitesse du son est fournie
  36 :           ICAS = 2
  37 : *
  37 :           CALL LIROBJ('MCHAML',IPCHA1,1,IRET1)
  38 :           IF (IRET1.EQ.0) RETURN
  39 : *
  39 :           CALL LIROBJ('MCHAML',IPCHA2,0,IRET2)
  40 : *         si les caractéristiques ne sont pas fournies
  41 : *         la taille l'est
  42 :           IF ( IRET2 .EQ. 0) IPCHA2 = IPCHA1
  43 :       ELSE IF (IVAL .EQ. 2 ) THEN
  44 :           ICAS = 3
  45 : *          champ de caractéristiques
  45 :           CALL LIROBJ('MCHAML',IPCHA1,1,IRET1)
  46 :           IF (IRET1.EQ.0) RETURN
  46 :           CALL LIROBJ('MCHAML',IPCHA3,1,IRET3)
  47 :           IF (IRET3.EQ.0) RETURN
  48 :       ENDIF
  49 : *
  50 : *  ipcha1 champ de caractéristiques
  51 : *  ipcha2 champ de vitesse du son composante 'cson'
  52 : *  ipcha3 champ de taille du maillage composante 'l' ( et 'l2h' facultatif)
  53 : *
  53 :       CALL CFL1(IPMODL,IPCHA1,IPCHA2,IPCHA3,IPCHA4,ICAS)
  54 : *
  55 : *  en retour on récupère le champ par élément de composante 'tcfl'
  56 : *
  57 :       IF ( IPCHA4 .EQ. 0) RETURN
  58 : *
  58 :       CALL ECROBJ('MCHAML',IPCHA4)
  59 :       RETURN
  60 :       END
  61 :  
  62 :  
  63 :