Télécharger la plaquette Cast3M


fr en

Télécharger atiso.eso

Retour à la liste

Numérotation des lignes : 

atiso

  1. C ATISO     SOURCE    CB215821  20/11/25    13:18:30     10792          
  2. C
  3.       SUBROUTINE ATISO(MELEME,ICPR,XPROJ,VCPCHA,VCHC,IVU,PTI,NISO,MCOUP,
  4.      >   mcham,BLOK)
  5. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  6. C                                                                     C
  7. C   APPEL DU TRACE DES ISOVALEURS D UN OBJET DE TYPE CHAMPOINT        C
  8. C     ==>EN DIMENSION 2 ET 3                                          C
  9. C     ==>AVEC VERIFICATIOM SI ICACHE=1 SI LES FACES SONT VISIBLES OU  C
  10. C        NON                                                          C
  11. C                                                                     C
  12. C     AOUT 85                                                         C
  13. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  14. C
  15.       IMPLICIT INTEGER(I-N)
  16.       SEGMENT ICPR(0)
  17.       SEGMENT VCPCHA(nbpts)
  18.       SEGMENT XPROJ(3,ICPR(/1))
  19.       SEGMENT IVU(ICPR(/1))
  20.       REAL BLOK
  21.       REAL VCHC
  22.       DIMENSION VCHC(*)
  23.       DIMENSION XX(4),YY(4),ZZ(4),VV(4)
  24.  
  25. -INC PPARAM
  26. -INC CCOPTIO
  27. -INC CCGEOME
  28. -INC SMELEME
  29. -INC SMCHPOI
  30. -INC SMCOORD
  31. -INC SMLREEL
  32. -INC SMCHAML
  33. C
  34. C  ON FAIT LE TRACER DES ISOVALEURS D UN CHAMPOINT ou d'un chamelem
  35. C  SI ON EST EN OPTION ISOV SURF ON VA APPELE FACED2
  36. C
  37. C SG Pourquoi separe-t-on le cas ISOTYP=0 et les autres ?
  38. C     Cela semble etre essentiellement du au fait que tciso ne gere pas
  39. C     bien ISOTYP=0. Egalement, il y a le trace des labels au bord
  40. C     gestion des enveloppes egalement ?
  41. *dbg      Write(ioimp,*) 'coucou atiso'
  42.       IF (ISOTYP.GT.0) THEN
  43.          CALL FACED2(MELEME,XPROJ,ICPR,VCHC,VCPCHA,PTI,NISO,IVU,MCOUP
  44.      $        ,mcham,BLOK)
  45.          RETURN
  46.       ENDIF
  47.       SEGACT ICPR
  48.       NELEM=0
  49.       IPT1=MELEME
  50.       SEGACT MELEME
  51.       DO 3101 IOB=1,MAX(1,LISOUS(/1))
  52.          IF (LISOUS(/1).NE.0) THEN
  53.           IPT1=LISOUS(IOB)
  54.           SEGACT IPT1
  55.          ENDIF
  56.          NELEM=NELEM+IPT1.NUM(/2)
  57. 3101  CONTINUE
  58. C
  59. C  ON TRAVAILLE PAR OBJET ELEMENTAIRE
  60. C
  61.       IPT1=MELEME
  62.       DO 20100 IOB=1,MAX(1,LISOUS(/1))
  63.          IF (LISOUS(/1).NE.0) IPT1=LISOUS(IOB)
  64.          if (mcham.ne.0) melval=lisref(iob)
  65.          IF (KSURF(IPT1.ITYPEL).EQ.0) GOTO 20150
  66.          NBELEM=IPT1.NUM(/2)
  67.          NBFAC=LTEL(1,IPT1.ITYPEL)
  68.          IAD=LTEL(2,IPT1.ITYPEL)-1
  69.          IF (NBFAC.EQ.0) GOTO 20150
  70.          DO 20200 IFAC=1,NBFAC
  71.             ITYP=LDEL(1,IAD+IFAC)
  72.             JAD=LDEL(2,IAD+IFAC)-1
  73.             IDEP=KDFAC(2,ITYP)
  74.             IFEP=IDEP+3*(KDFAC(3,ITYP)-1)
  75.             IF (ITYP.EQ.1.OR.ITYP.EQ.3.OR.ITYP.EQ.7.OR.ITYP.EQ.8) THEN
  76.                DO 20225 ITRIAN=IDEP,IFEP,3
  77.                   IAFA=LFAC(JAD+KFAC(ITRIAN))
  78.                   IBFA=LFAC(JAD+KFAC(ITRIAN+1))
  79.                   ICFA=LFAC(JAD+KFAC(ITRIAN+2))
  80.                   if(mcham.ne.0) then
  81.                      iafam=min(iafa,velche(/1))
  82.                      ibfam=min(ibfa,velche(/1))
  83.                      icfam=min(icfa,velche(/1))
  84.                   endif
  85.                   DO 20226 IEL=1,NBELEM
  86.                      ielm=iel
  87.                      if(mcham.ne.0) then
  88.                         ielm=min(iel,velche(/2))
  89.                      endif
  90.                      IA=ICPR(IPT1.NUM(IAFA,IEL))
  91.                      IB=ICPR(IPT1.NUM(IBFA,IEL))
  92.                      IC=ICPR(IPT1.NUM(ICFA,IEL))
  93.                      IF (IVU(IA).EQ.1.AND.IVU(IB).EQ.1.AND.IVU(IC).EQ.1)
  94.      $                    THEN
  95.                         XX(1)=XPROJ(1,IA)
  96.                         XX(2)=XPROJ(1,IB)
  97.                         XX(3)=XPROJ(1,IC)
  98.                         YY(1)=XPROJ(2,IA)
  99.                         YY(2)=XPROJ(2,IB)
  100.                         YY(3)=XPROJ(2,IC)
  101.                         ZZ(1)=XPROJ(3,IA)
  102.                         ZZ(2)=XPROJ(3,IB)
  103.                         ZZ(3)=XPROJ(3,IC)
  104.                         if (mcham.eq.0) then
  105.                            VV(1)=VCPCHA(IPT1.NUM(IAFA,IEL))
  106.                            VV(2)=VCPCHA(IPT1.NUM(IBFA,IEL))
  107.                            VV(3)=VCPCHA(IPT1.NUM(ICFA,IEL))
  108.                         else
  109.                            VV(1)=velche(IAFAm,IELm)
  110.                            VV(2)=velche(IBFAm,IELm)
  111.                            VV(3)=velche(ICFAm,IELm)
  112.                         endif
  113.                         CALL TRISO(VCHC,XX,YY,ZZ,VV,3,NISO)
  114.                      ENDIF
  115. 20226             continue
  116. 20225          CONTINUE
  117.             ELSEIF (ITYP.EQ.2) THEN
  118. *  cas des quadrangles on rajoute un point au milieu
  119.                DO 20235 IEL=1,NBELEM
  120.                   IAFA=LFAC(JAD+1)
  121.                   IBFA=LFAC(JAD+2)
  122.                   ICFA=LFAC(JAD+3)
  123.                   IDFA=LFAC(JAD+4)
  124.                   if(mcham.ne.0) then
  125.                      ielm=min(iel,velche(/2))
  126.                      iafam=min(iafa,velche(/1))
  127.                      ibfam=min(ibfa,velche(/1))
  128.                      icfam=min(icfa,velche(/1))
  129.                      idfam=min(idfa,velche(/1))
  130.                   endif
  131.                   ia=icpr(ipt1.num(iafa,iel))
  132.                   ib=icpr(ipt1.num(ibfa,iel))
  133.                   ic=icpr(ipt1.num(icfa,iel))
  134.                   id=icpr(ipt1.num(idfa,iel))
  135.                   IF (IVU(IA).EQ.1.AND.IVU(IB).EQ.1.AND.IVU(IC).EQ.1
  136.      $                 .and.IVU(Id).EQ.1) THEN
  137.                      xxm=(XPROJ(1,IA)+XPROJ(1,IB)
  138.      $                   +XPROJ(1,IC)+XPROJ(1,ID))/4
  139.                      yym=(XPROJ(2,IA)+XPROJ(2,IB)
  140.      $                   +XPROJ(2,IC)+XPROJ(2,ID))/4
  141.                      zzm=(XPROJ(3,IA)+XPROJ(3,IB)
  142.      $                   +XPROJ(3,IC)+XPROJ(3,ID))/4
  143.                      if (mcham.eq.0) then
  144.                         vvm=(VCPCHA(ipt1.num(iafa,iel))
  145.      $                      +VCPCHA(ipt1.num(ibfa,iel))
  146.      $                      +VCPCHA(ipt1.num(icfa,iel))
  147.      $                      +VCPCHA(ipt1.num(idfa,iel)))/4
  148.                      else
  149.                         vvm=(velche(iafam,ielm)+velche(ibfam,ielm)
  150.      *                      +velche(icfam,ielm)+velche(idfam,ielm))/4
  151.                      endif
  152.                      XX(1)=XPROJ(1,IA)
  153.                      XX(2)=XPROJ(1,IB)
  154.                      YY(1)=XPROJ(2,IA)
  155.                      YY(2)=XPROJ(2,IB)
  156.                      ZZ(1)=XPROJ(3,IA)
  157.                      ZZ(2)=XPROJ(3,IB)
  158.                      if (mcham.eq.0) then
  159.                         VV(1)=VCPCHA(IPT1.NUM(IAFA,IEL))
  160.                         VV(2)=VCPCHA(IPT1.NUM(IBFA,IEL))
  161.                      else
  162.                         VV(1)=velche(IAFAm,IELm)
  163.                         VV(2)=velche(IBFAm,IELm)
  164.                      endif
  165.                      XX(3)=XXM
  166.                      YY(3)=YYM
  167.                      ZZ(3)=ZZM
  168.                      VV(3)=VVM
  169.                      CALL TRISO(VCHC,XX,YY,ZZ,VV,3,NISO)
  170.                      XX(2)=XPROJ(1,ID)
  171.                      YY(2)=XPROJ(2,ID)
  172.                      ZZ(2)=XPROJ(3,ID)
  173.                      if (mcham.eq.0) then
  174.                         VV(2)=VCPCHA(IPT1.NUM(IDFA,IEL))
  175.                      else
  176.                         VV(2)=velche(IDFAm,IELm)
  177.                      endif
  178.                      CALL TRISO(VCHC,XX,YY,ZZ,VV,3,NISO)
  179.                      XX(1)=XPROJ(1,IC)
  180.                      YY(1)=XPROJ(2,IC)
  181.                      ZZ(1)=XPROJ(3,IC)
  182.                      if (mcham.eq.0) then
  183.                         VV(1)=VCPCHA(IPT1.NUM(ICFA,IEL))
  184.                      else
  185.                         VV(1)=velche(ICFAm,IELm)
  186.                      endif
  187.                      CALL TRISO(VCHC,XX,YY,ZZ,VV,3,NISO)
  188.                      XX(2)=XPROJ(1,IB)
  189.                      YY(2)=XPROJ(2,IB)
  190.                      ZZ(2)=XPROJ(3,IB)
  191.                      if (mcham.eq.0) then
  192.                         VV(2)=VCPCHA(IPT1.NUM(IBFA,IEL))
  193.                      else
  194.                         VV(2)=velche(IBFAm,IELm)
  195.                      endif
  196.                      CALL TRISO(VCHC,XX,YY,ZZ,VV,3,NISO)
  197.                   ENDIF
  198. 20235          CONTINUE
  199.             ELSEIF (ITYP.EQ.4) THEN
  200. *  cas des quadrangles a 8 noeuds on decoupe en 4 triangles et un
  201. *  quadrangle.
  202.                DO 20237 IEL=1,NBELEM
  203.                   IAFA=LFAC(JAD+1)
  204.                   IBFA=LFAC(JAD+2)
  205.                   ICFA=LFAC(JAD+3)
  206.                   IDFA=LFAC(JAD+4)
  207.                   IEFA=LFAC(JAD+5)
  208.                   IFFA=LFAC(JAD+6)
  209.                   IGFA=LFAC(JAD+7)
  210.                   IHFA=LFAC(JAD+8)
  211.                   if(mcham.ne.0) then
  212.                      ielm=min(iel,velche(/2))
  213.                      iafam=min(iafa,velche(/1))
  214.                      ibfam=min(ibfa,velche(/1))
  215.                      icfam=min(icfa,velche(/1))
  216.                      idfam=min(idfa,velche(/1))
  217.                      iefam=min(iefa,velche(/1))
  218.                      iffam=min(iffa,velche(/1))
  219.                      igfam=min(igfa,velche(/1))
  220.                      ihfam=min(ihfa,velche(/1))
  221.                   endif
  222.                   ia=icpr(ipt1.num(iafa,iel))
  223.                   ib=icpr(ipt1.num(ibfa,iel))
  224.                   ic=icpr(ipt1.num(icfa,iel))
  225.                   id=icpr(ipt1.num(idfa,iel))
  226.                   ie=icpr(ipt1.num(iefa,iel))
  227.                   if=icpr(ipt1.num(iffa,iel))
  228.                   ig=icpr(ipt1.num(igfa,iel))
  229.                   ih=icpr(ipt1.num(ihfa,iel))
  230.                   IF (IVU(IA).EQ.1.AND.IVU(IB).EQ.1.AND.IVU(IH).EQ.1)
  231.      $                 then
  232.                      XX(1)=XPROJ(1,IA)
  233.                      XX(2)=XPROJ(1,IB)
  234.                      XX(3)=XPROJ(1,IH)
  235.                      YY(1)=XPROJ(2,IA)
  236.                      YY(2)=XPROJ(2,IB)
  237.                      YY(3)=XPROJ(2,IH)
  238.                      ZZ(1)=XPROJ(3,IA)
  239.                      ZZ(2)=XPROJ(3,IB)
  240.                      ZZ(3)=XPROJ(3,IH)
  241.                      if (mcham.eq.0) then
  242.                         VV(1)=VCPCHA(IPT1.NUM(IAFA,IEL))
  243.                         VV(2)=VCPCHA(IPT1.NUM(IBFA,IEL))
  244.                         VV(3)=VCPCHA(IPT1.NUM(IHFA,IEL))
  245.                      else
  246.                         VV(1)=velche(IAFAm,IELm)
  247.                         VV(2)=velche(IBFAm,IELm)
  248.                         VV(3)=velche(IHFAm,IELm)
  249.                      endif
  250.                      CALL TRISO(VCHC,XX,YY,ZZ,VV,3,NISO)
  251.                   ENDIF
  252.                   IF (IVU(IB).EQ.1.AND.IVU(IC).EQ.1.AND.IVU(ID).EQ.1)
  253.      $                 then
  254.                      XX(1)=XPROJ(1,IB)
  255.                      XX(2)=XPROJ(1,IC)
  256.                      XX(3)=XPROJ(1,ID)
  257.                      YY(1)=XPROJ(2,IB)
  258.                      YY(2)=XPROJ(2,IC)
  259.                      YY(3)=XPROJ(2,ID)
  260.                      ZZ(1)=XPROJ(3,IB)
  261.                      ZZ(2)=XPROJ(3,IC)
  262.                      ZZ(3)=XPROJ(3,ID)
  263.                      if (mcham.eq.0) then
  264.                         VV(1)=VCPCHA(IPT1.NUM(IBFA,IEL))
  265.                         VV(2)=VCPCHA(IPT1.NUM(ICFA,IEL))
  266.                         VV(3)=VCPCHA(IPT1.NUM(IDFA,IEL))
  267.                      else
  268.                         VV(1)=velche(IBFAm,IELm)
  269.                         VV(2)=velche(ICFAm,IELm)
  270.                         VV(3)=velche(IDFAm,IELm)
  271.                      endif
  272.                      CALL TRISO(VCHC,XX,YY,ZZ,VV,3,NISO)
  273.                   ENDIF
  274.                   IF (IVU(ID).EQ.1.AND.IVU(IE).EQ.1.AND.IVU(IF).EQ.1)
  275.      $                 then
  276.                      XX(1)=XPROJ(1,ID)
  277.                      XX(2)=XPROJ(1,IE)
  278.                      XX(3)=XPROJ(1,IF)
  279.                      YY(1)=XPROJ(2,ID)
  280.                      YY(2)=XPROJ(2,IE)
  281.                      YY(3)=XPROJ(2,IF)
  282.                      ZZ(1)=XPROJ(3,ID)
  283.                      ZZ(2)=XPROJ(3,IE)
  284.                      ZZ(3)=XPROJ(3,IF)
  285.                      if (mcham.eq.0) then
  286.                         VV(1)=VCPCHA(IPT1.NUM(IDFA,IEL))
  287.                         VV(2)=VCPCHA(IPT1.NUM(IEFA,IEL))
  288.                         VV(3)=VCPCHA(IPT1.NUM(IFFA,IEL))
  289.                      else
  290.                         VV(1)=velche(IDFAm,IELm)
  291.                         VV(2)=velche(IEFAm,IELm)
  292.                         VV(3)=velche(IFFAm,IELm)
  293.                      endif
  294.                      CALL TRISO(VCHC,XX,YY,ZZ,VV,3,NISO)
  295.                   ENDIF
  296.                   IF (IVU(IF).EQ.1.AND.IVU(IG).EQ.1.AND.IVU(IH).EQ.1)
  297.      $                 then
  298.                      XX(1)=XPROJ(1,IF)
  299.                      XX(2)=XPROJ(1,IG)
  300.                      XX(3)=XPROJ(1,IH)
  301.                      YY(1)=XPROJ(2,IF)
  302.                      YY(2)=XPROJ(2,IG)
  303.                      YY(3)=XPROJ(2,IH)
  304.                      ZZ(1)=XPROJ(3,IF)
  305.                      ZZ(2)=XPROJ(3,IG)
  306.                      ZZ(3)=XPROJ(3,IH)
  307.                      if (mcham.eq.0) then
  308.                         VV(1)=VCPCHA(IPT1.NUM(IFFA,IEL))
  309.                         VV(2)=VCPCHA(IPT1.NUM(IGFA,IEL))
  310.                         VV(3)=VCPCHA(IPT1.NUM(IHFA,IEL))
  311.                      else
  312.                         VV(1)=velche(IFFAm,IELm)
  313.                         VV(2)=velche(IGFAm,IELm)
  314.                         VV(3)=velche(IHFAm,IELm)
  315.                      endif
  316.                      CALL TRISO(VCHC,XX,YY,ZZ,VV,3,NISO)
  317.                   ENDIF
  318. *  cas du quadrangle on rajoute un point central
  319.                   xxm=-0.25*XPROJ(1,IA)+0.5*XPROJ(1,IB)
  320.      *                -0.25*XPROJ(1,IC)+0.5*XPROJ(1,ID)
  321.      *                -0.25*XPROJ(1,IE)+0.5*XPROJ(1,IF)
  322.      *                -0.25*XPROJ(1,IG)+0.5*XPROJ(1,IH)
  323.                   yym=-0.25*XPROJ(2,IA)+0.5*XPROJ(2,IB)
  324.      *                -0.25*XPROJ(2,IC)+0.5*XPROJ(2,ID)
  325.      *                -0.25*XPROJ(2,IE)+0.5*XPROJ(2,IF)
  326.      *                -0.25*XPROJ(2,IG)+0.5*XPROJ(2,IH)
  327.                   zzm=-0.25*XPROJ(3,IA)+0.5*XPROJ(3,IB)
  328.      *                -0.25*XPROJ(3,IC)+0.5*XPROJ(3,ID)
  329.      *                -0.25*XPROJ(3,IE)+0.5*XPROJ(3,IF)
  330.      *                -0.25*XPROJ(3,IG)+0.5*XPROJ(3,IH)
  331.                   if (mcham.eq.0) then
  332.                      vvm=-0.25*VCPCHA(ipt1.num(iafa,iel))
  333.      $                    +0.5*VCPCHA(ipt1.num(ibfa,iel))
  334.      $                   -0.25*VCPCHA(ipt1.num(icfa,iel))
  335.      $                    +0.5*VCPCHA(ipt1.num(idfa,iel))
  336.      $                   -0.25*VCPCHA(ipt1.num(iefa,iel))
  337.      $                    +0.5*VCPCHA(ipt1.num(iffa,iel))
  338.      $                   -0.25*VCPCHA(ipt1.num(igfa,iel))
  339.      $                    +0.5*VCPCHA(ipt1.num(ihfa,iel))
  340.                   else
  341.                      vvm=-0.25*velche(iafam,ielm)
  342.      $                    +0.5*velche(ibfam,ielm)
  343.      *                   -0.25*velche(icfam,ielm)
  344.      $                    +0.5*velche(idfam,ielm)
  345.      $                   -0.25*velche(iefam,ielm)
  346.      $                    +0.5*velche(iffam,ielm)
  347.      $                   -0.25*velche(igfam,ielm)
  348.      $                    +0.5*velche(ihfam,ielm)
  349.                   endif
  350.                   IAFA=LFAC(JAD+2)
  351.                   IBFA=LFAC(JAD+4)
  352.                   ICFA=LFAC(JAD+6)
  353.                   IDFA=LFAC(JAD+8)
  354.                   if(mcham.ne.0) then
  355.                      iafam=min(iafa,velche(/1))
  356.                      ibfam=min(ibfa,velche(/1))
  357.                      icfam=min(icfa,velche(/1))
  358.                      idfam=min(idfa,velche(/1))
  359.                   endif
  360.                   ia=icpr(ipt1.num(iafa,iel))
  361.                   ib=icpr(ipt1.num(ibfa,iel))
  362.                   ic=icpr(ipt1.num(icfa,iel))
  363.                   id=icpr(ipt1.num(idfa,iel))
  364.                   IF (IVU(IA).EQ.1.AND.IVU(IB).EQ.1.AND.IVU(IC).EQ.1
  365.      $                 .and.IVU(Id).EQ.1) THEN
  366.                      XX(1)=XPROJ(1,IA)
  367.                      XX(2)=XPROJ(1,IB)
  368.                      YY(1)=XPROJ(2,IA)
  369.                      YY(2)=XPROJ(2,IB)
  370.                      ZZ(1)=XPROJ(3,IA)
  371.                      ZZ(2)=XPROJ(3,IB)
  372.                      if (mcham.eq.0) then
  373.                         VV(1)=VCPCHA(IPT1.NUM(IAFA,IEL))
  374.                         VV(2)=VCPCHA(IPT1.NUM(IBFA,IEL))
  375.                      else
  376.                         VV(1)=velche(IAFAm,IELm)
  377.                         VV(2)=velche(IBFAm,IELm)
  378.                      endif
  379.                      XX(3)=XXM
  380.                      YY(3)=YYM
  381.                      ZZ(3)=ZZM
  382.                      VV(3)=VVM
  383.                      CALL TRISO(VCHC,XX,YY,ZZ,VV,3,NISO)
  384.                      XX(2)=XPROJ(1,ID)
  385.                      YY(2)=XPROJ(2,ID)
  386.                      ZZ(2)=XPROJ(3,ID)
  387.                      if (mcham.eq.0) then
  388.                         VV(2)=VCPCHA(IPT1.NUM(IDFA,IEL))
  389.                      else
  390.                         VV(2)=velche(IDFAm,IELm)
  391.                      endif
  392.                      CALL TRISO(VCHC,XX,YY,ZZ,VV,3,NISO)
  393.                      XX(1)=XPROJ(1,IC)
  394.                      YY(1)=XPROJ(2,IC)
  395.                      ZZ(1)=XPROJ(3,IC)
  396.                      if (mcham.eq.0) then
  397.                         VV(1)=VCPCHA(IPT1.NUM(ICFA,IEL))
  398.                      else
  399.                         VV(1)=velche(ICFAm,IELm)
  400.                      endif
  401.                      CALL TRISO(VCHC,XX,YY,ZZ,VV,3,NISO)
  402.                      XX(2)=XPROJ(1,IB)
  403.                      YY(2)=XPROJ(2,IB)
  404.                      ZZ(2)=XPROJ(3,IB)
  405.                      if (mcham.eq.0) then
  406.                         VV(2)=VCPCHA(IPT1.NUM(IBFA,IEL))
  407.                      else
  408.                         VV(2)=velche(IBFAm,IELm)
  409.                      endif
  410.                      CALL TRISO(VCHC,XX,YY,ZZ,VV,3,NISO)
  411.                   ENDIF
  412. 20237          CONTINUE
  413.             ELSEIF (ITYP .EQ. 6) THEN
  414. *     
  415. *     Element polygone
  416. *
  417.                NBCOT = IPT1.NUM(/1)
  418.                DO 20300 IEL=1,NBELEM
  419. *
  420. *          Recherche des coordonnees et valeurs au centre du polygone
  421. *
  422.                   XXM = 0.
  423.                   YYM = 0.
  424.                   ZZM = 0.
  425.                   VVM = 0.
  426.                   IVUE = 1
  427.                   DO 20330 ICT = 1, NBCOT
  428. *     
  429.                      NUPT = IPT1.NUM(ICT, IEL)
  430.                      IDPT = ICPR(NUPT)
  431.                      XXM = XPROJ(1,IDPT) + XXM
  432.                      YYM = XPROJ(2,IDPT) + YYM
  433.                      ZZM = XPROJ(3,IDPT) + ZZM
  434.                      IF (MCHAM.EQ.0) THEN
  435.                         VVM = VCPCHA(NUPT) + VVM
  436.                      ELSE
  437.                         VVM = VELCHE(ICT, IEL) + VVM
  438.                      ENDIF
  439.                      IF (IVU(IDPT).NE.1)  IVUE = 0
  440. *     
  441. 20330             CONTINUE
  442. *     
  443.                   IF (IVUE.EQ.1) THEN
  444. *     
  445.                      XX(3)=XXM/NBCOT
  446.                      YY(3)=YYM/NBCOT
  447.                      ZZ(3)=ZZM/NBCOT
  448.                      VV(3)=VVM/NBCOT
  449. *     
  450.                      NUPT2 = IPT1.NUM(NBCOT, IEL)
  451.                      ID2 = ICPR(NUPT2)
  452.                      XX(2)=XPROJ(1,ID2)
  453.                      YY(2)=XPROJ(2,ID2)
  454.                      ZZ(2)=XPROJ(3,ID2)
  455.                      IF (MCHAM.EQ.0) THEN
  456.                         VV(2) = VCPCHA(NUPT2)
  457.                      ELSE
  458. C##   MC               VV(2) = VELCHE(ICT2, IEL)
  459.                         VV(2) = VELCHE(NBCOT, IEL)
  460.                      ENDIF
  461. *
  462. *            Boucle sur tous les triangles
  463. *
  464.                      DO 20360 , ICT = 1, NBCOT
  465. *     
  466.                         NUPT1 = IPT1.NUM(ICT, IEL)
  467.                         ID1 = ICPR(NUPT1)
  468. *     
  469.                         XX(1)=XPROJ(1,ID1)
  470.                         YY(1)=XPROJ(2,ID1)
  471.                         ZZ(1)=XPROJ(3,ID1)
  472.                         IF (MCHAM.EQ.0) THEN
  473.                            VV(1)= VCPCHA(NUPT1)
  474.                         ELSE
  475.                            VV(1) = VELCHE(ICT, IEL)
  476.                         ENDIF
  477.                         CALL TRISO(VCHC,XX,YY,ZZ,VV,3,NISO)
  478.                         XX(2)= XX(1)
  479.                         YY(2)= YY(1)
  480.                         ZZ(2)= ZZ(1)
  481.                         VV(2)= VV(1)
  482. *     
  483. 20360                CONTINUE
  484. *     
  485.                   ENDIF
  486. *     
  487. 20300          CONTINUE
  488.             ENDIF
  489. 20200    CONTINUE
  490. 20150    CONTINUE
  491.          IF (LISOUS(/1).NE.0) SEGDES IPT1
  492. 20100 CONTINUE
  493. *     METTRE LES LABELS SUR LES ISO QUI COUPENT LE CONTOUR
  494.       CALL TRCISO(MELEME,VCPCHA,ICPR,XPROJ,NISO,VCHC,IVU)
  495.       SEGDES MELEME
  496.       SEGDES ICPR
  497.       RETURN
  498.       END
  499.  
  500.  
  501.  
  502.  
  503.  
  504.  
  505.  
  506.  
  507.  
  508.  
  509.  
  510.  
  511.  
  512.  
  513.  
  514.  
  515.  
  516.  
  517.  
  518.  
  519.  
  520.  

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