55 |
C ligand balance in surface layer |
C ligand balance in surface layer |
56 |
C in surface layer |
C in surface layer |
57 |
|
|
58 |
DO j=jMin,jMax |
DO j=1-OLy,sNy+OLy |
59 |
DO i=iMin,iMax |
DO i=1-OLx,sNx+OLx |
60 |
DO k=1,nR |
DO k=1,nR |
61 |
IF(hFacC(i,j,k,bi,bj) .gt. 0.0)THEN |
IF(hFacC(i,j,k,bi,bj) .gt. 0.0)THEN |
62 |
|
|
70 |
|
|
71 |
FeL = ligand_tot-lig |
FeL = ligand_tot-lig |
72 |
freefe(i,j,k,bi,bj) = fe (i,j,k,bi,bj)-FeL |
freefe(i,j,k,bi,bj) = fe (i,j,k,bi,bj)-FeL |
73 |
|
#ifdef MINFE |
74 |
|
freefe(i,j,k,bi,bj) = min(freefe(i,j,k,bi,bj),freefemax) |
75 |
|
fe(i,j,k,bi,bj) = FeL+freefe(i,j,k,bi,bj) |
76 |
|
#endif |
77 |
END IF |
END IF |
78 |
ENDDO |
ENDDO |
79 |
ENDDO |
ENDDO |