----- Original Message -----

From:

"Liza Hoos" <>

To:

"Jason Roberts" <>

Cc:

"" <>

Sent:

Tue, 27 Aug 2013 16:34:39 -0400

Subject:

Re: [mget-help] Interpolate HYCOM GLBa0.08 Equatorial 4D Variables at Points

Hi Jason,

Thanks so much for taking the time to respond to this email while on your vacation! This issue is non-urgent, so please feel free to ignore this until you return.

I updated MGET to version 0.8a50 and am now receiving the below error when I specify the depth value as 20000. The same error occurs for other variables (e.g. salinity, u) at a depth value of 20000. If I specify the depth value as zero, the process completes without any errors but no points are interpolated.

Thanks!

Liza

Executing: HYCOMGLBa008Equatorial4DInterpolateAtArcGISPoints temperature bll_wgs84_full bottomdepthFloat RecordDate # 20000 Nearest false # # 60 120 # # 128 16 16 3 3
Start Time: Tue Aug 27 16:21:39 2013
Running script HYCOMGLBa008Equatorial4DInterpolateAtArcGISPoints...
Interpolating values for 1436 points of ArcGIS FeatureLayer "bll_wgs84_full" of FeatureClass "Z:\sharks\Models\Summer2013\summer2013.gdb\bll_wgs84_full".
Update stopped before all points were processed: 0:00:00 elapsed, 0 points updated, 0 deleted, 0 unchanged, 0:00:00 per point, 1436 points not processed.
IndexError: Indices must be integers or integer slices.

Traceback (most recent call last):
  File "C:\Python27\ArcGIS10.1\Lib\site-packages\GeoEco\ArcGISToolbox\Scripts\HYCOMGLBa008Equatorial4DInterpolateAtArcGISPoints.py", line 5, in <module>
    ExecuteMethodFromCommandLineAsArcGISTool('GeoEco.DataProducts.HYCOM', 'HYCOMGLBa008Equatorial4D', 'InterpolateAtArcGISPoints')
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\ArcGISScripts.py", line 210, in ExecuteMethodFromCommandLineAsArcGISTool
    exec sourceCode in globals(), locals()
  File "<string>", line 1, in <module>
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\DataProducts\HYCOM.py", line 1843, in InterpolateAtArcGISPoints
    Interpolator.InterpolateGridsValuesForTableOfPoints(grids, ArcGISTable(points), valueFields, zField=zField, tField=tField, zValue=zValue, where=where, orderBy=orderBy, method=method, noDataValue=noDataValue, gridsWrap=True, useAbsZ=True, seafloorZValue=20000., numBlocksToCacheInMemory=numBlocksToCacheInMemory, xBlockSize=xBlockSize, yBlockSize=yBlockSize, zBlockSize=zBlockSize, tBlockSize=tBlockSize)
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\SpatialAnalysis\Interpolation.py", line 545, in InterpolateGridsValuesForTableOfPoints
    dataForAllDepths = grids[i].Data.__getitem__(tuple(indices))
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\Datasets\__init__.py", line 3635, in __getitem__
    return getattr(self._Grid(), self._GetMethod)(key)
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\Datasets\__init__.py", line 3249, in _GetUnscaledDataAsArray
    flippedKey = self._ValidateAndFlipKey(key)
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\Datasets\__init__.py", line 3346, in _ValidateAndFlipKey
    raise IndexError(_(u'Indices must be integers or integer slices.'))
IndexError: Indices must be integers or integer slices.

Failed to execute (HYCOMGLBa008Equatorial4DInterpolateAtArcGISPoints).
Failed at Tue Aug 27 16:21:44 2013 (Elapsed Time: 5.00 seconds)

On Sun, Aug 25, 2013 at 1:29 PM, Jason Roberts <> wrote:

Hi Liza,

Sorry for the slow response. I am on vacation and it is hard to keep up on email. (I certainly don't keep up on ALL email while on vacation! but do try to respond to mget-help requests when I can).

I think this problem was fixed in a recent build of MGET. Please try the latest MGET 0.8a50 and let me know if it works.

Thanks,
Jason

----- Original Message -----

From:

"Liza Hoos" <>

To:

"" <>

Cc:

Sent:

Tue, 20 Aug 2013 12:30:06 -0400

Subject:

[mget-help] Interpolate HYCOM GLBa0.08 Equatorial 4D Variables at Points

Hello,

I am having trouble using the tool "Interpolate HYCOM GLBa0.08 Equatorial 4D Variables at Points". If I run the tool specifying 0 as the depth value, the tool works, but if I specify 20000 (bottom depth) as the depth value, the tool runs without any errors but no points are actually interpolated. Do you know what could be causing this problem?

Thank you!

Liza

Executing: HYCOMGLBa008Equatorial4DInterpolateAtArcGISPoints temperature bll_wgs84 BottomDepthHC RecordDate # 20000 Nearest false # # 60 120 # # 128 16 16 3 3
Start Time: Tue Aug 20 12:29:42 2013
Running script HYCOMGLBa008Equatorial4DInterpolateAtArcGISPoints..
Interpolating values for 1436 points of ArcGIS FeatureLayer "bll_wgs84" of FeatureClass "Z:\sharks\Models\Summer2013\summer2013.gdb\bll_wgs84".
Update complete: 0:00:02 elapsed, 1436 points updated, 0 deleted, 0 unchanged, 0:00:00.001605 per point.
Completed script HYCOMGLBa008Equatorial4DInterpolateAtArcGISPoints...
Succeeded at Tue Aug 20 12:29:48 2013 (Elapsed Time: 6.00 seconds)