Esope source files

   1 : C SUBDEN    SOURCE    CHAT      07/10/29    21:28:48     5940           
   2 : C   LECTURE DE LA DENSITE   ATTENTION ON VA CORRIGER LA VARIABLE &DENS
   3 : C
   4 : C PP 12/98: forcer la densite des point d'un maillage
   5 :       SUBROUTINE SUBDEN
   6 :       IMPLICIT INTEGER(I-N)
   7 : -INC CCGEOME
   8 : -INC CCOPTIO
   9 : C+PP
  10 : -INC SMCOORD
  11 : -INC SMELEME
  12 : C+PP
  13 :       REAL*8 XXX
  13 :       CALL MESLIR(-238)
  13 :       CALL LIRREE(XXX,1,IRETOU)
  14 :       IF (IERR.NE.0) RETURN
  15 : C+PP
  15 :       CALL LIROBJ('MAILLAGE',IPT1,0,IRETOU)
  16 :       IF(IRETOU.EQ.1)THEN
  17 :         SEGACT,IPT1,MCOORD*MOD
  18 :         NBSOUS=IPT1.LISOUS(/1)
  19 :         DO IE1=1,MAX(1,NBSOUS)
  20 :           IF(NBSOUS.EQ.0)THEN
  21 :             MELEME=IPT1
  22 :           ELSE
  23 :             MELEME=IPT1.LISOUS(IE1)
  24 :             SEGACT,MELEME
  25 :           ENDIF
  26 :           NBNN=NUM(/1)
  27 :           NBELEM=NUM(/2)
  28 :           DO IE2=1,NBELEM
  29 :             DO IE3=1,NBNN
  30 :               NUME=NUM(IE3,IE2)
  31 :               IREF=NUME*(IDIM+1)
  32 :               XCOOR(IREF)=ABS(XXX)
  33 :             ENDDO
  34 :           ENDDO
  35 :           SEGDES,MELEME
  36 :         ENDDO
  37 :         SEGDES,IPT1
  38 :         RETURN
  39 :       ENDIF
  40 : C+PP
  41 :       DENSIT=ABS(XXX)
  41 :       CALL NOMREE('&DENS   ',XXX)
  42 :       END
  43 : 
  44 : 
  45 : 
  46 :  
  47 :  
  48 :  
  49 :