Download modp.eso

Back to the list

   1 : C MODP      SOURCE    PV        05/09/22    21:20:54     5181           
   2 : C   DEPLACEMENT INTERACTIF D'UN NOEUD                                           
   3 : C                                                                               
   4 :       SUBROUTINE MODP(XPROJ,IVU,KON1,ICPR,IDCP,IBOUJ)                           
   5 :       IMPLICIT INTEGER(I-N)
   6 :       COMMON /CMODI/LIGMAX,XPREC,YPREC                                          
   7 :       DIMENSION XTR(2),YTR(2),ZTR(2)                                            
   8 :       SEGMENT /KON1/(KON(NBCON,NMAX)),KON2.KON1                                 
   9 :       SEGMENT XPROJ(3,ITE)                                                      
  10 :       SEGMENT IVU(0)                                                            
  11 :       SEGMENT ICPR(0)                                                           
  12 :       SEGMENT IDCP(0)                                                           
  13 :       SEGMENT IBOUJ(0)    
  14 :       ztr(1)=0
  15 :       ztr(2)=0                                                      
  16 :       NBCON=KON(/1)                                                             
  17 :       NBCONR=NBCON-1                                                            
  17 :       CALL TRMESS('POINTEZ LE POINT A DEPLACER')                                
  18 :   30  CONTINUE                                                                  
  18 :       CALL MOPF3                                                                
  18 :       CALL TRDIG(X,Y,INCLE)                                                     
  19 :       IF (INCLE.EQ.3) RETURN                                                    
  20 : *  CHERCHER LE POINT DU MAILLAGE                                                
  21 :       NBP=XPROJ(/2)                                                             
  22 :       CRIT=XPREC**2                                                             
  23 :       DO 10 IP=1,NBP                                                            
  24 :       IF (IVU(IP).NE.1) GOTO 10                                                 
  25 :       DIST=(XPROJ(1,IP)-X)**2+(XPROJ(2,IP)-Y)**2                                
  26 :       IF (DIST.LT.CRIT) GOTO 20                                                 
  27 :   10  CONTINUE                                                                  
  28 :       GOTO 30                                                                   
  29 :   20  CONTINUE                                                                  
  29 :       CALL TRMESS('POINTEZ LA NOUVELLE POSITION DU POINT')                      
  29 :       CALL TRDIG(X,Y,INCLE)                                                     
  29 :       CALL TRMESS(' ')                                                          
  30 :       IF (INCLE.EQ.3) RETURN                                                    
  31 :       XPROJ(1,IP)=X                                                             
  32 :       XPROJ(2,IP)=Y                                                             
  33 :       IOEE=IDCP(IP)                                                             
  33 :       CALL PROMOD(ICPR,XPROJ,IOEE,4,IBOUJ)                                      
  34 : *  REAFFICHER LES LIGNES MODIFIES 
  34 :       call insegt(3,iresu)                                            
  34 :       CALL CHCOUL(4)                                                            
  35 :       ICHAIN=IP                                                                 
  36 :   70  CONTINUE                                                                  
  37 :       DO 50 ICON=1,NBCONR                                                       
  38 :       IP=KON(ICON,ICHAIN)                                                       
  39 :       IF (IP.EQ.0) GOTO 60                                                      
  40 :       XTR(1)=X                                                                  
  41 :       YTR(1)=Y                                                                  
  42 :       XTR(2)=XPROJ(1,IP)                                                        
  43 :       YTR(2)=XPROJ(2,IP)                                                        
  43 :       CALL POLRL(2,XTR,YTR,ZTR)                                                     
  44 :   50  CONTINUE                                                                  
  45 :       ICHAIN=KON(NBCON,ICHAIN)                                                  
  46 :       IF (ICHAIN.NE.0) GOTO 70                                                  
  47 :   60  CONTINUE                                                                  
  48 :       END                                                                       
  49 :  
  50 :  

© Cast3M 2003 - All rights reserved.
Disclaimer