Procedures | Cast3M

   1 : * EXIC      PROCEDUR  GOUNAND   98/04/06    21:26:34     3135           
   2 : *X   EXIC (Procedure)
   3 : 
   4 : 
   5 : DEBPROC EXIC ;
   6 : ARGU RV*TABLE ;
   7 : 
   8 :  niter = rv.'NITER' ;
   9 :  itma  = rv.'ITMA' ;
  10 :  omeg  = rv.'OMEGA' ;
  11 :  epsm  = rv.'EPS' ;
  12 :  impr  = rv.'IMPR' ;
  13 : 
  14 :  si (exist rv 'DT') ;
  15 :  dt    = rv.'DT'  ;
  16 :  sinon ;
  17 :  itma=1 ;
  18 :  dt=0. ;
  19 :  finsi ;
  20 : 
  21 :  si (exist rv 'TPSI') ;
  22 :  tps=rv.'TPSI' ;
  23 :  sinon ;
  24 :  TPS=0. ;
  25 :  finsi ;
  26 : 
  27 :  repeter btemp itma  ;
  28 :  TPS=TPS+DT ;
  29 :  repeter bloc1 niter ;
  30 : nbop=dime (rv.listoper) ;
  31 : repeter bloc2 nbop ;
  32 : 
  33 : 
  34 :     nomper = extr &bloc2 rv.LISTOPER ;
  35 :     notable= mot (text (chai &bloc2 nomper) ) ;
  36 : *   mess ' exec oper ' notable ;
  37 :     (text nomper) (rv.notable) ;
  38 : 
  39 : 
  40 : *klop &bloc2 rv ;
  41 : * --------------------------------------------------- *
  42 : *nomper= text (extr rv.listoper &bloc2) '           ';
  43 : *nomper=nomper ;
  44 : * --------------------------------------------------- *
  45 : si ( ega &bloc2 1) ;
  46 : ma1= (rv.notable.'MATELM') ;
  47 : sinon ;
  48 : si (exist (rv.notable) 'MATELM' ) ;
  49 : ma1=ma1 et (rv.notable.'MATELM') ;
  50 : finsi ;
  51 : finsi ;
  52 : fin bloc2 ;
  53 : *             list ma1 1 ;
  54 : si ( exist rv 'SMBR');
  55 : s2 = (rv.'SMBR') ;
  56 : *list s2 ;
  57 :  res = kres ma1 'TYPI' (rv . 'METHINV')
  58 :                 'CLIM' (rv.'CLIM') 'SMBR' s2 'IMPR' 0 ;
  59 : rv.'SMBR'=1 ;
  60 : sinon ;
  61 :  res = kres ma1 'TYPI' (rv . 'METHINV')
  62 :                 'CLIM' (rv.'CLIM') 'IMPR' 0 ;
  63 : finsi ;
  64 :  ma1=0 ;
  65 : 
  66 :  eps = tcrr res omeg (rv.'INCO') 'IMPR' impr ;
  67 :  res=0. ;
  68 : 
  69 : * reactualistion linearisation
  70 : 
  71 : *- 1/ vitesse
  72 :  si (exist (rv.'INCO') 'UF');
  73 :  uf0=rv.'INCO'.'UF';
  74 :  un = rv.'INCO'.'UN' ;
  75 :  rv.'INCO'.'UF'= ksof $MT un;detr uf0 ;
  76 :  finsi ;
  77 : 
  78 :  si ( exist rv 'PPHYS') ;
  79 :  phye rv ;
  80 :  finsi ;
  81 : 
  82 :  menage ;
  83 : 
  84 :  si (( 0.) ;
  88 :    mess ' TEMPS ' TPS ' secondes ' ;
  89 :    tcnm rv ;
  90 :   finsi ;
  91 : 
  92 :  si ( exist rv 'VISU');
  93 :  visu rv ;
  94 :  finsi ;
  95 : 
  96 :  fin btemp ;
  97 :  rv.'TPSI'=tps ;
  98 : 
  99 : FINPROC ;
 100 :  
 101 :