FORSKAL模型结构与原理课件

1、FORSKAL模型结构与原理10/10/20221FORSKAL模型结构与原理10/9/20221程序流程图 程序流程图10/10/20222程序流程图 程序流程图10/9/20222LGTMOD模块功能：读取立地参数和运行参数实现： 从PARAM.DAT数据文件中读取下列参数KCO: 消光系数PAT: 斑块面积（hm2）NPAT: 斑块数TIM: 模拟步长NTIM: 模拟步长数VINC: 垂直积分步长（m）DBHO: 最小胸径（cm）RHO: 死亡函数的峭度参数TREF: 参照温度FCO2: 内外CO2浓度TIMRET: 干扰平均重复期GAMMA: 干扰率10/10/20223LGTMOD模

2、块功能：读取立地参数和运行参数10/9/202CGTSPE模块功能：读取树种参数实现：从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/10/20224CGT

3、SPE模块功能：读取树种参数10/9/20224LAC: 幼叶起始叶面积D2值（m2/cm2）LAF: 边材周转率（/yr）BCF: 树干生物量转换系数（kg/cm2/m）R: 边材材积维持成本（/yr）Q10: 呼吸增加速率TMIN: 同化作用所需最低温度TMAX: 同化作用所需最高温度CCP: 树种补偿点DRI: 最大耐受干旱指数MINGDD: 生长期度日最低值MINTC: 最冷月最低温度（）MAXTC: 最冷月的最高温（）MINTW: 最暖月的最低温（）DORE: 落叶或常绿树种标识，0落叶种，1常绿种10/10/20225LAC: 幼叶起始叶面积D2值（m2/cm2）10/9/2CGT

4、OUT模块功能：读取控制参数实现：从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等于总模拟时间。 SWCO210/10/20226CGTOUT模块功能：读取控制参数10/9/20226GETENV模块功能：读取和计算环境数据实现：打开数据文件CLIMATE.INP，读取运行期间各年度的气象数据，读取每个月的mtc（月均温）、mtr（月蒸发）、mpr（月降水量）、mcl（月平均云量），调用Envsub模块，得到下列参数每一年的gdd、GSDRI、M4DRI、GSINS、M4INS、TCOLD、TWARM、ADDCCP、tftmc、sapres10/10/20227GETENV模块

6、功能：读取和计算环境数据10/9/20227Envsub模块功能：计算各模拟年度的gdd、gsdri、m4dri、gsins、m4ins、sapres、tftmx实现：调用daily模块，根据mtc、mpr、mcl的月数据推算出日数据10/10/20228Envsub模块功能：计算各模拟年度的gdd、gsdri、mCLIMATE模块功能：计算影响树木更新和生长的乘数因子实现：根据GETENV模块得到的各年GDD、GSDRI、M4DRI、GSINS、M4INS、TCOLD、TWARM、ADDCCP来计算当前所在步长时间内的年平均GDD、GSDRI、M4DRI、GSINS、M4INS、TCOLD、

7、TWARM、TCO2，将xGDD与MINGDD比较 xTCOLD与MINTC和MAXTC比较 xWARM与MINTC比较10/10/20229CLIMATE模块功能：计算影响树木更新和生长的乘数因子10将DRMX与MRMX比较如果SWCO21，则否则SWCO21。10/10/202210将DRMX与MRMX比较10/9/202210CINIT模块功能：初始化每个树种的最大生长效率和最大起始径生长 实现：计算某个斑块上总叶面积TLA，并计算（1）计算斑块上LAI 生长地林下光照 4以上时期林地光照10/10/202211CINIT模块功能：初始化每个树种的最大生长效率和最大起始径（2）计算上述光

8、强下每个树种的生长量幼树叶面积评价：由于叶面叶必然为0，所以对叶面积的计算没有必要，叶面积和光衰减等不必计算10/10/20221210/9/202212CETBL模块功能：模拟某个斑块内数目的更新实现：1、计算叶面积指数和林冠下光强 2、CO2MX与ALP相结合 CO2ALP(I)=ALP(I)*CO2MX(I) 10/10/202213CETBL模块功能：模拟某个斑块内数目的更新10/9/2023、计算在林下和全光下的P值若为落叶树种，则若为常绿树种，则10/10/2022143、计算在林下和全光下的P值10/9/2022144、修改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,