c Program assumes A, BB1, BB2 defined in calling program
c as follows:
c        A = DEXP(-2.847d0)
c        BB1 = DEXP(-.693d0)
c        BB2 = DEXP(-3.101d0)
        SUBROUTINE FNMDA (VSTOR, OPEN, numcells, numcomps,
     &    MG, L, A, BB1, BB2) 
c Computes "open" for all compartments of cell # L

           integer L, numcells, numcomps, i
       REAL*8 VSTOR(numcomps,numcells), OPEN(numcomps)
       REAL*8 A, BB1, BB2, VM, A1, A2, B1, B2, MG
c modify so that potential is absolute and not relative to
c  "rest"
C  TO DETERMINE VOLTAGE-DEPENDENCE OF NMDA CHANNELS
           DO 1, I = 1, numcomps
           VM = VSTOR(I,L)
           A1 = DEXP(-.016d0*VM - 2.91d0)
           A2 = 1000.d0 * MG * DEXP (-.045d0 * VM - 6.97d0)
           B1 = DEXP(.009d0*VM + 1.22d0)
           B2 = DEXP(.017d0*VM + 0.96d0)
        OPEN(I)     = 1.d0/(1.d0 + (A1+A2)*(A1*BB1 + A2*BB2) /
     X   (A*A1*(B1+BB1) + A*A2*(B2+BB2))  )
C  FROM JAHR & STEVENS, EQ. 4A
C               DO 124, J = 1, 19
C          OPEN(J) = 1./(1.+.667* EXP(-0.07*(VSTOR(J)-60.)) )
C  FROM CHUCK STEVENS
1               CONTINUE
                        END