Télécharger la plaquette Cast3M


fr en

Télécharger nlin_decent2d.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : nlin_decent2d.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. 'OPTION' echo 0 ;
  5. *BEGINPROCEDUR gmass
  6. ************************************************************************
  7. * NOM         : GMASS
  8. * DESCRIPTION : Une matrice de masse
  9. *
  10. *
  11. *
  12. * LANGAGE     : GIBIANE-CAST3M
  13. * AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  14. *               mél : gounand@semt2.smts.cea.fr
  15. **********************************************************************
  16. * VERSION    : v2, 14/03/2006, mise à jour NLIN évolué
  17. * VERSION    : v1, 13/05/2004, version initiale
  18. * HISTORIQUE : v1, 13/05/2004, création
  19. * HISTORIQUE :
  20. * HISTORIQUE :
  21. ************************************************************************
  22. * Prière de PRENDRE LE TEMPS de compléter les commentaires
  23. * en cas de modification de ce sous-programme afin de faciliter
  24. * la maintenance !
  25. ************************************************************************
  26. *
  27. *
  28. 'DEBPROC' GMASS ;
  29. 'ARGUMENT'     _mt*'MAILLAGE' ;
  30. 'ARGUMENT'     _smt/'MAILLAGE' ;
  31. 'ARGUMENT'   gdisc*'MOT     ' ;
  32. 'ARGUMENT'    nomt*'MOT     ' ;
  33. 'ARGUMENT'   disct*'MOT     ' ;
  34. 'ARGUMENT'    nomq*'MOT     ' ;
  35. 'ARGUMENT'   discq*'MOT     ' ;
  36. 'ARGUMENT' coef/'FLOTTANT' ;
  37.    'SI' ('NON' ('EXISTE' coef)) ;
  38.       'ARGUMENT' coef2/'CHPOINT ' ;
  39.       'SI' ('NON' ('EXISTE' coef2)) ;
  40.          'ERREUR' 'Il faut donner un coef FLOTTANT ou CHPOINT' ;
  41.       'SINON' ;
  42.          coef = coef2 ;
  43.          'ARGUMENT' discc*'MOT     ' ;
  44.       'FINSI' ;                  
  45.    'SINON' ;
  46.       discc = 'CSTE' ;
  47.    'FINSI' ;
  48. 'ARGUMENT' methgau/'MOT     ' ;
  49.    'SI' ('NON' ('EXISTE' methgau)) ;
  50.       methgau = 'GAU7' ;
  51.    'FINSI' ;
  52. 'ARGUMENT' chpop/'CHPOINT' ;
  53. 'ARGUMENT' chpod/'CHPOINT' ;
  54. *
  55.    vdim = 'VALEUR' 'DIME' ;
  56.    vmod = 'VALEUR' 'MODE' ;
  57.    idim = 0 ;
  58.    'SI' ('ET' ('EGA' vdim 2) ('EGA' vmod 'PLANDEFO')) ;
  59.       idim = 2 ;
  60.       iaxi = FAUX ;
  61.    'FINSI' ;
  62.    'SI' ('ET' ('EGA' vdim 2) ('EGA' vmod 'AXIS')) ;
  63.       idim = 2 ;
  64.       iaxi = VRAI ;
  65.    'FINSI' ;
  66.    'SI' ('ET' ('EGA' vdim 3) ('EGA' vmod 'TRID')) ;
  67.       idim = 3 ;
  68.       iaxi = FAUX ;
  69.    'FINSI' ;
  70.    'SI' ('EGA' vdim 1) ;
  71.       idim = 1 ;
  72.       iaxi = FAUX ;
  73.    'FINSI' ;
  74.    'SI' ('EGA' idim 0) ;
  75.       'ERREUR' ('CHAINE' 'vmod=' vmod ' et vdim=' vdim ' non prevu') ;
  76.    'FINSI' ;
  77.    'SI' iaxi ;
  78.       dprmt = '*' ('COORDONNEE' 1 _mt) ('*' PI 2.D0) ;
  79.    'FINSI' ;      
  80.    numop  = 1 ;
  81.    numder = idim ;
  82.    mmt = 'MOTS' nomt ;
  83.    mmq = 'MOTS' nomq ;
  84.    numvar = 1 ;
  85.    numdat = 1 ;
  86.    numcof = 1 ;
  87. *
  88.    A = ININLIN numop numvar numdat numcof numder ;   
  89.    A . 'VAR' . 1 . 'NOMDDL' = mmt ;
  90.    A . 'VAR' . 1 . 'DISC'   = disct ;
  91.    'SI' ('EXISTE' chpop) ;
  92.       A . 'VAR' . 1 . 'VALEUR'   = chpop ;
  93.    'FINSI' ;   
  94.    A . 'DAT' . 1 . 'NOMDDL' = 'MOTS' 'SCAL' ;
  95.    A . 'DAT' . 1 . 'DISC'   = discc ;
  96.    A . 'DAT' . 1 . 'VALEUR' = coef ;
  97.    A . 'COF' . 1 . 'COMPOR' = 'IDEN' ;
  98.    A . 'COF' . 1 . 'LDAT'   = 'LECT' 1 ;
  99. *
  100.    A . 1 . 1 . 0 = 'LECT' 1 ;
  101. *
  102.    'SI' iaxi ;
  103.       numdat = 1 ;
  104.       numcof = 1 ;
  105.    'SINON' ;      
  106.       numdat = 0 ;
  107.       numcof = 0 ;
  108.    'FINSI' ;      
  109.    B = ININLIN numop numvar numdat numcof numder ;   
  110.    B . 'VAR' . 1 . 'NOMDDL' = mmq ;
  111.    B . 'VAR' . 1 . 'DISC'   = discq ;
  112.    'SI' ('EXISTE' chpod) ;
  113.       B . 'VAR' . 1 . 'VALEUR'   = chpod ;
  114.    'FINSI' ;   
  115. *
  116.    'SI' iaxi ;
  117.       B . 'DAT' . 1 . 'NOMDDL' = 'MOTS' 'SCAL' ;
  118.       B . 'DAT' . 1 . 'DISC'   = gdisc ;
  119.       B . 'DAT' . 1 . 'VALEUR' = dprmt ;
  120.       B . 'COF' . 1 . 'COMPOR' = 'IDEN' ;
  121.       B . 'COF' . 1 . 'LDAT'   = 'LECT' 1 ;
  122.    'FINSI' ;
  123.    'SI' iaxi ;      
  124.       B . 1 . 1 . 0 = 'LECT' 1 ;
  125.    'SINON' ;      
  126.       B . 1 . 1 . 0 = 'LECT' ;
  127.    'FINSI' ;      
  128. *
  129.    'SI' ('EXISTE' _smt) ;
  130.        mgmass = 'NLIN' gdisc _mt _smt A B methgau ;
  131.    'SINON' ;       
  132.        mgmass = 'NLIN' gdisc _mt A B methgau ;
  133.    'FINSI' ;       
  134. *
  135. 'RESPRO' mgmass ;
  136. 'FINPROC' ;
  137. *
  138. * End of procedure file GMASS
  139. *
  140. *ENDPROCEDUR gmass
  141. *BEGINPROCEDUR gugrad
  142. ************************************************************************
  143. * NOM         : GUGRAD
  144. * DESCRIPTION : U . grad
  145. *
  146. *
  147. *
  148. * LANGAGE     : GIBIANE-CAST3M
  149. * AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  150. *               mél : gounand@semt2.smts.cea.fr
  151. **********************************************************************
  152. * VERSION    : v1, 13/05/2004, version initiale
  153. * HISTORIQUE : v1, 13/05/2004, création
  154. * HISTORIQUE :
  155. * HISTORIQUE :
  156. ************************************************************************
  157. * Prière de PRENDRE LE TEMPS de compléter les commentaires
  158. * en cas de modification de ce sous-programme afin de faciliter
  159. * la maintenance !
  160. ************************************************************************
  161. *
  162. *
  163. 'DEBPROC' GUGRAD ;
  164. 'ARGUMENT'     _mt*'MAILLAGE' ;
  165. 'ARGUMENT'   gdisc*'MOT     ' ;
  166. 'ARGUMENT'    nomp*'MOT     ' ;
  167. 'ARGUMENT'   discp*'MOT     ' ;
  168. 'ARGUMENT'    nomd*'MOT     ' ;
  169. 'ARGUMENT'   discd*'MOT     ' ;
  170. 'ARGUMENT'    vtot*'CHPOINT ' ;
  171. 'ARGUMENT'    vcomp*'LISTMOTS' ;
  172. 'ARGUMENT'   discv*'MOT     ' ;
  173. 'ARGUMENT' coef/'FLOTTANT' ;
  174.    'SI' ('NON' ('EXISTE' coef)) ;
  175.       'ARGUMENT' coef2/'CHPOINT ' ;
  176.       'SI' ('NON' ('EXISTE' coef2)) ;
  177.          'ERREUR' 'Il faut donner un coef FLOTTANT ou CHPOINT' ;
  178.       'SINON' ;
  179.          coef = coef2 ;
  180.          'ARGUMENT' discc*'MOT     ' ;
  181.       'FINSI' ;                  
  182.    'SINON' ;
  183.       discc = 'CSTE' ;
  184.    'FINSI' ;            
  185. 'ARGUMENT' methgau/'MOT     ' ;
  186.    'SI' ('NON' ('EXISTE' methgau)) ;
  187.       methgau = 'GAU7' ;
  188.    'FINSI' ;      
  189. 'ARGUMENT' chpop/'CHPOINT' ;
  190. 'ARGUMENT' chpod/'CHPOINT' ;
  191. *
  192.    vdim = 'VALEUR' 'DIME' ;
  193.    vmod = 'VALEUR' 'MODE' ;
  194.    idim = 0 ;
  195.    'SI' ('ET' ('EGA' vdim 2) ('EGA' vmod 'PLANDEFO')) ;
  196.       idim = 2 ;
  197.       iaxi = FAUX ;
  198. *      Passé en argument désormais      
  199. *      vcomp = 'MOTS' 'UX' 'UY' ;
  200.    'FINSI' ;
  201.    'SI' ('ET' ('EGA' vdim 2) ('EGA' vmod 'AXIS')) ;
  202.       idim = 2 ;
  203.       iaxi = VRAI ;
  204. *      vcomp = 'MOTS' 'UR' 'UZ' ;
  205.    'FINSI' ;
  206.    'SI' ('ET' ('EGA' vdim 3) ('EGA' vmod 'TRID')) ;
  207.       idim = 3 ;
  208.       iaxi = FAUX ;
  209. *      vcomp = 'MOTS' 'UX' 'UY' 'UZ' ;
  210.    'FINSI' ;
  211.    'SI' ('EGA' vdim 1) ;
  212.       idim = 1 ;
  213.       iaxi = FAUX ;
  214.    'FINSI' ;
  215.    'SI' ('EGA' idim 0) ;
  216.       'ERREUR' ('CHAINE' 'vmod=' vmod ' et vdim=' vdim ' non prevu') ;
  217.    'FINSI' ;
  218.    'SI' iaxi ;
  219.       dprmt = '*' ('COORDONNEE' 1 _mt) ('*' PI 2.D0) ;
  220.    'FINSI' ;      
  221. *
  222.    numop  = 1 ;
  223.    numder = idim ;
  224.    mmp = 'MOTS' nomp ;
  225.    mmd = 'MOTS' nomd ;
  226.    numvar = 1 ;
  227.    numdat = idim ;
  228.    numcof = idim ;
  229. *
  230.    A = ININLIN numop numvar numdat numcof numder ;   
  231.    A . 'VAR' . 1 . 'NOMDDL' = mmp ;
  232.    A . 'VAR' . 1 . 'DISC'   = discp ;
  233.    'SI' ('EXISTE' chpop) ;
  234.       A . 'VAR' . 1 . 'VALEUR'   = chpop ;
  235.    'FINSI' ;   
  236.    'REPETER' iidim idim ;
  237.       nomco = 'EXTRAIRE' vcomp &iidim ;
  238.       A . 'DAT' . &iidim . 'NOMDDL' = 'MOTS' nomco ;
  239.       A . 'DAT' . &iidim . 'DISC'   = discv ;
  240. *      A . 'DAT' . &iidim . 'VALEUR' = 'EXCO' nomco 'NOID' vtot nomco ;
  241.       A . 'DAT' . &iidim . 'VALEUR' = 'EXCO' nomco vtot nomco ;
  242. *      
  243.       A . 'COF' . &iidim . 'COMPOR' = 'IDEN' ;
  244.       A . 'COF' . &iidim . 'LDAT'   = 'LECT' &iidim ;
  245.    'FIN' iidim ;      
  246. *
  247.    'REPETER' iidim idim ;
  248.       A . 1 . 1 . &iidim = 'LECT' &iidim ;
  249.    'FIN' iidim ;      
  250. *
  251.    numvar = 1 ;
  252.    'SI' iaxi ;
  253.       numdat = 2 ;
  254.       numcof = 2 ;
  255.    'SINON' ;      
  256.       numdat = 1 ;
  257.       numcof = 1 ;
  258.    'FINSI' ;      
  259. *
  260.    B = ININLIN numop numvar numdat numcof numder ;   
  261.    B . 'VAR' . 1 . 'NOMDDL' = mmd ;
  262.    B . 'VAR' . 1 . 'DISC'   = discd ;
  263.    'SI' ('EXISTE' chpod) ;
  264.       B . 'VAR' . 1 . 'VALEUR'   = chpod ;
  265.    'FINSI' ;   
  266. *   
  267.    B . 'DAT' . 1 . 'NOMDDL' = 'MOTS' 'SCAL' ;
  268.    B . 'DAT' . 1 . 'DISC'   = discc ;
  269.    B . 'DAT' . 1 . 'VALEUR' = coef ;
  270.    B . 'COF' . 1 . 'COMPOR' = 'IDEN' ;
  271.    B . 'COF' . 1 . 'LDAT'   = 'LECT' 1 ;
  272. *
  273.    'SI' iaxi ;
  274.       B . 'DAT' . 2 . 'NOMDDL' = 'MOTS' 'SCAL' ;
  275.       B . 'DAT' . 2 . 'DISC'   = gdisc ;
  276.       B . 'DAT' . 2 . 'VALEUR' = dprmt ;
  277.       B . 'COF' . 2 . 'COMPOR' = 'IDEN' ;
  278.       B . 'COF' . 2 . 'LDAT'   = 'LECT' 2 ;
  279.    'FINSI' ;
  280. *
  281.    'SI' iaxi ;      
  282.       B . 1 . 1 . 0 = 'LECT' 1 2 ;
  283.    'SINON' ;      
  284.       B . 1 . 1 . 0 = 'LECT' 1 ;
  285.    'FINSI' ;      
  286. *
  287.    mgugrad = 'NLIN' gdisc _mt A B methgau ;
  288. *
  289. 'RESPRO' mgugrad ;
  290. 'FINPROC' ;
  291. *
  292. * End of procedure file GUGRAD
  293. *
  294. *ENDPROCEDUR gugrad
  295. *BEGINPROCEDUR glapn
  296. ************************************************************************
  297. * NOM         : GLAPN
  298. * DESCRIPTION : Un laplacien scalaire 
  299. *
  300. *
  301. *
  302. * LANGAGE     : GIBIANE-CAST3M
  303. * AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  304. *               mél : gounand@semt2.smts.cea.fr
  305. **********************************************************************
  306. * VERSION    : v2, 14/03/2006, mise à jour NLIN évolué
  307. * VERSION    : v1, 13/05/2004, version initiale
  308. * HISTORIQUE : v1, 13/05/2004, création
  309. * HISTORIQUE :
  310. * HISTORIQUE :
  311. ************************************************************************
  312. * Prière de PRENDRE LE TEMPS de compléter les commentaires
  313. * en cas de modification de ce sous-programme afin de faciliter
  314. * la maintenance !
  315. ************************************************************************
  316. *
  317. *
  318. 'DEBPROC' GLAPN ;
  319. 'ARGUMENT'     _mt*'MAILLAGE' ;
  320. 'ARGUMENT'   gdisc*'MOT     ' ;
  321. 'ARGUMENT'    nomt*'MOT     ' ;
  322. 'ARGUMENT'   disct*'MOT     ' ;
  323. 'ARGUMENT'    nomq*'MOT     ' ;
  324. 'ARGUMENT'   discq*'MOT     ' ;
  325. 'ARGUMENT' coef/'FLOTTANT' ;
  326.    'SI' ('NON' ('EXISTE' coef)) ;
  327.       'ARGUMENT' coef2/'CHPOINT ' ;
  328.       'SI' ('NON' ('EXISTE' coef2)) ;
  329.          'ERREUR' 'Il faut donner un coef FLOTTANT ou CHPOINT' ;
  330.       'SINON' ;
  331.          coef = coef2 ;
  332.          'ARGUMENT' discc*'MOT     ' ;
  333.       'FINSI' ;                  
  334.    'SINON' ;
  335.       discc = 'CSTE' ;
  336.    'FINSI' ;            
  337.    'ARGUMENT' methgau/'MOT     ' ;
  338.    'SI' ('NON' ('EXISTE' methgau)) ;
  339.       methgau = 'GAU7' ;
  340.    'FINSI' ;      
  341. 'ARGUMENT' chpop/'CHPOINT' ;
  342. 'ARGUMENT' chpod/'CHPOINT' ;
  343. *
  344.    vdim = 'VALEUR' 'DIME' ;
  345.    vmod = 'VALEUR' 'MODE' ;
  346.    idim = 0 ;
  347.    'SI' ('ET' ('EGA' vdim 2) ('EGA' vmod 'PLANDEFO')) ;
  348.       idim = 2 ;
  349.       iaxi = FAUX ;
  350.    'FINSI' ;
  351.    'SI' ('ET' ('EGA' vdim 2) ('EGA' vmod 'AXIS')) ;
  352.       idim = 2 ;
  353.       iaxi = VRAI ;
  354.    'FINSI' ;
  355.    'SI' ('ET' ('EGA' vdim 3) ('EGA' vmod 'TRID')) ;
  356.       idim = 3 ;
  357.       iaxi = FAUX ;
  358.    'FINSI' ;
  359.    'SI' ('EGA' vdim 1) ;
  360.       idim = 1 ;
  361.       iaxi = FAUX ;
  362.    'FINSI' ;
  363.    'SI' ('EGA' idim 0) ;
  364.       'ERREUR' ('CHAINE' 'vmod=' vmod ' et vdim=' vdim ' non prevu') ;
  365.    'FINSI' ;
  366.    'SI' iaxi ;
  367.       dprmt = '*' ('COORDONNEE' 1 _mt) ('*' PI 2.D0) ;
  368.    'FINSI' ;      
  369. * Test bête...
  370.    'SI' ('EGA' ('TYPE' coef) 'CHPOINT ') ;
  371.       mincoef = 'MINIMUM' coef ;
  372.    'SINON' ;
  373.       mincoef = coef ;
  374.    'FINSI' ;            
  375.    'SI' ('<' mincoef 0.D0) ;
  376.       'ERREUR' 'Le coef (une viscosité) doit etre positive' ;
  377.    'FINSI' ;      
  378. *
  379.    'SI' iaxi ;
  380.       lcoef = 'MOTS' 'NURR' 'NUZZ' ;
  381.    'SINON' ;
  382.       'SI' ('EGA' idim 2) ;
  383.          lcoef = 'MOTS' 'NUXX' 'NUYY' ;
  384.       'SINON' ;
  385.          lcoef = 'MOTS' 'NUXX' 'NUYY' 'NUZZ' ;
  386.       'FINSI' ;
  387.    'FINSI' ;      
  388. *
  389.    ltens = FAUX ;
  390.    'SI' ('EGA' ('TYPE' coef) 'CHPOINT') ;
  391.       ncomp = 'DIME' ('EXTRAIRE' coef 'COMP') ;
  392.       'SI' ('EGA' ncomp idim) ;
  393.          ltens = VRAI ;
  394.       'FINSI' ;
  395.    'FINSI' ;
  396. *
  397.    tcoef = 'TABLE' ;
  398.    'SI' ltens ;
  399.       'REPETER' iidim idim ;
  400.          tcoef . &iidim = 'EXCO' ('EXTRAIRE' lcoef &iidim) coef ;
  401.       'FIN' iidim ;
  402.    'SINON' ;               
  403.       'REPETER' iidim idim ;
  404.          tcoef . &iidim = coef ;
  405.       'FIN' iidim ;
  406.    'FINSI' ;
  407. *      
  408.    numop  = idim ;
  409.    numder = idim ;
  410.    mmt = 'MOTS' nomt ;
  411.    mmq = 'MOTS' nomq ;
  412.    numvar = 1 ;
  413.    numdat = idim ;
  414.    numcof = idim ;
  415. *
  416.    A = ININLIN numop numvar numdat numcof numder ;   
  417.    A . 'VAR' . 1 . 'NOMDDL' = mmt ;
  418.    A . 'VAR' . 1 . 'DISC'   = disct ;
  419.    'SI' ('EXISTE' chpop) ;
  420.       A . 'VAR' . 1 . 'VALEUR'   = chpop ;
  421.    'FINSI' ;   
  422.    'REPETER' iidim idim ;
  423.       A . 'DAT' . &iidim . 'NOMDDL' = 'MOTS' 'SCAL' ;
  424.       A . 'DAT' . &iidim . 'DISC'   = discc ;
  425.       A . 'DAT' . &iidim . 'VALEUR' = tcoef . &iidim ;
  426.       A . 'COF' . &iidim . 'COMPOR' = 'IDEN' ;
  427.       A . 'COF' . &iidim . 'LDAT'   = 'LECT' &iidim ;
  428.    'FIN' iidim ;
  429. *
  430.    'REPETER' iidim idim ;
  431.       A . &iidim . 1 . &iidim = 'LECT' &iidim ;
  432.    'FIN' iidim ;      
  433. *
  434.    'SI' iaxi ;
  435.       numdat = 1 ;
  436.       numcof = 1 ;
  437.    'SINON' ;      
  438.       numdat = 0 ;
  439.       numcof = 0 ;
  440.    'FINSI' ;      
  441. *
  442.    B = ININLIN numop numvar numdat numcof numder ;   
  443.    B . 'VAR' . 1 . 'NOMDDL' = mmq ;
  444.    B . 'VAR' . 1 . 'DISC'   = discq ;
  445.    'SI' ('EXISTE' chpod) ;
  446.       B . 'VAR' . 1 . 'VALEUR'   = chpod ;
  447.    'FINSI' ;   
  448.    'SI' iaxi ;
  449.       B . 'DAT' . 1 . 'NOMDDL' = 'MOTS' 'SCAL' ;
  450.       B . 'DAT' . 1 . 'DISC'   = gdisc ;
  451.       B . 'DAT' . 1 . 'VALEUR' = dprmt ;
  452.       B . 'COF' . 1 . 'COMPOR' = 'IDEN' ;
  453.       B . 'COF' . 1 . 'LDAT'   = 'LECT' 1 ;
  454.    'FINSI' ;
  455. *
  456.    'REPETER' iidim idim ;
  457.       'SI' iaxi ;      
  458.          B . &iidim . 1 . &iidim = 'LECT' 1 ;
  459.       'SINON' ;      
  460.          B . &iidim . 1 . &iidim = 'LECT' ;
  461.       'FINSI' ;      
  462.    'FIN' iidim ;      
  463. *
  464.    mglapn = 'NLIN' gdisc _mt A B methgau ;
  465. * Integration par parties
  466.    mglapn = '*' mglapn -1.D0 ;
  467. *
  468. 'RESPRO' mglapn ;
  469. 'FINPROC' ;
  470. *
  471. * End of procedure file GLAPN
  472. *
  473. *ENDPROCEDUR glapn
  474. *BEGINPROCEDUR gdecent
  475. ************************************************************************
  476. * NOM         : GDECENT
  477. * DESCRIPTION : Décentrement
  478. *
  479. *
  480. *
  481. * LANGAGE     : GIBIANE-CAST3M
  482. * AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  483. *               mél : gounand@semt2.smts.cea.fr
  484. **********************************************************************
  485. * VERSION    : v1, 13/05/2004, version initiale
  486. * HISTORIQUE : v1, 13/05/2004, création
  487. * HISTORIQUE :
  488. * HISTORIQUE :
  489. ************************************************************************
  490. * Prière de PRENDRE LE TEMPS de compléter les commentaires
  491. * en cas de modification de ce sous-programme afin de faciliter
  492. * la maintenance !
  493. ************************************************************************
  494. *
  495. *
  496. 'DEBPROC' GDECENT ;
  497. 'ARGUMENT' imeth*'ENTIER' ;
  498. 'ARGUMENT'     _mt*'MAILLAGE' ;
  499. 'ARGUMENT'   gdisc*'MOT     ' ;
  500. 'ARGUMENT'    nomp*'MOT     ' ;
  501. 'ARGUMENT'   discp*'MOT     ' ;
  502. 'ARGUMENT'    nomd*'MOT     ' ;
  503. 'ARGUMENT'   discd*'MOT     ' ;
  504. 'ARGUMENT'    vtot*'CHPOINT ' ;
  505. 'ARGUMENT'    vcomp*'LISTMOTS' ;
  506. 'ARGUMENT'   discv*'MOT     ' ;
  507. 'ARGUMENT'     rho*'FLOTTANT' ;
  508. 'ARGUMENT'      nu*'FLOTTANT' ;
  509. 'ARGUMENT'     Pec*'FLOTTANT' ;
  510. 'ARGUMENT' methgau/'MOT     ' ;
  511.    'SI' ('NON' ('EXISTE' methgau)) ;
  512.       methgau = 'GAU7' ;
  513.    'FINSI' ;      
  514. 'ARGUMENT' chpop/'CHPOINT' ;
  515. 'ARGUMENT' chpod/'CHPOINT' ;
  516. *
  517.    vdim = 'VALEUR' 'DIME' ;
  518.    vmod = 'VALEUR' 'MODE' ;
  519.    idim = 0 ;
  520.    'SI' ('ET' ('EGA' vdim 2) ('EGA' vmod 'PLANDEFO')) ;
  521.       idim = 2 ;
  522.       iaxi = FAUX ;
  523. *      Passé en argument désormais      
  524. *      vcomp = 'MOTS' 'UX' 'UY' ;
  525.    'FINSI' ;
  526.    'SI' ('ET' ('EGA' vdim 2) ('EGA' vmod 'AXIS')) ;
  527.       idim = 2 ;
  528.       iaxi = VRAI ;
  529. *      vcomp = 'MOTS' 'UR' 'UZ' ;
  530.    'FINSI' ;
  531.    'SI' ('ET' ('EGA' vdim 3) ('EGA' vmod 'TRID')) ;
  532.       idim = 3 ;
  533.       iaxi = FAUX ;
  534. *      vcomp = 'MOTS' 'UX' 'UY' 'UZ' ;
  535.    'FINSI' ;
  536.    'SI' ('EGA' vdim 1) ;
  537.       idim = 1 ;
  538.       iaxi = FAUX ;
  539.    'FINSI' ;
  540.    'SI' ('EGA' idim 0) ;
  541.       'ERREUR' ('CHAINE' 'vmod=' vmod ' et vdim=' vdim ' non prevu') ;
  542.    'FINSI' ;
  543.    'SI' iaxi ;
  544.       dprmt = '*' ('COORDONNEE' 1 _mt) ('*' PI 2.D0) ;
  545.    'FINSI' ;      
  546. *
  547.    mdumm = 'MOTS' 'DUMM' ;
  548.    numop  = 1 ;
  549.    numder = idim ;
  550.    mmp = 'MOTS' nomp ;
  551.    mmd = 'MOTS' nomd ;
  552.    numvar = 1 ;
  553.    numdat = '+' idim 3 ;
  554.    numcof = idim ;
  555.    lisdat = 'LECT' 1 PAS 1 ('+' idim 3) ;
  556. *
  557.    A = ININLIN numop numvar numdat numcof numder ;   
  558.    A . 'VAR' . 1 . 'NOMDDL' = mmp ;
  559.    A . 'VAR' . 1 . 'DISC'   = discp ;
  560.    'SI' ('EXISTE' chpop) ;
  561.       A . 'VAR' . 1 . 'VALEUR'   = chpop ;
  562.    'FINSI' ;
  563.    idat = 0 ;
  564.    idat = '+' idat 1 ;
  565.    A . 'DAT' . idat . 'NOMDDL' = mdumm  ;
  566.    A . 'DAT' . idat . 'DISC'   = 'CSTE' ;
  567.    A . 'DAT' . idat . 'VALEUR' = rho    ;
  568.    idat = '+' idat 1 ;
  569.    A . 'DAT' . idat . 'NOMDDL' = mdumm  ;
  570.    A . 'DAT' . idat . 'DISC'   = 'CSTE' ;
  571.    A . 'DAT' . idat . 'VALEUR' = nu    ;
  572.    'REPETER' iidim idim ;
  573.       nomco = 'EXTRAIRE' vcomp &iidim ;
  574.       idat = '+' idat 1 ;
  575.       A . 'DAT' . idat . 'NOMDDL' = 'MOTS' nomco ;
  576.       A . 'DAT' . idat . 'DISC'   = discv ;
  577. *      A . 'DAT' . idat . 'VALEUR' = 'EXCO' nomco 'NOID' vtot nomco ;
  578.       A . 'DAT' . idat . 'VALEUR' = 'EXCO' nomco vtot nomco ;
  579.    'FIN' iidim ;
  580.    idat = '+' idat 1 ;
  581.    A . 'DAT' . idat . 'NOMDDL' = mdumm  ;
  582.    A . 'DAT' . idat . 'DISC'   = 'CSTE' ;
  583.    A . 'DAT' . idat . 'VALEUR' = Pec    ;
  584. *
  585.    icof = 0 ;
  586.    'REPETER' iidim idim ;
  587.       icof = '+' icof 1 ;
  588.       A . 'COF' . icof . 'COMPOR' = 'CHAINE' 'MUSTAB' imeth &iidim ;
  589.       A . 'COF' . icof . 'LDAT'   = lisdat ;
  590.    'FIN' iidim ;         
  591. *
  592.    'REPETER' iidim idim ;
  593.       A . 1 . 1 . &iidim = 'LECT' &iidim ;
  594.    'FIN' iidim ;      
  595. *
  596.    numvar = 1 ;
  597.    'SI' iaxi ;
  598.       numdat = '+' idim 4 ;
  599.       numcof = '+' idim 1 ;
  600.    'SINON' ;      
  601.       numdat = '+' idim 3 ;
  602.       numcof = idim ;
  603.    'FINSI' ;      
  604. *
  605.    B = ININLIN numop numvar numdat numcof numder ;   
  606.    B . 'VAR' . 1 . 'NOMDDL' = mmd ;
  607.    B . 'VAR' . 1 . 'DISC'   = discd ;
  608.    'SI' ('EXISTE' chpod) ;
  609.       B . 'VAR' . 1 . 'VALEUR'   = chpod ;
  610.    'FINSI' ;   
  611.    idat = 0 ;
  612.    idat = '+' idat 1 ;
  613.    B . 'DAT' . idat . 'NOMDDL' = mdumm  ;
  614.    B . 'DAT' . idat . 'DISC'   = 'CSTE' ;
  615.    B . 'DAT' . idat . 'VALEUR' = rho    ;
  616.    idat = '+' idat 1 ;
  617.    B . 'DAT' . idat . 'NOMDDL' = mdumm  ;
  618.    B . 'DAT' . idat . 'DISC'   = 'CSTE' ;
  619.    B . 'DAT' . idat . 'VALEUR' = nu    ;
  620.    'REPETER' iidim idim ;
  621.       nomco = 'EXTRAIRE' vcomp &iidim ;
  622.       idat = '+' idat 1 ;
  623.       B . 'DAT' . idat . 'NOMDDL' = 'MOTS' nomco ;
  624.       B . 'DAT' . idat . 'DISC'   = discv ;
  625. *      A . 'DAT' . idat . 'VALEUR' = 'EXCO' nomco 'NOID' vtot nomco ;
  626.       B . 'DAT' . idat . 'VALEUR' = 'EXCO' nomco vtot nomco ;
  627.    'FIN' iidim ;
  628.    idat = '+' idat 1 ;
  629.    B . 'DAT' . idat . 'NOMDDL' = mdumm  ;
  630.    B . 'DAT' . idat . 'DISC'   = 'CSTE' ;
  631.    B . 'DAT' . idat . 'VALEUR' = Pec    ;
  632. *
  633.    icof = 0 ;
  634.    'REPETER' iidim idim ;
  635.       icof = '+' icof 1 ;
  636.       B . 'COF' . icof . 'COMPOR' = 'CHAINE' 'MUSTAB' imeth &iidim ;
  637.       B . 'COF' . icof . 'LDAT'   = lisdat ;
  638.    'FIN' iidim ;         
  639. *
  640.    'SI' iaxi ;
  641.       idat = '+' idat 1 ;
  642.       B . 'DAT' . idat . 'NOMDDL' = 'MOTS' 'SCAL' ;
  643.       B . 'DAT' . idat . 'DISC'   = gdisc ;
  644.       B . 'DAT' . idat . 'VALEUR' = dprmt ;
  645.       icof = '+' icof 1 ;
  646.       B . 'COF' . icof . 'COMPOR' = 'IDEN' ;
  647.       B . 'COF' . icof . 'LDAT'   = 'LECT' idat ;
  648.    'FINSI' ;
  649. *
  650.    'SI' iaxi ;      
  651.       'REPETER' iidim idim ;
  652.          B . 1 . 1 . &iidim = 'LECT' &iidim icof ;
  653.       'FIN' iidim ;      
  654.    'SINON' ;      
  655.       'REPETER' iidim idim ;
  656.          B . 1 . 1 . &iidim = 'LECT' &iidim ;
  657.       'FIN' iidim ;      
  658.    'FINSI' ;      
  659. *
  660.    mgdecent = 'NLIN' gdisc _mt A B methgau ;
  661. *
  662. 'RESPRO' mgdecent ;
  663. 'FINPROC' ;
  664. *
  665. * End of procedure file GDECENT
  666. *
  667. *ENDPROCEDUR gdecent
  668. *BEGINPROCEDUR gmustab
  669. ************************************************************************
  670. * NOM         : GMUSTAB
  671. * DESCRIPTION : Décentrement
  672. *
  673. *
  674. *
  675. * LANGAGE     : GIBIANE-CAST3M
  676. * AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  677. *               mél : gounand@semt2.smts.cea.fr
  678. **********************************************************************
  679. * VERSION    : v1, 13/05/2004, version initiale
  680. * HISTORIQUE : v1, 13/05/2004, création
  681. * HISTORIQUE :
  682. * HISTORIQUE :
  683. ************************************************************************
  684. * Prière de PRENDRE LE TEMPS de compléter les commentaires
  685. * en cas de modification de ce sous-programme afin de faciliter
  686. * la maintenance !
  687. ************************************************************************
  688. *
  689. *
  690. 'DEBPROC' GMUSTAB ;
  691. 'ARGUMENT' imeth*'ENTIER' ; 
  692. 'ARGUMENT'     _mt*'MAILLAGE' ;
  693. 'ARGUMENT'   gdisc*'MOT     ' ;
  694. 'ARGUMENT'    vtot*'CHPOINT ' ;
  695. 'ARGUMENT'    vcomp*'LISTMOTS' ;
  696. 'ARGUMENT'   discv*'MOT     ' ;
  697. 'ARGUMENT'     rho*'FLOTTANT' ;
  698. 'ARGUMENT'      nu*'FLOTTANT' ;
  699. 'ARGUMENT'     Pec*'FLOTTANT' ;
  700. 'ARGUMENT' methgau/'MOT     ' ;
  701.    'SI' ('NON' ('EXISTE' methgau)) ;
  702.       methgau = 'GAU7' ;
  703.    'FINSI' ;      
  704. *
  705.    vdim = 'VALEUR' 'DIME' ;
  706. *
  707.    mdumm = 'MOTS' 'DUMM' ;
  708.    numop  = 1 ;
  709.    numder = vdim ;
  710.    numvar = 1 ;
  711.    numdat = '+' vdim 3 ;
  712.    numcof = 1 ;
  713.    lisdat = 'LECT' 1 PAS 1 ('+' vdim 3) ;
  714. *
  715.    A = ININLIN numop numvar numdat numcof numder ;   
  716.    A . 'VAR' . 1 . 'NOMDDL' = mdumm  ;
  717.    A . 'VAR' . 1 . 'DISC'   = 'CSTE' ;
  718.    A . 'VAR' . 1 . 'VALEUR' = 1.D0   ;
  719.    idat = 0 ;
  720.    idat = '+' idat 1 ;
  721.    A . 'DAT' . idat . 'NOMDDL' = mdumm  ;
  722.    A . 'DAT' . idat . 'DISC'   = 'CSTE' ;
  723.    A . 'DAT' . idat . 'VALEUR' = rho    ;
  724.    idat = '+' idat 1 ;
  725.    A . 'DAT' . idat . 'NOMDDL' = mdumm  ;
  726.    A . 'DAT' . idat . 'DISC'   = 'CSTE' ;
  727.    A . 'DAT' . idat . 'VALEUR' = nu    ;
  728.    'REPETER' ivdim vdim ;
  729.       nomco = 'EXTRAIRE' vcomp &ivdim ;
  730.       idat = '+' idat 1 ;
  731.       A . 'DAT' . idat . 'NOMDDL' = 'MOTS' nomco ;
  732.       A . 'DAT' . idat . 'DISC'   = discv ;
  733. *      A . 'DAT' . idat . 'VALEUR' = 'EXCO' nomco 'NOID' vtot nomco ;
  734.       A . 'DAT' . idat . 'VALEUR' = 'EXCO' nomco vtot nomco ;
  735.    'FIN' ivdim ;
  736.    idat = '+' idat 1 ;
  737.    A . 'DAT' . idat . 'NOMDDL' = mdumm  ;
  738.    A . 'DAT' . idat . 'DISC'   = 'CSTE' ;
  739.    A . 'DAT' . idat . 'VALEUR' = Pec    ;
  740. *
  741.    A . 'COF' . 1 . 'COMPOR' = 'CHAINE' 'MUSTAB' imeth '0' ;
  742.    A . 'COF' . 1 . 'LDAT'   = lisdat ;
  743. *
  744.    A . 1 . 1 . 0 = 'LECT' 1 ;
  745. *
  746.    numvar = 1 ;
  747.    numdat = 0 ;
  748.    numcof = 0 ;
  749. *
  750.    B = ININLIN numop numvar numdat numcof numder ;   
  751.    B . 'VAR' . 1 . 'NOMDDL' = mdumm ;
  752.    B . 'VAR' . 1 . 'DISC'   = 'CSTE' ;
  753.    B . 'VAR' . 1 . 'VALEUR' = 1.D0 ;
  754. *
  755.    B . 1 . 1 . 0 = 'LECT' ;
  756. *
  757.    mgmustab = 'NLIN' gdisc _mt A B 'ERF1' methgau ;
  758. *
  759. 'RESPRO' mgmustab ;
  760. 'FINPROC' ;
  761. *
  762. * End of procedure file GMUSTAB
  763. *
  764. *ENDPROCEDUR gmustab
  765. *BEGINPROCEDUR formar
  766. ************************************************************************
  767. * NOM         : FORMAR
  768. * DESCRIPTION : formate un réel de facon courte
  769. *               pratique pour les noms de
  770. *               sauvegarde
  771. *               Exemples :
  772. * 'MESSAGE' ('CHAINE' (formar 2.9e5   1)) ;
  773. * 2.9E5
  774. * 'MESSAGE' ('CHAINE' (formar -2.9e5  1)) ;
  775. * -2.9E5
  776. * 'MESSAGE' ('CHAINE' (formar 2.9e-5  1)) ;
  777. * 2.9E-5
  778. * 'MESSAGE' ('CHAINE' (formar -2.9e-5 1)) ;
  779. * -2.9E-5
  780. * 'MESSAGE' ('CHAINE' (formar 2.9     1)) ;
  781. * 2.9
  782. * 'MESSAGE' ('CHAINE' (formar -2.9    1)) ;
  783. * -2.9
  784. * 'MESSAGE' ('CHAINE' (formar 0       1)) ;
  785. * 0
  786. * 'MESSAGE' ('CHAINE' (formar 0       1)) ;
  787. * 0
  788. * 'MESSAGE' ('CHAINE' (formar 2.9e5   0)) ;
  789. * 3E5
  790. * 'MESSAGE' ('CHAINE' (formar -2.9e5  0)) ;
  791. * -3E5
  792. * 'MESSAGE' ('CHAINE' (formar 2.9e-5  0)) ;
  793. * 3E-5
  794. * 'MESSAGE' ('CHAINE' (formar -2.9e-5 0)) ;
  795. * -3E-5
  796. * 'MESSAGE' ('CHAINE' (formar 2.9     0)) ;
  797. * 3
  798. * 'MESSAGE' ('CHAINE' (formar -2.9    0)) ;
  799. * -3
  800. * 'MESSAGE' ('CHAINE' (formar 0       0)) ;
  801. * 0
  802. * 'MESSAGE' ('CHAINE' (formar 0       0)) ;
  803. * 0
  804. *               
  805. *
  806. *
  807. * LANGAGE     : GIBIANE-CAST3M
  808. * AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  809. *               mél : gounand@semt2.smts.cea.fr
  810. **********************************************************************
  811. * VERSION    : v1, 18/02/2003, version initiale
  812. * HISTORIQUE : v1, 18/02/2003, création
  813. * HISTORIQUE :
  814. * HISTORIQUE :
  815. ************************************************************************
  816. * Prière de PRENDRE LE TEMPS de compléter les commentaires
  817. * en cas de modification de ce sous-programme afin de faciliter
  818. * la maintenance !
  819. ************************************************************************
  820. *
  821. *
  822. 'DEBPROC' FORMAR ;
  823. 'ARGUMENT' fl*'FLOTTANT' ;
  824. 'ARGUMENT' vir/'ENTIER  ' ;
  825. 'SI' ('NON' ('EXISTE' vir)) ;
  826.    vir = 1 ;
  827. 'SINON' ;   
  828.    'SI' ('<' vir 0) ;
  829.       'ERREUR' 'fournir un entier positif' ;
  830.    'FINSI' ;
  831. 'FINSI' ;         
  832. 'SI' ('<' ('ABS' fl) 10.D-100) ;
  833.    chfl = 'CHAINE' '0' ;
  834. 'SINON' ;
  835. *! sans le 1.D-10, ca ne fonctionne pas
  836. *! qd on entre pile poil une puissance de 10
  837.    lfl = LOG10 ('ABS' fl) ;
  838. *   lfl = '+' (LOG10 ('ABS' fl)) 1.D-10  ;
  839.    slfl = 'SIGNE' ('ENTIER' lfl) ;
  840.    'SI' ('EGA' slfl 1) ;
  841.       elfl = 'ENTIER' lfl ;
  842.    'SINON' ;
  843.       elfl = '-' ('ENTIER' lfl) 1 ;
  844.    'FINSI' ;
  845.    man = '/' fl ('**' 10.D0 elfl) ;
  846. *
  847. * Une verrue pour des histoires de précision...
  848. *   
  849.    'SI' ('EGA' man 10.D0 ('**' 10.D0 ('*' vir -1.D0))) ;
  850.       man  = '/' man 10.D0 ;
  851.       elfl = '+' elfl 1 ;
  852.    'FINSI' ;
  853. *   
  854.    sman = 'SIGNE' man ;
  855.    'SI' ('EGA' sman 1) ;   
  856.       fman = 'CHAINE' '(F' ('+' vir 2) '.0' vir ')' ;
  857.    'SINON' ;
  858.       fman = 'CHAINE' '(F' ('+' vir 3) '.0' vir ')' ;
  859.    'FINSI' ;
  860.    'SI' ('NEG' vir 0) ;   
  861.       'SI' ('NEG' elfl 0) ;
  862.          chfl = 'CHAINE' 'FORMAT' fman man 'E' elfl ;
  863.       'SINON' ;
  864.          chfl = 'CHAINE' 'FORMAT' fman man ;
  865.       'FINSI' ;
  866.    'SINON' ;
  867.       man2 = 'ENTIER' ('+' man ('*' 0.5D0 sman)) ;
  868.       'SI' ('NEG' elfl 0) ;
  869.          chfl = 'CHAINE' man2 'E' elfl ;
  870.       'SINON' ;
  871.          chfl = 'CHAINE' man2 ;
  872.       'FINSI' ;
  873.    'FINSI' ;            
  874. 'FINSI' ;   
  875. 'RESPRO' chfl ;
  876. *
  877. * End of procedure file FORMAR
  878. *
  879. 'FINPROC' ;
  880. *ENDPROCEDUR formar
  881. *BEGINPROCEDUR log10
  882. ************************************************************************
  883. * NOM         : LOG10
  884. * DESCRIPTION : Log_10
  885. *
  886. *
  887. *
  888. * LANGAGE     : GIBIANE-CAST3M
  889. * AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  890. *               mél : gounand@semt2.smts.cea.fr
  891. **********************************************************************
  892. * VERSION    : v1, 18/02/2003, version initiale
  893. * HISTORIQUE : v1, 18/02/2003, création
  894. * HISTORIQUE :
  895. * HISTORIQUE :
  896. ************************************************************************
  897. * Prière de PRENDRE LE TEMPS de compléter les commentaires
  898. * en cas de modification de ce sous-programme afin de faciliter
  899. * la maintenance !
  900. ************************************************************************
  901. *
  902. *
  903. 'DEBPROC' LOG10 ;
  904. 'REPETER' bouc ;
  905.    ok = FAUX ;
  906.    'ARGUMENT' fl/'FLOTTANT' ;
  907.    'ARGUMENT' lr/'LISTREEL' ;
  908.    'ARGUMENT' cp/'CHPOINT ' ;
  909.    'ARGUMENT' cm/'MCHAML  ' ;
  910.    'SI' ('EXISTE' fl) ;
  911.       ok = VRAI ;
  912.       'RESPRO' ('/' ('LOG' fl) ('LOG' 10.D0)) ;
  913.    'FINSI' ;
  914.    'SI' ('EXISTE' lr) ;
  915.       ok = VRAI ;
  916.       'RESPRO' ('/' ('LOG' lr) ('LOG' 10.D0)) ;
  917.    'FINSI' ;
  918.    'SI' ('EXISTE' cp) ;
  919.       ok = VRAI ;
  920.       'RESPRO' ('/' ('LOG' cp) ('LOG' 10.D0)) ;
  921.    'FINSI' ;
  922.    'SI' ('EXISTE' cm) ;
  923.       ok = VRAI ;
  924.       'RESPRO' ('/' ('LOG' cm) ('LOG' 10.D0)) ;
  925.    'FINSI' ;
  926.    'SI' ('NON' ok) ;
  927.       'QUITTER' bouc ;
  928.    'FINSI' ;
  929. 'FIN' bouc ;         
  930. *
  931. * End of procedure file LOG10
  932. *
  933. 'FINPROC' ;
  934. *ENDPROCEDUR log10
  935. *BEGINPROCEDUR dessevol
  936. ************************************************************************
  937. * NOM         : DESSEVOL
  938. * DESCRIPTION : Dessine des évolutions : choisit automatiquement
  939. *               les options, marqueurs, couleurs...
  940. *
  941. *
  942. * LANGAGE     : GIBIANE-CAST3M
  943. * AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  944. *               mél : gounand@semt2.smts.cea.fr
  945. **********************************************************************
  946. * VERSION    : v1, 16/11/2004, version initiale
  947. * HISTORIQUE : v1, 16/11/2004, création
  948. * HISTORIQUE :
  949. * HISTORIQUE :
  950. ************************************************************************
  951. * Prière de PRENDRE LE TEMPS de compléter les commentaires
  952. * en cas de modification de ce sous-programme afin de faciliter
  953. * la maintenance !
  954. ************************************************************************
  955. *
  956. *
  957. 'DEBPROC' DESSEVOL ;
  958. 'ARGUMENT' evtot*'EVOLUTION' ;
  959. 'ARGUMENT' tabt*'TABLE' ;
  960. 'ARGUMENT' tit*'MOT' ;
  961. 'ARGUMENT' tix*'MOT' ;
  962. 'ARGUMENT' tiy*'MOT' ;
  963. 'ARGUMENT' lnclk/'LOGIQUE' ;
  964. 'ARGUMENT' nb/'ENTIER' ;
  965. *
  966. 'SI' ('NON' ('EXISTE' lnclk)) ;
  967.    lnclk = FAUX ;
  968. 'FINSI' ;   
  969. *
  970. 'SI' ('NON' ('EXISTE' nb)) ;
  971.    nb = 3 ;
  972. 'FINSI' ;   
  973. *
  974. nt  = 'DIME' tabt ;
  975. nev = 'DIME' evtot ;
  976. *
  977. * Attention, dans evtot, il y a une évolution avec des noms de points ?
  978. *
  979. *'SI' ('NEG' nev nt) ;
  980. *   cherr = 'CHAINE' 'Evolution and title table : not same dim.' ;
  981. *   'ERREUR' cherr ;
  982. *'FINSI' ;
  983. *
  984. tev = 'TABLE' ;
  985. tev . 'TITRE' = tabt ;
  986. *
  987. toto = 'TABLE' ;
  988. *
  989. lcoul = 'MOTS' 'TURQ' 'VERT' 'JAUN' 'ROSE' 'ROUG' 'BLEU' ;
  990. lmarq = 'MOTS' 'TRIB' 'TRIA' 'LOSA' 'CARR' 'ETOI' 'PLUS' 'CROI' ;
  991. ltirr = 'MOTS' 'TIRR' 'TIRC' 'TIRL' 'TIRM' ;
  992. *
  993. 'SI' ('EGA' nb 0) ;
  994.    ev2 = evtot ;
  995. 'SINON' ;
  996.    icou = 0 ;
  997.    'REPETER' iev nev ;
  998.       ii = &iev ;
  999.       evi = 'EXTRAIRE' evtot 'COUR' ii ;
  1000.       'SI' ('NEG' ('TYPE' ('EXTRAIRE' evi 'ORDO')) 'LISTMOTS') ;
  1001.          icou = '+' icou 1 ;
  1002.       'FINSI' ;         
  1003. *      ii2 = '/' ('+' ii 1) 2 ;
  1004. *      ci  = EXMOMOD lcoul ii2 ;
  1005. *      ci  = EXMOMOD lcoul ii ;
  1006.       ci  = EXMOMOD lcoul icou ;
  1007.       APPEND toto 'EVOLUTION' ('COULEUR' evi ci) ;
  1008.    'FIN' iev ;
  1009.    ev2 = toto . 'EVOLUTION' ;
  1010. 'FINSI' ;
  1011. *      
  1012. 'REPETER' iev nev ;
  1013.    ii = &iev ;
  1014.    mi = EXMOMOD lmarq ii ;
  1015.    ti = EXMOMOD ltirr ii ;
  1016.    'SI' ('>' nb 2) ;
  1017.       tev . ii = 'CHAINE' 'MARQ ' mi ' ' ti ;
  1018.    'FINSI' ;      
  1019.    'SI' ('>' nb 1) ;
  1020.       tev . ii = 'CHAINE' 'MARQ ' mi ;
  1021.    'FINSI' ;      
  1022. 'FIN' iev ;   
  1023. *
  1024. 'SI' lnclk ;
  1025.    'DESSIN' ev2 'TITR' tit 'TITX' tix 'TITY' tiy 'LEGE' 'MIMA' tev
  1026.    'NCLK' ;
  1027. 'SINON' ;   
  1028.    'DESSIN' ev2 'TITR' tit 'TITX' tix 'TITY' tiy 'LEGE' 'MIMA' tev ;
  1029. 'FINSI' ;   
  1030. *
  1031. * End of procedure file DESSEVOL
  1032. *
  1033. 'FINPROC' ;
  1034. *ENDPROCEDUR dessevol
  1035. *BEGINPROCEDUR append
  1036. ************************************************************************
  1037. * NOM         : APPEND
  1038. * DESCRIPTION : Rajoute :
  1039. *               - un entier à un listentier
  1040. *               - un réel   à un listreel
  1041. *               - un objet (liste, evolution, matrice ou chpoint)
  1042. *                 à un indice de table ('MOT' ou 'ENTIER')
  1043. *                    * si l'indice n'existe pas
  1044. *                    * 'ET' si l'indice existe 
  1045. *                
  1046. *
  1047. *
  1048. * LANGAGE     : GIBIANE-CAST3M
  1049. * AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  1050. *               mél : gounand@semt2.smts.cea.fr
  1051. **********************************************************************
  1052. * VERSION    : v1, 10/09/2004, version initiale
  1053. * HISTORIQUE : v1, 10/09/2004, création
  1054. * HISTORIQUE :
  1055. * HISTORIQUE :
  1056. ************************************************************************
  1057. * Prière de PRENDRE LE TEMPS de compléter les commentaires
  1058. * en cas de modification de ce sous-programme afin de faciliter
  1059. * la maintenance !
  1060. ************************************************************************
  1061. *
  1062. *
  1063. 'DEBPROC' APPEND ;
  1064. 'ARGUMENT' tab/'TABLE' ;
  1065. 'SI' ('EXISTE' tab) ;
  1066.    'ARGUMENT' itab/'MOT' ;
  1067.    'SI' ('NON' ('EXISTE' itab)) ;
  1068.       'ARGUMENT' itab*'ENTIER' ;
  1069.    'FINSI' ;      
  1070.    lobj = FAUX ;
  1071.    'SI' ('NON' lobj) ;
  1072.       'ARGUMENT' lr/'LISTREEL' ;
  1073.       'SI' ('EXISTE' lr) ;
  1074.          obj = lr ; lobj = VRAI ;
  1075.       'FINSI' ;
  1076.    'FINSI' ;               
  1077.    'SI' ('NON' lobj) ;
  1078.       'ARGUMENT' le/'LISTENTI' ;
  1079.       'SI' ('EXISTE' le) ;
  1080.          obj = le ; lobj = VRAI ;
  1081.       'FINSI' ;
  1082.    'FINSI' ;               
  1083.    'SI' ('NON' lobj) ;
  1084.       'ARGUMENT' lev/'EVOLUTION' ;
  1085.       'SI' ('EXISTE' lev) ;
  1086.          obj = lev ; lobj = VRAI ;
  1087.       'FINSI' ;
  1088.    'FINSI' ;               
  1089.    'SI' ('NON' lobj) ;
  1090.       'ARGUMENT' lm/'MAILLAGE' ;
  1091.       'SI' ('EXISTE' lm) ;
  1092.          obj = lm ; lobj = VRAI ;
  1093.       'FINSI' ;
  1094.    'FINSI' ;               
  1095.    'SI' ('NON' lobj) ;
  1096.       'ARGUMENT' chpo/'CHPOINT' ;
  1097.       'SI' ('EXISTE' chpo) ;
  1098.          obj = chpo ; lobj = VRAI ;
  1099.       'FINSI' ;
  1100.    'FINSI' ;               
  1101.    'SI' ('NON' lobj) ;
  1102.       'ARGUMENT' rig/'RIGIDITE' ;
  1103.       'SI' ('EXISTE' rig) ;
  1104.          obj = rig ; lobj = VRAI ;
  1105.       'FINSI' ;
  1106.    'FINSI' ;               
  1107.    'SI' ('NON' lobj) ;
  1108.       'ARGUMENT' matk/'MATRIK' ;
  1109.       'SI' ('EXISTE' matk) ;
  1110.          obj = matk ; lobj = VRAI ;
  1111.       'FINSI' ;
  1112.    'FINSI' ;
  1113.    'SI' ('NON' lobj) ;
  1114.       cherr = 'CHAINE'
  1115.        'Il faut fournir un objet liste, evolution, matrice ou chpoint.'
  1116.         ; 
  1117.       'ERREUR' cherr ;
  1118.    'FINSI' ;      
  1119.    'SI' ('EXISTE' tab itab) ;
  1120.       'SI' ('EGA' ('TYPE' obj) 'CHPOINT') ;
  1121.          tab . itab = '+' (tab . itab) obj ;
  1122.       'SINON' ;         
  1123.          tab . itab = 'ET' (tab . itab) obj ;
  1124.       'FINSI' ;         
  1125.    'SINON' ;
  1126.       tab . itab = obj ;
  1127.    'FINSI' ;
  1128.    'RESPRO' tab ;      
  1129. 'FINSI' ;   
  1130. 'ARGUMENT' lenti/'LISTENTI' ;
  1131. 'ARGUMENT' lreel/'LISTREEL' ;
  1132. 'SI' ('EXISTE' lenti) ;
  1133.    'ARGUMENT' enti*'ENTIER' ;
  1134.    lenti = 'ET' lenti ('LECT' enti) ;
  1135.    'RESPRO' lenti ;
  1136. 'FINSI' ;   
  1137. 'SI' ('EXISTE' lreel) ;
  1138.    'ARGUMENT' reel*'FLOTTANT' ;
  1139.    lreel = 'ET' lreel ('PROG' reel) ;
  1140.    'RESPRO' lreel ;
  1141. 'FINSI' ;   
  1142. *
  1143. * End of procedure file APPEND
  1144. *
  1145. 'FINPROC' ;
  1146. *ENDPROCEDUR append
  1147. *BEGINPROCEDUR exmomod
  1148. ************************************************************************
  1149. * NOM         : EXMOMOD
  1150. * DESCRIPTION : Extraction d'un mot d'un listmots
  1151. *
  1152. *
  1153. *
  1154. * LANGAGE     : GIBIANE-CAST3M
  1155. * AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  1156. *               mél : gounand@semt2.smts.cea.fr
  1157. **********************************************************************
  1158. * VERSION    : v1, 23/06/2003, version initiale
  1159. * HISTORIQUE : v1, 23/06/2003, création
  1160. * HISTORIQUE :
  1161. * HISTORIQUE :
  1162. ************************************************************************
  1163. * Prière de PRENDRE LE TEMPS de compléter les commentaires
  1164. * en cas de modification de ce sous-programme afin de faciliter
  1165. * la maintenance !
  1166. ************************************************************************
  1167. *
  1168. *
  1169. 'DEBPROC' EXMOMOD ;
  1170. 'ARGUMENT' lm*'LISTMOTS' i*'ENTIER' ;
  1171. j = 'DIME' lm ;
  1172. k = '+' (MODULO ('-' i 1) j) 1 ;
  1173. lemot = 'EXTRAIRE' lm k ;
  1174. * Usage de l'opérateur text pour éviter que lemot
  1175. * ne soit interprété comme un opérateur
  1176. 'RESPRO' 'TEXTE' lemot ;
  1177. *
  1178. * End of procedure file EXMOMOD
  1179. *
  1180. 'FINPROC' ;
  1181. *ENDPROCEDUR exmomod
  1182. *BEGINPROCEDUR modulo
  1183. ************************************************************************
  1184. * NOM         : MODULO
  1185. * DESCRIPTION : Calcule un entier modulo un autre...
  1186. *
  1187. *
  1188. *
  1189. * LANGAGE     : GIBIANE-CAST3M
  1190. * AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  1191. *               mél : gounand@semt2.smts.cea.fr
  1192. **********************************************************************
  1193. * VERSION    : v1, 15/10/2002, version initiale
  1194. * HISTORIQUE : v1, 15/10/2002, création
  1195. * HISTORIQUE :
  1196. * HISTORIQUE :
  1197. ************************************************************************
  1198. * Prière de PRENDRE LE TEMPS de compléter les commentaires
  1199. * en cas de modification de ce sous-programme afin de faciliter
  1200. * la maintenance !
  1201. ************************************************************************
  1202. *
  1203. *
  1204. 'DEBPROC' MODULO ;
  1205. 'ARGUMENT' i*'ENTIER' j*'ENTIER' ;
  1206. 'SI' ('EGA' j 0) ;
  1207.   'MESSAGE' 'Impossible de faire modulo 0' ;
  1208.   'ERREUR' 5 ;
  1209. 'SINON' ;
  1210.   k=i '/' j ;
  1211.   mod=i '-' ( k '*'j ) ;
  1212.   'RESPRO' mod ;
  1213. 'FINSI' ;
  1214. *
  1215. * End of procedure file MODULO
  1216. *
  1217. 'FINPROC' ;
  1218. *ENDPROCEDUR modulo
  1219. *
  1220. 'DEBPROC' SOLEX ;
  1221. 'ARGUMENT' nu*'FLOTTANT' ;
  1222. 'ARGUMENT' mesh*'MAILLAGE' ;
  1223. *
  1224.    ct = '/' 1.D0 nu ;
  1225.    deno = '-' 1. ('EXP' ct) ;
  1226.    nume = 'EXP' ('*' ('COORDONNEE' 1 mesh) ct) ;
  1227.    nume = '*' nume -1. ;
  1228.    nume = '+' nume 1. ;
  1229.    vfon = '*' nume ('/' 1. deno) ;
  1230. *   
  1231. 'RESPRO' vfon ;
  1232. 'FINPROC' ;
  1233. *
  1234. 'DEBPROC' SOLAPP ;
  1235. 'ARGUMENT' _mt*'MAILLAGE' ;
  1236. 'ARGUMENT' gau*'MAILLAGE' ;
  1237. 'ARGUMENT' dro*'MAILLAGE' ;
  1238. 'ARGUMENT' disc*'MOT     ' ; 
  1239. 'ARGUMENT' nu*'FLOTTANT' ;
  1240. 'ARGUMENT' pec*'FLOTTANT' ;
  1241. 'ARGUMENT' imeth*'ENTIER' ;
  1242. *
  1243. lux = 'MOTS' 'UX' 'UY' ; lval = 'PROG' 1. 0. ;
  1244. cvit = 'MANUEL' 'CHPO' _mt lux lval ;
  1245. rho = 1.D0 ;
  1246. m1 = GUGRAD _mt disc 'T' disc 'Q' disc cvit lux disc 1. 'GAU7' ;
  1247. m2 = GLAPN  _mt disc 'T' disc 'Q' disc nu 'GAU7' ;
  1248. m2 = '*' m2 -1.D0 ;
  1249. 'SI' ('EGA' imeth 0) ;
  1250.    musmax = 0. ;
  1251. 'SINON' ;   
  1252.    m3 = GDECENT imeth _mt disc 'T' disc 'Q' disc cvit lux disc
  1253.                  rho nu pec 'GAU7' ;
  1254.    musmax = 'MAXIMUM' (GMUSTAB imeth
  1255.                _mt disc cvit lux disc rho nu pec 'GAU7') ;
  1256. 'FINSI' ;               
  1257. mblo = 'BLOQUE' 'T' (gau 'ET' dro) ;
  1258. fblo = 'DEPIMPOSE' mblo ('+' ('MANUEL' 'CHPO' gau 1 'T' 0.D0)
  1259.                              ('MANUEL' 'CHPO' dro 1 'T' 1.D0)) ;
  1260. 'SI' ('EGA' imeth 0) ;                             
  1261.    mtot = m1 'ET' m2 'ET' mblo ;
  1262. 'SINON' ;   
  1263.    mtot = m1 'ET' m2 'ET' m3 'ET' mblo ;
  1264. 'FINSI' ;   
  1265. ftot = fblo ;
  1266. sola = 'RESO' mtot ftot ;
  1267. 'RESPRO' sola musmax ;
  1268. 'FINPROC' ;                             
  1269. *
  1270. *
  1271. ************************************************************************
  1272. * NOM         : NLIN_DECENT2D
  1273. * DESCRIPTION : Cas-test du décentrement :
  1274. *               couche limite 1D en 2D, on regarde les
  1275. *               ordres de convergences           
  1276. *
  1277. *
  1278. * LANGAGE     : GIBIANE-CAST3M
  1279. * AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  1280. *               mél : gounand@semt2.smts.cea.fr
  1281. **********************************************************************
  1282. * VERSION    : v1, 31/08/2006, version initiale
  1283. * HISTORIQUE : v1, 31/08/2006, création
  1284. * HISTORIQUE :
  1285. * HISTORIQUE :
  1286. ************************************************************************
  1287. * Prière de PRENDRE LE TEMPS de compléter les commentaires
  1288. * en cas de modification de ce sous-programme afin de faciliter
  1289. * la maintenance !
  1290. ************************************************************************
  1291. *
  1292. *
  1293. interact= FAUX ;
  1294. graph   = FAUX ;
  1295. verbose = FAUX ;
  1296. 'SI' ('NON' interact) ;
  1297.   'OPTION' 'TRAC' 'PS' ;
  1298. 'SINON' ;
  1299.   'OPTION' 'TRAC' 'X' ;
  1300. *  'OPTION' 'TRAC' 'OPEN' ;
  1301. 'FINSI' ;  
  1302. *
  1303. 'OPTION' 'DIME' 2 ;
  1304. lnx = 'LECT' 4 10 20 30 ;
  1305. dlnx = 'DIME' lnx ;
  1306. nxfin = '*' ('EXTRAIRE' lnx ('DIME' lnx)) 2 ;
  1307. nu = 0.1 ; 
  1308. Pec = 2.D0 ;
  1309. *
  1310. moscal1 = mots 'SCAL';
  1311. *
  1312. idx = 0 ;
  1313. *
  1314. * Table des ordres espérés pour chaque méthode
  1315. * 'LINE' : éléments linéaires
  1316. * 'QUAD' : éléments quadratiques
  1317. * imeth = 0 : Pas de décentrement
  1318. * imeth = 1 : Décentrement full upwind
  1319. * imeth = 2 : Décentrement type SUPG avec une fonction de coupure
  1320. *             "critical approximation"
  1321. * imeth = 3 : Décentrement type SUPG avec une fonction de coupure
  1322. *             "doubly asymptotic"
  1323.  
  1324. tabordv = 'TABLE' ;
  1325. tabordv . 'LINE' = 'TABLE' ;
  1326. tabordv . 'LINE' . 0 = 2. ;
  1327. tabordv . 'LINE' . 1 = 1. ;
  1328. tabordv . 'LINE' . 2 = 2. ;
  1329. tabordv . 'LINE' . 3 = 2. ;
  1330. tabordv . 'QUAD' = 'TABLE' ;
  1331. tabordv . 'QUAD' . 0 = 3. ;
  1332. tabordv . 'QUAD' . 1 = 1. ;
  1333. tabordv . 'QUAD' . 2 = 3. ;
  1334. tabordv . 'QUAD' . 3 = 2. ;
  1335. *
  1336. tabordc = 'TABLE' ;
  1337. tabordc . 'LINE' = 'TABLE' ;
  1338. tabordc . 'QUAD' = 'TABLE' ;
  1339. tabev = 'TABLE' ;
  1340. tabt  = 'TABLE' ;
  1341. 'REPETER' iiquad 2 ;
  1342.    iquad = &iiquad ;
  1343.    'SI' ('EGA' iquad 1) ;
  1344.       disc = 'LINE' ;
  1345.    'SINON' ;      
  1346.       disc = 'QUAD' ;
  1347.    'FINSI' ;
  1348.    'REPETER' iimeth 4 ;
  1349.       imeth = '-' &iimeth 1 ;
  1350.       'SI' verbose ;
  1351.          tit   = 'CHAINE' ' disc = ' disc ' imeth=' imeth ;
  1352.          'MESSAGE' tit ;
  1353.       'FINSI' ;      
  1354.       lh = 'PROG' ;
  1355.       lerr = 'PROG' ;
  1356.       'REPETER' iinx ('DIME' lnx) ;
  1357.          inx = &iinx ;
  1358.          nx = 'EXTRAIRE' lnx inx ;
  1359.          'SI' ('EGA' iquad 1) ;
  1360.             'OPTION' 'ELEM' 'TRI3' ;
  1361.             nx2    = nx ;
  1362.             nxfin2 = nxfin ;
  1363.          'SINON' ;            
  1364.             'OPTION' 'ELEM' 'TRI6' ;
  1365.             nx2 = '/' nx 2 ;
  1366.             nxfin2 = '/' nxfin 2 ;
  1367.          'FINSI' ;            
  1368.          vh = '/' 1. ('FLOTTANT' nx2) ;
  1369.          lh = 'ET' lh ('PROG' vh) ;
  1370.       p1 = 0. 0. ; p2 = 1. 0. ; p3 = 1. 1. ; p4 = 0. 1. ;
  1371.       bas = 'DROIT' nx2 p1 p2 ; dro = 'DROIT' nx2 p2 p3 ;
  1372.       hau = 'DROIT' nx2 p3 p4 ; gau = 'DROIT' nx2 p4 p1 ;
  1373.       cmt = bas 'ET' dro 'ET' hau 'ET' gau ;
  1374.       mt = 'SURFACE' cmt 'PLAN' ;
  1375.       _mt = 'CHANGER' mt 'QUAF' ;
  1376.       basfin = 'DROIT' nxfin2 p1 p2 ; drofin = 'DROIT' nxfin2 p2 p3 ;
  1377.       haufin = 'DROIT' nxfin2 p3 p4 ; gaufin = 'DROIT' nxfin2 p4 p1 ;
  1378.       cmtfin = basfin 'ET' drofin 'ET' haufin 'ET' gaufin ;
  1379.       mtfin = 'SURFACE' cmtfin 'PLAN' ;
  1380.       _mtfin = 'CHANGER' mtfin 'QUAF' ;
  1381. *         p0 = 'POIN' 0. ; p1 = 'POIN' 1. ;
  1382. *         mt = 'DROIT' nx2 p0 p1 ;
  1383. *         _mt ='CHANGER' mt 'QUAF' ;
  1384. *         mtfin  = 'DROIT' nxfin2 p0 p1 ;
  1385. *         _mtfin ='CHANGER' mtfin 'QUAF' ;
  1386.          mas = GMASS _mtfin disc 'SCAL' disc 'SCAL' disc 1. ;         
  1387.          sol = SOLEX nu mtfin ;
  1388.          sola musmax = SOLAPP _mt gau dro disc nu pec imeth ;
  1389.          sola = 'EXCO' 'T' sola ;
  1390.          sola2 = proi mtfin ('CHANGER' 'CHAM' sola mt) 1.D-9 ;
  1391.          err = '-' sol sola2 ;
  1392. *         verr = '**' ('XTMX' err mas) 0.5D0 ;
  1393.          verr = '**' (XTY err (mas * err) moscal1 moscal1) 0.5D0 ;
  1394.          lerr = 'ET' lerr ('PROG' verr) ;
  1395.          'SI' verbose ;
  1396.             ch = 'CHAINE' 'h=' (formar vh 2)
  1397.                           ' => err=' (formar verr 2) ;
  1398.             'MESSAGE' ch ;
  1399.          'FINSI' ;            
  1400.          'SI' (graph 'ET' verbose) ;
  1401.             montagne sola sola mt 'TITR' ('CHAINE' tit ' ' ch) ;
  1402.          'FINSI' ;
  1403.          'OPTION' 'DIME' 2 ;
  1404.       'FIN' iinx ;      
  1405. * Calcul des ordres
  1406.       lh lerr = LOG10 lh lerr ;
  1407.       everr = 'EVOL' 'VERT' 'MANU' lh lerr ;
  1408.       di = 'DIME' lh ;
  1409.       lh2   = 'EXTRAIRE' lh   ('LECT' ('-' di 1) di) ;
  1410.       lerr2 = 'EXTRAIRE' lerr ('LECT' ('-' di 1) di) ;
  1411.       everr2 = 'EVOL' 'MANU' lh2 lerr2 ;
  1412.       cpl dummy = @POMI everr2 1 'IDEM' ;
  1413.       ord = cpl . 1 ;
  1414.       idx = '+' idx 1 ;
  1415.       tabordc . disc . imeth = ord ;
  1416.       tabev . idx = everr ;
  1417.       'SI' verbose ;
  1418.          tito = 'CHAINE' disc imeth
  1419.                         ';ord=' (formar ord 1) ;
  1420.          'MESSAGE' tito ;
  1421.       'FINSI' ;         
  1422.       tabt . idx = tito ;               
  1423.    'FIN' iimeth ;   
  1424. 'FIN' iiquad ;
  1425. tit  = 'CHAINE' 'Log Err. (Log h)' ;
  1426. tix  = 'CHAINE' 'Log h' ;
  1427. tiy  = 'CHAINE' 'Log Err' ;
  1428. evt = @STBL tabev ;
  1429. 'SI' graph ;
  1430.    dessevol evt tabt tit tix tiy FAUX ;
  1431. 'FINSI' ;
  1432. 'SAUTER' 2 'LIGNE' ;
  1433. * Tests des ordres
  1434. ok = VRAI ;
  1435. 'REPETER' iiquad 2 ;
  1436.    iquad = &iiquad ;
  1437.    'SI' ('EGA' iquad 1) ;
  1438.       disc = 'LINE' ;
  1439.    'SINON' ;      
  1440.       disc = 'QUAD' ;
  1441.    'FINSI' ;
  1442.    'REPETER' iimeth 4 ;
  1443.       imeth = '-' &iimeth 1 ;
  1444.       'SI' ('EGA' imeth 0) ;
  1445.          met = 'CHAINE' 'rien            ' ; 'FINSI' ;
  1446.       'SI' ('EGA' imeth 1) ;
  1447.          met = 'CHAINE' 'full upwind     ' ; 'FINSI' ;
  1448.       'SI' ('EGA' imeth 2) ;
  1449.          met = 'CHAINE' 'supg crit approx' ; 'FINSI' ;
  1450.       'SI' ('EGA' imeth 3) ;
  1451.          met = 'CHAINE' 'supg doub asymp ' ; 'FINSI' ;
  1452.       valvoul = tabordv . disc . imeth ;
  1453.       valobt  = tabordc . disc . imeth ;
  1454.       test = ('EGA' valobt valvoul ('*' valvoul 0.2D0)) ;
  1455.       ok = ok 'ET' test ;
  1456.       mess1 = 'CHAINE' 'Eléments ' disc ' Décentrement ' met
  1457.               ' ordre = ' valobt ;
  1458.       'MESSAGE' mess1 ;              
  1459.       'SI' ('NON' test) ;
  1460.          'MESSAGE' ('CHAINE' '!!!! On aurait voulu avoir :' valvoul) ;
  1461.       'FINSI' ;
  1462.    'FIN' iimeth ;
  1463. 'FIN' iiquad ;                  
  1464. *
  1465. 'SI' ('NON' ok) ;
  1466.    'MESSAGE' ('CHAINE' 'Il y a eu des erreurs') ;
  1467.    'ERREUR' 5 ;
  1468. 'SINON' ;
  1469.    'MESSAGE' ('CHAINE' 'Tout sest bien passe !') ;
  1470. 'FINSI' ;
  1471. *
  1472. 'SI' interact ;
  1473.    'OPTION' 'ECHO' 1 ;
  1474.    'OPTION' 'DONN' 5 ;
  1475. 'FINSI' ;
  1476. *
  1477. * End of dgibi file NLIN_DECENT2D
  1478. *
  1479. 'FIN' ;
  1480.  
  1481.  
  1482.  
  1483.  
  1484.  
  1485.  
  1486.  
  1487.  
  1488.  
  1489.  

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