VolumeLibrary/fia_vol_carib.for at master · FMSC-Measurements/VolumeLibrary · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
! This is from FIA package code NIMS_VOL_FIA_CARIB
! Created by YW 2018/09/18
! ---------------------------------------------------------------------
! Brandeis, T.J. etal 2005
! Equations for Merchantable Volume for Subtropical Moist and Wet Forests of Puerto Rico
! USDA Forest Service Southern Research Station Research Paper SRS-39
! CU000122 CV4 FUNCTION CARIB_BRANDEIS_2005(SPN IN NUMBER, DBH IN NUMBER, HT IN NUMBER)
! Brandeis, T.J. etal 2006
! Development of equations for predicting Puerto Rican subtropical dry forest biomass and volume
! Forest Ecology and Management 233 (2006) 133–142
! CU000123 CV4 FUNCTION CARIB_BRANDEIS_2006(SPN IN NUMBER, DBH IN NUMBER, HT IN NUMBER)
! Valid species code:
! NVEL EQUATION NUMBER:
! S99BRA****
      SUBROUTINE BRANDEIS_VOL(VOLEQ,DBHOB,HTTOT,VOL,ERRFLG)
      CHARACTER*10 VOLEQ
      INTEGER SPN, ERRFLG
      REAL DBHOB,HTTOT,DBH_CM,HT_MT,CVT_M,CV4_M,VOL(15)
      REAL A0,A1,B0,B1,B2,B3,C0,C1
      ERRFLG = 0
      VOL = 0.0
      IF(DBHOB.LT.5.0) RETURN
      READ(VOLEQ(7:10),'(I4)')SPN
      DBH_CM = DBHOB*2.54
      HT_MT = HTTOT*0.3048
      IF(SPN.EQ.854.OR.SPN.EQ.6008.OR.SPN.EQ.6284)THEN
!  --   Dry forests
!  --BUBU
        IF(SPN.EQ.6284)THEN
          A0 = 0.0001
          A1 = 0.000047
          C0 = -0.0697
          C1 = 0.000046
        ELSEIF(SPN.EQ.854)THEN
!  --BUSI
          A0 = 0.0370
          A1 = 0.000028
          C0 = 0.02890
! THE 0.2890 ON THE PUBLICATION PAPER SEEMS NOT RIGHT
          C1 = 0.000025
        ELSEIF(SPN.EQ.6008)THEN
!  --ACMA
          A0 = -0.0173
          A1 = 0.000047
          C0 = -0.0593
          C1 = 0.000046
        ENDIF
        CVT_M = A0 + A1 * (DBH_CM ** 2) * HT_MT
        CVT = CVT_M*35.314667
        CV4_M = C0 + C1 * (DBH_CM ** 2) * HT_MT
        CV4 = CV4_M*35.314667
        VOL(1) = CVT
        VOL(4) = CV4
      ELSE
!  -- Wet/moist/lower montane wet/rain forests
        IF(SPN.EQ.6114)THEN
!--ANIN
          B0 = -10.7353102290492
          B1 = 2.33132884613613
          B2 = 0.774680339670667
          B3 = -0.0074807863759593
        ELSEIF(SPN.EQ.6443)THEN
!--CESC9
          B0 = -13.5633594255738
          B1 = 3.5507061416754
          B2 = 0.787970783049223
          B3 = -0.0469524275398158
        ELSEIF(SPN.EQ.6728)THEN
!--COAL
          B0 = -11.9506641352761
          B1 = 3.06860415362186
          B2 = 0.576660209391705
          B3 = -0.029577504262509
        ELSEIF(SPN.EQ.7011)THEN
!--ERPO5
          B0 = -14.2743202553787
          B1 = 3.63718110738798
          B2 = 0.576660209391705
          B3 = -0.0326995133584038
        ELSEIF(SPN.EQ.7290)THEN
!--GUGU
          B0 = -12.5614979281204
          B1 = 3.05691245763579
          B2 = 0.748845391883614
          B3 = -0.0261745352292469
        ELSEIF(SPN.EQ.7470)THEN
!--INLA
          B0 = -12.7021831143844
          B1 = 2.83592579314056
          B2 = 1.01688921838352
          B3 = -0.0206949734897774
        ELSEIF(SPN.EQ.7474)THEN
!--INVE
          B0 = -10.8303561289069
          B1 = 2.50262939251883
          B2 = 0.731086851638329
          B3 = -0.015600143438074
        ELSEIF(SPN.EQ.885)THEN
!--MAIN3
          B0 = -10.1626005016519
          B1 = 1.92057572622235
          B2 = 1.08296645946787
          B3 = -0.00293548732354752
        ELSEIF(SPN.EQ.7994)THEN
!--OCLE
          B0 = -11.1068755264823
          B1 = 2.19819049034765
          B2 = 1.05004167825012
          B3 = -0.00256209197493238
        ELSEIF(SPN.EQ.8558)THEN
!--SCMO10
          B0 = -10.7075669936423
          B1 = 2.1613616649571
          B2 = 0.979394658095212
          B3 = -0.00322527623337223
        ELSEIF(SPN.EQ.8644)THEN
!--SPCA2
          B0 = -11.8315998001912
          B1 = 2.7089366261757
          B2 = 0.869934254242142
          B3 = -0.0166886450211204
        ELSEIF(SPN.EQ.8701)THEN
!--SYJA
          B0 = -10.5433370090313
          B1 = 2.08250274792873
          B2 = 0.893677401952679
          B3 = 0.0000773127731972417
        ELSEIF(SPN.EQ.8713)THEN
!--TAHE
          B0 = -11.5743019513462
          B1 = 2.49204956936753
          B2 = 0.89277350958636
          B3 = -0.00702587297575205
        ELSEIF(SPN.EQ.6001)THEN
!--ACAN4
          B0 = -12.0234809372091
          B1 = 2.75276356832771
          B2 = 0.88961180433259
          B3 = -0.0207220137274765
        ELSE
          ERRFLG = 6
          RETURN
        ENDIF
        CV4_M = EXP(B0+B1*LOG(DBH_CM)+B2*LOG(HT_MT)+B3*DBH_CM)
        CV4 = CV4_M*35.314667
        VOL(4) = CV4
      ENDIF
      RETURN
      END SUBROUTINE BRANDEIS_VOL