      SUBROUTINE ETR3P1D(coor,car,iopt,AE)
      -------------------------------------
      in

      coor(npo,ndim) : coordones
      car            : caracteristiques des materiaux
                       car(1) = e car(2) = nu    si iopt = 2
                       car(1:9) mat orthotrope   si iopt = 9
      iopt           : option de calcul ici scalaire mais peut devenir 
                       vecteur si necessaire
      out

      ae            : matrice triangulaire sup

      SUBROUTINE ETM3P1D(coor,car,iopt,AE)
      -------------------------------------
      in

      coor(npo,ndim) : coordones
      car            : caracteristiques des materiaux
                       car(1) = ro    si iopt = 1
      iopt           : ouvert si masse lumping ou autre ds futur

      out

      ae            : matrice triangulaire sup

      SUBROUTINE ETS3P1D(coor,volf,sur,pres,noref,BE)
      -------------------------------------
      in

      coor(npo,ndim) : coordones
      volf           : efforts volumiques (aux noeuds) 
                       volf(3,noe)
      surf           : efforts surfaciues (aux noeuds)
                       surf(ndim,noe,nbface) 
      noref(nbface,2)          : noref(nbface,1) =  0 pas int surface
                                                 =/ 0 integrale de surface
                                 noref(nbface,2) =  0 pas de pression
                                                 =/ 0 pression
                       attention : si noref(j)=0 sur(.,.,j)=0
      out

      be            : second membre

      SUBROUTINE ETC3P1D(coor,car,iopt,SIGMA,U)
      -----------------------------------------
      in

      coor(npo,ndim) : coordones
      car            : caracteristiques des materiaux
                       car(1) = e car(2) = nu    si iopt = 1
                       car(1:9) mat orthotrope   si iopt = 2
      u              : solution aux noeudsu(noe)
      out

      sigma            : tenseur des contraintes (aux noeuds)

     Remarque : si necessaire il est possible de prevoir plusieurs u,
                alors (u(noe,ndsm) sigma(6,ndsm)) et iopt(2) = ndsm

