WRF模拟出的数据，用grads怎么画出某地区的边界线？

造字仓颉

我用WRF模拟出来的数据，wlmq003.ctl和wlmq003.dat，
其中，ctl文件内容为：
dset ^wlmq003.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM WRF V3.3.1 MODEL
pdef  135 117 lcc  43.355   86.498   68.000   59.000  60.00000  30.00000   87.65000   3333.333   3333.333
xdef  412 linear   83.36471   0.01501502
ydef  264 linear   41.34251   0.01501502
zdef   27 linear 1 1  
tdef    5 linear 18Z31OCT2013     360MN      
VARS  117
LU_INDEX       1  0  LAND USE CATEGORY (-)
U             27  0  x-wind component (m s-1)
V             27  0  y-wind component (m s-1)
W             27  0  z-wind component (m s-1)
PH            27  0  perturbation geopotential (m2 s-2)
PHB           27  0  base-state geopotential (m2 s-2)
T             27  0  perturbation potential temperature (theta-t0) (K)
MU             1  0  perturbation dry air mass in column (Pa)
MUB            1  0  base state dry air mass in column (Pa)
NEST_POS       1  0  - (-)
P             27  0  perturbation pressure (Pa)
PB            27  0  BASE STATE PRESSURE (Pa)
P_HYD         27  0  hydrostatic pressure (Pa)
Q2             1  0  QV at 2 M (kg kg-1)
T2             1  0  TEMP at 2 M (K)
TH2            1  0  POT TEMP at 2 M (K)
PSFC           1  0  SFC PRESSURE (Pa)
U10            1  0  U at 10 M (m s-1)
V10            1  0  V at 10 M (m s-1)
QVAPOR        27  0  Water vapor mixing ratio (kg kg-1)
QCLOUD        27  0  Cloud water mixing ratio (kg kg-1)
QRAIN         27  0  Rain water mixing ratio (kg kg-1)
LANDMASK       1  0  LAND MASK (1 FOR LAND, 0 FOR WATER) (-)
TSLB           4  0  SOIL TEMPERATURE (K)
SMOIS          4  0  SOIL MOISTURE (m3 m-3)
SH2O           4  0  SOIL LIQUID WATER (m3 m-3)
SMCREL         4  0  RELATIVE SOIL MOISTURE (-)
SEAICE         1  0  SEA ICE FLAG (-)
XICEM          1  0  SEA ICE FLAG (PREVIOUS STEP) (-)
SFROFF         1  0  SURFACE RUNOFF (mm)
UDROFF         1  0  UNDERGROUND RUNOFF (mm)
IVGTYP         1  0  DOMINANT VEGETATION CATEGORY (-)
ISLTYP         1  0  DOMINANT SOIL CATEGORY (-)
VEGFRA         1  0  VEGETATION FRACTION (-)
GRDFLX         1  0  GROUND HEAT FLUX (W m-2)
ACGRDFLX       1  0  ACCUMULATED GROUND HEAT FLUX (J m-2)
SNOW           1  0  SNOW WATER EQUIVALENT (kg m-2)
SNOWH          1  0  PHYSICAL SNOW DEPTH (m)
CANWAT         1  0  CANOPY WATER (kg m-2)
SST            1  0  SEA SURFACE TEMPERATURE (K)
SSTSK          1  0  SKIN SEA SURFACE TEMPERATURE (K)
LAI            1  0  Leaf area index (area/area)
H_DIABATIC    27  0  MICROPHYSICS LATENT HEATING (K s-1)
MAPFAC_M       1  0  Map scale factor on mass grid (-)
MAPFAC_MX      1  0  Map scale factor on mass grid, x direction (-)
MAPFAC_MY      1  0  Map scale factor on mass grid, y direction (-)
MF_VX_INV      1  0  Inverse map scale factor on v-grid, x direction (-)
F              1  0  Coriolis sine latitude term (s-1)
E              1  0  Coriolis cosine latitude term (s-1)
SINALPHA       1  0  Local sine of map rotation (-)
COSALPHA       1  0  Local cosine of map rotation (-)
HGT            1  0  Terrain Height (m)
TSK            1  0  SURFACE SKIN TEMPERATURE (K)
RAINC          1  0  ACCUMULATED TOTAL CUMULUS PRECIPITATION (mm)
RAINSH         1  0  ACCUMULATED SHALLOW CUMULUS PRECIPITATION (mm)
RAINNC         1  0  ACCUMULATED TOTAL GRID SCALE PRECIPITATION (mm)
SNOWNC         1  0  ACCUMULATED TOTAL GRID SCALE SNOW AND ICE (mm)
GRAUPELNC      1  0  ACCUMULATED TOTAL GRID SCALE GRAUPEL (mm)
HAILNC         1  0  ACCUMULATED TOTAL GRID SCALE HAIL (mm)
CLDFRA        27  0  CLOUD FRACTION (-)
SWDOWN         1  0  DOWNWARD SHORT WAVE FLUX AT GROUND SURFACE (W m-2)
GLW            1  0  DOWNWARD LONG WAVE FLUX AT GROUND SURFACE (W m-2)
SWNORM         1  0  NORMAL SHORT WAVE FLUX AT GROUND SURFACE (SLOPE-DEPENDENT) (W m-2)
OLR            1  0  TOA OUTGOING LONG WAVE (W m-2)
XLAT           1  0  LATITUDE, SOUTH IS NEGATIVE (degree_north)
XLONG          1  0  LONGITUDE, WEST IS NEGATIVE (degree_east)
ALBEDO         1  0  ALBEDO (-)
ALBBCK         1  0  BACKGROUND ALBEDO (-)
EMISS          1  0  SURFACE EMISSIVITY (-)
NOAHRES        1  0  RESIDUAL OF THE NOAH SURFACE ENERGY BUDGET (W m{-2})
TMN            1  0  SOIL TEMPERATURE AT LOWER BOUNDARY (K)
XLAND          1  0  LAND MASK (1 FOR LAND, 2 FOR WATER) (-)
UST            1  0  U* IN SIMILARITY THEORY (m s-1)
PBLH           1  0  PBL HEIGHT (m)
HFX            1  0  UPWARD HEAT FLUX AT THE SURFACE (W m-2)
QFX            1  0  UPWARD MOISTURE FLUX AT THE SURFACE (kg m-2 s-1)
LH             1  0  LATENT HEAT FLUX AT THE SURFACE (W m-2)
ACHFX          1  0  ACCUMULATED UPWARD HEAT FLUX AT THE SURFACE (J m-2)
ACLHF          1  0  ACCUMULATED UPWARD LATENT HEAT FLUX AT THE SURFACE (J m-2)
SNOWC          1  0  FLAG INDICATING SNOW COVERAGE (1 FOR SNOW COVER) (-)
SR             1  0  fraction of frozen precipitation (-)
POTEVP         1  0  accumulated potential evaporation (W m-2)
SNOPCX         1  0  snow phase change heat flux (W m-2)
SOILTB         1  0  bottom soil temperature (K)
AVG_FUEL_F     1  0  fuel remaining averaged to atmospheric grid (1)
UAH            1  0  wind at fire_wind_height (m/s)
VAH            1  0  wind at fire_wind_height (m/s)
pressure      27  0  Model pressure (hPa)
geopt         27  0  Geopotential (m2/s2)
height        27  0  Model height (km)
tk            27  0  Temperature (K)
tc            27  0  Temperature (C)
theta         27  0  Potential Temperature (K)
td            27  0  Dewpoint Temperature (C)
td2            1  0  Dewpoint Temperature at 2m (C)
rh            27  0  Relative Humidity (%)
clflo          1  0  Low Cloud Fraction (%)
clfmi          1  0  Mid Cloud Fraction (%)
clfhi          1  0  High Cloud Fraction (%)
rh2            1  0  Relative Humidity at 2m (%)
wspd          27  0  Wind Speed (m s-1)
wdir          27  0  Wind Direction (Degrees)
ws10           1  0  Wind Speed at 10 M (m s-1)
wd10           1  0  Wind Direction at 10 M (Degrees)
umet          27  0  Rotated wind component (m s-1)
vmet          27  0  Rotated wind component (m s-1)
u10m           1  0  Rotated wind component (m s-1)
v10m           1  0  Rotated wind component (m s-1)
slp            1  0  Sea Levelp Pressure (hPa)
dbz           27  0  Reflectivity (-)
max_dbz        1  0  Max Reflectivity (-)
cape          27  0  CAPE (J/kg)
cin           27  0  CIN (J/kg)
mcape          1  0  MCAPE (J/kg)
mcin           1  0  MCIN (J/kg)
lcl            1  0  LCL (meters AGL)
lfc            1  0  LFC (meters AGL)
ENDVARS
@ global String comment TITLE =  OUTPUT FROM WRF V3.3.1 MODEL
@ global String comment START_DATE = 2013-10-31_18:00:00
@ global String comment SIMULATION_START_DATE = 2013-10-31_18:00:00
@ global String comment WEST-EAST_GRID_DIMENSION =   136
@ global String comment SOUTH-NORTH_GRID_DIMENSION =   118
@ global String comment BOTTOM-TOP_GRID_DIMENSION =    28
@ global String comment DX =      3333.33
@ global String comment DY =      3333.33
@ global String comment GRIDTYPE = C
@ global String comment DIFF_OPT =     1
@ global String comment KM_OPT =     4
@ global String comment DAMP_OPT =     0
@ global String comment DAMPCOEF =         0.20
@ global String comment KHDIF =         0.00
@ global String comment KVDIF =         0.00
@ global String comment MP_PHYSICS =     3
@ global String comment RA_LW_PHYSICS =     1
@ global String comment RA_SW_PHYSICS =     1
@ global String comment SF_SFCLAY_PHYSICS =     1
@ global String comment SF_SURFACE_PHYSICS =     2
@ global String comment BL_PBL_PHYSICS =     1
@ global String comment CU_PHYSICS =     0
@ global String comment SURFACE_INPUT_SOURCE =     1
@ global String comment SST_UPDATE =     0
@ global String comment GRID_FDDA =     0
@ global String comment GFDDA_INTERVAL_M =     0
@ global String comment GFDDA_END_H =     0
@ global String comment GRID_SFDDA =     0
@ global String comment SGFDDA_INTERVAL_M =     0
@ global String comment SGFDDA_END_H =     0
@ global String comment SF_URBAN_PHYSICS =     0
@ global String comment FEEDBACK =     1
@ global String comment SMOOTH_OPTION =     0
@ global String comment SWRAD_SCAT =         1.00
@ global String comment W_DAMPING =     0
@ global String comment MOIST_ADV_OPT =     1
@ global String comment SCALAR_ADV_OPT =     1
@ global String comment TKE_ADV_OPT =     1
@ global String comment DIFF_6TH_OPT =     0
@ global String comment DIFF_6TH_FACTOR =         0.12
@ global String comment OBS_NUDGE_OPT =     0
@ global String comment BUCKET_MM =        -1.00
@ global String comment BUCKET_J =        -1.00
@ global String comment PREC_ACC_DT =         0.00
@ global String comment OMLCALL =     0
@ global String comment ISFTCFLX =     0
@ global String comment ISHALLOW =     0
@ global String comment DFI_OPT =     0
@ global String comment SHCU_PHYSICS =     0
@ global String comment WEST-EAST_PATCH_START_UNSTAG =     1
@ global String comment WEST-EAST_PATCH_END_UNSTAG =   135
@ global String comment WEST-EAST_PATCH_START_STAG =     1
@ global String comment WEST-EAST_PATCH_END_STAG =   136
@ global String comment SOUTH-NORTH_PATCH_START_UNSTAG =     1
@ global String comment SOUTH-NORTH_PATCH_END_UNSTAG =   117
@ global String comment SOUTH-NORTH_PATCH_START_STAG =     1
@ global String comment SOUTH-NORTH_PATCH_END_STAG =   118
@ global String comment BOTTOM-TOP_PATCH_START_UNSTAG =     1
@ global String comment BOTTOM-TOP_PATCH_END_UNSTAG =    27
@ global String comment BOTTOM-TOP_PATCH_START_STAG =     1
@ global String comment BOTTOM-TOP_PATCH_END_STAG =    28
@ global String comment GRID_ID =     3
@ global String comment PARENT_ID =     2
@ global String comment I_PARENT_START =    46
@ global String comment J_PARENT_START =    31
@ global String comment PARENT_GRID_RATIO =     3
@ global String comment DT =        20.00
@ global String comment CEN_LAT =        43.36
@ global String comment CEN_LON =        86.50
@ global String comment TRUELAT1 =        30.00
@ global String comment TRUELAT2 =        60.00
@ global String comment MOAD_CEN_LAT =        43.78
@ global String comment STAND_LON =        87.65
@ global String comment POLE_LAT =        90.00
@ global String comment POLE_LON =         0.00
@ global String comment GMT =        18.00
@ global String comment JULYR =  2013
@ global String comment JULDAY =   304
@ global String comment MAP_PROJ =     1
@ global String comment MMINLU = USGS
@ global String comment NUM_LAND_CAT =    24
@ global String comment ISWATER =    16
@ global String comment ISLAKE =    -1
@ global String comment ISICE =    24
@ global String comment ISURBAN =     1
@ global String comment ISOILWATER =    14

编写的Grads作图脚本为：
'reinit'
'open D:/Firsttest/wlmq003.ctl'
'set  lat 42.9 44.5'
'set  lon 86.9 89.1'
'set t 2'
'set  gxout  shaded'
*'set mpdraw on'
*'set mpdraw off'
*'set map auto'
*'set mpdset hres'
'set grads off'
*'d slp'
'set  gxout  shaded'
*'d HGTsfc'
*'d slp'
*'d tk'
'cbarn 1 1'
'set  gxout  shaded'
'd height*1000'
*'print'
'set gxout vector'
'd skip(u10m,5,5);skip(v10m,5,5)'
*'set gxout barb'
*'d skip(u10m,5,5);skip(v10m,5,5)'
*'set z 11'
'print D:/Firsttest/a02shi.eps'
'disable print'
;

请问怎样添加某地的边界线？比如乌鲁木齐市？用shp文件能添加吗？