A vgf file is created using the interactive product generation tools within NMAP2. A closed line is drawn enclosing the region of interest, and a text label is grouped with the line assigning a value of 1 to the contour.
Once the region of interest is defined, a grid file can be created using the GRPHGD standalone program (It can also be created using the graph-to-grid option in NMAP2).
Create grid from VGF file with GRPHGD
GEMPAK-GRPHGD> GDOUTF = 2007071900_cmask.grd GUESS = PROJ = MER GRDAREA = KXKY = 10;10 MAXGRD = 200 CPYFIL = #A218 ANLYSS = CNTRFL = cmask.info GDATTIM = 070719/0000f000 GFUNC = cmask GLEVEL = 0 GVCORD = none KEYCOL = KEYLINE = OLKDAY = GGLIMS = HISTGRD = NO BOUNDS = TYPE = C GAMMA = 0.3 SEARCH = 20 NPASS = 2 QCNTL = GUESFUN = CATMAP = DISCRETE = DLINES = yes;no|-0.5 GGVGF = conusoutline_cmask.vgf EDGEOPTS = GEMPAK-GRPHGD>r
Check with GDINFO
GEMPAK-GDINFO> GDFILE = 2007071900_cmask.grd LSTALL = YES OUTPUT = T GDATTIM = all GLEVEL = 0 GVCORD = none GFUNC = cmask GEMPAK-GDINFO>r GRID FILE: cmask.gem GRID NAVIGATION: PROJECTION: LCC ANGLES: 25.0 -95.0 25.0 GRID SIZE: 614 428 LL CORNER: 12.19 -133.46 UR CORNER: 57.33 -49.42 Number of grids in file: 1 NUM TIME1 TIME2 LEVL1 LEVL2 VCORD PARM 1 070719/0000F000 0 NONE CMASK GEMPAK-GDINFO>
The parameter CMASK is created where values greater than 1 are enclosed by the contour. DLINES defines the epsilon of (here -0.5 is used since the contour was drawn counterclockwise) to add to vlaues on either side of the single contour to be define greater / less than the contour value of 1.0.
The grid point values of SGT(cmask,1) (greater than 1.0) as shown here:
Plot with GDPLOT2
The resultant grid can ge used with the MASK() function and logical operators to define mask or clipping regions of interest. As an example, the 24 hour precipitation in the top panel is masked by the region of interest in the lower panel so that only data within the region of interest will be considered.
GDFILE = nam12 + 2007071900_cmask.grd GDATTIM = f030 GLEVEL = 0 GVCORD = none PANEL = t ! b SKIP = 0 SCALE = 0 GDPFUN = p24i ! mask(p24i,sgt(cmask^070719/0000f000+2,1)) TYPE = f CONTUR = 3/3 CINT = 0 LINE = 2/1/2 FINT = .01;.1;.25;.5;.75;1;1.25;1.5;1.75;2;2.5;3;4;5;6;7;8;9 FLINE = 0;21-30;14-20;5 HILO = HLSYM = CLRBAR = 1 WIND = BM1 REFVEC = TITLE = 1 TEXT = 0.7/2/SW CLEAR = YES GAREA = us IJSKIP = PROJ = STR/90;-100;0
Using GDCSV to find local maximas
The HIGH() function can be used to obtain local maximas from the grid region of interest. The GWFS() gaussian weighted smoothing function can be used to reduce higher frequency features in the grid and focus on broader areas of interest. By using the GDCSV program, the locations of HIGH() output are output to a file for use in determining local mesoscale model domain centers. By masking the region, the model domains are ensured to be within the desired region.
GEMPAK-GDCSV>l GDATTIM = f030 GDFILE = nam12 + 2007071900_cmask.grd GLEVEL = 0 GAREA = grid PROJ = def GVCORD = none GFUNC = high(mask(gwfs(p24i,40),sgt(cmask^070719/0000f000+2,1)),30) SCALE = 0 OUTPUT = f/p24i_highs.dat GEMPAK-GDCSV>r
OUTPUT is to the text file
p24i_highs.dat. You can use the
sort command to return the top two local maxima:
sort -t, -k 5bnr p24i_highs.dat | head -2 508, 244, 42.3889, -72.2385, 2.53370 338, 236, 42.9560, -96.3721, 0.58850
These two values are now scriptable (use
|head -1 for the first and
| head -2 | tail -1 for the second) both for plotting regional areas of interest and for WRF domains. A WRF lesson is outside the reach of this document, but as an example of what is run at Unidata, here is are the primary and secondary regions
Plotting again with GDPLOT2 and GPANOT
The locations of the 2 greatest local maximas (42.3889,-72.2385 and 42.9560,-96.3721) are shown along with the precipitation forecast field. The primary and secondary domains can be drawn using GPANOT to overlay the boxes.
First plot masked precipitation, same as above, but change
PANEL = 0 and
CLEAR = n
GDFILE = nam12 + 2007071900_cmask.grd GDATTIM = f030 GLEVEL = 0 GVCORD = none PANEL = 0 SKIP = 0 SCALE = 0 GDPFUN = mask(p24i,sgt(cmask^070719/0000f000+2,1)) TYPE = f CONTUR = 3/3 CINT = 0 LINE = 2/1/2 FINT = .01;.1;.25;.5;.75;1;1.25;1.5;1.75;2;2.5;3;4;5;6;7;8;9 FLINE = 0;21-30;14-20;5 HILO = HLSYM = CLRBAR = 1 WIND = BM1 REFVEC = TITLE = 1 TEXT = 0.7/2/SW CLEAR = YES GAREA = us IJSKIP = PROJ = STR/90;-100;0 CLEAR = n r