74 |
ELSE |
ELSE |
75 |
C Base of domain is in this cell |
C Base of domain is in this cell |
76 |
C Set hFac to the fraction of the cell that is open. |
C Set hFac to the fraction of the cell that is open. |
77 |
C hFacC(I,J,K,bi,bj) = (rF(K)*rkFac-H(I,J,bi,bj)*rkFac)*recip_drF(K) |
hFacC(I,J,K,bi,bj) = (rF(K)*rkFac-H(I,J,bi,bj)*rkFac)*recip_drF(K) |
|
CcnhDebugStarts |
|
|
C Impose full-step |
|
|
hFacC(I,J,K,bi,bj) = 1. |
|
|
CCnhDebugEnds |
|
78 |
ENDIF |
ENDIF |
79 |
C Impose minimum fraction |
C Impose minimum fraction |
80 |
IF (hFacC(I,J,K,bi,bj).LT.hFacMin) THEN |
IF (hFacC(I,J,K,bi,bj).LT.hFacMin) THEN |