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

    2011 Land Use Land Cover (nlcd) grids, 30 m

    http://www.mrlc.gov/nlcd11_leg.php