$$$$ PPRE
* PPRE      PROCEDUR  MAGN      96/09/28    21:28:39     2303           

  debproc ppre  rx*table                                            ;

* Operateur PPRE
* -------------
*
* Description: The PPRE operator computes the pressure gradient term
*              for the particle qdm equation.
*
*              It computes the term: grad(P)/rho_p
*
* Syntax (EQEX) :
*
*         PPRE  VN rho_p
*
* rho_p   FLOTTANT Particle density
* VN      MOT      CHPOINT (VECT SOMMET) in table INCO for the
*                  particle velocity

  uname   = rx.arg1                                                 ;
  coef    = rx.arg2                                                 ;
  rv      = rx.eqex                                                 ;
  rvp     = rx.eqex.'PRESSION'                                      ;
  tdom    = rv.domaine                                              ;

* modifica pressione

  si (exist rvp gradp)                                              ;
     gpp             = (rvp.'GRADP') *  coef                        ;
     fp              = rv.'KIZG'.uname                              ;
     rv.'KIZG'.uname = kops fp '-' gpp                              ;
  finsi                                                             ;

  finproc                                                           ;