《FORSKAL模型结构与原理课件.ppt》由会员分享，可在线阅读，更多相关《FORSKAL模型结构与原理课件.ppt（17页珍藏版）》请在三一办公上搜索。
1、FORSKAL模型结构与原理,10/9/2022,1,FORSKAL模型结构与原理10/9/20221,程序流程图,程序流程图,10/9/2022,2,程序流程图 程序流程图10/9/20222,LGTMOD模块,功能：读取立地参数和运行参数实现： 从PARAM.DAT数据文件中读取下列参数KCO: 消光系数PAT: 斑块面积（hm2）NPAT: 斑块数TIM: 模拟步长NTIM: 模拟步长数VINC: 垂直积分步长（m）DBHO: 最小胸径（cm）RHO: 死亡函数的峭度参数TREF: 参照温度FCO2: 内外CO2浓度TIMRET: 干扰平均重复期GAMMA: 干扰率,10/9/2022,
2、3,LGTMOD模块功能：读取立地参数和运行参数10/9/202,CGTSPE模块,功能：读取树种参数实现：从PARAM.DAT中读取要模拟的树种数（NTAX），以及每个树种的下列参数NAM: 树种名HMX: 最大高（m）HDS: 胸径树高曲线起始坡度（m/cm）ALP: 半饱和点（umol/m2/s）LCP: 补偿点（umol/m2/s）GSC: 生长常数（cm2/m/yr）EST: 幼树更新速率（hm2/yr）TDI: 导致死亡增加的相对生长效率阈值UMN: 内禀死亡率（/yr）UMX: 受压死亡率（yr）SPR: 每株死树能够萌生幼树数量SMN: 具有萌生能力的死树最小胸径（cm）,10
3、/9/2022,4,CGTSPE模块功能：读取树种参数10/9/20224,LAC: 幼叶起始叶面积D2值（m2/cm2）LAF: 边材周转率（/yr）BCF: 树干生物量转换系数（kg/cm2/m）R: 边材材积维持成本（/yr）Q10: 呼吸增加速率TMIN: 同化作用所需最低温度TMAX: 同化作用所需最高温度CCP: 树种补偿点DRI: 最大耐受干旱指数MINGDD: 生长期度日最低值MINTC: 最冷月最低温度（）MAXTC: 最冷月的最高温（）MINTW: 最暖月的最低温（）DORE: 落叶或常绿树种标识，0落叶种，1常绿种,10/9/2022,5,LAC: 幼叶起始叶面积D2值（
4、m2/cm2）10/9/2,CGTOUT模块,功能：读取控制参数实现：从PARAM.DAT中读取INIT、IOUT、LPR、WPAT、START、SWCO2INIT：控制模型从裸地开始模拟，还是以某个阶段开始 若INIT=0，裸地 若INIT0，从STAND.INI中读取各斑块上所有树木变量值IOUT：控制存储树木变量文件的参数 若IOUT0，则将最后一年各斑块模拟的结果存入NEWSTAND.INI中LPR：运行结果输出年度间隔WPAT：控制程序输出内容的参数 若WPAT0，程序将每一个斑块的模拟结果都输入到数据文件中 若WPAT1，程序只将所有斑块的模拟结果之和存入到数据文件中。 START
5、：控制模拟结果输出的起始年代的参数 若START10，则只将10年以后的模拟数据输出。 若START大于总的模拟时间，则令START等于总模拟时间。 SWCO2,10/9/2022,6,CGTOUT模块功能：读取控制参数10/9/20226,GETENV模块,功能：读取和计算环境数据实现：打开数据文件CLIMATE.INP，读取运行期间各年度的气象数据，读取每个月的mtc（月均温）、mtr（月蒸发）、mpr（月降水量）、mcl（月平均云量），调用Envsub模块，得到下列参数每一年的gdd、GSDRI、M4DRI、GSINS、M4INS、TCOLD、TWARM、ADDCCP、tftmc、sap
6、res,10/9/2022,7,GETENV模块功能：读取和计算环境数据10/9/20227,Envsub模块,功能：计算各模拟年度的gdd、gsdri、m4dri、gsins、m4ins、sapres、tftmx实现：调用daily模块，根据mtc、mpr、mcl的月数据推算出日数据,10/9/2022,8,Envsub模块功能：计算各模拟年度的gdd、gsdri、m,CLIMATE模块,功能：计算影响树木更新和生长的乘数因子实现：根据GETENV模块得到的各年GDD、GSDRI、M4DRI、GSINS、M4INS、TCOLD、TWARM、ADDCCP来计算当前所在步长时间内的年平均GDD、
7、GSDRI、M4DRI、GSINS、M4INS、TCOLD、TWARM、TCO2，将xGDD与MINGDD比较 xTCOLD与MINTC和MAXTC比较 xWARM与MINTC比较,10/9/2022,9,CLIMATE模块功能：计算影响树木更新和生长的乘数因子10,将DRMX与MRMX比较如果SWCO21，则否则SWCO21。,10/9/2022,10,将DRMX与MRMX比较10/9/202210,CINIT模块,功能：初始化每个树种的最大生长效率和最大起始径生长 实现：计算某个斑块上总叶面积TLA，并计算（1）计算斑块上LAI 生长地林下光照 4以上时期林地光照,10/9/2022,11
8、,CINIT模块功能：初始化每个树种的最大生长效率和最大起始径,（2）计算上述光强下每个树种的生长量幼树叶面积评价：由于叶面叶必然为0，所以对叶面积的计算没有必要，叶面积和光衰减等不必计算,10/9/2022,12,10/9/202212,CETBL模块,功能：模拟某个斑块内数目的更新实现：1、计算叶面积指数和林冠下光强 2、CO2MX与ALP相结合 CO2ALP(I)=ALP(I)*CO2MX(I),10/9/2022,13,CETBL模块功能：模拟某个斑块内数目的更新10/9/202,3、计算在林下和全光下的P值若为落叶树种，则若为常绿树种，则,10/9/2022,14,3、计算在林下和全
9、光下的P值10/9/202214,4、修改ESTAMDEST(I)=EST(I)*GDDMX(I)*DRMX(I)*TCMX(I)*TWMX(I)*PMX(I)*XTFTMX(I)*TWARMX(I)*CO2MX(I) 5、调用PRAND，得到幼树形成数量NSAP，将这部分数据增加到数组DBH(J, K)，HBC(J, K)，AGE(J, K)，LEA(J, K)，并用数组TAX(J, K)来标明第K个个体的树种。将总个体数存到NTREES中。,10/9/2022,15,4、修改EST10/9/202215,结果,10/9/2022,16,结果 10/9/202216,讨论,10/9/2022,17,讨论 10/9/202217,
