Download repec2.eso

Back to the list

   1 : C REPEC2    SOURCE    KICH      08/06/26    21:16:50     6133           
   2 :       SUBROUTINE REPEC2(ITAB1,imo1,ICHP1)                                                        
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)
   5 : -INC CCOPTIO                                                                    
   6 : -INC CCREEL                                                                     
   7 : -INC SMCHPOI                                                                                                                                    
   8 : -INC SMLREEL                                                                    
   9 :       LOGICAL L0,L1
  10 : *                  1      2      3      4      5                                    
  11 : *      DATA MOOPT/'DEPL','CONT','VITE','ACCE','REAC'/                            
  12 : 
  13 :          ICHP1=0
  13 :          call lirree(val,0,ire1)
  13 :          CALL ACCTAB(ITAB1,'MOT',I0,X0,'TEMPS',L0,IP0,                
  14 :      &                'TABLE',I1,X1,' ',L1,ITABT)
  15 : 
  15 :          call dimen7(itabt,NDIM)
  15 :          if (imo1.eq.1) CALL ACCTAB(ITAB1,'MOT',I0,X0,
  16 :      &'DEPLACEMENTS',L0,IP0,'TABLE',I1,X1,' ',L1,ITAB2)
  16 :           if (imo1.eq.2) CALL ACCTAB(ITAB1,'MOT',I0,X0,
  17 :      &'CONTRAINTES',L0,IP0,'TABLE',I1,X1,' ',L1,ITAB2)
  17 :          if (imo1.eq.3) CALL ACCTAB(ITAB1,'MOT',I0,X0,
  18 :      &'VITESSES',L0,IP0,'TABLE',I1,X1,' ',L1,ITAB2)
  18 :          if (imo1.eq.4) CALL ACCTAB(ITAB1,'MOT',I0,X0,
  19 :      &'ACCELERATIONS',L0,IP0,'TABLE',I1,X1,' ',L1,ITAB2)
  19 :          if (imo1.eq.5) CALL ACCTAB(ITAB1,'MOT',I0,X0,
  20 :      &'REACTIONS',L0,IP0,'TABLE',I1,X1,' ',L1,ITAB2)
  21 :            if (ierr.ne.0) return
  22 :                       
  23 :         ISUCCES = 0                                                      
  24 :        if (ire1.gt.0) then                                             
  25 :         ICOMP = -1                                                               
  25 :          CALL ACCTAB(ITABT,'ENTIER',(NDIM - 1),X0,' ',L0,IP0,                
  26 :      &                'FLOTTANT',I1,X1,' ',L1,IP1)
  27 :         VALREF = X1 * 1.E-6 
  28 :         DO IBOUC = 1,NDIM                                                                                            
  29 :          ICOMP = ICOMP + 1                                                                                                            
  29 :          CALL ACCTAB(ITABT,'ENTIER',ICOMP,X0,' ',L0,IP0,                
  30 :      &                'FLOTTANT',I1,VAL1,' ',L1,IP1)
  31 :           if (dabs(val - val1).le. valref) then                                       
  31 :          CALL ACCTAB(ITAB2,'ENTIER',ICOMP,X0,' ',L0,IP0,                
  32 :      &                'CHPOINT',I1,VAL1,' ',L1,ICHP1)                         
  33 :            ISUCCES = 1
  34 : c         write(6,*)'ON A TROUVE UN RESULTAT AU TEMPS ', VAL1                                                                   
  35 :           RETURN                                                     
  36 :           endif                                                              
  37 :         ENDDO                                                               
  38 :        else                                                          
  38 :          CALL ACCTAB(ITAB2,'ENTIER',(NDIM-1),X0,' ',L0,IP0,                
  39 :      &                'CHPOINT',I1,VAL1,' ',L1,ICHP1)                                   
  40 :          ISUCCES = 1                           
  41 :        endif   
  42 : 
  43 :        if (ISUCCES.EQ.0) then
  43 :         call erreur(-10)
  44 : *       'ERREUR' '*** On n a pas trouve le résultat cherché.'                  
  45 :        endif
  46 :        return
  47 :        end
  48 :  

© Cast3M 2003 - All rights reserved.
Disclaimer