Esope source files

   1 : C FOBESS    SOURCE    CHAT      05/01/13    00:05:34     5004           
   2 :       SUBROUTINE FOBESS                                                         
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
   5 : C                                                                               
   6 : C  OPERATEUR FONC                                                               
   7 : C  CAS 'BESSEL'                                                                 
   8 : C  CALCUL DES FONCTIONS DE BESSEL J0 J1 Y0 Y1                                   
   9 : C                                                                               
  10 : C E. DE LANGRE                                                                  
  11 : C                                                                               
  12 :       CHARACTER*40  CHA1                                                        
  13 : C                                                                               
  14 : -INC CCOPTIO                                                                    
  15 : -INC SMLREEL                                                                    
  16 : C                                                                               
  17 : C                                                                               
  18 : C                                                                               
  18 :       CALL LIRCHA (CHA1,1,IRETOU)                                               
  19 :       IF (IERR.NE.0) RETURN                                                     
  19 :       CALL LIROBJ ('LISTREEL',MLREEL,1,IRETOU)                                  
  20 :       IF (IERR.NE.0) RETURN                                                     
  21 : C                                                                               
  22 : C                                                                               
  23 :       SEGACT MLREEL                                                             
  24 :       JG  = PROG(/1)                                                            
  25 :       SEGINI MLREE1                                                             
  26 : C                                                                               
  27 : C                                                                               
  28 :       IF     (CHA1(1:2) .EQ.'J0') THEN                                          
  28 :           CALL FOBEJ0(PROG,JG,MLREE1.PROG)                                      
  29 :       ELSEIF (CHA1(1:2) .EQ.'J1') THEN                                          
  29 :           CALL FOBEJ1(PROG,JG,MLREE1.PROG)                                      
  30 :       ELSEIF (CHA1(1:2) .EQ.'Y0') THEN                                          
  30 :           CALL FOBEY0(PROG,JG,MLREE1.PROG)                                      
  31 :       ELSEIF (CHA1(1:2) .EQ.'Y1') THEN                                          
  31 :           CALL FOBEY1(PROG,JG,MLREE1.PROG)                                      
  32 :       ENDIF                                                                     
  33 : C                                                                               
  33 :       CALL ECROBJ ('LISTREEL',MLREE1)                                           
  34 :       IF (IERR.EQ.0) RETURN                                                     
  35 : C                                                                               
  36 :       RETURN                                                                    
  37 :       END                                                                       
  38 :