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