pro grpresp1, weight, parm, temparr, resparr ;Calculates respiration and rates for given temperature, weight ;and species- size-specific parameters according to Hewett & Johnson 1992 ; for use with Tyler version of parameter array ; adapted from work of Horne, Luo, Mason et al. ; parameter array (parm) differs from the parameter array (param) used in ; the Luo/Horne versions of the bioenergetics routines. Req=parm(9) ; respiration equation set Ralpha=parm(10) ; resp alpha or RA Rbeta=parm(11) ; resp beta or RB Rthet=parm(12) ; resp theta or RQ Rtopt=parm(13) ; resp temp optimal or RTO Rtmax=parm(14) ; resp temp max or RTM Rtl=parm(15) ; resp RK1=parm(16) ; K1 RK4=parm(17) ; K4 Ract=parm(18) ; ACT Rbact=parm(19) ; BACT preded=parm(29) ; predator energy density ;calpgm=1200. ;Ralpha=Ralpha*3280./calpgm ;convert calories to grams.... kept from horne Ralpha=Ralpha*3280./preded ;convert calories to grams rsp= Ralpha * weight ^Rbeta ; Ft and activity computed differently in the two eq sets if (Req eq 1) then begin Ft= exp(Rthet * temparr) tv=temparr * 0.0 ; temp var for activity computation tp1=where(temparr gt Rtl,np1) tp2=where(temparr le Rtl,np2) if (np1 gt 0) then tv(tp1)= RK1 * weight(tp1) ^ RK4 if (np2 gt 0) then tv(tp2)= Ract * exp(Rbact * temparr) * weight(tp2)^RK4 activ=exp((Rtopt-(Rtmax*temparr))*tv) endif if (Req eq 2) then begin activ=Ract ty=alog(Rthet)*(Rtmax-Rtopt+2.0) tz=alog(Rthet)*(Rtmax-Rtopt) tx=(tz^2 * (1+(1+40/ty)^0.5)^2.)/400 tv=(Rtmax-temparr)/(Rtmax-Rtopt) Ft=tv^tx * exp(tx*(1-tv)) endif resparr=rsp * Ft * activ ;read,'in grpresp1 ',junk return end