!! Event Monitor Program to calculate Spruce Beetle Risk Rating, based upon Schmid and Frye 1976
!! (Stand Ratings for Spruce Beetles.  RMF&RES, USDA FS, Research Note RM309)
!!
!! Physiographic location is one of the risk rating system's criteria.  This location information is
!! not available in the Event Monitor.  Thus, users need to identify these stands and then place
!! them into their own "groups" in FVS/Suppose and attach the appropriate Event Monitor addfile to
!! the appropriate groups.
!!
!! This Event Monitor application is for use with creek bottom stands.
!! A separate Event Monitor addfile should be used for non-creek bottom stands.
!!
!! The system:
!!
!! RISK FACTOR
!!    VALUE   PHYSIOGRAPHIC LOCATION   DBH(of trees>10")  STAND BA    %SPRUCE
!!      3        Creek bottoms           >16"              >150         >65%
!!      2        SI(100)>= 80            >12-16"           >100-150     >50-65%
!!      1        SI(100) < 80             10-12"           <=100        <=50%
!!
!! Total risk is then determined by summing the risk factor values for the 4 criteria.
!!
!! Note: SI in table above is on 100-year basis.
!! Note: since this addfile is for use with creek bottom stands, site index is not used!!
!!
!! ************************************************************************************************
!!
!! This EM Addfile written 11/01 by A J McMahan (INTECS International, Inc.)
!! for the Forest Health Technology Enterprise Team, USDA Forest Service,
!! 2150 Centre Ave., Bldg. A, Ft. Collins, CO  80526-1891.
!!
!! Contact:   Eric Smith     970 295-5841 (email:  elsmith@fs.fed.us)
!!
!! Slightly revised 8/02 AJM
!!
!!*************************************************************************************************
!!
!! Variable Names:
!!
!! SB_HAZRD = Spruce Beetle Hazard Rating
!!   Note: The Spruce Beetle Hazard Rating equals the sum of the 4 component rating factors:
!!   (SBF1, SBF2, SBF3, SBF4), ** UNLESS ** : (1) there is no spruce in the stand with a dbh >= 10"
!!   (in which case the rating will be "1"); or (2) less than 5% of stand BA is spruce, in which
!!   case the overall rating is zero.
!!
!! QMD_SP10 = QMD OF SPRUCE GREATER THAN (or =) 10"
!! TOT_BA   = TOTAL STAND BA
!! PER_SPRC = PERCENT OF STND_BA THAT IS SPRUCE (of all trees >2.6"dbh)
!!
!! Note that some FVS variants (EM, KT, NI) recognize "S"-not "ES"-for Engelmann spruce; hence
!! if you are using any of these variants, you will need to edit the SPMCDBH functions below!
!!
***************************************************************************************************



!!!!!!!!!!!!!!!!!!!!
!!!  INITIALIZE  !!!
!!!!!!!!!!!!!!!!!!!!

COMPUTE           0
SBF1=-1
SBF2=-1
SBF3=-1
SBF4=-1
SB_HAZRD=-1

QMD_SP10 = SPMCDBH(5,ES,0,10,999)
TOT_BA   = SPMCDBH(2,0,0)
PER_SPRC = SPMCDBH(2,ES,0,2.6,999) / SPMCDBH(2,0,0,2.6,999)

END

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!!!!!!!!!!!!!!!!!!!!
!!!  FIRST FACTOR  !!!
!!!!!!!!!!!!!!!!!!!!!!

!! assume this IS a creek bottom stand !
!! (if it is a not creek bottom stand, a different *.kcp file should be attached, not this one!).
!! Since we assume a creek bottom, Site Index won't be used.

COMPUTE           0
SBF1=3
END


!!!!!!!!!!!!!!!!!!!!!!!
!!!  SECOND FACTOR  !!!
!!!!!!!!!!!!!!!!!!!!!!!

IF
QMD_SP10 GT 16
THEN
COMPUTE           0
SBF2 = 3
END
ENDIF

IF
QMD_SP10 GT 12 AND QMD_SP10 LE 16
THEN
COMPUTE           0
SBF2 = 2
END
ENDIF

IF
QMD_SP10 LE 12 AND QMD_SP10 GE 10
THEN
COMPUTE           0
SBF2 = 1
END
ENDIF

IF
QMD_SP10 LT 10
THEN
COMPUTE           0
SBF2 = 0
END
ENDIF

!!!!!!!!!!!!!!!!!!!!!!
!!!  THIRD FACTOR  !!!
!!!!!!!!!!!!!!!!!!!!!!

IF
TOT_BA GT 150
THEN
COMPUTE           0
SBF3 = 3
END
ENDIF

IF
TOT_BA GT 100 AND TOT_BA LE 150
THEN
COMPUTE           0
SBF3 = 2
END
ENDIF

IF
TOT_BA LE 100
THEN
COMPUTE           0
SBF3 = 1
END
ENDIF

!!!!!!!!!!!!!!!!!!!!!!!
!!!  FOURTH FACTOR  !!!
!!!!!!!!!!!!!!!!!!!!!!!

IF
PER_SPRC GT 0.65
THEN
COMPUTE           0
SBF4 = 3
END
ENDIF

IF
PER_SPRC GT 0.50 AND PER_SPRC LE 0.65
THEN
COMPUTE           0
SBF4 = 2
END
ENDIF

IF
PER_SPRC LE 0.50 AND PER_SPRC GE 0.05
THEN
COMPUTE           0
SBF4 = 1
END
ENDIF

IF
PER_SPRC LT 0.05
THEN
COMPUTE
SBF4 = 0
END
ENDIF

!!!!!!!!!!!!!!!!!!!!
!!! RISK RATING  !!!
!!!!!!!!!!!!!!!!!!!!

IF
SBF1 NE -1 AND SBF2 NE -1 AND SBF3 NE -1 AND SBF4 NE -1 AND &
     SB_HAZRD EQ -1
THEN
COMPUTE           0
SB_HAZRD = SBF1 + SBF2 + SBF3 + SBF4
END
ENDIF


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!  CONDITION: IF NO SPRUCE GREATER THAN 10 INCHES   !!!
!!!!             THEN RISK IS VERY LOW                 !!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

IF
SB_HAZRD GT -1 AND SBF2 EQ 0 AND SBF4 NE 0
THEN
COMPUTE
SB_HAZRD = 1
END
ENDIF

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!  CONDITION: IF SPRUCE BA IS VERY LOW   !!!
!!!!             THEN RISK VALUE IS ZERO    !!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

IF
SB_HAZRD GT -1 AND SBF4 EQ 0
THEN
COMPUTE
SB_HAZRD = 0
END
ENDIF