1! { dg-do compile }
2! { dg-options "-O2" }
3      SUBROUTINE UMPSE(AIBJ,NOC,NDIM,NOCA,NVIRA,NOCCA,E2)
4      DIMENSION AIBJ(NOC,NDIM,*)
5      DO 20 MA=1,NVIRA
6      DO 20 MI=1,NOCA
7         DO 10 MB=1,MA
8         MBI = MI+NOCA*(MB-1)
9         DO 10 MJ=1,NOCCA
10            DUM = AIBJ(MJ,MAI,MB)-AIBJ(MJ,MBI,MA)
11            E2A = E2A-DUM
12   10    CONTINUE
13   20 CONTINUE
14      E2  = E2+E2A
15      END
16