Procedures | Cast3M

   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 :