      SUBROUTINE BMCM (TOV,P1,P2,H1,H2,DEL,CI,BM,CM,DT,SV,AM)
      DIMENSION BM( 7 , 7 ), CM( 7 , 7 ),
     1  TOV( 7 ), P1( 7 ), P2( 7 ), H1( 7 ), H2( 7 ),
     2  DEL( 7 ), CI( 8 ), SV( 7 ), AM( 7 , 7 )
      RK =  287.05  /  1005.
      DO 1 LE=1, 6
1     H1(LE) = P1(LE) * TOV(LE+1) - TOV(LE)
      H1( 7 ) = 0.
      H2( 1 ) = 0.
      DO 2 LE=2, 7
2     H2(LE) = TOV(LE) - P2(LE) * TOV(LE-1)
      DO 3 K=1, 7
      KK = K
      SAV1 = RK*TOV(K) + (CI(K+1)*H1(K)+CI(K)*H2(K)) / (2.*DEL(K))
      SAV2 = (H1(K)+H2(K)) / (2.*DEL(K))
      DO 33 J=1, 7
33    BM(K,J) = -SAV1*DEL(J)
      DO 34 J=1,KK
34    BM(K,J) = BM(K,J) + SAV2*DEL(J)
      BM(K,K) = BM(K,K) - H2(K) / 2.
3     CONTINUE
      RAA =  287.05  / ( 6370000. * 6370000. )
      DO 4 I=1, 7
      DO 4 J=1, 7
      CM(I,J) = 0.
      DO 44 K=1, 7
44    CM(I,J) = CM(I,J) + AM(I,K) * BM(K,J)
      CM(I,J) = (CM(I,J)+RAA*TOV(I)*SV(J))*DT*DT
4     CONTINUE
C(DEBUGLINE) VERSION CREATED JULY 8 1982 BY J.SELA,NMC.
      RETURN
      END
