Procedures | Cast3M

   1 : * PPRE      PROCEDUR  MAGN      96/09/28    21:28:39     2303           
   2 : 
   3 :   debproc ppre  rx*table                                            ;
   4 : 
   5 : * Operateur PPRE
   6 : * -------------
   7 : *
   8 : * Description: The PPRE operator computes the pressure gradient term
   9 : *              for the particle qdm equation.
  10 : *
  11 : *              It computes the term: grad(P)/rho_p
  12 : *
  13 : * Syntax (EQEX) :
  14 : *
  15 : *         PPRE  VN rho_p
  16 : *
  17 : * rho_p   FLOTTANT Particle density
  18 : * VN      MOT      CHPOINT (VECT SOMMET) in table INCO for the
  19 : *                  particle velocity
  20 : 
  21 :   uname   = rx.arg1                                                 ;
  22 :   coef    = rx.arg2                                                 ;
  23 :   rv      = rx.eqex                                                 ;
  24 :   rvp     = rx.eqex.'PRESSION'                                      ;
  25 :   tdom    = rv.domaine                                              ;
  26 : 
  27 : * modifica pressione
  28 : 
  29 :   si (exist rvp gradp)                                              ;
  30 :      gpp             = (rvp.'GRADP') *  coef                        ;
  31 :      fp              = rv.'KIZG'.uname                              ;
  32 :      rv.'KIZG'.uname = kops fp '-' gpp                              ;
  33 :   finsi                                                             ;
  34 : 
  35 :   finproc                                                           ;
  36 :