DATE     94/01/10

     CETTE PROCEDURE A ETE MISE GRACIEUSEMENT
    A DISPOSITION DE LA COMMUNAUTE  CASTEM2000
  PAR  DELERUYELLE F.   (SOCOTEC-INDUSTRIE à L IPSN/DES)

                                                                         
      Procédure @CDG
      --------------
                                                                         
      XG1 XG2 (XG3)  =  @CDG  MODL1 (CAR1) ;
                                                                         

     Description :
     _____________

     This procedure calculates the coordinates of the center of gravity
     of a geometry contained in a model.

     Comments :
     __________

      MODL1  : model object (MMODEL type).
                                                                         
      CAR1   : geometrical properties field (not necessary)
               (MCHAML type, CARACTERISTIQUES subtype)
               If one wants to account for a non uniform cubic mass
               distribution, it is necessary to supply in CAR1 the
               field of material properties (to give RHO) and the
               field of geometrical properties (to give EPAI,
               SECT, ...).
                                                                         
      XG1    : 1-st coordinate of the center of gravity (FLOTTANT type).
                                                                         
      XG2    : 2-nd coordinate of the center of gravity (FLOTTANT type).
                                                                         
      XG3    : 3-rd coordinate of the center of gravity in a 3D case
               (FLOTTANT type).
                                                                         
      Example :
      ________

      mo1 = mode s3d mecanique elastique isotrope coq3 coq4 ;
      ca1 = cara mo1 epai 0.5 ;
      xg1 xg2 xg3 = @cdg mo1 ca1 ;
      list xg1 ; list xg2 ; list xg3 ;

      Notes :
      ________

      1) This procedure has only been tested on mechanical models.

      2) For shells, the integration domain is the mean surface of the
         shell and for beams, the integration domain is the mean line of
         the beam. If an integration over the volume is wanted, the field
         of geometrical properties CAR1 of the elements must be given.

      3) For shells, even when CAR1 is given, the excentricity of the
         shell is not taken into account.

      4) A mechanical model containing several types of elements can be
         supplied.