Procedures | Cast3M

   1 : * T_IPOL    PROCEDUR  MAGN      00/10/12    21:17:53     3973           
   2 : DEBPROC T_IPOL t*FLOTTANT rxt*TABLE ;
   3 : 
   4 : *------------------------------------------------------------
   5 : * Procedure interpolation des CLIM imposee en temperature
   6 : * Fonction du temps
   7 : *
   8 : * ENTREE : t : FLOTTANT : Temps auquel on veut faire interpolation
   9 : *          rxt : TABLE contenant les informations necessaires
  10 : * SORTIE : TBT.'cTp0' CHAMPOINT SOMMET actualise avec interpolation
  11 : *------------------------------------------------------------
  12 :   GEO=rxt.'GEO' ;
  13 :   TBT=rxt.'TBT' ;
  14 :   $mtpi = GEO.'$mtpi' ;
  15 :   NCLIM = TBT.'NCLIM' ;
  16 :   T_MCLIM = TBT.'T_MCLIM' ;
  17 :   T_CTP0 = TBT.'T_CTP0' ;
  18 : *--- Creation des champs points sur MTPI : cTp0
  19 :   cTp0    = 'KCHT' $mtpi SCAL SOMMET 0.0 ;
  20 : 
  21 :  'REPETER' BOUC3 NCLIM ;
  22 :   maili = T_MCLIM . &BOUC3 ;
  23 : *--- Calcul de la temperature a partir du nom de la table
  24 :   ts       = rxt.'scenario'.'t'                        ;
  25 :   ctp1     = 'IPOL' t ts  (rxt.'scenario'. (T_CTP0 . &BOUC3) ) ;
  26 :   tempo = 'MANU' 'CHPO' maili 1 'SCAL' ctp1
  27 :                       'NATURE' 'DISCRET' ;
  28 :   cTp0    = 'KCHT' $mtpi SCAL SOMMET cTp0 tempo ;
  29 :  'FIN' BOUC3 ;
  30 :   ctp0c   = NOEL $mtpi cTp0    ;
  31 :   cTp0    = ELNO $mtpi ctp0c   ;
  32 :   'SI' (('EGA' TBT.'discr' 'LINE' ) 'OU' ('EGA' ('VALE' 'DIME') 2)) ;
  33 :        FACT = 1.0 ;
  34 :   'SINON' ;
  35 :        FACT = 4.0/ 3.0 ;
  36 :   'FINSI' ;
  37 :   cTp0    = 'KCHT' $mtpi SCAL SOMMET (cTp0 * FACT ) ;
  38 :   TBT.'cTp0' = cTp0 ;
  39 : FINPROC                ;
  40 :