Esope source files

   1 : C CYLSYM    SOURCE    CHAT      05/01/12    22:34:10     5004           
   2 :       SUBROUTINE CYLSYM(ITYPE,EDICON)                                           
   3 : C                                                                               
   4 : C     SET UP SYMMETRY CONDITIONS FOR CYLINDER COORDINATES CASE                  
   5 : C-----------------------------------------------------------------------        
   6 : C                                                                               
   7 :       IMPLICIT INTEGER(I-N)
   8 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   9 : C     PROBLEM CONSTANTS                                                         
  10 : C                                                                               
  11 : C                                                                               
  12 : C     CONSTANTS DESCRIBING SYMMETRY CONDITIONS                                  
  13 : C                                                                               
  14 :       segment EDICON                                                            
  15 :         integer  KSTRT, KSTEP, NMIR, is                                         
  16 :         real*8  crot, srot, symfct                                              
  17 :         LOGICAL LREAL, LIMAG                                                    
  18 :       endsegment                                                                
  19 : C-----------------------------------------------------------------------        
  20 : C                                                                               
  21 :       ISCAL= 2                                                                  
  22 :       IS = ISCAL - 1                                                            
  23 : C                                                                               
  24 :       IF(ITYPE) 10, 20, 30                                                      
  25 : C                                                                               
  26 : C     DIRICHLET CONDITION ON R-AXIS --- USE ALL ODD ORDERS                      
  27 : C                                                                               
  28 :    10 KSTRT = 2                                                                 
  29 :       KSTEP = 2                                                                 
  30 :       RETURN                                                                    
  31 : C                                                                               
  32 : C     NO SYMMETRY ON R-AXIS --- USE ALL ORDERS                                  
  33 : C                                                                               
  34 :    20 KSTRT = 1                                                                 
  35 :       KSTEP = 1                                                                 
  36 :       RETURN                                                                    
  37 : C                                                                               
  38 : C     NEUMANN CONDITION ON R-AXIS -- USE ALL EVEN ORDERS                        
  39 : C                                                                               
  40 :    30 KSTRT = 1                                                                 
  41 :       KSTEP = 2                                                                 
  42 :       RETURN                                                                    
  43 : C***********************************************************************        
  44 :       END                                                                       
  45 :