def get_glacier_poly():
<EOL>rgi_fn = os.path.join(datadir, '<STR_LIT>')<EOL>if not os.path.exists(rgi_fn):<EOL><INDENT>cmd = ['<STR_LIT>',]<EOL>sys.exit("<STR_LIT>" % cmd[<NUM_LIT:0>])<EOL><DEDENT>return rgi_fn<EOL>
Calls external shell script `get_rgi.sh` to fetch:

    Randolph Glacier Inventory (RGI) glacier outline shapefiles 

    Full RGI database: rgi50.zip is 410 MB

    The shell script will unzip and merge regional shp into single global shp

    http://www.glims.org/RGI/