Télécharger la plaquette Cast3M


fr en

Télécharger exic.procedur

Retour à la liste

Numérotation des lignes : 

  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 ((<eg eps epsm)) ; quitter bloc1 ; finsi ;
  85.  fin bloc1 ;
  86.  
  87.   si ( dt > 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.  
  102.  

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