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