summaryrefslogtreecommitdiff
path: root/gcc/testsuite/gfortran.dg/graphite/id-21.f
blob: 4fa047ed6f283aca4eb8dc7c9c54028b846bea73 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
      MODULE LES3D_DATA
      DOUBLE PRECISION,ALLOCATABLE,DIMENSION(:,:,:) ::
     >             P, T, H
      DOUBLE PRECISION,ALLOCATABLE,DIMENSION(:,:,:,:) ::
     >             HF
      DOUBLE PRECISION,ALLOCATABLE,DIMENSION(:,:,:,:,:) ::
     >             Q
      END MODULE LES3D_DATA
      USE LES3D_DATA
      DO K = 1, KMAX - 1
         DO J = 1, JMAX - 1
            DO I = 1, I2
               T(I,J,K) = (EI - HF(I,J,K,1)) / HF(I,J,K,3)
            ENDDO
            P(1:I2,J,K) = Q(1:I2,J,K,1,M) * HF(1:I2,J,K,4) * T(1:I2,J,K)
            IF(ISGSK .EQ. 1) H(1:I2,J,K) =
     >                   (Q(1:I2,J,K,5,M) + P(1:I2,J,K))
         END DO
      ENDDO
      END