20250312
12542000
1.0
FALSE
RasterCalculator " "DEM_2022_1ft" - "Detrending_Hydrosurface_2022_1ft_1"" P:\MPWEO_jmcsloy\WettedAreaModels_2022\2022_FocalReach_DEM.gdb\DtDEM_2022_1ft_1
Rename P:\MPWEO_jmcsloy\DtDEM_2022\2022_FocalReach_DEM.gdb\DtDEM_2022_1ft_1 P:\MPWEO_jmcsloy\DtDEM_2022\2022_FocalReach_DEM.gdb\DtDEM_2022_1ft RasterDataset
ProjectRaster V:\Users\jmcsloy\4_Eric\2022_FocalReach_DtDEM.gdb\DtDEM_2022_1ft C:\Users\jmcsloy\Documents\ArcGIS\Projects\MyProject_transfer2M\MyProject_transfer2M.gdb\DtDEM_2022_1ft_WebMerc PROJCS["WGS_1984_Web_Mercator_Auxiliary_Sphere",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mercator_Auxiliary_Sphere"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",0.0],PARAMETER["Standard_Parallel_1",0.0],PARAMETER["Auxiliary_Sphere_Type",0.0],UNIT["Meter",1.0]] "Nearest neighbor" "0.304800609601219 0.304800609601219" WGS_1984_(ITRF00)_To_NAD_1983 # PROJCS["NAD_1983_StatePlane_California_I_FIPS_0401_Feet",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["False_Easting",6561666.666666666],PARAMETER["False_Northing",1640416.666666667],PARAMETER["Central_Meridian",-122.0],PARAMETER["Standard_Parallel_1",40.0],PARAMETER["Standard_Parallel_2",41.66666666666666],PARAMETER["Latitude_Of_Origin",39.33333333333334],UNIT["Foot_US",0.3048006096012192]],VERTCS["NAVD88_height_(ftUS)",VDATUM["North_American_Vertical_Datum_1988"],PARAMETER["Vertical_Shift",0.0],PARAMETER["Direction",1.0],UNIT["Foot_US",0.3048006096012192]] NO_VERTICAL
002
DtDEM_2022_1ft_WebMerc
-13706923.620521
-13668166.089406
4959074.293589
4979513.308067
1
file://\\IBRSACGISHV009\C$\Users\jmcsloy\Documents\ArcGIS\Projects\MyProject_transfer2M\MyProject_transfer2M.gdb
Local Area Network
Projected
GCS_WGS_1984
Linear Unit: Meter (1.000000)
WGS_1984_Web_Mercator_Auxiliary_Sphere
<ProjectedCoordinateSystem xsi:type='typens:ProjectedCoordinateSystem' xmlns:xsi='http://www.w3.org/2001/XMLSchema-instance' xmlns:xs='http://www.w3.org/2001/XMLSchema' xmlns:typens='http://www.esri.com/schemas/ArcGIS/3.3.0'><WKT>PROJCS["WGS_1984_Web_Mercator_Auxiliary_Sphere",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mercator_Auxiliary_Sphere"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",0.0],PARAMETER["Standard_Parallel_1",0.0],PARAMETER["Auxiliary_Sphere_Type",0.0],UNIT["Meter",1.0],AUTHORITY["EPSG",3857]]</WKT><XOrigin>-20037700</XOrigin><YOrigin>-30241100</YOrigin><XYScale>148923141.92838538</XYScale><ZOrigin>-100000</ZOrigin><ZScale>10000</ZScale><MOrigin>-100000</MOrigin><MScale>10000</MScale><XYTolerance>0.001</XYTolerance><ZTolerance>0.001</ZTolerance><MTolerance>0.001</MTolerance><HighPrecision>true</HighPrecision><WKID>102100</WKID><LatestWKID>3857</LatestWKID></ProjectedCoordinateSystem>
150000000
5000
ISO19139
20250501
18102300
20250501
18102300
Raster-based data viewing and analysis.
<DIV STYLE="text-align:Left;"><DIV><DIV><P><SPAN>This is a raster with 1-ft cell resolution representing a detrended version of the ArcGIS Terrain dataset, "2022 Trinity River Focal Reach Digital Terrain" (YT TSP 2023). For further details of the source data can be found in: Pryor, C., T. Grey, and B. Hocker. 2023. 2022 Trinity River focal reach digital terrain model survey report. Report for the Trinity River Restoration Program (TRRP). Yurok Tribal Fisheries Program, Arcata, California. Available: https://www.trrp.net/library/document?id=2627.</SPAN></P><P><SPAN>Detrending used a mass point file derived from wetted cells (water depth > 0) in the 2022 SRH2D modeled flow output for 450cfs, along with the 2022 DTM's "Terrain_Extent" polygon data plus a regional buffer to create a "hydro-planar" Terrain surface that represents the water surface from the source data. The regional buffer was created by buffering the "Terrain_Extent" polygon by 1,000,000 feet. Elevation for the regional buffer was set to a uniform 1596.69 ft* as previously used for the elevation mid-point of the focal reach. Both the original DTM and this hydro-planar surface were converted to 1 ft resolution rasters. The hydro-planer raster was then subtracted from the DTM-derived raster (a DEM) to yield this Detrended Digital Elevation Model (DtDEM)</SPAN></P><P><SPAN>* = This elevation is the midpoint between the water surface elevation at the lower end of the DTM and the water surface elevation immediately below Lewiston Dam. While this uniform elevation theoretically could make water flow out of the river near the dam (river is higher than the value of 1596.69) and into the river near the North Fork, the slope created by this is minuscule. At the upper and lower ends of the river where this potential slope would be greatest, it could account for only about 0.025 ft elevation per 100 feet lateral to the river.</SPAN></P><P><SPAN>Citation: Yurok Tribe Fisheries Department Design and Technical Services Program (YT TSP), 2023. Data Package for 2022 Trinity River Focal Reach Digital Terrain. </SPAN></P></DIV></DIV></DIV>
Conversion by Trinity River Restoration Program from data provided by YT TSP.
DtDEM_2022_1ft
TRRP (Trinity River Restoration Program) and YT TSP (Yurok Tribe Fisheries Department Design and Technical Services Program ) 2023. Detrended Digital Elevation Model (DtDEM) Derived from the 2022 Trinity River Focal Reach Digital Terrain. Data Package. TRRP, Weaverville, California.
2023-10-09T00:00:00
2023-10-10T00:00:00
Eric Peterson
Trinity River Restoration Program
Science Coordinator
Weaverville
CA
96093
US
530-623-1800
<DIV STYLE="text-align:Left;"><DIV><DIV><P><SPAN>Data provided as-is, with no warranty or guarantee. This data is for planning purposes only and should not be used for legal or cadastral purposes. Please contact the Trinity River Restoration Program (TRRP ) for information regarding the use of this data.</SPAN></P></DIV></DIV></DIV>
-123.131376
-122.783212
40.777054
40.637871
1
Microsoft Windows 10 Enterprise N, version 22H2 (Build 19045.5487), ESRI ArcGIS Pro 3.3.1
Eric Peterson
Trinity River Restoration Program
Science Coordinator
Eric Peterson
Jeanne McSloy
Trinity River Restoration Program
Natural Resource Specialist
Jeanne McSloy
Weaverville
CA
96093
US
530-623-1800
Cort Pryor
Yurok Tribal Fisheries Program
Geomatics and Survey Branch Manager
180 Klamath Blvd
Klamath
CA
95548
cpryor@yuroktribe.nsn.us
US
707-482-1350
Detrended DEM
20250501
dataset
Eric Peterson
Trinity River Restoration Program
Science Coordinator
Eric Peterson
Jeanne McSloy
Trinity River Restoration Program
Natural Resource Specialist
<DIV STYLE="text-align:Left;"><DIV><P><SPAN>Data provided as-is, with no warranty or guarantee.</SPAN></P></DIV></DIV>
File Geodatabase Raster Dataset
na
LZW
File Geodatabase Raster Dataset
3.3.1
Eric Peterson
Trinity River Restoration Program
Science Coordinator
Eric Peterson
EPSG
5.3(9.0.0)
Elevation relative to the water surface as represented by water surface elevation (WSE) modeled by SRH2D at 450cfs.
A brief description of the lineage is also provided in the metadata abstract. Intermediate data layers mentioned in the process steps are available from TRRP upon request.
import arcpy
import re
# Define the list of mass points files
mass_points_files = [
"P:\\MPWEO_jmcsloy\\WettedAreaModels_2022\\WettedAreaModelComponents.gdb\\AllFLows_FDS\\F23D_F2P_s450cfs"
]
# Define the paths for hard line and soft clip
hard_line = "P:\\MPWEO_jmcsloy\\WettedAreaModels_2022\\WettedAreaModelComponents.gdb\\AllFLows_FDS\\DtDEM_beyond_OG"
soft_clip = "P:\\MPWEO_jmcsloy\\WettedAreaModels_2022\\WettedAreaModelComponents.gdb\\AllFLows_FDS\\T2022_Terrain_Extent_polygon"
# Output feature dataset to store the terrains
output_fds = "P:\\MPWEO_jmcsloy\\WettedAreaModels_2022\\WettedAreaModelComponents.gdb\\AllFLows_FDS"
# Check out the 3D Analyst extension
arcpy.CheckOutExtension("3D")
# Function to extract the number before "cfs"
def extract_number(file_path):
match = re.search(r'(\d+)cfs', file_path)
return match.group(1) if match else ""
# Iterate through each mass points file to create terrains
for mass_points in mass_points_files:
try:
# Extract the relevant part of the mass points file name
number = extract_number(mass_points)
terrain_name = f"Ter_{number}cfs"
# Create a new Terrain Dataset in the output feature dataset
terrain_path = arcpy.ddd.CreateTerrain(in_feature_dataset=output_fds, out_terrain_name=terrain_name, average_point_spacing=50, pyramid_type="ZTOLERANCE")[0]
# Add terrain pyramid levels
arcpy.ddd.AddTerrainPyramidLevel(in_terrain=terrain_path, pyramid_level_definition=["1 1000"])
# Add the hard line, soft clip, and mass points file to the Terrain
arcpy.ddd.AddFeatureClassToTerrain(in_terrain=terrain_path, in_features=[
[hard_line, "SHAPE", "Hard_Clip", 3, 0, 1, True, False, "DtDEM_beyond_OG_embed", "<None>", False],
[soft_clip, "<None>", "Soft_Clip", 4, 0, 1, True, False, "T2022_Terrain_Extent_polygon_embed", "<None>", False],
[mass_points, "Shape", "Mass_Points", 5, 0, 1, True, False, f"{number}cfs_embed", "<None>", False]
])
# Build the Terrain Dataset
arcpy.ddd.BuildTerrain(in_terrain=terrain_path)
print(f"Created terrain: {terrain_name}")
except Exception as e:
print(f"Error processing {mass_points}: {e}")
# Check in the 3D Analyst extension
arcpy.CheckInExtension("3D")
Creation of "DtDEM_beyond_OG". Created from 2016TERRAIN "DtDEM_beyond": the "Focal_Reach_Terrain_Extents_Polygon" was buffered by 1,000,000 feet. A new 3 dimensional polygon Feature Class was created and the buffer polygon (2D) was copied and pasted into it. The new feature class was then edited and all vertices were given an elevation of 1596.69 ft (the midpoint between the upper and lower elevations within the reach used previously).
Detrending requires extending the water surface laterally from the river. The large buffer allows the convenience of a single elevation for the lateral perimeter without introducing substantive error. See DtDEM_1ft from 2016 TERRAIN/2016_FocalReach_DTM.gdb for more details.
Jeanne McSloy
Trinity River Restoration Program
Natural Resource Specialist
Jeanne McSloy
Weaverville
CA
96093
US
530-623-1800
Jeanne McSloy
Creation of DtDEM_HydroSurface (Terrain format). A Terrain format layer was created using SRH2D model output for 450cfs wetted (water depth>0) cell centroids with water surface elevation (WSE) assigned to Z (mass points), DtDEM_beyond_OG (hard lines), and T2022_Terrain_Extents_Polygon" (soft clip).
Provides the laterally extended water surface elevation, clipped to match the 2022 DTM.
Creation of Detrending_HydroSurface_2022_1ft (raster). The Terrain to Raster tool was used with 1 foot resolution and processing extent set to match the DEM__2022_1ft raster.
Rasterization of the extended water surface, to use for final calculation of the DtDEM.
Creation of DtDEM_2022_1ft (raster). The Spatial Analyst tool Raster Calculator was used to subtract the Detrending_Hydrosurfase_2022_1ft from DEM_2022_1ft. Positive values are relative elevations above the water surface; negative values are below the water surface (bathymetric).
Provides the final raster DtDEM.
Pixel
67057
127157
0.304801
0.304801
32
1
Upper Left
FALSE
LZ77
1
pixel codes
FGDBR
TRUE
row and column
0.304801
0.304801
/9j/4AAQSkZJRgABAQEAAQABAAD/2wBDAAgGBgcGBQgHBwcJCQgKDBQNDAsLDBkSEw8UHRofHh0a
HBwgJC4nICIsIxwcKDcpLDAxNDQ0Hyc5PTgyPC4zNDL/2wBDAQkJCQwLDBgNDRgyIRwhMjIyMjIy
MjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjL/wAARCADIASwDASIA
AhEBAxEB/8QAHwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUFBAQA
AAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkKFhcYGRolJicoKSo0NTY3
ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWm
p6ipqrKztLW2t7i5usLDxMXGx8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/8QAHwEA
AwEBAQEBAQEBAQAAAAAAAAECAwQFBgcICQoL/8QAtREAAgECBAQDBAcFBAQAAQJ3AAECAxEEBSEx
BhJBUQdhcRMiMoEIFEKRobHBCSMzUvAVYnLRChYkNOEl8RcYGRomJygpKjU2Nzg5OkNERUZHSElK
U1RVVldYWVpjZGVmZ2hpanN0dXZ3eHl6goOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3
uLm6wsPExcbHyMnK0tPU1dbX2Nna4uPk5ebn6Onq8vP09fb3+Pn6/9oADAMBAAIRAxEAPwDs6KKK
5jsCiiigBVJVgw6jmt2ykCPuX7p5FYNamnyAx7WPTj6VpTs7xls0eXmtOTpKpHeLOhyHwR6Vi6mh
tpjcKMM5AFXoZzA+1+V9amvrdLq2ZeDuHB968WvRdCaT2ezNcqxsXLmfo0c9qKDdFKFAEibiR3Pe
jTplSR4XBKzLsGOxPSmXDyJCtu4+5wCRzVeNtkqP12sDXVSjzUOVns1laaQjoY5GQ9QcU2rmpAfa
RIg/dyjcD71TropT54KRjUpuD9QooorQzCiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooA
6fwn924+orpa5rwn924+orpa3hsc1T4goooqiAooooAKKKKAPM6KKK5jsCiiigAqSGQxyA54PB+l
R0UCaTVmdDG2YVD8qfuv2NWLacofLY5XsfSs3TLxWha0mYYPKE/yq3PF5QWSMkoeuexoi4106VVW
d/6Z81i8NPCVeek9P07fIl1OzSeMyMMbRkkVzssPlqGDBlPcV1NnOJ4fvA44INULzTxbxSeXgxuc
4PavJjKWGqOnLofTYDExxFJRlv0KQVZ9LCsfmQ8VmjmrVqzbZYx34/Gq0qmNixHHQ10UZck2ujPU
nh+eHL1WwlFA5GR0oruTueZKLi7MKKKKBBRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQBsaJqs
WmiXzVZt/TFa3/CU2n/PKSuRoqlJohwTdzrH8V24HyQOx9M4pn/CWx/8+j/99/8A1q5ainzsXs4n
VL4shLYa1dR67s/0rUs9Ws77Aikw/wDcbg1wNOR2jcOhIYHIIoU2DproelUVQ0e+N/YLI3+sX5W+
tX62TuYNW0PM6KKK5jrCiiigAooooAOldHp1wtzZ7GUscYOO1c5VzTr1rKfcBuVuCM4rKsm4NrdG
Nah7Zcq3NBkazuN6Ag+9aqSJdQAnlG+8KjlWO+thInJqjBK1jIc5aNjg1M0swoc6X7yP9f8ADHhx
c8DX9nJ2T/r/AIcfNpixOWhXAPPFZl/btEPmzhh+tdEtzET1ZfqKg1KH7RZ/KAxXnivKtVpTXtE7
H0+DzKFSSXMmcnEcqVPUU6kmQwy5IwDSghhkV61KaaOvGUHf2kdmFFLikrY88KKKKACiiigAoooo
AKKKKACiiigAooooAKKKKACiiigAooooAKKKKAOp8JE+RcjPAZf610dc54S/1N1/vL/WujreGxzT
+JnmdFFFYHSFFFFABRRRQAUUUUAbOh3YWVoHbG77vpVy7iCSMCPlbkVzQJByCQfUVqW1+00RglJL
fwN/SsaUHRr88fhe/wDmcOZ4X6zRc0/eX4mhbv8AMIn5B6GrKjypcHpVG0kD89xxWlKhO1s5OOlX
jXGFTklopfmeDhYylFyS96Ji6nZN5hwAUblT6ViyQyQ5DDArtGiM1qVGNw6EjODXM3V3MVMMyLvV
sZ215lGUlJwXQ+8y/FyrUldFVTkZoPWkVgc4GM1ehspyqzrGHQj8675VVFXZnUopSaeiKNFPuIpL
c4kTBIzjNQmTjgfnVqrBq6ZEcJVlsh9FR+d2IqQfMMiqU0xVMLVp6yQUUUVRzhRRRQAUUUUAFFFF
ABRRRQAUUUUAFFFFABRRRQB1PhL/AFN1/vL/AFro65jwmx3XCZ+UgGunreGxzVPiPM6KKKwOkKKK
KACiiigAooooAKVWKMGHBByKSigDUtpxHcbv4ZPmFbiljtHBHUGuWjlzGI24wcqfQ1q6dqA/1E5w
B0PpV4qn9YoqaV3HdfqfPzoyw1ZtOyez/Q2NxVugGaqXumx3Y3hcSeoq5glT0YHoRREeNrdRXzs0
/jjuj08NiJUZpLr/AFY5eXSUh3eZcCNuwIODRaanPbZjLCRVOMGuplhjlXDqCPesO8s4YXYGLg9C
K0pYhS0qK57ka8cRFwnqyrdX0F2cyW3zAYBDEVXtbe2uJikjmIY+U9s08afM6l4wHUHoOtVWUoxB
BBHY13wp05RtTdiFWqUfcWy7dCxJoz+biK4gdc9d4GKgki8iUp5iuB3XpSZPrQASeB1rSFJxd5SG
8XVqLlWo/YrJkP8AN/dxUdKysjYYEEdjTpJd8YyvzDuO4rRNx80YWU32YyiiitDFpp2YUUUUAFFF
FABRRRQAUUUUAFFFFABRRRQBoaZqr6Y0jJGrl+OTitH/AISyf/n1j/76Nc9RVKTRLinuFFFFSUFF
FFABRRRQAUUUUAFFFFABUsbr918gf3h1FRUVUZSi7xdmRUpxqR5ZK6NO11Ce2fCsJYv5VqpqaOBu
QowrmUdo23KcGrsV9GcLOpH+0v8AhRUWHqvmqw17r/I8urgcRHTDu67Pf5HSR3ccuApzTpEWT5XU
Ee9YYCrh423KehFaNreK0flu3I715mKyvkj7XDu67CoZlafs8QuVrqNn+0xOFhtwYvVTzmo5LaK7
OZY9rfTnNaauGXKEGqsl2FZvNVdo6EHmvPpxqylywWq+89j63CMeZtW733/zOevrZLaZMbijdRmq
zspcGJWUehbJrUvQuoLm0JZozlk71lqTHKD3U5xXr0lJR5anxLozeFVVEp03o+w8Q3Fw/CO7H2qT
+z7rOGiZR6mtG/v3e0iktnwD98gdD6VmpfTgsHkZ1YYKsTUQnWnG8UkKyunPqVmXDEZzg9RVqOxl
mtzNFhgv3l7ioZIwCTGdy9fpRDPLbyB4nKn271u1JxvB6jqTbdpL0I+hwaKvM9ve5LjyZz3H3Wqm
6NG5VhgiqhU5tGrMylG2qG0UUVoSFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUA
FFFFABRRRQAUUUUAWrK5EMoVz8jHr6H1rZa2SdgG+ST+F1/irnK2NLuXZQHYssTDHrg1HNUi/c/p
/wDBOHMMNTlT9vf3luu6/wCAWZLe7gXhi6+ootoUmOZc59DWoGBwex54prQpJz9054IrnWaS5XGo
uXzX6nmSy9c3NSd/J/oZtxCLSQXMAKlThhjqKpanaiW4+0Rj5HXIx61sur7WjYb1PGarQmOMNbTj
C5+UntVT56tLnT5pR6rrH/gHTgcWsJiOVpxi+j6MxvIubWFZgMI45H+IoKRXKjyV2zd0HRvpWhLH
cWM6vs86I9x2FVJ7dnT7TCm0YJYA/drmhX5nd6Pv/mfS+zTV07plSKRraU7k9mRhVuaySa3+1WgJ
T+JOpU09ZoL+MJckJN2kxgfjUBNzpk5CsR7jo1bNylK8dJfgzDVe6ymRjrVtLtGjEVxHvUDCsPvD
8auNPY36gznyJv7wGQaiktbEj5btQR/snn9KJVk9JxaZMY31RnMAGO05FJV+1iinmFqi7gxP7wjm
qc0TQzNGw5U4reFVSfL1CpDl1GUUUVqZhRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAF
FFFABRRRQAUUUUAFX9LYC5IPQrVCp7Q7byE54B5qKjtG41S9rGUO6Zq30klvIlzESVX5WXsRVuz1
KK5UYO2QdVP9Ki/dTIyBg8Z+XI7VhyxS20205BB4I71hUowls99Uefg/31P2VRWnDTzt0OpFzg8r
VC6Inu/LXgsOKsWUEz2qtO2WYZHGCPrVS43WupxuwOzHJrXAex537LSVnoeXjqVeMuWrLminuW4L
h4v3VwmB0z2NMvdPJjaSyfDH+HPBq4xinTcCCPaotpgAZHBXrtNcM1GtK8Vyz6p7M7cNiKuD1+KH
c5tleaSQkBZEyWzx0q9bX0U0ZgvgW/utiptQt45l3W8yo7HLA96ozQvIoEjw7x/EG5Nawp1Jqzi1
8tj3KmLwtaCamvvsyw+lwSDME/PYGs+4g+zy+WXDEdSvakZHiIzkehBphJJyetdMYyvqzKEeXVvm
X9dVuTw3LWx3QDa2PvHmk3G7uSZGw7nr2zTYreecERRljSzM28b0COvBI71i+WM7R+LudsY+0XNb
QjkjaJyrjBFMq9G6XiLC42yD7retVJYmhkKOMEVrSq83uy3RzVaXJqthlFFFbGQUUUUAFFFFABRR
RQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFLyehxSUUmrlQk4u6LNnP5Em0k7W4OP51tqc
H97GGCnglcj8K5zqc1YtnlMqxb2IbgDPrUrT3ZK6Zz5hhI1/39OXLNde/kdOLnAwUcHtxUF1GbyE
xMhGejHqKtQS+ZArY6DBHeq2q3LWtrvjA3Mcc9q8yNWNKv7sLST7nDTw9XEQjBTun5GfHpE0YAW+
YEds4q19gl/imLe9YAlkU58xsk881estUlt3CyuXi755xXoVMVi/iun8tToq5LSatGV2af8AZwPV
yaibRo2bduJIq3Ff2s7YjmXd6HirKsQ2TWCzHErVM86eAjSlyzTRgLZNBclJfnhb+H0PtU8dhZRn
LI7H0Y1d1AZTf0IPWm288d3ERkFl6111JqvTjXmvJ2MoYnE0HKhSnZdCWF4oYxHGAienSo7y3juV
K/IM9z608ZjbpkVVmuJIXJdNyE8GvNq0YRkpK9vU9DAV8RWb5Wudd76/iYksUtlchXUbl5B9amuX
+0RK5H7xeDx1Fapntb/KMdpAGM0kdhCoIDqV75NZ+3akpNao+idWLhaqrM56ip723ks7gxuBhuVI
qurBuO9epCrGaTRzTw8ormWqFooorQwCiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKK
KKACiiigAqxb3SWp3eQHfqGLYxVemyHAFZVY3R1YVKcvZtaPc6Gwu1mJaPcrdwTRqcczRiWNvMjH
34z/ADFZ+jOA7g45HHrW3JJFGAVfOe1YOKrNfzen5njYlvK8W/Zr3NHa/wCRzbRP5nyxPtPK8dqj
b5chgQfQjFdTFcROMAgH6VQ1C6WGUJNbpIDyOxqXUrU5ezlE9PDY2ninzQXyuYq9cAflWraX8trE
PPfcp6KetVGu4pAyJAIyBlSvrT7Oxa5zNKSIx+Zon+80mrfmdNfljTc5rTrc3La9t76MqCfdW61m
XEbabqCvH/q3qwbeAjCKEdO8bcj606dTdWTRzMvmr9x+ma2wXNRqaq8JaPr8z5rGfV5xvB8slrZ6
O3l5GhGVkQOOQRTHg3AgDIPY1l6dqAj/AHMvGDj6Vto6sAVYZPpWGLpVMJUaa0e3YxoctaPZr7/k
ZEmnxF9yExuO3aojbzxyqACefwrbeJX+8MH1FReU8bcHcKwbp1FeLs+z2+89GlmGKovlq+/Hv1Mv
UrBrxg6ygFOCGqu1nEQiKvCjG4963pIQ68AeuKrPbqBliE+prB1m3aNz08JiXGkozktDnJoGib1X
1qGukez3qflDL6ism6sGjJKDj0ruoYxP3Z7lSpxnrTKNFFFd97mDTWjCiiigQUUUUAFFFFABRRRQ
AUUUUAFFFFABRRRQAUUUUAFFFFABQwDDBFFFJq5UJuDuiS1I80BRhga21v4o49s6q34VgAlG3qOa
Q5lO5+g7Vwyoy59HY9KrGhiYc1TVHRQyadIxYSBCexOKlmtY5gM7ZV7c8j8a5jYg6ClB2njj6Vvy
Vn8UrrzPFqYDCwlzUG4vyNtNLtlkzIsmM8A1Yup0tYA5G1RxGnqawlu7lRhZ5APZjUbs0rh5ZGdh
3Y5pezle9te92/uKjhnVa+sVbxWtrWv62Jobp4bgy9dx+Yeta8srvaCW2cfQjNYFWLW7ktm+XlT1
WrnS920Sq0VUqKq0m132a7f1sXYlS6XbcKIpv4ZB0P1p+bqzbB3YHQjoasIIruPzYiBnqOnNIkl7
G2yONio/vDg1lDHVqN4zXNHs/wBH+hzVsvw2K96k+SS3XYdFrTLxIoIqwNVhkXKgg+9ReV9o5mt4
1PfHWq9xBZQr87bT7GrTy+TUpQcX5f8AAOB4fGNunTlzW+dvmSyXdzMdsWOfSqw0t533XkrE9gtQ
xy2yyYE3Hbcp/nWgkr5ASTd/uNmu11IwXLhZKPy1OV0K9F82Jpya7kUUNxpkm6OYyQnqrcitWOSC
6ToM9x6VGJImiKyEtnqT2rPmD2j+ZEd4/wBmuOdH67dVFaoutrJlQxcsPJSg7xfS+qH3miCQloTh
j2rEms7m3YiSI4HcCt6DW4pPkf5H9GqZ76JxyA30rnhhsdRfK43PZjnNHl5auq/E5WitKa3t2kZs
SLk5wuKpvCu75XOPcV6caNWS+FkrHYZ7T/MhooorM6wooooAKKKKACiiigAooooAKKKKACiiigAo
oooAKKKKACiiigAooooAKKKKACiiigBfPkhwY3ZPcGtKC9n8rc844rLZd64pQD0HWuapScnvY7k6
Lormim/QuzajNJ8qsQO59aqFyzZYlqlt7YzN8zqi9y1WXlsIAoihMxxyXOKScIPlirv+upg1LorL
t0KgRpmxHGevAFO+zXaFRsMZJ4ZjjNWF1WWNv3UUSL6bagmvZ51KO5Kk5x6VLhUl0sjaniXF2dmi
6L+d4QrtHGw4LcZoTVQqMrEsRxn1rNZAEVtwOeo9KZSjhYvW4pzpv3XBGk2ow3DbZbZQv97uKsRJ
xiF0x2Df41jBWPRSfoKVXdDlWIrpg6lPSlLTs9UedXwmGrr95Gz7o6eG0XaDIyMe+elWRBb44jT8
q5mK/kVSr/MPY4NTpqrBcEn8a83ELFSm25M1oZZR5LQs7eRU+w3X/PCT8qPsN1/zwk/KiivZ5SOZ
iGyuQMmB8fSoKKKlqw07hRRRSKCiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACi
iigApQSrBgeRRRSaTVmVCTjJNCyMXYsByat272Mcf7xXZ2GCOw+lFFZOkjplUlOHLt6DnOnMvyiR
fp/9enCTTAM+U+fTJooo9ku7OZp9yqJIUuC6xFkxwrHv+FSnUDs2i3gUey/40UVXsYdRyqSk7sWL
U5o8/JGR6bcfyqtNJ50pfYqZ7L0oopqnFO6RC0d+5JbWct1nygOOuTSTW8kEmx1w1FFcsq0vacpv
QadVwtpa5//Z