! print out many of the biomass pool components (used for debugging purposes) ! These files use the following columnar format. Some are filled with zeros ! to make it easier to select specific columns for comparisons between the ! crop and individual biomass pools (not all pools have the same variables) subroutine bpools (cd,cm,cy,isr) integer cd,cm,cy,isr real total, saitotal !added by Simon include 'p1werm.inc' include 'file.fi' include 'm1flag.inc' include 'b1glob.inc' include 'd1glob.inc' include 'c1glob.inc' include 'c1db1.inc' include 'decomp/decomp.inc' ! statements below added by Simon include 'w1clig.inc' ! include 's1layr.inc' ! include 'p1const.inc' ! include 'm1dbug.inc' ! These hydrology common blocks provide soil temp, moisture and irrigation ! include 'h1temp.inc' ! include 'h1db1.inc' ! include 'h1hydro.inc' 1234 format (3i3,16f10.3) 2345 format (3i3,3f10.5,13f10.3) ! day, month, year ! flat crop cover, standing crop cover, total crop cover ! crop cover fract, crop SAI, crop LAI ! total crop biomass, 0.0, standing crop mass ! crop root mass, 0.0, crop yield mass ! number crop stems per area, crop height, crop root depth, repr stem dia if ((am0dfl .eq. 1).or.(am0dfl.eq.3)) then ! day, month, year ! flat crop cover, standing crop cover, total crop cover ! crop cover fract, crop SAI, crop LAI ! total crop biomass, 0.0, standing crop mass ! (no "flat crop biomass") ! crop root mass, 0.0, crop yield mass ! (no "buried crop biomass") ! qty crop stems per area, crop height, crop root depth, repr stem dia ! Crop Pool if ((cd .eq. 1) .and. (cm .eq. 1) .and. (cy .eq. 1)) then write(luocrp1,*) 'dd mm yy tmin tmax tavg tempc', & & ' precip waterc waterc dd dd', & & ' stand1 stand2 stand3 stand', & & ' flat1 flat2 flat3 flat total', & & ' belo1 belo2 belo3 belo ', & & ' root1 root2 root3 root ', & & ' scov1 scov2 scov3 scov fcov1 fcov2', & & ' fcov3 fcov tcov1 tcov2 tcov3 tcov ', & & ' SAI1 SAI2 SAI3 SAI SAI SAI', & & ' LAI1 LAI2 LAI3 LAI LAI ', & & ' silho', & & ' stem1 stem2 stem3 stem hght1 hght2 hght3' write(luocrp1,*) ' ', & & ' stand flat stand flat ', & & ' ', & & ' ', & & ' ', & & ' ', & & ' ', & & ' dead crop bmass', & & ' crop bmass', & & ' bmass' write(luocrp1,*) ' C C C ', & & ' mm kg/m2 kg/m2 kg/m2', & & ' kg/m2 kg/m2 kg/m2 kg/m2 kg/m2 kg/m2', & & ' kg/m2 kg/m2 kg/m2 kg/m2', & & ' kg/m2 kg/m2 kg/m2 kg/m2', & & ' m2/m2 m2/m2 m2/m2 m2/m2 m2/m2 m2/m2', & & ' m2/m2 m2/m2 m2/m2 m2/m2 m2/m2 m2/m2', & & ' m2/m2 m2/m2 m2/m2 m2/m2 m2/m2 m2/m2', & & ' m2/m2 m2/m2 m2/m2 m2/m2 m2/m2', & & ' m2/m2', & & ' no/m2 no/m2 no/m2 no/m2 m m m' endif total = abmf(1) + abmst(1) !sum of standing and flat residue mass, all pools saitotal = adrsai(1,1) + adrsai(2,1) + adrsai(3,1)!stem area index, sum of all residue pools write(luocrp1,2222) cd, cm, cy, & !day, month, year & awtdmn, awtdmx, awtdav, ditca, & !tmin, tmax, tavg, tf & aqua, diwcs, diwcf, didds, diddf, & !precip, wf standing, wf flat, dd standing, dd flat & admst(1,1), admst(2,1), admst(3,1), abmst(1), & !mass, standing & admf(1,1), admf(2,1), admf(3,1), abmf(1), & !mass, flat & total, & !sum of standing and flat residue mass, all pools & admbg(1,1), admbg(2,1), admbg(3,1), abmbg(1), & !mass, below ground & admrt(1,1), admrt(2,1), admrt(3,1), abmrt(1), & !mass, roots & adfscv(1,1), adfscv(2,1), adfscv(3,1),abfscv(1), & !cover provided by standing residue (fraction) & adffcv(1,1), adffcv(2,1), adffcv(3,1),abffcv(1), & !cover provided by flat residue (fraction) & adftcv(1,1), adftcv(2,1), adftcv(3,1),abftcv(1), & !cover provided by standing+flat residue (fraction) & adrsai(1,1), adrsai(2,1), adrsai(3,1), & !stem area index for each residue pool & saitotal, & !stem area index, sum of all residue pools & acrsai(1), & !stem area index, crop & abrsai(1), & !stem area index, total biomass & adrlai(1,1), adrlai(2,1), adrlai(3,1), & !leaf area index for each residue pool & acrlai(1), & !leaf area index, crop & abrlai(1), & !leaf area index, total biomass & abrcd(1), & !biomass silhouette area & addstm(1,1), addstm(2,1), addstm(3,1), abdstm(1), & !stems (no/m2) & adzht(1,1), adzht(2,1), adzht(3,1) !stubble for each residue pool ! tf=temperature factor, wf=water factor, dd=decomposition day 2222 format (2i3,i4,3f7.1,f7.3,f7.2,4f7.3,17f7.4,24f7.4,4f7.1,3f7.2) !added by Simon ! & acffcv(isr), acfscv(isr), acftcv(isr), ! & accovfact(isr), acrsai(isr), acrlai(isr), ! & acm(isr), 0.0, acmst(isr), ! & acmrt(isr), 0.0, acmyld(isr), ! & acdstm(isr), aczht(isr), aczrtd(isr), acxstmrep(isr) ! day, month, year ! flat residue cover, standing residue cover, total residue cover ! residue cover fract, residue SAI, residue LAI ! total residue biomass, flat residue mass, standing residue mass ! residue root mass, below gnd residue mass ! qty residue stems per area, "ave" residue height, 0.0, 0.0 ! (no "ave" root depth or stem dia computed across residue pools) ! All Residue Pools Combined write(luobio1,2345) cd, cm, cy, & & abffcv(isr), abfscv(isr), abftcv(isr), & & 0.0, abrsai(isr), abrlai(isr), & & abm(isr), abmf(isr), abmst(isr), & & abmrt(isr), abmbg(isr), & & abdstm(isr), abzht(isr), 0.0, 0.0 endif if ((am0dfl .eq. 2).or.(am0dfl.eq.3)) then ! day, month, year ! flat residue cover, standing residue cover, total residue cover ! residue cover fract, residue SAI, residue LAI ! total residue biomass, flat residue mass, standing residue mass ! residue root mass, below gnd residue mass, 0.0 ! (no residue yield mass) ! qty residue stems per area, residue height, 0.0, rep stem dia ! (no root depth for residue pools) ! Residue Pool #1 write(luodec1,2345) cd, cm, cy, & & adffcv(1,isr), adfscv(1,isr), adftcv(1,isr), & & covfact(1,isr), adrsai(1,isr), adrlai(1,isr), & & adm(1,isr), admf(1,isr), admst(1,isr), & & admrt(1,isr), admbg(1,isr), 0.0, & & addstm(1,isr), adzht(1,isr), 0.0, adxstmrep(1,isr) ! Residue Pool #2 write(luodec2,2345) cd, cm, cy, & & adffcv(2,isr), adfscv(2,isr), adftcv(2,isr), & & covfact(2,isr), adrsai(2,isr), adrlai(2,isr), & & adm(2,isr), admf(2,isr), admst(2,isr), & & admrt(2,isr), admbg(2,isr), 0.0, & & addstm(2,isr), adzht(2,isr), 0.0, adxstmrep(2,isr) ! Residue Pool #3 write(luodec3,2345) cd, cm, cy, & & adffcv(3,isr), adfscv(3,isr), adftcv(3,isr), & & covfact(3,isr), adrsai(3,isr), adrlai(3,isr), & & adm(3,isr), admf(3,isr), admst(3,isr), & & admrt(3,isr), admbg(3,isr), 0.0, & & addstm(3,isr), adzht(3,isr), 0.0, adxstmrep(3,isr) endif end