c file name: nuts.for subroutine nuts i (y1, y2, m uu) c + + + PURPOSE + + + c This subroutine calculates a nutrient stress factor caused by limited c supply of N or P. c + + + KEYWORDS + + + c Nutrient stress c + + + COMMON BLOCKS + + + *$noereference include 'crop/cparm.inc' *$reference c + + + LOCAL VARIABLES + + + real y1, y2, yy, uu c + + + LOCAL VARIABLE DEFINITIONS + + + c y1 - cummulative amount of N or P taken by the plant (kg/ha) - supply c y2 - potential amount of N or P needed by the plant (kg/ha) - demand c yy - scaled ratio of supply over demand c uu - N or P stress factor c yy replaces uu where appropriate to minimize confusion c a_s8,b_s8 are used instead of scrp(8,1) and scrp(8,2) c + + + END OF SPECIFICATIONS + + + if (y2 .eq. 0.) goto 3 yy = 200. * (y1 / y2 - .5) if (yy .gt. 0.) go to 2 uu = 0. go to 3 2 continue uu = yy / (yy + exp(a_s8 - b_s8 * yy)) 3 continue return end