|
This file can be copied as a text document to be used
in part to repeat this study.
A description of the methods
used to create this is on this site.
/* Coding the Critical
Lands into Parcels
&echo &off
&s work [response 'What is the name of this analsysis']
cw %work%
w %work%
/************************************************************************
&echo &off
&s agland [response
'Enter the Value for Agricultural Lands']
&s primeimp [response 'Enter the Value for Prime and Important
Agricultural Soils']
&s greenway [response 'Enter the Value for Greenways']
&s histdist [response 'Enter the Value for Historic Districts']
&s scenic [response 'Enter the Value for Scenic Landscapes']
&s gwreserv [response 'Enter the Value for Ground Water Reservoir']
&s cwellprt [response 'Enter the Value for Community Wellhead
Protection Areas']
&s gwrech [response 'Enter the Value for Groundwater Recharge
Areas']
&s ncwelprt [response 'Enter the Value for Non-Community Wellhead
Protection Areas']
&s briver [response 'Enter the Value for 300ft River Buffers']
&s bwetland [response 'Enter the Value for 300ft Wetland Buffers']
&s forbrush [response 'Enter the Value for Forest and Brushland
Areas']
&s openbuf [response 'Enter the Value for 300ft Openspace Buffers']
&s rarespec [response 'Enter the Value for Rare Species Habitat']
&s wetland [response 'Enter the Value for Wetlands']
&s build0 [response
'Enter the Value for Parcels That Will Have 1-2 Houses At Buildout']
&s build2 [response 'Enter the Value for Parcels That Will Have
3-10 Houses At Buildout']
&s build3 [response 'Enter the Value for Parcels That Will Have
10-20 Houses At Buildout']
&s build4 [response 'Enter the Value for Parcels That Will Have
Over 20 Houses At Buildout']
&s size1 [response 'Enter the Value for Parcels That are between
5 acres']
&s size2 [response 'Enter the Value for Parcels That are between
15 acres']
&s size3 [response 'Enter the Value for Parcels That are over
25 acres']
&echo &on
/***********************************************************************
copy /QUAHOG/USERS/UNDERGRAD/MATTA/POINTS/parcel
parcel
copy /QUAHOG/USERS/UNDERGRAD/MATTA/POINTS/critical
additem parcel.pat parcel.pat parcel_area 8 18 F 5
ap
res critical poly agland_point
> 0
calc critical poly agland_point = %agland%
clearselect
res critical poly primeimp_point
> 0
calc critical poly primeimp_point = %primeimp%
clearselect
res critical poly greenway_point
> 0
calc critical poly greenway_point = %greenway%
clearselect
res critical poly histdist_point
> 0
calc critical poly histdist_point = %histdist%
clearselect
res critical poly scenic_point
> 0
calc critical poly scenic_point = %scenic%
clearselect
res critical poly gwreserv_point
> 0
calc critical poly gwreserv_point = %gwreserv%
clearselect
res critical poly cwellprt_point
> 0
calc critical poly cwellprt_point = %cwellprt%
clearselect
res critical poly ncwelprt_point
> 0
calc critical poly ncwelprt_point = %ncwelprt%
clearselect
res critical poly briver_point
> 0
calc critical poly briver_point = %briver%
clearselect
res critical poly bwetland_point
> 0
calc critical poly bwetland_point = %bwetland%
clearselect
res critical poly forbrush_point
> 0
calc critical poly forbrush_point = %forbrush%
clearselect
res critical poly openbuf_point
> 0
calc critical poly openbuf_point = %openbuf%
clearselect
res critical poly gwrech_point
> 0
calc critical poly gwrech_point = %gwrech%
clearselect
res critical poly wetland_point
> 0
calc critical poly wetland_point = %wetland%
clearselect
res critical poly rarespec_point
> 0
calc critical poly rarespec_point = %rarespec%
clearselect
/************************************************************************
calc critical poly total_point
= agland_point + primeimp_point + greenway_point + histdist_point
+ scenic_point + gwreserv_point + cwellprt_point + gwrech_point
+ ncwelprt_point + briver_point + bwetland_point + forbrush_point
+ openbuf_point + rarespec_point + wetland_point
calc parcel poly parcel_area = area
quit
union parcel critical
union
additem parcel.pat parcel.pat point 8 18 F 3
additem union.pat union.pat avg_point 8 18 F 5
ap
calc union poly avg_point = ( area / parcel_area ) * total_point
/************************************************************************
&echo &off
&do count = 2 &to
7120 &by 1
&type *********UP TO %count% of 7120 *****************
clearselect
reselect union poly parcel# = %count%
statistics union poly
sum avg_point
end
&sv sumpoint [show
statistic 1 1]
&if [null %sumpoint%] &then
&sv sumpoint 0
clearselect
reselect parcel poly parcel# = %count%
calc parcel poly point = %sumpoint%
&end
quit
/************************************************************************
&echo &off
additem parcel.pat parcel.pat buildout_point 8 18 I
ap
res parcel poly unit_max
> 0
calc parcel poly buildout_point = %build0%
clearselect
res parcel poly unit_max
> 2
calc parcel poly buildout_point = %build2%
clearselect
res parcel poly unit_max
> 10
calc parcel poly buildout_point = %build3%
clearselect
res parcel poly unit_max
> 20
calc parcel poly buildout_point = %build4%
clearselect
quit
/************************************************************************
additem parcel.pat parcel.pat
area_acres 8 18 F 5
additem parcel.pat parcel.pat area_point 8 18 I
ap
calc parcel poly area_acres = area / 43514
res parcel poly area_acres
> 5
calc parcel poly area_point = %size1%
clearselect
res parcel poly area_acres
> 15
calc parcel poly area_point = %size2%
clearselect
res parcel poly area_acres > 25
calc parcel poly area_point = %size3%
clearselect
quit
/************************************************************************
additem parcel.pat parcel.pat
final_point 8 18 F 2
ap
calc parcel poly final_point = area_point + buildout_point + point
quit
/***********************************************************************
rename parcel %work%
kill critical all
kill union all
&return
|
