Create location from command line using "r.in.gdal" and similar

[mankoff]

I find that I am continually creating some new dummy location just so I can then create the actual location I want using r.in.gdal. I could maintain a global dummy location and use that, but it makes the code less portable when I share it with others. Therefore, I often do this:

rm -fR grass
mkdir grass
grass70 -c EPSG:4326 ./grass/latlon
r.in.gdal -c input=/path/to/somefile.TIF location=foo output=tmp
exit
grass70 ./grass/foo/PERMANENT

I find it difficult to know what EPSG code I should use for somefile.TIF, because it reports many (see below), therefore it would be nice to create the location from outside of GRASS, but using GRASS. For example,

grass70 -c input=/path/to/somefile importer=r.in.gdal

Coordinate System is:
PROJCS["WGS 84 / UTM zone 33N",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",15],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH],
    AUTHORITY["EPSG","32633"]]
Origin = (444892.500000000000000,8690407.500000000000000)
Pixel Size = (15.000000000000000,-15.000000000000000)

[sbl]

See: ​http://grass.osgeo.org/grass70/manuals/grass7.html

There is
grass70 -c geofile