# Get PSCF layer
mf = miapp.getMapDocument().getActiveMapFrame()
layerName = 'winterPSCF.shp'
layer = mf.getMapView().getLayer(layerName)
# Get grid x/y data
sn = layer.getShapeNum()
xn = 0
oldx = 0.0
for shp in layer.getShapes():
x = shp.getPoints().get(0).X
if xn == 0:
oldx = x
else:
if x < oldx:
break
else:
oldx = x
xn += 1
yn = sn / xn
print 'Total number: ' + str(sn)
print 'X number: ' + str(xn)
print 'Y number: ' + str(yn)
idx = sn - xn
shp = layer.getShapes().get(idx)
point = shp.getExtent().getCenterPoint()
minx = point.X
miny = point.Y
rpoint = layer.getShapes().get(idx + 1).getExtent().getCenterPoint()
xdelta = rpoint.X - minx
tpoint = layer.getShapes().get(idx - xn).getExtent().getCenterPoint()
ydelta = tpoint.Y - miny
# Create GridData and set value
print 'Create GridData and set value...'
gData = GridData(minx, xdelta, xn, miny, ydelta, yn)
fieldName = 'WPSCF'
for i in range(0, yn):
for j in range(0, xn):
idx = (yn - i - 1) * xn + j
v = layer.getCellValue(fieldName, idx)
gData.data[i][j] = v
# Output GridData file
print 'Output GridData file...'
outfn = 'D:/0303HYSPLITDoc/Meteoinfo_out/cs20132019/cluster/winter/pscf.dat'
gData.saveAsSurferASCIIFile(outfn)
print 'Finished!'