* fichier : gravite.dgibi 'OPTI' 'ECHO' 0 ; ************************************************************************ * NOM : GRAVITE * DESCRIPTION : Cas-test gravité servant à tester la méthode * méthode de projection incrémentale. * On doit trouver vitesse nulle et pression * hydrostatique * * LANGAGE : GIBIANE-CAST3M * AUTEUR : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF) * mél : gounand@semt2.smts.cea.fr ********************************************************************** * VERSION : v1, 21/11/2007, version initiale * HISTORIQUE : v1, 21/11/2007, création * HISTORIQUE : * HISTORIQUE : ************************************************************************ * Prière de PRENDRE LE TEMPS de compléter les commentaires * en cas de modification de ce sous-programme afin de faciliter * la maintenance ! ************************************************************************ 'SAUTER' 2 'LIGNE' ; 'MESSAGE' ' Execution de gravite.dgibi' ; 'SAUTER' 2 'LIGNE' ; * interact = FAUX ; graph = FAUX ; * * * Procédure pour afficher la résultante d'un chpoint et de sa valeur * absolue (utile pour debugger sans lister le CHPOINT en question) * *'DEBPROC' HCHPO ; *'ARGUMENT' chp*'CHPOINT' ; *'ARGUMENT' tit*'MOT' ; ** *achp = 'ABS' chp ; *nachp = 'MAXIMUM' achp ; **'SI' ('' vv 1.D-10) ; errud = '/' errud vv ; 'FINSI' ; tit = 'CHAINE' (ttit . &imet) '/non prec' ; 'SINON' ; errud = 'MAXIMUM' und 'ABS' ; tit = 'CHAINE' (ttit . &imet) ; 'FINSI' ; chmes = 'CHAINE' tit ' max err. vitesse =' errud ; 'MESSAGE' chmes ; erruda = 'EXTRAIRE' lerru &imet ; tst = ('<' errud erruda) ; 'SI' ('NON' tst) ; chmes = 'CHAINE' '!!! Erreur, pas inferieure a ' erruda ; 'MESSAGE' chmes ; 'FINSI' ; lok = lok 'ET' tst ; * pnd = tresp . &imet ; 'SI' ('OU' ('EGA' &imet 3) ('EGA' &imet 5)) ; pndd = tresp . ('-' &imet 1) ; errpd = 'MAXIMUM' ('-' pnd pndd) 'ABS' ; vv = 'MAXIMUM' pndd 'ABS' ; 'SI' ('>' vv 1.D-10) ; errpd = '/' errpd vv ; 'FINSI' ; tit = 'CHAINE' (ttit . &imet) '/non prec' ; 'SINON' ; errpd = 'MAXIMUM' pnd 'ABS' ; tit = 'CHAINE' (ttit . &imet) ; 'FINSI' ; chmes = 'CHAINE' tit ' max err. pression =' errpd ; 'MESSAGE' chmes ; errpda = 'EXTRAIRE' lerrp &imet ; tst = ('<' errpd errpda) ; 'SI' ('NON' tst) ; chmes = 'CHAINE' '!!! Erreur, pas inferieure a ' errpda ; 'MESSAGE' chmes ; 'FINSI' ; lok = lok 'ET' tst ; * dnd = tresd . &imet ; errdd = 'MAXIMUM' dnd 'ABS' ; chmes = 'CHAINE' (ttit . &imet) ' max err. divergence =' errdd ; 'MESSAGE' chmes ; errdda = 'EXTRAIRE' lerrd &imet ; tst = ('<' errdd errdda) ; 'SI' ('NON' tst) ; chmes = 'CHAINE' '!!! Erreur, pas inferieure a ' errdda ; 'MESSAGE' chmes ; 'FINSI' ; lok = lok 'ET' tst ; 'FIN' imet ; * 'FIN' iaxi ; 'FIN' imesh ; * 'SAUTER' 2 'LIGNE' ; 'SI' lok ; 'MESSAGE' 'Tout sest bien passe' ; 'SINON' ; 'MESSAGE' 'Il y a eu des erreurs' ; 'FINSI' ; 'SAUTER' 2 'LIGNE' ; * 'SI' interact ; 'OPTION' 'DONN' 5 'ECHO' 1 ; 'FINSI' ; 'SI' ('NON' lok) ; 'ERREUR' 5 ; 'FINSI' ; * * End of dgibi file GRAVITE * 'FIN' ;