Marten
; Habitat model for marten (Martes americana)
; In RMLANDS workspace, use 10 800 RUNLOOP 'marten.txt'
; Full version includes roads, but for HRV there aren’t any roads, so we’ve simplified this model
gridpath = 'z:\RMLANDS\SJNF\hrv_a1\run1\elk' ; Location of zipped input grids (uses same grid as elk)
BLOCK 140 140 36 ; 1100x1100 blocks, with 69 cell edge buffer
covcond = READ 'covcond.asc' ; read grids
lri1 = (GET 'mcover.txt') CATEGORY covcond ; lri1 = cover type
SETWINDOW 50 10 ; logistic function determines window size
lri2 = (GET 'medge.txt') 'LOGISTIC' (50 10 .5) (50 10 .6) (50 10 .8) BADEDGE covcond ; lri2 = adverse edge effects
lra = lri1 * lri2 ; lra = product of lris
hrc = 892 ('MEAN' 'L') COMBINE lra ; hrc1 = lra over homerange
SAVE 'hrc'
END
Three-toed Woodpecker
; Habitat model for Three-toed Woodpecker (Picoides tridactylus)
; 13 May 2003, 2 Mar 2004
; In RMLANDS workspace, use 10 800 RUNLOOP 'TTWO.TXT'
; 8 March 2004, Modifications to include insect disturbances
gridpath = 'd:\habit@\input\ttwo'
BLOCK 1100 1100 14 ; 1100x1100 blocks, with 14 cell edge buffer
covercond = READ 'covcond.asc' ; read grids
covercondt1 = -10 READ 'covcond.asc' ; covercond for step t-1
agewfire = READ 'wfashm.asc'
agesbeetleh = READ 'sbashm.asc'
agesbeetlel = READ 'sbaslm.asc'
agefbeetleh = READ 'fbashm.asc'
agefbeetlel = READ 'fbaslm.asc'
nodist = (GET 'ttwocover.txt') CATEGORY covercond ; undisturbed based on cover type
withdist = (GET 'ttwocover.txt') CATEGORY covercondt1 ; burned based on cover type in step t-1
dist1 = 1.33 * withdist ; it's always better after a disturbance
sbh = (-9999 0, 0 1, 10 .75, 20 0) ARBFN agesbeetleh
sbl = (-9999 0, 0 .75, 10 .25, 20 0) ARBFN agesbeetlel
fbh = (-9999 0, 0 1, 10 .75, 20 0) ARBFN agefbeetleh
fbl = (-9999 0, 0 .75, 10 .25, 20 0) ARBFN agefbeetlel
fire = (-9999 0, 0 1, 10 1, 20 0) ARBFN agewfire ; <10 y since burn or bugs is optimal;
; >20 y has no additional value
dist = MAX (dist1 * fire) (dist1 * sbh) (dist1 * sbl) (dist1 * fbh) (dist1 * fbl)
hrc1 = 125 COMBINE dist ; burned/bug-killed hrc
hrc2 = 350 COMBINE nodist ; undisturbed hrc (larger home range)
hrc = MAX hrc1 hrc2 ; hrc is best of disturbed/undisturbed
SAVE 'hrc' ; save files with year attached to names
END
Olive-sided Flycatcher
; Habitat model for Olive-Sided Flycatcher (Contopus cooperi)
; "Quick, three beers!"
; In RMLANDS workspace, use 10 800 RUNLOOP 'osfl.txt'
;test = 1 ; uncomment this line to do a one-block time test
gridpath = 'z:\RMLANDS\SJNF\hrv_a1\run1\osfl' ; Location of zipped input grids
BLOCK 1100 1100 6 ; 1100x1100 blocks, with 6 cell edge buffer
covcond = READ 'covcond.asc' ; read grids
cond = READ 'cond.asc'
f1 = (GET 'osflforage.txt') CATEGORY covcond
edgemax = THRESHOLD 100 -20
f2 = (GET 'osfledge.txt') (1 4 1, 1 3 .75, 1 2 .05, 2 4 .75, 2 3 .5, 3 4 .1) 'LOGISTIC' (100 -20) GOODEDGE cond
forage = MIN f1 f2
nest = (GET 'osflnest.txt') CATEGORY covcond
hrcf = 150 ('MEAN' 'S') COMBINE forage ; food based on squared distance-weighted mean
hrcn = 150 ('AREA' .5 .1 .6 .75 .1 1) COMBINE nest ; at least 10% of nest > .75 for prime nesting;
; or 10% better than .5 for value of 0.6
hrc = GEOMEAN hrcf hrcn
SAVE 'hrc'
END
Elk
; Habitat model for Elk (Cervus elaphus)
; In RMLANDS workspace, use 10 800 RUNLOOP 'elk.txt'
gridpath = 'z:\RMLANDS\SJNF\hrv_a1\run1\elk' ; Location of zipped input grids
BLOCK 300 300 69 ; 1100x1100 blocks, with 69 cell edge buffer
covcond = READ 'covcond.asc' ; read grids
slope = 1 READ 'slope.asc'
forage1 = (GET 'elkforage.txt') CATEGORY covcond ; cover types for foraging
steep = 25 -3 .8 LOGISTIC slope ; like shallow slopes
forage = forage1 * steep
cov1 = (GET 'elkcover.txt') CATEGORY covcond ; cover type for cov
cov = cov1 * steep
SETWINDOW 300 -40
forage_r = forage * 300 -40 LOGISTIC JUXTAPOSE cov >= .75 ; realized forage must be near cover
SETWINDOW 500 -70
cov_r = cov * 500 -70 LOGISTIC JUXTAPOSE forage >= .5 ; realized cover must be near forage
hrc1 = 1725 ('AREA' .5 .66 1) COMBINE forage_r ; at least 2/3 of forage must be > 0.5
hrc2 = 1725 ('AREA' .75 .5 1 .85 .25 1) COMBINE cov_r ; at least half of cov must be > 0.75, or
; a quarter better than 0.85.
hrc = GEOMEAN hrc1 hrc2 ; geometric mean
SAVE 'hrc'
END
10> |