draw topo geogrid boxes
Plotting WRF nested domains in a single plot using line-boxes to represent the extension of the nested domains in the parent one represented using its topography
Plotting different geo_em.d[nn].nc topography from WPS files
draw_topo_geogrid_boxes(ncfiles, values)
  ncfiles= ',' list of geo_em.d[nn].nc files to use (fisrt as topographyc reference)
  values= [minTopo],[maxTopo]:[lonlatL]:[title]:[graphic_kind]:[mapvalues]:[labels]:[legvals]:[close]
    [min/max]Topo: minimum and maximum values of topography to draw
    lonlatL: limits of longitudes and latitudes [lonmin, latmin, lonmax, latmax] or None
    title: title of the graph ('!' for spaces)
    graphic_kind: kind of figure (jpg, pdf, png)
    mapvalues: map characteristics [proj],[res]
      see full documentation: http://matplotlib.org/basemap/
      [proj]: projection
        * 'cyl', cilindric
        * 'lcc', lambert conformal
      [res]: resolution:
        * 'c', crude
        * 'l', low
        * 'i', intermediate
        * 'h', high
        * 'f', full
    legvals: [locleg]|[fontsize]:
      [locleg]: location of the legend (0, autmoatic)
        1: 'upper right', 2: 'upper left', 3: 'lower left', 4: 'lower right',
        5: 'right', 6: 'center left', 7: 'center right', 8: 'lower center',
        9: 'upper center', 10: 'center'
      [fontsize]: font size for the legend (auto for 12)
    labels: labels to write in the graph ('!' for spaces)
    close: Whether figure should be finished or not
$ python ${pyHOME}/drawing.py -f geo_em.d01.nc,geo_em.d02.nc -o draw_topo_geogrid_boxes -S '0.,1500.:None:WRF!domain!centered!at!SIRTA:png:cyl,i:d01$_{15k}$,d02$_{3k}$:0|10:True'