 = node_val(KK,i)*node_val(KK,i) / ( node_val(eps,i)*node_val(eps,i) )             an = tau0 * node_val(NN0,i)             ! clip an at minimum value             an = max(an,anLimitFact*anMin)              ! compute the equilibrium value of as             tmp0  = -d0 - (d0 + nt0)*an - (d0 + nt0)*an*an             tmp0  = -d0 + n0 + 