AME9.6standardWedJun16133522GMT2010node00001node00003node00009node00011node00013node00015node00017node00021node00023node00025node00027node00029node00031node00033node00035node00037node00039node00041node00043node00045node00047node00051node00053node00055node00057node00061node00063node00065node00067node00069node00071node00073node00075node00077node00079node00083node00085node00087node00089node00091node00093node00095node00097node00099node00101node00103node00105node00107node00109node00113node00115node00117node00206node00207node00208node00209node00301node00303node00305node00311node00312node00313node00314node00315node00316node00317node00318node00323node00324node00327node00328node00329node00330node00331node00332node00333node00334node00341node00342node00343node00344node00345node00346node00351node00352node00353node00354node00355node00356node00357node00358node00359node00360node00361node00362node00363node00364node00365node00366node00367node00368node00369node00370node00371node00372node00373node00374node00375node00376node00377node00378node00379node00380node00381node00382node00387node00388node00389node00390node00391node00392node00393node00394node00395node00396node00403node00404node00411node00412node00413node00414node00419node00420node00501node00502node00503node00504node00505node00506node00509node00510node00515node00516node00517node00518node00523node00524node00525node00526node00527node00528node00529node00530node00531node00532node00535node00536node00537node00538node00539node00540node00541node00542node00543node00544node00545node00546node00547node00548node00549node00550node00551node00553node00555node00557node00558node00601node00602node00603node00604node00606node00607node00701node00703node00704node00705node00706node00711node00712node00713node00714node00718node00719node00720node00721node00726node00727node00729node00730node00731node00737node00738node00745node00746node00747node00748node00749node00750node00751node00752node00772node00773node0078025/1/2006 RM All labels, descriptions and units for variables are exactly as given in the technical documentation for the VB, spreadsheet version. Sometimes the descriptions and units are rather dubious: e.g. Monthly transpiration rate in current period (mm) probably should be Transpiration rate (mm.month-1) (since month is the basic unit of time). 16/6/2010 JM The timestep MUST be 1 and the time unit is 1 month.true3PG model: derived from excel spreadsheet versionnonenonenone3pg_v10E. globus = 1.0trueRatio of foliage stem partitioning at B=2pFS2 file1E. globus = 0.15trueRatio of foliage:stem partitioning at B = 20 cmpFS20file1E.globus = 0.8trueMaximum fraction of NPP to rootspRxfile1E. globus = 0.25trueMinimum fraction of NPP to rootspRnfile1E. globus = 0.027trueLitterfall rate at t = 0, month-1gammaF0file1E. globus = 0.001trueLitterfall rate for mature stands Litterfall rate for mature stands, month-1gammaF1file1E. globus = 12trueAge at which litterfall rate has median value, monthtgammaFfile1E. globus = 8.5trueMinimum temperature for growth, degreeCTmin1file1E. globus = 16trueOptimum temperature for growth, degreeCTopt1file1E. globus = 0trueNumber of days production lost for each frost day , dayskFfile1E. globus = 0trueValue of m when FR = 0, dimensionlessm0file1E. globus = 1trueValue of fN when FR = 0fN0file1E. globus = 0truePower of (1-FR) in fN Power of (1-FR) in fN, dimensionlessfNnfile1E. globus = 0.05trueDefines stomatal response to VPD, mbarCoeffCondfile1E. globus, sandy = 0.7 E. globus, sandy loam = 0.6 E. globus, Clay loam = 0.5 E. globus, Clay = 0.4trueMoisture ratio deficit which gives f? = 0.5, dimensionlessSWconstfileanyE. globus, sandy = 9 E. globus, sandy loam = 7 E. globus, clay loam = 5 E. globus, clay = 3truePower of moisture ratio deficit in f?, dimensionlessSWpowerfileanyE. globus = 50trueMaximum stand age used to compute relative age, yearMaxAgefile1E. globus = 4trueParameter in age modifier for fAge, dimensionlessnAgefile1E. globus = 0.95trueRelative age to give fage = 0.5, dimensionlessrAgefile1E. globus = 0.02trueMaximum canopy conductance, m.s-1MaxCondfile1E. globus = 3.33trueCanopy LAI for maximum canopy conductance, m2.m-2LAIgcxfile1E. globus = 0trueSeedling mortality rate at t=0, year-1gammaN0file1E. globus = 0trueMortality rate for older stands at large t, year-1gammaN1file1E. globu s= 2trueAge at which ?N = \xbd\(?N0+?N1), yeartgammaNfile1E. globus = 1trueShape of mortality response, dimensionlessngammaNfile1E. globus = 0.0trueFraction of foliage biomass pool per tree on each dying tree, dimensionlessmFfile1E. globus = 0.2trueFraction of root biomass pool per tree on each dying tree, dimensionlessmRfile1E. globus = 0.2trueFraction of stem biomass pool per tree on each dying tree, dimensionlessmSfile1E. globus = 11trueSpecific leaf area at stand age 0, m2.kg-1SLA0file1E. globus = 4trueSpecific leaf area for mature aged stands, m2.kg-1SLA1file1E. globus = 2.5trueAge at which specific leaf area = 12(?0+?1), yeartSLAfile1E. globus = 0.15trueMaximum fraction of rainfall intercepted by canopy, dimensionlessMaxIntcptnfile1E. globus = 0trueLAI for maximum rainfall interception, m2.m-2LAImaxIntcptnfile1E. globus = 0.5trueExtinction coefficient for PAR absorption by canopy, dimensionlesskfile1E. globus = 0 (?)trueAge at full canopy cover, yearfullCanAgefile1E. globus = 0.47trueRatio NPP/GPP, dimensionlessyfile1E. globus = 0.75trueBranch and bark fraction at stand age 0, dimensionlessfracBB0file1E. globus = 0.15trueBranch and bark fraction for mature aged stands, dimensionlessfracBB1file1E. globus = 2trueAge at which pBB = (p?0+p?1)/2, yeartBBfile1E. globus = 0.5trueMinimum basic density XI0c for older trees, t.m-3rho1file1E. globus = 4trueAge at which ? = \xbd\ density of old and young trees, yeartRhofile1E. globus = 0.0trueConstant in stem height relationship, dimensionlessaHfile1E. globus = 0.0truePower of DBH in stem height relationship, dimensionlessnHBfile1E. globus = 0.0truePower of stocking in stem height relationship, dimensionlessnHNfile1E. globus = 0.0trueConstant in stem volume relationship, dimensionlessaVfile1E. globus = 0.0truePower of DBH in stem volume relationship, dimensionlessnVBfile1E. globus = 0.0truePower of stocking in stem volume relationship, dimensionlessnVNfile1E. globus = -90trueIntercept of net radiation vs solar radiation relationship, W.m-2Qafile1E. globus = 0.8trueSlope of net radiation v solar radiation relationship, dimensionlessQbfile1E. globus = 2.3trueConversion of solar radiation to PAR, mol.MJ-1molPAR_MJfile1E. globus = 40trueMaximum temperature for growth, degreeCTmax1file1E. globus = 0.5trueMinimum basic density XI0c for young trees, t.m-3rho0file1trueTmintruefn3element([Tmin],month)1elementTminmonthtrueTmaxtruefn4element([Tmax],month)1elementTmaxmonth1= sand 2 = sandy-loam 3 = clay-loam 4 = claytrueSoil class, dimensionlessSoilClassfile1Seems to lie between 0 and 1.trueFertility rating, dimensionlessFRfile1Site Forcett: 190 mmtrueMaximum available soil water, mmmaxASWfile1trueMean number of frost days per month, day.month-1FrostDaystruefn6element([FrostDays],month)1elementFrostDaysmonthtrueMean daily incident solar radiation, MJ.m-2.day-1SolarRadtruefn7element([SolarRad],month)1elementSolarRadmonthtrueMean daily temperature, degreeCTavtruefn8(Tmax+Tmin)/21TmaxTmin2trueMean day-time VPD, mbarVPDtruefn9element([VPD],month)1elementVPDmonthtrueStand age, yearStandAgetruefn12time(1)/12.01time112.0trueStand basal area, m2.ha-1BasalAreatruefn14(((avDBH/200)^2)*3.141592654)*StemNo111avDBH20023.141592654StemNotrueStand volume excluding branch and bark, m3.ha-1StandVoltruefn15if aV>0 then aV*avDBH^nVB*StemNo^nVN else WS*(1-fracBB)/Density1aVavDBHnVBStemNonVNgtaV0WS1fracBBDensitytrueStand-basedmean DBH, cmavDBHtruefn16(AvStemMass/aWs)^(1.0/nWs)1AvStemMassaWs1.0nWstrueMean annual volume increement, m3.ha-1.year-1MAItruefn17if StandAge>0 then StandVol/StandAge else 01StandVolStandAgegtStandAge00trueSpecific leaf area, m2.kg-1SLAtruefn21if tSLA != 0 then SLA1 + (SLA0 - SLA1) * exp(-0.693147181 * (StandAge/ tSLA ^ 2.0)) else tSLA1SLA1SLA0SLA1exp0.693147181StandAgetSLA2.0neqtSLA0tSLAtrueFraction of ground area covered by canopy, dimensionlessCanCovertruefn22if fullCanAge>0 and StandAge<fullCanAge then (StandAge+0.01)/fullCanAge else 11StandAge0.01fullCanAgegtfullCanAge0ltStandAgefullCanAge1trueCanopy leaf area index, m2.m-2LAItruefn23WF*SLA*0.11WFSLA0.1trueTotal biomass, t.ha-1TotalWtruefn26WF+WS+WR1WFWSWRtrueMean stem biomass per tree, kg.tree-1AvStemMasstruefn27WS*1000/StemNo1WS1000StemNotrueBasic density, t.m-3Densitytruefn28if tRho != 0 then rho1+ (rho0 - rho1) * exp(-0.693147181 * (StandAge / tRho) ) else rho11rho1rho0rho1exp0.693147181StandAgetRhoneqtRho0rho1trueFraction of stem biomass as branch and bark, dimensionlessfracBBtruefn29if tBB != 0 then fracBB1 + (fracBB0 - fracBB1) * exp(-0.693147181 * (StandAge/ tBB) ) else fracBB11fracBB1fracBB0fracBB1exp0.693147181StandAgetBBneqtBB0fracBB1trueAge dependent modifier, dimensionlessfAgetruefn30if nAge==0 then 1 else (1/(1+(RelAge/rAge)^nAge))11eqnAge011RelAgerAgenAgetrueVPD dependent modifier, dimensionlessfVPDtruefn31exp(-1*CoeffCond*VPD)1exp1CoeffCondVPDtrueTemperature dependent modifier, dimensionlessfTtruefn3210leqTavTmin1geqTavTmax1TavTmin1Topt1Tmin1Tmax1TavTmax1Topt1Tmax1Topt1Topt1Tmin1trueFrost dependent modifier, dimensionlessfFrosttruefn331-kF*(FrostDays/30)11kFFrostDays30trueSoil water dependent modifier, dimensionlessfSWtruefn341/(1+((1-MoistRatio)/SWconst)^SWpower)1111MoistRatioSWconstSWpowertrueNutrition dependent modifier, dimensionlessfNutrtruefn35if fNn==0 then 1 else 1-(1-fN0)*(1-FR)^fNn11eqfNn011fN01FRfNntruePhysiological modifier of canopy conductance, dimensionlessPhysModtruefn36min(fVPD,fSW)*fAge1minfVPDfSWfAgetrueGross primary production in current period, t(DM).ha-1GPPtruefn37TranspScaleFactor*epsilon*RADint/100.01TranspScaleFactorepsilonRADint100.0trueNet primary production in current period, t(DM).ha-1NPPtruefn38GPP*y1GPPytrueTotal solar radiation intercepted by canopy, MJ.m-2.month-1RADinttruefn39RAD*lightIntcptn*CanCover1RADlightIntcptnCanCovertrueCanopy quantum efficiency after modifiers, mol.mol-1alphaCtruefn40alphaCx*fNutr*fT*fFrost*PhysMod1alphaCxfNutrfTfFrostPhysModtrueLight utilisation efficiency based on total biomass, g(DM).MJ-1epsilontruefn41gDM_mol*molPAR_MJ*alphaC1gDM_molmolPAR_MJalphaCtrueFR modifier of root biomass allocation, dimensionlessmtruefn44m0+(1-m0)*FR1m01m0FRtrueFraction of NPP allocated to roots, dimensionlesspRtruefn45pRx*pRn/(pRn+(pRx-pRn)*PhysMod*m)1pRxpRnpRnpRxpRnPhysModmtrueFraction of NPP allocated to stem, dimensionlesspStruefn46(1-pR)/(1+pFS)11pR1pFStrueFraction of NPP allocated to foliage, dimensionlesspFtruefn471-pR-pS11pRpStrueRatio of foliage to stem biomass allocation, dimensionlesspFStruefn48pfsConst*avDBH^pfsPower1pfsConstavDBHpfsPowertrueDensity independent mortality rate, month-1gammaNtruefn52if tgammaN != 0 then gammaN1 + (gammaN0 - gammaN1) * Exp(-0.693147181* (StandAge / tgammaN) ^ ngammaN) else gammaN11gammaN1gammaN0gammaN1Exp0.693147181StandAgetgammaNngammaNneqtgammaN0gammaN1trueFraction of rainfall intercepted by canopy, dimensionlessfracRainIntcptntruefn56if LAImaxIntcptn>0 then MaxIntcptn*min(1,LAI/LAImaxIntcptn)else MaxIntcptn1MaxIntcptnmin1LAILAImaxIntcptngtLAImaxIntcptn0MaxIntcptntrueCanopy conductance, m.s-1CanCondtruefn57MaxCond*PhysMod*min(1,LAI/LAIgcx)1MaxCondPhysModmin1LAILAIgcxtrueMonthly transpiration rate in current period, mmRainIntcptntruefn60Rain*fracRainIntcptn1RainfracRainIntcptntruefn1010.8-0.1*SoilClass10.80.1SoilClasstruefn10211-2*SoilClass1112SoilClasstruepfsPowertruefn103log(pFS20/pFS2)/log(20/2)1logpFS20pFS2log202truepfsConsttruefn104pFS2/2^pfsPower1pFS22pfsPowertruenetRadtruefn106Qa+Qb*(SolarRad*10^6/DayLength)1QaQbSolarRad106DayLengthtrueDayLengthtruefn109element([DayLength],month)1elementDayLengthmonthtruedefTermtruefn1101.2*2460000*(0.000622*VPD)*BLcond11.224600000.000622VPDBLcondtruedivtruefn113CanCond*(1+2.2)+BLcond1CanCond12.2BLcondtrueEtransptruefn114CanCond*(2.2*netRad+defTerm)/div1CanCond2.2netRaddefTermdivtrueMoistRatiotruefn115ASW/maxASW1ASWmaxASWtrueTransptruefn1161Etransp2460000DayLengthtruegammaFtruefn117if tgammaF*gammaF1==0 then gammaF1 else gammaF1*gammaF0/(gammaF0+(gammaF1-gammaF0)*exp(-1*kgammaF*StandAge))1gammaF1eqtgammaFgammaF10gammaF1gammaF0gammaF0gammaF1gammaF0exp1kgammaFStandAgetrueWFtruefn119414trueWStruefn120414trueWRtruefn121515trueWLtruefn122010truecd1truefn123NPP*pF1NPPpFtruecd2truefn124NPP*pS1NPPpStruecd3truefn125NPP*pR1NPPpRtrueStemNotruefn126143011430E. globus = 0.06trueCanopy quantum efficiencyalphaCxfile1E. globus = 0.095trueStem allometric constant, dimensionlessaWsfile1E. globus = 2.4trueStem allometric constant, dimensionlessnWsfile1trueHeighttruefn130aH*avDBH^nHB*StemNo^nHN1aHavDBHnHBStemNonHNtrueASWtruefn131010trueRelAgetruefn132StandAge/MaxAge1StandAgeMaxAgetruelightIntcptntruefn201(1-(exp(-1*k*LAI/CanCover)))11exp1kLAICanCoverE. globus = 24trueConversion of mol to gDMgDM_molfile1trueRADtruefn302SolarRad*daysInMonth1SolarRaddaysInMonthtruedaysInMonthtruefn303element([daysInMonth],month)1elementdaysInMonthmonthtruecd101truefn306Rain1Raintruecd102truefn307min(ASW,Transp+RainIntcptn)1minASWTranspRainIntcptntrueRaintruefn309element([Rain],month)1elementRainmonthtrueTranspScaleFactortruefn310EvapTransp/(Transp+RainIntcptn)1EvapTranspTranspRainIntcptntruecd106truefn313gammaR*WR1gammaRWRtruekgammaFtruefn31412*log(1+gammaF1/gammaF0)/tgammaF112log1gammaF1gammaF0tgammaFE. globus = 0.015trueRoot turnover rate per month, 1/monthgammaRfile1truedelStemstruefn318gammaN*StemNo/12/1001gammaNStemNo12100truecd108truefn322mF*delStems*(WF/StemNo)1mFdelStemsWFStemNotruecd109truefn323mS*delStems*(WS/StemNo)1mSdelStemsWSStemNotruecd110truefn324mR*delStems*(WR/StemNo)1mRdelStemsWRStemNotruefn325gammaF*WF1gammaFWFnode00201node00203node00210node00211node00753node00757node00759node00761node00763node00765node00766node00770node00771node00774node00775node00777truewhitenone12Monthly data0trueTminfile1trueTmaxfile1var6var7true1000SolarRadfile1true1000VPD1true1000FrostDaysfile1true240DayLengthfile1true1000Rainfile1var22true400daysInMonthfile1var25var26var28var29var31truemonthtruefn334floor(fmod(time(1),1)+1)intfloorfmodtime111E. globus = 0.2trueCanopy boundary layer conductance, assumed constantBLcondfile1i5i606i6i7i8truei613truei614i9i12i14i15truei191truei194truei195truei182truei192i16i17truei199truei200i21truei158truei159truei160truei161i22truei307truei308i23truei189i26truei178i27truei198i28truei166truei167truei168truei169truei197i29truei162truei163truei164truei165i30truei301truei302i31truei205truei206i32truei201truei202truei203truei204i33truei217truei218i34truei212truei213i35truei214truei215truei216i36truei304truei305truei306i37i38truei421truei422truei419i39truei416i40truei402truei403truei404truei405truei420i41truei406truei409i44truei441truei442i45truei446truei447truei448truei449i46truei450i47truei452truei453truei451i48truei445i52truei476truei477truei478truei479truei506i56truei429truei432truei433i57truei113truei114truei115truei116i60truei430i101truei102i103truei104truei444i105truei107truei110truei443i106truei109truei111i117truei122truei123truei539truei121i120i125truei128i130truei545i133truei134truei135truei136truei548truei211i138truei208i139truei140truei141i142truei143truei144truei145truei171truei155truei190i148truei156truei174truei196i149truei157i150truelossWFi151trueincrWFi152truei458truei461trueincrWSi153truei459truei462trueincrWRi154truei460truei463truei175truei183truei193i172truei401truei179truei180i181truei184truei185truei186truei187truei188truei210i209truei303i219truei220truei221truei415i309truei310truei311truei312truei410truei414i408truei412truei413i411trueRainInflowi423trueEvapTranspi424truei434truei435truei538truei428truei431i427truei440i436truei437truei438truei439truelossWRi456truei473truei468i464truei465truei466truei467truei472i480truei481truei507truemort_WFi491truei494truei497truei500truei505truemort_WSi492truei495truei498truei501truei504truefmort_WRi493truei496truei499truei502truei503i508truei509truei510truei518truei522truei524truei530truei533truei541truei608truei610truei607truei609truei517truei521truei523truei529truei532truei540truei526truei527truei528truei531truei534truei542truei611truei612i525truei546truei547