Télécharger la plaquette Cast3M


fr en

Télécharger rdiv.procedur

Retour à la liste

Numérotation des lignes : 

  1. * RDIV      PROCEDUR  MAGN      96/07/05    21:17:04     2215           
  2. *X RDIV (Procedure)
  3. *X Remise a divergence nulle
  4.  
  5. *   Procedure RDIV   Objet : Remise a divergence nulle
  6.  
  7.   DEBPROC RDIV ;
  8.   ARGU RX*TABLE ;
  9.  
  10.  rv=rx.'EQEX' ;
  11.  rvp=rv.'PRESSION' ;
  12.  iarg=rx.'IARG' ;
  13.  si( non ( ega iarg 1)) ;
  14.  mess 'Procedure RDIV : nombre d arguments incorrect ' iarg ;
  15.  quitter rdiv ;
  16.  finsi ;
  17.  si ( ega ('TYPE' rx.'ARG1') 'MOT     ') ;
  18.  U1=rv.'INCO'.(rx.'ARG1') ;
  19.  sinon ;
  20.  si ( ega ('TYPE' (rx.'ARG1')) 'CHPOINT ') ;
  21.  U1=rx.'ARG1' ;
  22.  sinon ;
  23.  mess 'Procedure RDIV : type argument invalide ' ;
  24.  quitter rdiv ;
  25.  finsi ;
  26.  finsi ;
  27.  
  28.  si (exist rv 'CLIM') ;
  29.  rvp.'CLIM'=rv.'CLIM' ;
  30.  finsi ;
  31.  
  32. NOMUT=chai ( text (extr rx.'LISTINCO' 1) ) ;
  33.  
  34. *si (non(exist rv kizd)) ;
  35. *mess 'On cree LA DIAGONALE ' ;
  36. *kdia rv ;
  37. *finsi ;
  38.  
  39.  Du =rv.'KIZD'.NOMUT;
  40.  test= exist rvp 'MATC' ;
  41.    si (test) ;
  42.    C=rvp.'MATC' ;
  43.    sinon ;
  44.    C=kmac rvp ;
  45.    rvp.'MATC'=C ;
  46.      si (non(exist rvp diagv)) ;
  47.      rvp.diagv=Du ;
  48.      finsi ;
  49.    finsi ;
  50.  
  51.    dt=1. ;
  52.  si (exist rvp 'DELTAT') ;
  53.  dt= rvp.deltat ;
  54.  finsi ;
  55.  
  56.  rvp.deltat =1. ;
  57.  P=kmf C u1 ;
  58.  
  59. *? peu conseillé
  60.    kpimp=rvp.'KPIMP' ;
  61.    pimp=rvp.'PIMP' ;
  62.    kbeta=rvp.'KBETA' ;
  63.    beta=rvp.'BETA' ;
  64.    beta=1.e-5  ;
  65.    kres rvp P 'BETA' kbeta beta 'PIMP' kpimp pimp ;
  66. *?
  67. *  kres rvp P   ;
  68.    gp= kops (( 1.) * (kmtp C P)) '/' Du ;
  69.    UR= kops gp '-' u1 ;
  70.    si (exist rv 'CLIM') ;
  71.    UR= kops ur 'CLIM' (rv.'CLIM') ;
  72.    finsi ;
  73.  
  74.    RV.'INCO'.NOMUT=UR ;
  75.  
  76.    rvp.deltat =dt;
  77.  
  78. FINPROC ;
  79.  
  80.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales