我需要将这段代码（fortran77）与openmp 4.5并行化。
DO 700 K=1,KMM1
K1 = KB1(K)
K2 = KB2(K)
DO 700 J=2,JM
JSB= JSBLK(J)
J1 = JB1(J)
J2 = JB2(J)
DO 700 I=1,IMM1
I1 = IB1(I)
I2 = IB2(I)
DELTA = STORE(I,J,K)
B1CHG(I1,J1,K1) = B1CHG(I1,J1,K1) + DELTA
B2CHG(I2,J2,K2) = B2CHG(I2,J2,K2) + DELTA
SBCHG(JSB) = SBCHG(JSB) + DELTA
700 CONTINUE
我在代码和检查结果的前后都使用了简单的语句！OMP PARALLEL DO，！OMP END PARALLEL DO，但是它们并不相同。我已经使用随机值初始化了向量KB1，KB2，JSBLK，JB1，JB2，IB1，IB2，STORE，B1CHG，B2CHG，SBCHG，并且在并行循环中，我对变量使用了不同的名称（例如，使用STOREP而不是STORE ），然后检查STOREP == STORE，B1CHGP == B2CHG等。如果我删除openmp编译指示，结果是相同的。这是我用于比较和测试串行和并行版本的代码：
PARAMETER(KMM1=10, JM=5 , IMM1=5 )
COMMON/VEC/ KB1(KMM1),KB2(KMM1),JSBLK(JM),JB1(JM),JB2(JM),
& IB1(IMM1),IB2(IMM1),STORE(IMM1,JM,KMM1),
& B1CHG(IMM1,JM,KMM1),B2CHG(IMM1,JM,KMM1),SBCHG(JM)
COMMON/PARVEC/ STOREP(IMM1,JM,KMM1),
& B1CHGP(IMM1,JM,KMM1),B2CHGP(IMM1,JM,KMM1),SBCHGP(JM)
include 'omp_lib.h'
C---------------------------------------------------------------------
C INITIALIZING DUMMY ARRAYS
C RANDOM_NUMBER E' UNA SUBROUTINE DEL FORTRAN 95
CALL RANDOM_SEED()
DO 701 K=1,KMM1
CALL RANDOM_NUMBER(U)
KB1(K)=FLOOR((KMM1+1)*U)
CALL RANDOM_NUMBER(U)
JH = 0 + FLOOR((KMM1+1)*u)
KB2(K)=JH
DO 701 J=2,JM
CALL RANDOM_NUMBER(U)
JH = 0 + FLOOR((JM+1)*u)
JSBLK(J)=JH
CALL RANDOM_NUMBER(U)
JH = 0 + FLOOR((JM+1)*U)
JB1(J)=JH
CALL RANDOM_NUMBER(U)
JH = 0 + FLOOR((JM+1)*U)
JB2(J)= JH
DO 701 I=1,IMM1
CALL RANDOM_NUMBER(U)
JH = 0 + FLOOR((IMM1+1)*U)
IB1(I)=JH
CALL RANDOM_NUMBER(U)
JH = 0 + FLOOR((IMM1+1)*U)
IB2(I)=JH
CALL RANDOM_NUMBER(STORE(I,J,K))
CALL RANDOM_NUMBER(B1CHG(I,J,K))
CALL RANDOM_NUMBER(B2CHG(I,J,K))
CALL RANDOM_NUMBER(SBCHG(J))
STOREP(I,J,K) = STORE(I,J,K)
B1CHGP(I,J,K) = B1CHG(I,J,K)
B2CHGP(I,J,K) = B2CHG(I,J,K)
SBCHGP(J)= SBCHG(J)
701 CONTINUE
C---------------------------------------------------------------------
C PLOTTING THE INITIALIZED ARRAYS FOR VISUAL INSPECTION
DO 702 K=1,KMM1
WRITE(6,*) KB1(K),KB2(K)
DO 702 J=2,JM
WRITE(6,*) JSBLK(J),JB1(J),JB2(J)
DO 702 I=1,IMM1
WRITE(6,*) IB1(I),IB2(I),STORE(I,J,K)
WRITE(6,*) B1CHG(I,J,K)
WRITE(6,*) B2CHG(I,J,K)
WRITE(6,*) SBCHG(J)
702 CONTINUE
C
C
C NOW THE SERIAL/PARALLEL CODE COMPARING
C
c---------------------------------------------------------------------
C ENTERING SERIAL REGION
SEC1 = omp_get_wtime ( )
DO 700 K=1,KMM1
K1 = KB1(K)
K2 = KB2(K)
DO 700 J=2,JM
JSB= JSBLK(J)
J1 = JB1(J)
J2 = JB2(J)
DO 700 I=1,IMM1
I1 = IB1(I)
I2 = IB2(I)
DELTA = STORE(I,J,K)
B1CHG(I1,J1,K1) = B1CHG(I1,J1,K1) + DELTA
B2CHG(I2,J2,K2) = B2CHG(I2,J2,K2) + DELTA
SBCHG(JSB) = SBCHG(JSB) + DELTA
700 CONTINUE
SEC = omp_get_wtime ( ) - SEC1
WRITE(6,*) '#CPU SERIAL=', SEC, 'SECONDS.'
C---------------------------------------------------------------------
C LEAVING SERIAL REGION
C---------------------------------------------------------------------
C ENTERING PARALLEL REGION
SEC1 = omp_get_wtime ( )
!$OMP PARALLEL DO
DO K=1,KMM1
K1 = KB1(K)
K2 = KB2(K)
DO J=2,JM
JSB= JSBLK(J)
J1 = JB1(J)
J2 = JB2(J)
DO I=1,IMM1
I1 = IB1(I)
I2 = IB2(I)
DELTAP = STOREP(I,J,K)
B1CHGP(I1,J1,K1) = B1CHGP(I1,J1,K1) + DELTAP
B2CHGP(I2,J2,K2) = B2CHGP(I2,J2,K2) + DELTAP
SBCHGP(JSB) = SBCHGP(JSB) + DELTAP
END DO
END DO
END DO
!$OMP END PARALLEL DO
SEC = omp_get_wtime ( ) - SEC1
WRITE(6,*) '#CPU PARALLEL=', SEC, 'SECONDS.'
C---------------------------------------------------------------------
C LEAVING PARALLEL REGION
C----------------------------------------------------------------------
C VERIFYING DATA CORRECTNESS
DO K=1,KMM1
DO J=2,JM
DO I=1,IMM1
IF(B1CHG(I,J,K).NE.B1CHGP(I,J,K)) THEN
T=B1CHG(I,J,K)-B1CHGP(I,J,K)
WRITE(6,*) 'ERROR ', T
END IF
IF(B2CHG(I,J,K).NE.B2CHGP(I,J,K)) THEN
WRITE(6,*) 'ERROR'
END IF
IF(SBCHG(J).NE.SBCHGP(J)) THEN
WRITE(6,*) 'ERROR'
END IF
END DO
END DO
END DO
RETURN
END