PRO altitude, targetName=targetName if not keyword_set(targetName) then targetName='lutetia' verfile='/astrodata/SHAPE_MODELS/VER/' + targetName + '.ver' read_ver_shape_model, verFile, $ objVertex, objPolygons, $ nbVertex=nbVertex, $ nbFacet=nbFacet COMMON vertexCoord, objVertex result = altitude_func(50) END