FORSKAL模型结构与原理课件.ppt

1、1/1/20231l程序流程图1/1/20232l功能：读取立地参数和运行参数l实现：从PARAM.DAT数据文件中读取下列参数lKCO:消光系数lPAT:斑块面积（hm2）lNPAT:斑块数lTIM:模拟步长lNTIM:模拟步长数lVINC:垂直积分步长（m）lDBHO:最小胸径（cm）lRHO:死亡函数的峭度参数lTREF:参照温度lFCO2:内外CO2浓度lTIMRET:干扰平均重复期lGAMMA:干扰率1/1/20233l功能：读取树种参数l实现：从PARAM.DAT中读取要模拟的树种数（NTAX），以及每个树种的下列参数lNAM:树种名lHMX:最大高（m）lHDS:胸径树高曲线起始

2、坡度（m/cm）lALP:半饱和点（umol/m2/s）lLCP:补偿点（umol/m2/s）lGSC:生长常数（cm2/m/yr）lEST:幼树更新速率（hm2/yr）lTDI:导致死亡增加的相对生长效率阈值lUMN:内禀死亡率（/yr）lUMX:受压死亡率（yr）lSPR:每株死树能够萌生幼树数量lSMN:具有萌生能力的死树最小胸径（cm）1/1/20234lLAC:幼叶起始叶面积D2值（m2/cm2）lLAF:边材周转率（/yr）lBCF:树干生物量转换系数（kg/cm2/m）lR:边材材积维持成本（/yr）lQ10:呼吸增加速率lTMIN:同化作用所需最低温度lTMAX:同化作用所需最

3、高温度lCCP:树种补偿点lDRI:最大耐受干旱指数lMINGDD:生长期度日最低值lMINTC:最冷月最低温度（）lMAXTC:最冷月的最高温（）lMINTW:最暖月的最低温（）lDORE:落叶或常绿树种标识，0落叶种，1常绿种1/1/20235l功能：读取控制参数l实现：从PARAM.DAT中读取INIT、IOUT、LPR、WPAT、START、SWCO2lINIT：控制模型从裸地开始模拟，还是以某个阶段开始l 若INIT=0，裸地l 若INIT0，从STAND.INI中读取各斑块上所有树木变量值lIOUT：控制存储树木变量文件的参数l 若IOUT0，则将最后一年各斑块模拟的结果存入NEW

4、STAND.INI中lLPR：运行结果输出年度间隔lWPAT：控制程序输出内容的参数l 若WPAT0，程序将每一个斑块的模拟结果都输入到数据文件中l 若WPAT1，程序只将所有斑块的模拟结果之和存入到数据文件中。lSTART：控制模拟结果输出的起始年代的参数l 若START10，则只将10年以后的模拟数据输出。l 若START大于总的模拟时间，则令START等于总模拟时间。lSWCO21/1/20236l功能：读取和计算环境数据l实现：打开数据文件CLIMATE.INP，读取运行期间各年度的气象数据，读取每个月的mtc（月均温）、mtr（月蒸发）、mpr（月降水量）、mcl（月平均云量），调用

5、Envsub模块，得到下列参数l每一年的gdd、GSDRI、M4DRI、GSINS、M4INS、TCOLD、TWARM、ADDCCP、tftmc、sapres1/1/20237l功能：计算各模拟年度的gdd、gsdri、m4dri、gsins、m4ins、sapres、tftmxl实现：调用daily模块，根据mtc、mpr、mcl的月数据推算出日数据1/1/20238l功能：计算影响树木更新和生长的乘数因子l实现：根据GETENV模块得到的各年GDD、GSDRI、M4DRI、GSINS、M4INS、TCOLD、TWARM、ADDCCP来计算当前所在步长时间内的年平均GDD、GSDRI、M4D

6、RI、GSINS、M4INS、TCOLD、TWARM、TCO2，将xGDD与MINGDD比较l xTCOLD与MINTC和MAXTC比较l xWARM与MINTC比较2)(1)(IDRIXGSDRTIDRMX1/1/20239l将DRMX与MRMX比较l如果SWCO21，则l否则SWCO21。)(3302221)(2ICCPFCOXTCOFCOIMXCO1/1/202310l功能：初始化每个树种的最大生长效率和最大起始径生长 l实现：计算某个斑块上总叶面积TLA，并计算l（1）计算斑块上LAIl l 生长地林下光照l l 4以上时期林地光照10000PATTLALAILAIKCOeXGSINS

7、GSLITELAIKCOeINSXMLITEM441/1/202311l（2）计算上述光强下每个树种的生长量l幼树叶面积l评价：由于叶面叶必然为0，所以对叶面积的计算没有必要，叶面积和光衰减等不必计算LCP(I)-ALP(I)GSLITEKCOLCP(I)-GSLITEKCOGSC(I)P2DBHOLAC(I)FA1/1/202312l功能：模拟某个斑块内数目的更新l实现：1、计算叶面积指数和林冠下光强l l l 2、CO2MX与ALP相结合l CO2ALP(I)=ALP(I)*CO2MX(I)10000PATTLALAILAIKCOeXGSINSGSLITELAIKCOeINSXMLITEM

8、441/1/202313l3、计算在林下和全光下的P值l若为落叶树种，则ll若为常绿树种，则LCP(I)-CO2ALP(I)GSLITEKCOLCP(I)-GSLITEKCOPFORET(I)LCP(I)-CO2ALP(I)XGSINSKCOLCP(I)-XGSINSKCOPFULLLCP(I)-CO2ALP(I)M4LITEKCOLCP(I)-M4LITEKCOPFORETLCP(I)-CO2ALP(I)M4INSKCOLCP(I)-XM4INSKCOPFULL1/1/202314l4、修改ESTlAMDEST(I)=EST(I)*GDDMX(I)*DRMX(I)*TCMX(I)*TWMX(I)*PMX(I)*XTFTMX(I)*TWARMX(I)*CO2MX(I)l5、调用PRAND，得到幼树形成数量NSAP，将这部分数据增加到数组DBH(J,K)，HBC(J,K)，AGE(J,K)，LEA(J,K)，并用数组TAX(J,K)来标明第K个个体的树种。将总个体数存到NTREES中。1/1/2023151/1/2023161/1/202317