Download bgmo_bcn.dgibi

Back to the list

* fichier :  bgmo_bcn.dgibi
****************************************************
** TESTING FILE FOR THE OPERATOR BGMO EVALUATING   *        
** THE FUNCIONS INVOLVED IN THE BRUNO GERARD MODEL *
****************************************************
GRAPH='N';                                                                     

nn1 = 100;
tkk = PROG;
cokk = PROG;
t=0.;
repeter kk1 nn1;
  con1 = BGMO 'COND' t 20. .8e-3;
  tkk = tkk INSE &kk1 t;
  cokk = cokk INSE &kk1 con1;
  t= t + (.5 / 50.);
fin kk1;
num1 = EXTR cokk 50 ;
flag1 = ABS(num1 - 1.6577);
evol1 = EVOL MANU 'temps' tkk 'COND' cokk;
SI (NEG GRAPH 'N');
dess evol1;
FINSI;

nn1 = 100;
tkk = PROG;
cokk = PROG;
t=0.;
repeter kk1 nn1;
  con1 = BGMO 'DCON' t 20. .8e-3;
  tkk = tkk INSE &kk1 t;
  cokk = cokk INSE &kk1 con1;
  t= t + (.5 / 50.);
fin kk1;
num1 = EXTR cokk 50 ;
flag2 =  ABS(num1 + 2.4476);
evol1 = EVOL MANU 'temps' tkk 'COND' cokk;
SI (NEG GRAPH 'N');
dess evol1;
FINSI;

nn1 = 100;
tkk = PROG;
cokk = PROG;
t=0.;
repeter kk1 nn1;
  con1 = BGMO 'CAPA' t 20. 700.;
  tkk = tkk INSE &kk1 t;
  cokk = cokk INSE &kk1 con1;
  t= t + (.5 / 50.);
fin kk1;
num1 = EXTR cokk 50 ;
flag3 = ABS(num1 - 0.33457);
evol1 = EVOL MANU 'temps' tkk 'COND' cokk;
SI (NEG GRAPH 'N');
dess evol1;
FINSI;

nn1 = 100;
tkk = PROG;
cokk = PROG;
t=0.;
repeter kk1 nn1;
  con1 = BGMO 'DCAP' t 20. 700.;
  tkk = tkk INSE &kk1 t;
  cokk = cokk INSE &kk1 con1;
  t= t + (.5 / 50.);
fin kk1;
num1 = EXTR cokk 50 ;
flag4 =  ABS(num1 + 0.65976);
evol1 = EVOL MANU 'temps' tkk 'COND' cokk;
SI (NEG GRAPH 'N');
dess evol1;
FINSI;

flag5 = (flag1+flag2+flag3+flag4);
SI (flag5 < 1.e-4);
ERRE 0;
SINON;
ERRE 5;
FINSI;


FIN;
 

© Cast3M 2003 - All rights reserved.
Disclaimer