Numérotation des lignes :

carsym
C CARSYM    SOURCE    CHAT      05/01/12    21:49:26     5004      SUBROUTINE CARSYM (ITYPE,EDICON)CC     SET UP SYMMETRY CONDITIONS FOR THE CARTESIAN (X,Y) COORDINATESC-----------------------------------------------------------------------CC     PROBLEM CONSTANTS      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8(A-H,O-Z)-INC CCREELCCC     CONSTANTS DESCRIBING SYMMETRY CONDITIONSC      segment EDICON        integer  KSTRT, KSTEP, NMIR, is        real*8  crot, srot,  symfct        LOGICAL LREAL, LIMAG      endsegmentC-----------------------------------------------------------------------C      IF(ITYPE .GE. 10) GO TO 100      IF(ITYPE .LE. 0) GO TO 10CC     "ITYPE" LESS THAN 10 --- PREDEFINED SYMMETRY CODEC      GO TO (10, 20, 30, 40, 50, 60, 70, 80, 90), ITYPEC   10 KSTRT = 1      KSTEP = 1      LREAL = .TRUE.      LIMAG = .TRUE.      GO TO 200C   20 KSTRT = 1      KSTEP = 1      LREAL = .TRUE.      LIMAG = .FALSE.      GO TO 200C   30 KSTRT = 2      KSTEP = 2      LREAL = .TRUE.      LIMAG = .FALSE.      GO TO 200C   40 KSTRT = 2      KSTEP = 4      LREAL = .TRUE.      LIMAG = .FALSE.      GO TO 200C   50 KSTRT = 2      KSTEP = 4      LREAL = .FALSE.      LIMAG = .TRUE.      GO TO 200C   60 KSTRT = 1      KSTEP = 2      LREAL = .TRUE.      LIMAG = .FALSE.      GO TO 200C   70 KSTRT = 3      KSTEP = 6      LREAL = .TRUE.      LIMAG = .FALSE.      GO TO 200C   80 KSTRT = 2      KSTEP = 2      LREAL = .TRUE.      LIMAG = .FALSE.      GO TO 200C   90 KSTRT = 4      KSTEP = 8      LREAL = .TRUE.      LIMAG = .FALSE.      GO TO 200CC     "ITYPE" MUST BE DECODEDC  100 KSTRT = ITYPE / 100      KSTEP = (ITYPE - 100 * KSTRT) / 10      MEDIAN = ITYPE - 100 * KSTRT - 10 * KSTEP      LREAL = MEDIAN .GT. 0      LIMAG = MEDIAN .NE. 1CC     COEFFICIENTS FOR ROTATIONAL SYMMETRYC  200 PHI = 2. * XPI / FLOAT(KSTEP)      CROT = COS(PHI)      SROT = SIN(PHI)      NMIR = 2      IF(LREAL .AND. LIMAG) NMIR = 1      SYMFCT = 1.0      IF(KSTEP .GT. KSTRT) SYMFCT = - 1.0      RETURN      END    

© Cast3M 2003 - Tous droits réservés.
Mentions légales