@@ -290,15 +290,13 @@ def _add_transects_to_moc(mesh, mocMask, southernBoundaryEdges,
290290 mocMask ['transectVertexGlobalIDs' ] = \
291291 (('nTransects' , 'maxVerticesInTransect' ), transectVertexGlobalIDs )
292292
293- mocMask ['transectNames' ] = mocMask .regionNames .rename (
294- {'nRegions' : 'nTransects' })
295-
296293 if 'nRegionsInGroup' not in mocMask :
297294 nRegions = mocMask .sizes ['nRegions' ]
298295 nRegionGroups = 2
299296 nRegionsInGroup = nRegions * numpy .ones (nRegionGroups , dtype = int )
300297 regionsInGroup = numpy .zeros ((nRegionGroups , nRegions ), dtype = int )
301298 regionGroupNames = ['MOCBasinRegionsGroup' , 'all' ]
299+ regionNames = mocMask .regionNames .values
302300 nChar = 64
303301 for index in range (nRegionGroups ):
304302 regionsInGroup [index , :] = numpy .arange (1 , nRegions + 1 )
@@ -315,6 +313,17 @@ def _add_transects_to_moc(mesh, mocMask, southernBoundaryEdges,
315313 for index in range (nRegionGroups ):
316314 mocMask ['regionGroupNames' ][index ] = regionGroupNames [index ]
317315
316+ # we need to make sure the region names use the same string length
317+ mocMask ['regionNames' ] = \
318+ (('nRegions' ,), numpy .zeros ((nRegions ,),
319+ dtype = '|S{}' .format (nChar )))
320+
321+ for index in range (nRegions ):
322+ mocMask ['regionNames' ][index ] = regionNames [index ]
323+
324+ mocMask ['transectNames' ] = mocMask .regionNames .rename (
325+ {'nRegions' : 'nTransects' })
326+
318327 mocMask ['nTransectsInGroup' ] = mocMask .nRegionsInGroup .rename (
319328 {'nRegionGroups' : 'nTransectGroups' })
320329
0 commit comments