Download dyne25.eso

Back to the list

   1 : C DYNE25    SOURCE    CHAT      05/01/12    23:16:03     5004           
   2 :       SUBROUTINE DYNE25(IBAS,KCPR,IK)                                           
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 : *--------------------------------------------------------------------*          
   6 : *                                                                    *          
   7 : *     Op{rateur DYNE : algorithme de Fu - de Vogelaere               *          
   8 : *     ________________________________________________               *          
   9 : *                                                                    *          
  10 : *     Recherche des points de r{f{rence d{finis par les points       *          
  11 : *     rep}res de la base modale.                                     *          
  12 : *                                                                    *          
  13 : *     Param}tres:                                                    *          
  14 : *                                                                    *          
  15 : * e   IBAS    Table repr{sentant une base modale                     *          
  16 : *  s  KCPR    Segment des points                                     *          
  17 : * es  IK      Compteur                                               *          
  18 : *                                                                    *          
  19 : *     Auteur, date de cr{ation:                                      *          
  20 : *                                                                    *          
  21 : *     Lionel VIVAN, le 24 octobre 1989.                              *          
  22 : *                                                                    *          
  23 : *--------------------------------------------------------------------*          
  24 : -INC CCOPTIO                                                                    
  25 : -INC SMCOORD                                                                    
  26 : *                                                                               
  27 :       SEGMENT,ICPR(XCOOR(/1)/(IDIM+1))                                          
  28 : *                                                                               
  29 :       LOGICAL L0,L1                                                             
  30 :       CHARACTER*4  MO2                                                          
  31 :       CHARACTER*8  TYPRET,CHARRE                                                       
  32 : *                                                                               
  33 :       ICPR = KCPR                                                               
  34 : *                                                                               
  35 :       IMODE = 0                                                                 
  36 :  10   CONTINUE                                                                  
  37 :       ITBMOD = 0                                                                
  38 :       IMODE = IMODE + 1                                                         
  39 :       TYPRET = ' '                                                              
  39 :       CALL ACCTAB(IBAS,'ENTIER',IMODE,X0,' ',L0,IP0,                            
  40 :      &                  TYPRET,I1,X1,CHARRE,L1,ITBMOD)                             
  41 :       IF (IERR.NE.0) RETURN                                                     
  42 :       IF (ITBMOD.NE.0) THEN                                                     
  42 :          CALL ACCTAB(ITBMOD,'MOT',I0,X0,'SOUSTYPE',L0,IP0,                      
  43 :      &               'MOT',I1,X1,MO2,L1,IP1)                                    
  44 :          IF (IERR.NE.0) RETURN                                                  
  45 :          IF (MO2.EQ.'MODE') THEN                                                
  46 : *                                                                               
  47 : *           Le point rep}re nous {choit                                         
  48 : *                                                                               
  48 :             CALL ACCTAB(ITBMOD,'MOT',I0,X0,'POINT_REPERE',L0,IP0,               
  49 :      &                         'POINT',I1,X1,' ',L1,KNOE)                       
  50 :             IF (IERR.NE.0) RETURN                                               
  51 :             IF (KNOE.NE.0) THEN                                                 
  52 :                IF (ICPR(KNOE).EQ.0) THEN                                        
  53 :                   IK = IK + 1                                                   
  54 :                   ICPR(KNOE) = IK                                               
  55 :       IF (IIMPI.EQ.333) THEN                                                    
  56 :       WRITE(IOIMP,*)'DYNE25 : basemo. ICPR(',KNOE,')=',ICPR(KNOE)               
  57 :       ENDIF                                                                     
  58 :                ENDIF                                                            
  59 :             ENDIF                                                               
  60 :          ELSE                                                                   
  60 :             CALL ERREUR(314)                                                    
  61 :             RETURN                                                              
  62 :          ENDIF                                                                  
  63 :          GOTO 10                                                                
  64 :       ENDIF                                                                     
  65 : *                                                                               
  66 :       END                                                                       
  67 :  
  68 :  

© Cast3M 2003 - All rights reserved.
Disclaimer