def METHOD_NAME(self,k):
    """
    REQUIRED INPUT:
        - all Jarvis stress functions
        - windspeed at 2m height (m/s)
        - sgamma (slope of saturation vapour function)
    PARAMETERS:
        - rstmin (minimal stomatal resistence, depending on land use, s/m)
        - gamma (psychrometric constant)
    """
    JC_all = self.JC_rad * self.JC_vpd * self.JC_temp * self.JC_sm
    stomRes1 = ifthenelse(JC_all == 0, 50000, self.rst_lai[k] / JC_all)
    stomRes2 = ifthenelse(stomRes1 > 50000, 50000, stomRes1)
    self.stomRes = stomRes2 / (3600*24)
    self.k = 1 / (1 + self.stomRes * self.gamma / (self.aeroRes * (self.sgamma + self.gamma)))
    self.JC_k_[k] = self.k