fvs-python 0.2.3__py3-none-any.whl

pyfvs/cfg/sn_diameter_growth_coefficients.json

@@ -0,0 +1,300 @@
+{
+  "metadata": {
+    "title": "Large Tree Diameter Growth Coefficients for FVS Southern Variant",
+    "description": "Diameter growth model (DDS - change in squared diameter inside bark) coefficients extracted from official USDA FVS source code",
+    "source": "https://github.com/USDAForestService/ForestVegetationSimulator/blob/master/sn/dgf.f",
+    "extracted_date": "2025-11-29",
+    "equation_reference": "FVS-SN Variant Overview Section 4.7",
+    "threshold_diameter": 3.0,
+    "threshold_description": "Trees >= 3.0 inches DBH use large tree diameter growth model"
+  },
+  "equation": {
+    "main_equation": "ln(DDS) = CONSPP + INTERC + LDBH*ln(D) + DBH2*D^2 + LCRWN*ln(CR) + HREL*RELHT + PLTB*BA + PNTBL*PBAL + [forest_type_terms] + [eco_unit_terms]",
+    "conspp_equation": "CONSPP = DGCON + COR, where DGCON = ISIO*SI + TANS*SLOPE + FCOS*SLOPE*cos(ASPECT) + FSIN*SLOPE*sin(ASPECT) + [eco_unit_terms]",
+    "variable_definitions": {
+      "DDS": "Change in squared diameter (inside bark), inches^2",
+      "D": "DBH (diameter at breast height), inches",
+      "CR": "Crown ratio as integer percentage (e.g., 50 for 50%)",
+      "RELHT": "Relative height = HT/AVH, capped at 1.5",
+      "BA": "Stand basal area, sq ft/acre",
+      "PBAL": "Point basal area in larger trees = PBA * (1 - PCT/100)",
+      "SI": "Site index, feet (base age 50)",
+      "SLOPE": "Slope as tangent (rise/run)",
+      "ASPECT": "Aspect in radians"
+    },
+    "bounds": {
+      "DDS_minimum": -9.21,
+      "RELHT_maximum": 1.5,
+      "BA_minimum": 25.0,
+      "CR_minimum": 25
+    }
+  },
+  "species_order": {
+    "description": "FVS species index order (1-based Fortran arrays)",
+    "indices": {
+      "1": "FR", "2": "JU", "3": "PI", "4": "PU", "5": "SP", "6": "SA",
+      "7": "SR", "8": "LL", "9": "TM", "10": "PP", "11": "PD", "12": "WP",
+      "13": "LP", "14": "VP", "15": "BY", "16": "PC", "17": "HM", "18": "FM",
+      "19": "BE", "20": "RM", "21": "SV", "22": "SM", "23": "BU", "24": "BB",
+      "25": "SB", "26": "AH", "27": "HI", "28": "CA", "29": "HB", "30": "RD"
+    }
+  },
+  "core_coefficients": {
+    "SP": {
+      "species_index": 5,
+      "common_name": "Shortleaf Pine",
+      "INTERC": -0.008942,
+      "LDBH": 1.238170,
+      "DBH2": -0.001170,
+      "LCRWN": 0.053076,
+      "HREL": 0.040334,
+      "ISIO": 0.004723,
+      "PLTB": -0.004394,
+      "PNTBL": -0.003271,
+      "TANS": -0.704687,
+      "FCOS": 0.127667,
+      "FSIN": 0.028391,
+      "observations": 26008
+    },
+    "SA": {
+      "species_index": 6,
+      "common_name": "Slash Pine",
+      "INTERC": -1.641698,
+      "LDBH": 1.461093,
+      "DBH2": -0.002530,
+      "LCRWN": 0.265872,
+      "HREL": 0.069104,
+      "ISIO": 0.006851,
+      "PLTB": -0.002939,
+      "PNTBL": -0.004873,
+      "TANS": -0.018479,
+      "FCOS": -0.193157,
+      "FSIN": -0.251016,
+      "observations": 18736
+    },
+    "LL": {
+      "species_index": 8,
+      "common_name": "Longleaf Pine",
+      "INTERC": -1.331052,
+      "LDBH": 1.098112,
+      "DBH2": -0.001834,
+      "LCRWN": 0.184512,
+      "HREL": 0.388018,
+      "ISIO": 0.008774,
+      "PLTB": -0.002182,
+      "PNTBL": -0.002898,
+      "TANS": 0.225213,
+      "FCOS": 0.086883,
+      "FSIN": 0.107445,
+      "observations": 9190
+    },
+    "LP": {
+      "species_index": 13,
+      "common_name": "Loblolly Pine",
+      "INTERC": 0.222214,
+      "LDBH": 1.163040,
+      "DBH2": -0.000863,
+      "LCRWN": 0.028483,
+      "HREL": 0.006935,
+      "ISIO": 0.005018,
+      "PLTB": -0.003408,
+      "PNTBL": -0.004184,
+      "TANS": -0.759347,
+      "FCOS": 0.185360,
+      "FSIN": -0.072842,
+      "observations": 92313
+    }
+  },
+  "forest_type_coefficients": {
+    "description": "Categorical variables applied based on forest type classification (IFORTP)",
+    "forest_type_groups": {
+      "LOHD": {
+        "name": "Lowland Hardwood",
+        "forest_types": [168, 508, 601, 602, 605, 606, 607, 608, 702, 703, 704, 705, 706, 708],
+        "coefficients": {
+          "SP": 0.106418,
+          "SA": 0.325861,
+          "LL": 0.048216,
+          "LP": 0.126441
+        }
+      },
+      "NOHD": {
+        "name": "No Hardwood (Pure Softwood)",
+        "forest_types": [701, 801, 805],
+        "coefficients": {
+          "SP": 0.455020,
+          "SA": 0.000000,
+          "LL": 0.000000,
+          "LP": -0.122163
+        }
+      },
+      "OKPN": {
+        "name": "Oak-Pine",
+        "forest_types": [165, 403, 404, 405, 406, 407, 409],
+        "coefficients": {
+          "SP": 0.017518,
+          "SA": 0.116235,
+          "LL": 0.088872,
+          "LP": 0.050835
+        }
+      },
+      "SFHP": {
+        "name": "Spruce-Fir/Hemlock-Pine",
+        "forest_types": [104, 105, 121, 124],
+        "coefficients": {
+          "SP": 0.000000,
+          "SA": 0.000000,
+          "LL": 0.000000,
+          "LP": 0.000000
+        }
+      },
+      "UPHD": {
+        "name": "Upland Hardwood",
+        "forest_types": [103, 167, 181, 401, 402, 506, 507, 511, 512, 513, 519, 520, 802, 803, 807, 809],
+        "coefficients": {
+          "SP": 0.066811,
+          "SA": 0.162020,
+          "LL": 0.086720,
+          "LP": 0.063669
+        }
+      },
+      "UPOK": {
+        "name": "Upland Oak",
+        "forest_types": [501, 502, 503, 504, 505, 510, 514, 515],
+        "coefficients": {
+          "SP": -0.040181,
+          "SA": 0.410684,
+          "LL": 0.106061,
+          "LP": -0.016885
+        }
+      },
+      "YLPN": {
+        "name": "Yellow Pine",
+        "forest_types": [141, 142, 161, 162, 163, 164, 166],
+        "coefficients": {
+          "SP": 0.000000,
+          "SA": 0.000000,
+          "LL": 0.000000,
+          "LP": 0.000000
+        }
+      }
+    },
+    "plant_effect": {
+      "description": "Applied when MANAGD keyword is set (plantation effect)",
+      "coefficients": {
+        "SP": 0.000000,
+        "SA": 0.227572,
+        "LL": 0.110751,
+        "LP": 0.245669
+      }
+    }
+  },
+  "ecological_unit_coefficients": {
+    "description": "Categorical variables based on ecological unit classification (Bailey's ecoregion codes)",
+    "provinces": {
+      "PM221": {
+        "name": "Mountain Province 221",
+        "description": "Mountainous terrain in province 221",
+        "coefficients": {"SP": -0.569409, "SA": 0.000000, "LL": 0.000000, "LP": -0.069716}
+      },
+      "PM222": {
+        "name": "Mountain Province 222/223",
+        "description": "Mountainous terrain in provinces 222 or 223",
+        "coefficients": {"SP": -0.252741, "SA": 0.000000, "LL": 0.000000, "LP": 0.581967}
+      },
+      "PM231": {
+        "name": "Mountain Province 231",
+        "description": "Mountainous terrain in province 231",
+        "coefficients": {"SP": -0.265699, "SA": 0.000000, "LL": 0.000000, "LP": 0.790149}
+      },
+      "P221": {
+        "name": "Province 221 (non-mountain)",
+        "description": "Coastal plains/piedmont in province 221",
+        "coefficients": {"SP": -0.694484, "SA": 0.000000, "LL": 0.000000, "LP": -0.584818}
+      },
+      "P222": {
+        "name": "Province 222/223",
+        "description": "Province 222 or 223",
+        "coefficients": {"SP": -0.285112, "SA": 0.000000, "LL": -0.034553, "LP": -0.364073}
+      },
+      "S231T": {
+        "name": "Section 231 (Transition)",
+        "description": "Province 231 transition sections (A,C,D,I)",
+        "coefficients": {"SP": -0.504565, "SA": -0.025549, "LL": -0.175073, "LP": -0.183317}
+      },
+      "S231L": {
+        "name": "Section 231 (Lowland)",
+        "description": "Province 231 lowland sections (B,E,F,G)",
+        "coefficients": {"SP": 0.000000, "SA": 0.324111, "LL": -0.067793, "LP": 0.256273}
+      },
+      "P232": {
+        "name": "Province 232",
+        "description": "Outer Coastal Plain Mixed Forest Province",
+        "coefficients": {"SP": -0.113258, "SA": 0.000000, "LL": 0.000000, "LP": 0.000000}
+      },
+      "P234": {
+        "name": "Province 234",
+        "description": "Lower Mississippi Riverine Forest Province",
+        "coefficients": {"SP": 0.114097, "SA": 0.306793, "LL": 0.123262, "LP": 0.281790}
+      },
+      "P255": {
+        "name": "Province 255",
+        "description": "Prairie Parkland Province",
+        "coefficients": {"SP": 0.092458, "SA": 0.000000, "LL": 0.000000, "LP": 0.274618}
+      },
+      "P411": {
+        "name": "Province 411",
+        "description": "Everglades Province",
+        "coefficients": {"SP": 0.000000, "SA": -0.342293, "LL": 0.000000, "LP": 0.000000}
+      }
+    }
+  },
+  "fort_bragg_equations": {
+    "description": "Special diameter growth equations for Fort Bragg forest (IFOR=20)",
+    "applies_to_species": ["LL", "LP"],
+    "LL": {
+      "species_index": 8,
+      "equation": "DG5 = (D*BARK) * ((-0.4553*(0.09737-exp(-0.2428*D))) + 0.05574*(CR/100) - 0.0002965*BA - 0.00002481*PBA - 0.001192*((PCT/100)^(-0.9663)) + 0.0010110*SI - 0.007711*RELHT)",
+      "coefficients": {
+        "dbh_exp_a": -0.4553,
+        "dbh_exp_b": 0.09737,
+        "dbh_exp_c": -0.2428,
+        "crown_ratio": 0.05574,
+        "basal_area": -0.0002965,
+        "point_basal_area": -0.00002481,
+        "percentile_a": -0.001192,
+        "percentile_b": -0.9663,
+        "site_index": 0.0010110,
+        "relative_height": -0.007711
+      },
+      "observations": 2056
+    },
+    "LP": {
+      "species_index": 13,
+      "equation": "DG5 = (D*BARK) * ((-0.3428*(-0.1741-exp(-0.1328*D))) + 0.1145*(CR/100) - 0.0001682*BA - 0.00003978*PBA - 0.159400*((PCT/100)^(-0.1299)) + 0.0006204*SI + 0.02474*RELHT)",
+      "coefficients": {
+        "dbh_exp_a": -0.3428,
+        "dbh_exp_b": -0.1741,
+        "dbh_exp_c": -0.1328,
+        "crown_ratio": 0.1145,
+        "basal_area": -0.0001682,
+        "point_basal_area": -0.00003978,
+        "percentile_a": -0.159400,
+        "percentile_b": -0.1299,
+        "site_index": 0.0006204,
+        "relative_height": 0.02474
+      },
+      "observations": 689
+    },
+    "conversion_note": "DDS = ln((DG5*(2.0*D*BARK+DG5))) where DG5 >= 0.01"
+  },
+  "implementation_notes": {
+    "array_indexing": "FVS uses 1-based Fortran arrays. Subtract 1 for Python 0-based indexing.",
+    "crown_ratio_units": "Crown ratio (ICR) is stored as integer percentage (0-100), not proportion",
+    "aspect_units": "Aspect is in radians for cos/sin calculations",
+    "slope_units": "Slope is tangent (rise/run), not percent or degrees",
+    "dds_to_diameter_growth": "After calculating DDS, convert back: DG = sqrt(D^2 + exp(DDS)) - D",
+    "bark_ratio": "All diameter growth is inside bark. Apply bark ratio to convert to outside bark.",
+    "calibration": "DGCON includes site-specific calibration factor COR from DGDRIV"
+  }
+}

pyfvs/cfg/sn_height_diameter_coefficients.json

@@ -0,0 +1,97 @@
+{
+  "FR": {"P2": 2163.946776, "P3": 6.26880851, "P4": -0.2161439, "Dbw": 0.1, "Wykoff_B1": 4.5084, "Wykoff_B2": -6.0116},
+  "JU": {"P2": 212.7932729, "P3": 3.47154903, "P4": -0.3258523, "Dbw": 0.3, "Wykoff_B1": 4.0374, "Wykoff_B2": -4.2964},
+  "PI": {"P2": 2163.946776, "P3": 6.26880851, "P4": -0.2161439, "Dbw": 0.2, "Wykoff_B1": 4.5084, "Wykoff_B2": -6.0116},
+  "PU": {"P2": 3919.995225, "P3": 6.87312726, "P4": -0.19063343, "Dbw": 0.5, "Wykoff_B1": 4.2899, "Wykoff_B2": -4.1019},
+  "SP": {"P2": 444.0921666, "P3": 4.11876312, "P4": -0.30617043, "Dbw": 0.5, "Wykoff_B1": 4.6271, "Wykoff_B2": -6.4095},
+  "SP*": {"P2": 444.0921666, "P3": 4.11876312, "P4": -0.30617043, "Dbw": 0.5, "Wykoff_B1": 4.705, "Wykoff_B2": -7.904},
+  "SA": {"P2": 1087.101439, "P3": 5.10450596, "P4": -0.24284896, "Dbw": 0.5, "Wykoff_B1": 4.6561, "Wykoff_B2": -6.2258},
+  "SA*": {"P2": 110.3, "P3": 7.0670, "P4": -1.0420, "Dbw": 0.5, "Wykoff_B1": 4.787, "Wykoff_B2": -8.015},
+  "SR": {"P2": 333.3145742, "P3": 4.13108244, "P4": -0.37092539, "Dbw": 0.5, "Wykoff_B1": 4.7258, "Wykoff_B2": -6.7703},
+  "LL": {"P2": 98.56082813, "P3": 3.89930709, "P4": -0.86730393, "Dbw": 0.5, "Wykoff_B1": 4.5991, "Wykoff_B2": -5.9111},
+  "LL*": {"P2": 114.6, "P3": 4.1840, "P4": -0.6940, "Dbw": 0.5, "Wykoff_B1": 4.562, "Wykoff_B2": -7.314},
+  "TM": {"P2": 691.5411919, "P3": 4.19801014, "P4": -0.1856823, "Dbw": 0.5, "Wykoff_B1": 4.2139, "Wykoff_B2": -4.5419},
+  "PP": {"P2": 208.7773185, "P3": 3.72806565, "P4": -0.410875, "Dbw": 0.5, "Wykoff_B1": 4.3898, "Wykoff_B2": -5.7183},
+  "PD": {"P2": 142.7468108, "P3": 3.97260802, "P4": -0.5870983, "Dbw": 0.5, "Wykoff_B1": 4.5457, "Wykoff_B2": -6.8000},
+  "PD*": {"P2": 623.9, "P3": 4.7396, "P4": -0.2763, "Dbw": 0.5, "Wykoff_B1": 4.806, "Wykoff_B2": -9.573},
+  "WP": {"P2": 2108.844224, "P3": 5.65948135, "P4": -0.18563136, "Dbw": 0.4, "Wykoff_B1": 4.6090, "Wykoff_B2": -6.1896},
+  "LP": {"P2": 243.860648, "P3": 4.28460566, "P4": -0.47130185, "Dbw": 0.5, "Wykoff_B1": 4.6897, "Wykoff_B2": -6.8801},
+  "LP*": {"P2": 184.3, "P3": 4.2660, "P4": -0.5496, "Dbw": 0.5, "Wykoff_B1": 4.79, "Wykoff_B2": -8.5},
+  "VP": {"P2": 926.1802712, "P3": 4.46209203, "P4": -0.20053974, "Dbw": 0.5, "Wykoff_B1": 4.4718, "Wykoff_B2": -5.0078},
+  "BY": {"P2": 119.5749091, "P3": 4.13535453, "P4": -0.79625456, "Dbw": 0.2, "Wykoff_B1": 4.6171, "Wykoff_B2": -6.2684},
+  "PC": {"P2": 162.6505825, "P3": 3.20796415, "P4": -0.47880203, "Dbw": 0.2, "Wykoff_B1": 4.4603, "Wykoff_B2": -5.0577},
+  "HM": {"P2": 266.4562239, "P3": 3.99313675, "P4": -0.38600287, "Dbw": 0.1, "Wykoff_B1": 4.5084, "Wykoff_B2": -6.0116},
+  "FM": {"P2": 603.6736175, "P3": 3.9896005, "P4": -0.21651785, "Dbw": 0.2, "Wykoff_B1": 4.3164, "Wykoff_B2": -4.0582},
+  "BE": {"P2": 287.9445676, "P3": 3.27674704, "P4": -0.26617485, "Dbw": 0.2, "Wykoff_B1": 4.2378, "Wykoff_B2": -4.1080},
+  "RM": {"P2": 268.5564351, "P3": 3.11432843, "P4": -0.29411156, "Dbw": 0.2, "Wykoff_B1": 4.3379, "Wykoff_B2": -3.8214},
+  "SV": {"P2": 80.51179925, "P3": 26.98331005, "P4": -2.02202808, "Dbw": 0.2, "Wykoff_B1": 4.5991, "Wykoff_B2": -6.6706},
+  "SM": {"P2": 209.8555358, "P3": 2.95281334, "P4": -0.36787496, "Dbw": 0.2, "Wykoff_B1": 4.4834, "Wykoff_B2": -4.5431},
+  "BU": {"P2": 630.9504602, "P3": 4.51086779, "P4": -0.26826208, "Dbw": 0.3, "Wykoff_B1": 4.5697, "Wykoff_B2": -5.7172},
+  "BB": {"P2": 170.5253403, "P3": 2.68833651, "P4": -0.40080716, "Dbw": 0.1, "Wykoff_B1": 4.4388, "Wykoff_B2": -4.0872},
+  "SB": {"P2": 68.92234069, "P3": 43.33832185, "P4": -2.44448482, "Dbw": 0.1, "Wykoff_B1": 4.4522, "Wykoff_B2": -4.5758},
+  "AH": {"P2": 628.0209077, "P3": 3.88103963, "P4": -0.15387585, "Dbw": 0.2, "Wykoff_B1": 3.8550, "Wykoff_B2": -2.6623},
+  "HI": {"P2": 337.6684758, "P3": 3.62726466, "P4": -0.32083172, "Dbw": 0.3, "Wykoff_B1": 4.5128, "Wykoff_B2": -4.9918},
+  "CA": {"P2": 190.9797059, "P3": 3.69278884, "P4": -0.52730469, "Dbw": 0.3, "Wykoff_B1": 4.9396, "Wykoff_B2": -8.1838},
+  "HB": {"P2": 484.7529797, "P3": 3.93933286, "P4": -0.25998833, "Dbw": 0.1, "Wykoff_B1": 4.4207, "Wykoff_B2": -5.1435},
+  "RD": {"P2": 103.1767713, "P3": 2.21695491, "P4": -0.3596216, "Dbw": 0.2, "Wykoff_B1": 3.7512, "Wykoff_B2": -2.5539},
+  "DW": {"P2": 863.0501053, "P3": 4.38560239, "P4": -0.14812185, "Dbw": 0.1, "Wykoff_B1": 3.7301, "Wykoff_B2": -2.7758},
+  "PS": {"P2": 488.9349192, "P3": 4.06503751, "P4": -0.27180547, "Dbw": 0.2, "Wykoff_B1": 4.4091, "Wykoff_B2": -4.8464},
+  "AB": {"P2": 526.1392688, "P3": 3.89232121, "P4": -0.22587084, "Dbw": 0.1, "Wykoff_B1": 4.4772, "Wykoff_B2": -4.7206},
+  "AS": {"P2": 251.4042514, "P3": 3.26919806, "P4": -0.35905996, "Dbw": 0.2, "Wykoff_B1": 4.4819, "Wykoff_B2": -4.5314},
+  "WA": {"P2": 91.35276617, "P3": 6.99605268, "P4": -1.22937669, "Dbw": 0.2, "Wykoff_B1": 4.5959, "Wykoff_B2": -6.4497},
+  "BA": {"P2": 178.9307637, "P3": 4.92861465, "P4": -0.63777014, "Dbw": 0.2, "Wykoff_B1": 4.6155, "Wykoff_B2": -6.2945},
+  "GA": {"P2": 404.9692122, "P3": 3.39019741, "P4": -0.255096, "Dbw": 0.2, "Wykoff_B1": 4.6155, "Wykoff_B2": -6.2945},
+  "HL": {"P2": 778.9356784, "P3": 4.20756452, "P4": -0.18734197, "Dbw": 0.1, "Wykoff_B1": 4.3734, "Wykoff_B2": -5.3135},
+  "LB": {"P2": 265.7422693, "P3": 3.59041788, "P4": -0.35232417, "Dbw": 0.2, "Wykoff_B1": 4.4009, "Wykoff_B2": -5.0560},
+  "HA": {"P2": 2620.585492, "P3": 5.84993689, "P4": -0.18030935, "Dbw": 0.2, "Wykoff_B1": 4.4931, "Wykoff_B2": -4.6501},
+  "HY": {"P2": 1467.643523, "P3": 5.33438509, "P4": -0.17395792, "Dbw": 0.1, "Wykoff_B1": 4.0151, "Wykoff_B2": -4.3314},
+  "BN": {"P2": 285.8797853, "P3": 3.52138815, "P4": -0.3193688, "Dbw": 0.3, "Wykoff_B1": 4.5018, "Wykoff_B2": -5.6123},
+  "WN": {"P2": 93.71042027, "P3": 3.6575094, "P4": -0.88246833, "Dbw": 0.4, "Wykoff_B1": 4.5018, "Wykoff_B2": -5.6123},
+  "SU": {"P2": 290.90548, "P3": 3.6239536, "P4": -0.3720123, "Dbw": 0.2, "Wykoff_B1": 4.5920, "Wykoff_B2": -5.1719},
+  "YP": {"P2": 625.7696614, "P3": 3.87320571, "P4": -0.23349496, "Dbw": 0.2, "Wykoff_B1": 4.6892, "Wykoff_B2": -4.9605},
+  "MG": {"P2": 585.6609078, "P3": 3.41972033, "P4": -0.17661706, "Dbw": 0.2, "Wykoff_B1": 4.4004, "Wykoff_B2": -4.7519},
+  "CT": {"P2": 660.1996521, "P3": 3.92077102, "P4": -0.21124354, "Dbw": 0.2, "Wykoff_B1": 4.6067, "Wykoff_B2": -5.2030},
+  "MS": {"P2": 139.3315132, "P3": 2.89981329, "P4": -0.48514023, "Dbw": 0.2, "Wykoff_B1": 4.4004, "Wykoff_B2": -4.7519},
+  "MV": {"P2": 184.1931837, "P3": 2.84569124, "P4": -0.36952511, "Dbw": 0.2, "Wykoff_B1": 4.3609, "Wykoff_B2": -4.1423},
+  "ML": {"P2": 366.4744742, "P3": 2.8733336, "P4": -0.1819814, "Dbw": 0.2, "Wykoff_B1": 4.4004, "Wykoff_B2": -4.7519},
+  "AP": {"P2": 574.0200612, "P3": 3.86373895, "P4": -0.16318776, "Dbw": 0.2, "Wykoff_B1": 3.9678, "Wykoff_B2": -3.2510},
+  "MB": {"P2": 750.1823388, "P3": 4.14262749, "P4": -0.15940723, "Dbw": 0.2, "Wykoff_B1": 3.9613, "Wykoff_B2": -3.1993},
+  "WT": {"P2": 163.9728054, "P3": 2.76819717, "P4": -0.44098009, "Dbw": 0.2, "Wykoff_B1": 4.4330, "Wykoff_B2": -4.5383},
+  "BG": {"P2": 319.9788466, "P3": 3.67313408, "P4": -0.30651323, "Dbw": 0.2, "Wykoff_B1": 4.3802, "Wykoff_B2": -4.7903},
+  "TS": {"P2": 252.3566527, "P3": 3.24398683, "P4": -0.33343129, "Dbw": 0.2, "Wykoff_B1": 4.4334, "Wykoff_B2": -4.5709},
+  "HH": {"P2": 109.7324294, "P3": 2.25025802, "P4": -0.41297463, "Dbw": 0.2, "Wykoff_B1": 4.0322, "Wykoff_B2": -3.0833},
+  "SD": {"P2": 690.4917743, "P3": 4.15983216, "P4": -0.18613455, "Dbw": 0.2, "Wykoff_B1": 4.1352, "Wykoff_B2": -3.7450},
+  "RA": {"P2": 257.0532628, "P3": 3.4047448, "P4": -0.30291274, "Dbw": 0.2, "Wykoff_B1": 4.0965, "Wykoff_B2": -3.9250},
+  "SY": {"P2": 644.3567687, "P3": 3.92045786, "P4": -0.21444786, "Dbw": 0.1, "Wykoff_B1": 4.6355, "Wykoff_B2": -5.2776},
+  "CW": {"P2": 190.9797059, "P3": 3.69278884, "P4": -0.52730469, "Dbw": 0.1, "Wykoff_B1": 4.9396, "Wykoff_B2": -8.1838},
+  "BT": {"P2": 66.6488871, "P3": 135.4825559, "P4": -2.88622709, "Dbw": 0.2, "Wykoff_B1": 4.9396, "Wykoff_B2": -8.1838},
+  "BC": {"P2": 364.0247807, "P3": 3.55987361, "P4": -0.27263121, "Dbw": 0.1, "Wykoff_B1": 4.3286, "Wykoff_B2": -4.0922},
+  "WO": {"P2": 170.1330787, "P3": 3.27815866, "P4": -0.48744214, "Dbw": 0.2, "Wykoff_B1": 4.5463, "Wykoff_B2": -5.2287},
+  "SO": {"P2": 196.0564703, "P3": 3.0067167, "P4": -0.38499624, "Dbw": 0.2, "Wykoff_B1": 4.5225, "Wykoff_B2": -4.9401},
+  "SK": {"P2": 150.4300023, "P3": 3.13270999, "P4": -0.49925872, "Dbw": 0.1, "Wykoff_B1": 4.5142, "Wykoff_B2": -5.2205},
+  "CB": {"P2": 182.6306309, "P3": 3.12897883, "P4": -0.46391125, "Dbw": 0.1, "Wykoff_B1": 4.7342, "Wykoff_B2": -6.2674},
+  "TO": {"P2": 2137.575644, "P3": 5.80907868, "P4": -0.15590506, "Dbw": 0.2, "Wykoff_B1": 3.9365, "Wykoff_B2": -4.4599},
+  "LK": {"P2": 208.2300233, "P3": 3.13834277, "P4": -0.37158262, "Dbw": 0.1, "Wykoff_B1": 4.4375, "Wykoff_B2": -4.6654},
+  "OV": {"P2": 184.0856396, "P3": 3.49535241, "P4": -0.46211544, "Dbw": 0.2, "Wykoff_B1": 4.5710, "Wykoff_B2": -6.0922},
+  "BJ": {"P2": 157.4828626, "P3": 3.38919504, "P4": -0.39151499, "Dbw": 0.2, "Wykoff_B1": 3.9191, "Wykoff_B2": -4.3503},
+  "SN": {"P2": 281.3413276, "P3": 3.51695826, "P4": -0.3336282, "Dbw": 0.2, "Wykoff_B1": 4.6135, "Wykoff_B2": -5.7613},
+  "CK": {"P2": 72.7907469, "P3": 3.67065539, "P4": -1.09878979, "Dbw": 0.1, "Wykoff_B1": 4.3420, "Wykoff_B2": -5.1193},
+  "WK": {"P2": 470.0617193, "P3": 3.78892643, "P4": -0.25123824, "Dbw": 0.1, "Wykoff_B1": 4.5577, "Wykoff_B2": -4.9595},
+  "CO": {"P2": 94.54465221, "P3": 3.42034111, "P4": -0.818759, "Dbw": 0.2, "Wykoff_B1": 4.4618, "Wykoff_B2": -4.8786},
+  "RO": {"P2": 700.0636452, "P3": 4.10607389, "P4": -0.21392785, "Dbw": 0.2, "Wykoff_B1": 4.5202, "Wykoff_B2": -4.8896},
+  "QS": {"P2": 215.0009406, "P3": 3.14204012, "P4": -0.39067352, "Dbw": 0.1, "Wykoff_B1": 4.6106, "Wykoff_B2": -5.4380},
+  "PO": {"P2": 765.2907525, "P3": 4.22375114, "P4": -0.18974706, "Dbw": 0.1, "Wykoff_B1": 4.2496, "Wykoff_B2": -4.8061},
+  "BO": {"P2": 224.716279, "P3": 3.11648501, "P4": -0.35982064, "Dbw": 0.2, "Wykoff_B1": 4.4747, "Wykoff_B2": -4.8698},
+  "LO": {"P2": 153.9588254, "P3": 3.11348786, "P4": -0.38947124, "Dbw": 0.2, "Wykoff_B1": 4.2959, "Wykoff_B2": -5.3332},
+  "BK": {"P2": 880.2844971, "P3": 4.59642097, "P4": -0.21824277, "Dbw": 0.1, "Wykoff_B1": 4.4299, "Wykoff_B2": -4.9920},
+  "WI": {"P2": 408.2772475, "P3": 3.81808285, "P4": -0.27210505, "Dbw": 0.1, "Wykoff_B1": 4.4911, "Wykoff_B2": -5.7928},
+  "SS": {"P2": 755.1038099, "P3": 4.39496421, "P4": -0.21778831, "Dbw": 0.1, "Wykoff_B1": 4.3383, "Wykoff_B2": -4.5018},
+  "BD": {"P2": 293.5715132, "P3": 3.52261899, "P4": -0.35122247, "Dbw": 0.1, "Wykoff_B1": 4.582, "Wykoff_B2": -5.0903},
+  "EL": {"P2": 1005.80672, "P3": 4.6473994, "P4": -0.20336143, "Dbw": 0.1, "Wykoff_B1": 4.3744, "Wykoff_B2": -4.5257},
+  "WE": {"P2": 1001.672885, "P3": 4.57310438, "P4": -0.18898217, "Dbw": 0.1, "Wykoff_B1": 4.5992, "Wykoff_B2": -7.7428},
+  "AE": {"P2": 418.5941897, "P3": 3.17038578, "P4": -0.18964025, "Dbw": 0.1, "Wykoff_B1": 4.6008, "Wykoff_B2": -7.2732},
+  "RL": {"P2": 1337.547184, "P3": 4.48953501, "P4": -0.14749529, "Dbw": 0.1, "Wykoff_B1": 4.6238, "Wykoff_B2": -7.4847},
+  "OS": {"P2": 212.7932729, "P3": 3.47154903, "P4": -0.3258523, "Dbw": 0.3, "Wykoff_B1": 4.3898, "Wykoff_B2": -5.7183},
+  "OH": {"P2": 109.7324294, "P3": 2.25025802, "P4": -0.41297463, "Dbw": 0.2, "Wykoff_B1": 3.9392, "Wykoff_B2": -3.4279},
+  "OT": {"P2": 31021.35552, "P3": 8.3958757, "P4": -0.10372075, "Dbw": 0.2, "Wykoff_B1": 3.9089, "Wykoff_B2": -3.0149}
+} 

pyfvs/cfg/sn_large_tree_diameter_growth.json

@@ -0,0 +1,191 @@
+{
+  "metadata": {
+    "title": "Large Tree Height Growth for FVS Southern Variant",
+    "description": "Large tree height growth model following Wensel et al. (1987) approach with crown ratio and relative height modifiers",
+    "source": "FVS Southern Variant Documentation - Section 4.7.2",
+    "threshold_diameter": 3.0,
+    "threshold_description": "Trees equal to or larger than 3.0 inches DBH are considered 'large trees'",
+    "model_type": "Diameter growth driven (diameter growth predicted first, then height growth)",
+    "equations": {
+      "4.7.2.1": "HTG = POTHTG * (0.25 * HGMDCR + 0.75 * HGMDRH)",
+      "4.7.2.2": "HGMDCR = 100 * CR^3.0 * exp(-5.0*CR)",
+      "4.7.2.3": "FCTRKX = ((RHK / RHYXS)^(RHM - 1)) - 1",
+      "4.7.2.4": "FCTRRB = (-1.0 * RHR) / (1 - RHB)",
+      "4.7.2.5": "FCTRXB = RELHT^(1 - RHB) - RHXS^(1 - RHB)",
+      "4.7.2.6": "FCTRM = 1 / (1 - RHM)",
+      "4.7.2.7": "HGMDRH = RHK * (1 + FCTRKX * exp(FCTRRB * FCTRXB))^FCTRM"
+    },
+    "variable_definitions": {
+      "HTG": "Periodic height growth",
+      "POTHTG": "Potential periodic height growth (see section 4.6.1)",
+      "HGMDCR": "Crown ratio modifier (bounded to HGMDCR < 1.0)",
+      "HGMDRH": "Relative height modifier",
+      "CR": "Crown ratio expressed as a proportion",
+      "RELHT": "Subject tree's height relative to the 40 tallest trees in the stand",
+      "RHK": "Shade tolerance coefficient",
+      "RHYXS": "Shade tolerance coefficient",
+      "RHXS": "Shade tolerance coefficient",
+      "RHR": "Shade tolerance coefficient",
+      "RHB": "Shade tolerance coefficient",
+      "RHM": "Shade tolerance coefficient",
+      "FCTRKX": "Intermediate calculation factor",
+      "FCTRRB": "Intermediate calculation factor",
+      "FCTRXB": "Intermediate calculation factor",
+      "FCTRM": "Intermediate calculation factor"
+    }
+  },
+  "model_methodology": {
+    "approach": "Follows Wensel et al. (1987) where potential height growth is calculated for every tree and modified based on individual tree crown ratio and relative height in the stand",
+    "potential_height_growth": "Calculated using the methodology described in the small-tree height increment model (section 4.6.1)",
+    "crown_ratio_modifier": {
+      "function_type": "Hoerl's Special Function (HSF) form (Cuthbert and Wood 1971, p. 23)",
+      "range": "0.0 to 1.0",
+      "optimization": "Parameters chosen so that height growth is maximized for crown ratios between 45 and 75%"
+    },
+    "relative_height_modifier": {
+      "function_type": "Generalized Chapman-Richards function (Donnelly and Betters 1991, Donnelly and others 1992, and Pienaar and Turnbull 1973)",
+      "purpose": "Attenuate height growth based on relative height and shade tolerance",
+      "range": "0.0 to 1.0",
+      "behavior": "Decreases with decreasing relative height and species intolerance",
+      "asymptote": "Height growth reaches upper asymptote of 1.0 at relative height of 1.0 for intolerant species and 0.7 for tolerant species"
+    },
+    "weighting": {
+      "crown_ratio_weight": 0.25,
+      "relative_height_weight": 0.75
+    }
+  },
+  "shade_tolerance_coefficients": {
+    "very_tolerant": {
+      "RHR": 20,
+      "RHYXS": 0.20,
+      "RHM": 1.1,
+      "RHB": -1.10,
+      "RHXS": 0,
+      "RHK": 1
+    },
+    "tolerant": {
+      "RHR": 16,
+      "RHYXS": 0.15,
+      "RHM": 1.1,
+      "RHB": -1.20,
+      "RHXS": 0,
+      "RHK": 1
+    },
+    "intermediate": {
+      "RHR": 15,
+      "RHYXS": 0.10,
+      "RHM": 1.1,
+      "RHB": -1.45,
+      "RHXS": 0,
+      "RHK": 1
+    },
+    "intolerant": {
+      "RHR": 13,
+      "RHYXS": 0.05,
+      "RHM": 1.1,
+      "RHB": -1.60,
+      "RHXS": 0,
+      "RHK": 1
+    },
+    "very_intolerant": {
+      "RHR": 12,
+      "RHYXS": 0.01,
+      "RHM": 1.1,
+      "RHB": -1.60,
+      "RHXS": 0,
+      "RHK": 1
+    }
+  },
+  "species_shade_tolerance": {
+    "FR": "very_tolerant",
+    "JU": "intolerant",
+    "PI": "tolerant",
+    "PU": "intolerant",
+    "SP": "intolerant",
+    "SA": "intolerant",
+    "SR": "very_tolerant",
+    "LL": "intolerant",
+    "TM": "intolerant",
+    "PP": "intolerant",
+    "PD": "intolerant",
+    "WP": "intermediate",
+    "LP": "intolerant",
+    "VP": "intolerant",
+    "BY": "intermediate",
+    "PC": "intermediate",
+    "HM": "very_tolerant",
+    "FM": "tolerant",
+    "BE": "tolerant",
+    "RM": "tolerant",
+    "SV": "tolerant",
+    "SM": "very_tolerant",
+    "BU": "tolerant",
+    "BB": "intolerant",
+    "SB": "intolerant",
+    "AH": "very_tolerant",
+    "HI": "intermediate",
+    "CA": "intolerant",
+    "HB": "intermediate",
+    "RD": "tolerant",
+    "DW": "very_tolerant",
+    "PS": "very_tolerant",
+    "AB": "very_tolerant",
+    "AS": "tolerant",
+    "WA": "intolerant",
+    "BA": "intolerant",
+    "GA": "tolerant",
+    "HL": "intolerant",
+    "LB": "tolerant",
+    "HA": "tolerant",
+    "HY": "very_tolerant",
+    "BN": "intolerant",
+    "WN": "intolerant",
+    "SU": "intolerant",
+    "YP": "intolerant",
+    "MG": "tolerant",
+    "CT": "intermediate",
+    "MS": "tolerant",
+    "MV": "intermediate",
+    "ML": "tolerant",
+    "AP": "intolerant",
+    "MB": "tolerant",
+    "WT": "intolerant",
+    "BG": "tolerant",
+    "TS": "intolerant",
+    "HH": "tolerant",
+    "SD": "tolerant",
+    "RA": "tolerant",
+    "SY": "intermediate",
+    "CW": "very_intolerant",
+    "BT": "very_intolerant",
+    "BC": "intolerant",
+    "WO": "intermediate",
+    "SO": "very_intolerant",
+    "SK": "intermediate",
+    "CB": "intolerant",
+    "TO": "intolerant",
+    "LK": "tolerant",
+    "OV": "intermediate",
+    "BJ": "intolerant",
+    "SN": "intolerant",
+    "CK": "intolerant",
+    "WK": "intolerant",
+    "CO": "intermediate",
+    "RO": "intermediate",
+    "QS": "intolerant",
+    "PO": "intolerant",
+    "BO": "intermediate",
+    "LO": "intermediate",
+    "BK": "very_intolerant",
+    "WI": "very_intolerant",
+    "SS": "intolerant",
+    "BD": "tolerant",
+    "EL": "intermediate",
+    "WE": "tolerant",
+    "AE": "intermediate",
+    "RL": "tolerant",
+    "OS": "intermediate",
+    "OH": "intermediate",
+    "OT": "intermediate"
+  }
+}