Esope source files

   1 : C THCOND2   SOURCE    CHAT      09/10/09    21:25:24     6519           
   2 : 
   3 : C=======================================================================
   4 : C=                            T H C O N D 2                            =
   5 : C=                            -----------                              =
   6 : C=                      (TCOND2 dans le cas de la thermique)           =
   7 : C=  Fonction :                                                         =
   8 : C=  ----------                                                         =
   9 : C=  Ce sous-programme remplit le segment DESCRipteur de la matrice de  =
  10 : C=  CAPACITE (CONDUCTIVITE)dans le cas d'un modele THERMOHYDRIQUE.     =
  11 : C=                                                                     =
  12 : C=  Parametres :  (E)=Entree   (S)=Sortie                              =
  13 : C=  ------------                                                       =
  14 : C=   NBNN    (E)  Nombre de NOEUDS dans le type d'element considere    =
  15 : C=   IDESCR  (S)  Pointeur sur le segment DESCR                        =
  16 : C=                                                                     =
  17 : C=                                                                     =
  18 : C=   Zakaria HABIBI le 30 juin 2008.                                   =
  19 : C=======================================================================
  20 : 
  21 :       SUBROUTINE THCOND2 (NBNN,IDESCR)
  22 : 
  23 :       IMPLICIT INTEGER(I-N)
  24 :       IMPLICIT REAL*8 (A-H,O-Z)
  25 : 
  26 : -INC SMRIGID
  27 : 
  28 :       CHARACTER*4 MOTHER,MOCHAL,MOTGAZ,MOCGAZ,MOTCAP,MOCCAP
  29 : 
  30 :       PARAMETER ( MOTHER='T   ' )
  31 :       PARAMETER ( MOCHAL='Q   ' )
  32 :       PARAMETER ( MOTGAZ='PG  ' )
  33 :       PARAMETER ( MOCGAZ='QG  ' )
  34 :       PARAMETER ( MOTCAP='PC  ' )
  35 :       PARAMETER ( MOCCAP='QC  ' )
  36 : 
  37 : 
  38 : C  1 - Cas des elements MASSIFS
  39 : C ==============================
  40 : 
  41 :         NLIGRP=3*NBNN
  42 :         NLIGRD=3*NBNN
  43 :         SEGINI,DESCR
  44 :         DO i=1,NLIGRP
  45 :             IF (i.LE.NBNN) THEN
  46 :               LISINC(i)=MOTGAZ
  47 :               LISDUA(i)=MOCGAZ
  48 :               NOELEP(i)=i
  49 :               NOELED(i)=i
  50 :             ELSE IF (i.GT.NBNN.AND.i.LE.(2*NBNN)) THEN
  51 :               LISINC(i)=MOTCAP
  52 :               LISDUA(i)=MOCCAP
  53 :               NOELEP(i)=i-NBNN
  54 :               NOELED(i)=i-NBNN
  55 :             ELSE
  56 :               LISINC(i)=MOTHER
  57 :               LISDUA(i)=MOCHAL
  58 :               j=i-(2*NBNN)
  59 :               NOELEP(i)=j
  60 :               NOELED(i)=j
  61 :             ENDIF
  62 :         ENDDO
  63 : 
  64 :       IDESCR=DESCR
  65 :       SEGDES,DESCR
  66 : 
  67 :       RETURN
  68 :       END
  69 : 
  70 :  
  71 :  
  72 :  
  73 :