<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 12 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:"Malgun Gothic";
        panose-1:2 11 5 3 2 0 0 2 0 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:Tahoma;
        panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
        {font-family:Verdana;
        panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
        {font-family:"\@Malgun Gothic";
        panose-1:2 11 5 3 2 0 0 2 0 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
span.apple-tab-span
        {mso-style-name:apple-tab-span;}
span.EmailStyle18
        {mso-style-type:personal-reply;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
.MsoChpDefault
        {mso-style-type:export-only;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-GB" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Hi,
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">I don’t know exactly what reason is cause that problem.
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">I recommend to check coordinate variables of rain_tot_mask(using printVarSummary() function).<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Maybe those are absent.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Or use this command before plotting script.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"> Copy_varCoords(rain_tot_tend, rain_tot_mask)
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"> or<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"> Copy_varMeta(rain_tot_tend, rain_tot_mask)<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Cheers!<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span lang="EN-US" style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> ncl-talk-bounces@ucar.edu [mailto:ncl-talk-bounces@ucar.edu]
<b>On Behalf Of </b>grace<br>
<b>Sent:</b> 05 August 2016 10:14<br>
<b>To:</b> ncl-talk<br>
<b>Subject:</b> [ncl-talk] data did not match with original plot after mask data with shapefile<o:p></o:p></span></p>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">Hi:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> All,I am trying to mask data with shapefile according to the example shapefiles_14_mask.ncl on the NCL website.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">But the masked plot did not match the original plot,Ican not figure out what's wrong with it.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Left plot is the masked plot,right is the original one.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><img width="751" height="295" id="_x0000_i1025" src="cid:image001.jpg@01D1EF04.A60EF200"><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">This is my script:<o:p></o:p></p>
</div>
<div>
<div>
<p class="MsoNormal">; Example script to produce plots for a WRF real-data run,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; with the ARW coordinate dynamics option.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" <o:p>
</o:p></p>
</div>
<div>
<p class="MsoNormal">load"./shapefile_mask_data.ncl"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">begin<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; The WRF ARW input file. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; This needs to have a ".nc" appended, so just do it.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> a = addfile("/public/home/huanglei/data/20160724/oldwrf_cu5/wrfout_d03_2016-07-24_00:00:00"+".nc","r")<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">; We generate plots, but what kind do we prefer?<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ; type = "x11"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> type = "png"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; type = "ps"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; type = "ncgm"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> wks = gsn_open_wks(type,"mask_plt_Precip_oldwrf5km_zoom_cu5_from24")<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> plot = new(3,graphic)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">; Set some basic resources<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> res = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> res@MainTitle = "REAL-TIME WRF"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> mpres = True ; Map resources<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> mpres@mpOutlineOn = False ; Turn off map outlines<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> mpres@mpFillOn = False ; Turn off map fill<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> mpres@mpGridAndLimbOn = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ;res@mpProjection = "Lambert"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> pltres = True ; Plot resources<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> pltres@PanelPlot = True ; Tells wrf_map_overlays not to remove overlays<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">; What times and how many time steps are in the data set?<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> FirstTime = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> times = wrf_user_getvar(a,"times",-1) ; get all times in the file<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ntimes = dimsizes(times) ; number of times in the file<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ; print(times)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ; exit<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> it_start = 10<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> it_end = 17<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ; print("Working on time: " + times(it) )<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ; res@TimeLabel = times(it) ; Set Valid time to use on plots<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; First get the variables we will need <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> ; Get non-convective, convective and total precipitation of 5km <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_exp = wrf_user_getvar(a,"RAINNC",it_end)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_con = wrf_user_getvar(a,"RAINC",it_end)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_tot = rain_exp + rain_con<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_tot@description = "Total Precipitation"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><span class="apple-tab-span"> </span>;calculate the precipitation
<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_exp_save = wrf_user_getvar(a,"RAINNC",it_start)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_con_save = wrf_user_getvar(a,"RAINC",it_start)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_tot_save = rain_exp_save + rain_con_save<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> times_sav = times(it_start)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_tot_tend = rain_tot - rain_tot_save<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rainc_tend = rain_con - rain_con_save ; CUMULUS PRECIPITATION<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rainnc_tend= rain_exp - rain_exp_save ; SCALE PRECIPITATION<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_tot_tend@description = "Precipitation of 5km(old wrf)"<o:p></o:p></p>
</div>
<div>
<div>
<p class="MsoNormal"> rain_tot_tend@lat2d = a->XLAT(1,:,:)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_tot_tend@lon2d = a->XLONG(1,:,:)<o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> rainc_tend@description = "RAINC of 5km(old wrf)"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rainnc_tend@description = "RAINNC of 5km(old wrf)"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;mask data with shapefile;;;;;;;;;;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">shp_filename = "/public/home/huanglei/map/xian.shp"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">rain_tot_mask = shapefile_mask_data(rain_tot_tend,shp_filename,True)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> ; Plotting options for Precipitation<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r = res <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@UnitLabel = "mm"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@cnLevelSelectionMode = "ExplicitLevels"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@cnLevels = (/ .1, .2, .4, .8, 1.6, 3.2, 6.4, \<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> 12.8, 25.6, 51.2/)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@cnFillColors = (/"White","DarkOliveGreen1", \<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> "DarkOliveGreen3","Chartreuse", \<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> "Chartreuse3","Green","ForestGreen", \<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> "Yellow","Orange","Red","Violet"/)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@cnInfoLabelOn = False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@cnConstFLabelOn = False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@cnFillOn = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ; opts_r@vpHeightF = 0.1<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ; opts_r@vpWidthF = 0.9 <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@gsnDraw = False <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@gsnFrame = False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@lbLabelBarOn = False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@Footer = False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@NoHeaderFooter =True <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ;;;;;;set zoom ;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> tes = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> tes@returnInt = False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> loc1=wrf_user_ll_to_ij(a,107.50,33.5,tes)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> print("X/Y location is: "+ loc1)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> loc2=wrf_user_ll_to_ij(a,110.0,35.2,tes)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> print("X/Y location is: "+ loc2)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ; exit<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> x_start = 35<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> x_end = 81<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> y_start = 47<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> y_end = 84<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> mpres1 = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> mpres1@ZoomIn = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> mpres1@Xstart = x_start<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> mpres1@Ystart = y_start<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> mpres1@Xend = x_end<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> mpres1@Yend = y_end<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_tot_zoom = rain_tot_tend(y_start:y_end,x_start:x_end)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rainc_zoom = rainc_tend(y_start:y_end,x_start:x_end)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rainnc_zoom = rainnc_tend(y_start:y_end,x_start:x_end)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ; Precipitation Tendencies<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ; opts_r@SubFieldTitle = "from " + times(it_start) + " to " + times(it_end)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> contour_tend = wrf_contour(a,wks, rain_tot_mask,opts_r) ; total (color)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> contour_rainc_tend = wrf_contour(a,wks, rainc_zoom,opts_r) ; total cumulus precipitation (color)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> contour_rainnc_tend = wrf_contour(a,wks, rainnc_zoom,opts_r) ; total scale precipitation(color)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> delete(opts_r)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> ; MAKE PLOTS <o:p>
</o:p></p>
</div>
<div>
<p class="MsoNormal"> plot(0) = wrf_map_overlays(a,wks,contour_tend,pltres,mpres1)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> plot(1) = wrf_map_overlays(a,wks,contour_rainc_tend,pltres,mpres1)
<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> plot(2) = wrf_map_overlays(a,wks,contour_rainnc_tend,pltres,mpres1)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">;>============================================================<<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; add China map<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">;>------------------------------------------------------------<<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> shp_name1 = "/public/home/huanglei/map/xian.shp"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> lnres = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> lnres@gsLineColor = "gray25"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> lnres@gsLineThicknessF = 0.5 <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> id = new(3,graphic)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> id(0) = gsn_add_shapefile_polylines(wks,plot(0),shp_name1,lnres)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> id(1) = gsn_add_shapefile_polylines(wks,plot(1),shp_name1,lnres)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> id(2) = gsn_add_shapefile_polylines(wks,plot(2),shp_name1,lnres)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> shp_name2 = "/public/home/huanglei/map/shaanxi_city_l.shp"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> prres=True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> prres@gsLineThicknessF = 1.0<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> prres@gsLineColor = "black"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> id2 = new(3,graphic)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> id2(0) = gsn_add_shapefile_polylines(wks,plot(0),shp_name2,prres)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> id2(1) = gsn_add_shapefile_polylines(wks,plot(1),shp_name2,prres)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> id2(2) = gsn_add_shapefile_polylines(wks,plot(2),shp_name2,prres)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;creat panel;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@gsnMaximize = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@gsnFrame =False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@gsnPanelLabelBar = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@gsnPanelBottom = 0.05<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@gsnPanelTop = 0.95 ; Make sure not too close to
<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@gsnPanelBottom = 0.001 ; edge, so it maximizes better.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@gsnPanelLabelBar = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@lbLabelFontHeightF = 0.006<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@gsnPanelFigureStrings= (/"Total precipitation","RAINC","RAINNC"/) ; add strings to panel<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@amJust = "TopRight" <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@gsnPanelFigureStringsFontHeightF = 0.005<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@lbLabelAutoStride = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@gsnPaperOrientation = "landscape"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> gsn_panel(wks,plot,(/1,3/),resP)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> ;psres = True
<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; maximize_output(wks,psres) ; calls draw and frame for you
<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">; draw(plot) ; This will draw the map and the shapefile outlines.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> frame(wks)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">end<o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Verdana","sans-serif""> How can I slove the problem?</span><o:p></o:p></p>
</div>
</div>
</body>
</html>