<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:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D">Hi,
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D"><o:p>&nbsp;</o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D">I don&#8217;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:&quot;Calibri&quot;,&quot;sans-serif&quot;;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:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D">Maybe those are absent.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D"><o:p>&nbsp;</o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:&quot;Calibri&quot;,&quot;sans-serif&quot;;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:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D"><o:p>&nbsp;</o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D">&nbsp; &nbsp;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:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;or<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D">&nbsp;&nbsp;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:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D"><o:p>&nbsp;</o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D">Cheers!<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D"><o:p>&nbsp;</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:&quot;Tahoma&quot;,&quot;sans-serif&quot;">From:</span></b><span lang="EN-US" style="font-size:10.0pt;font-family:&quot;Tahoma&quot;,&quot;sans-serif&quot;"> 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>&nbsp;</o:p></p>
<div>
<p class="MsoNormal">Hi:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; All,I am trying to mask data with shapefile according &nbsp;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">; &nbsp; Example script to produce plots for a WRF real-data run,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; &nbsp; with the ARW coordinate dynamics option.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl&quot;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot; <o:p>
</o:p></p>
</div>
<div>
<p class="MsoNormal">load&quot;./shapefile_mask_data.ncl&quot;<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. &nbsp;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; This needs to have a &quot;.nc&quot; appended, so just do it.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; a = addfile(&quot;/public/home/huanglei/data/20160724/oldwrf_cu5/wrfout_d03_2016-07-24_00:00:00&quot;&#43;&quot;.nc&quot;,&quot;r&quot;)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</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">&nbsp;; type = &quot;x11&quot;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; type = &quot;png&quot;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; type = &quot;ps&quot;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; type = &quot;ncgm&quot;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; wks = gsn_open_wks(type,&quot;mask_plt_Precip_oldwrf5km_zoom_cu5_from24&quot;)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; plot = new(3,graphic)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">; Set some basic resources<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; res = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; res@MainTitle = &quot;REAL-TIME WRF&quot;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; mpres &nbsp;= True &nbsp;; Map resources<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; mpres@mpOutlineOn = False &nbsp;; Turn off map outlines<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; mpres@mpFillOn &nbsp; &nbsp;= False &nbsp;; Turn off map fill<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; mpres@mpGridAndLimbOn = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp;;res@mpProjection &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;Lambert&quot;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; pltres = True ; Plot resources<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; pltres@PanelPlot &nbsp;= True &nbsp; ; Tells wrf_map_overlays not to remove overlays<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</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>&nbsp;</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">&nbsp; FirstTime = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; times = wrf_user_getvar(a,&quot;times&quot;,-1) &nbsp;; get all times in the file<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; ntimes = dimsizes(times) &nbsp; &nbsp; &nbsp; &nbsp; ; number of times in the file<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp;; print(times)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp;; exit<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp;it_start = 10<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp;it_end = 17<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp;; print(&quot;Working on time: &quot; &#43; times(it) )<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp;; res@TimeLabel = times(it) &nbsp; ; Set Valid time to use on plots<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</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 &nbsp; &nbsp; &nbsp; &nbsp;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; ; Get non-convective, convective and total precipitation of 5km &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; rain_exp = wrf_user_getvar(a,&quot;RAINNC&quot;,it_end)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; rain_con = wrf_user_getvar(a,&quot;RAINC&quot;,it_end)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; rain_tot = rain_exp &#43; rain_con<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; rain_tot@description = &quot;Total Precipitation&quot;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal"><span class="apple-tab-span">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>;calculate the precipitation
<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; rain_exp_save = wrf_user_getvar(a,&quot;RAINNC&quot;,it_start)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; rain_con_save = wrf_user_getvar(a,&quot;RAINC&quot;,it_start)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; rain_tot_save = rain_exp_save &#43; rain_con_save<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; times_sav = times(it_start)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; rain_tot_tend = rain_tot - rain_tot_save<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; rainc_tend = rain_con - rain_con_save &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; CUMULUS PRECIPITATION<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; rainnc_tend= rain_exp - rain_exp_save &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; SCALE PRECIPITATION<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; rain_tot_tend@description = &quot;Precipitation of 5km(old wrf)&quot;<o:p></o:p></p>
</div>
<div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp;&nbsp;rain_tot_tend@lat2d = a-&gt;XLAT(1,:,:)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp;&nbsp;rain_tot_tend@lon2d = a-&gt;XLONG(1,:,:)<o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; rainc_tend@description = &quot;RAINC of 5km(old wrf)&quot;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; rainnc_tend@description = &quot;RAINNC of 5km(old wrf)&quot;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;mask data with shapefile;;;;;;;;;;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">shp_filename = &quot;/public/home/huanglei/map/xian.shp&quot;<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>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; ; Plotting options for Precipitation<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; opts_r = res &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; opts_r@UnitLabel &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;mm&quot;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; opts_r@cnLevelSelectionMode = &quot;ExplicitLevels&quot;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; opts_r@cnLevels &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = (/ .1, .2, .4, .8, 1.6, 3.2, 6.4, \<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 12.8, 25.6, 51.2/)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; opts_r@cnFillColors &nbsp; &nbsp; &nbsp; &nbsp; = (/&quot;White&quot;,&quot;DarkOliveGreen1&quot;, \<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &quot;DarkOliveGreen3&quot;,&quot;Chartreuse&quot;, \<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &quot;Chartreuse3&quot;,&quot;Green&quot;,&quot;ForestGreen&quot;, \<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &quot;Yellow&quot;,&quot;Orange&quot;,&quot;Red&quot;,&quot;Violet&quot;/)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; opts_r@cnInfoLabelOn &nbsp; &nbsp; &nbsp; &nbsp;= False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; opts_r@cnConstFLabelOn &nbsp; &nbsp; &nbsp;= False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; opts_r@cnFillOn &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; ; &nbsp;opts_r@vpHeightF &nbsp; = 0.1<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; ; &nbsp;opts_r@vpWidthF &nbsp; &nbsp;= 0.9 <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; opts_r@gsnDraw &nbsp; &nbsp; &nbsp;= &nbsp;False &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; opts_r@gsnFrame &nbsp; &nbsp; = &nbsp;False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; opts_r@lbLabelBarOn &nbsp;= False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; opts_r@Footer &nbsp; &nbsp; = False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; opts_r@NoHeaderFooter =True <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; ;;;;;;set zoom ;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;tes = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;tes@returnInt = False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;loc1=wrf_user_ll_to_ij(a,107.50,33.5,tes)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;print(&quot;X/Y location is: &quot;&#43; loc1)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;loc2=wrf_user_ll_to_ij(a,110.0,35.2,tes)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;print(&quot;X/Y location is: &quot;&#43; loc2)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; &nbsp;exit<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;x_start = 35<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;x_end &nbsp; = 81<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;y_start = 47<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;y_end &nbsp; = 84<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; mpres1 = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; mpres1@ZoomIn = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; mpres1@Xstart = x_start<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; mpres1@Ystart = y_start<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; mpres1@Xend &nbsp; = x_end<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; mpres1@Yend &nbsp; = y_end<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; rain_tot_zoom = rain_tot_tend(y_start:y_end,x_start:x_end)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; rainc_zoom = rainc_tend(y_start:y_end,x_start:x_end)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; rainnc_zoom = rainnc_tend(y_start:y_end,x_start:x_end)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; ; Precipitation Tendencies<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp;; &nbsp; opts_r@SubFieldTitle = &quot;from &quot; &#43; times(it_start) &#43; &quot; to &quot; &#43; times(it_end)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; contour_tend = wrf_contour(a,wks, rain_tot_mask,opts_r) ; total (color)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; 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">&nbsp; &nbsp; &nbsp; &nbsp; 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>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; delete(opts_r)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; ; MAKE PLOTS &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; <o:p>
</o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;plot(0) = wrf_map_overlays(a,wks,contour_tend,pltres,mpres1)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;plot(1) = wrf_map_overlays(a,wks,contour_rainc_tend,pltres,mpres1)
<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;plot(2) = wrf_map_overlays(a,wks,contour_rainnc_tend,pltres,mpres1)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">;&gt;============================================================&lt;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;add China map<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">;&gt;------------------------------------------------------------&lt;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; &nbsp; &nbsp;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; shp_name1 &nbsp; &nbsp;= &quot;/public/home/huanglei/map/xian.shp&quot;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; lnres &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; lnres@gsLineColor &nbsp; &nbsp; &nbsp;= &quot;gray25&quot;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; lnres@gsLineThicknessF = 0.5 &nbsp; <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp;id = new(3,graphic)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp;id(0) = gsn_add_shapefile_polylines(wks,plot(0),shp_name1,lnres)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp;id(1) = gsn_add_shapefile_polylines(wks,plot(1),shp_name1,lnres)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp;id(2) = gsn_add_shapefile_polylines(wks,plot(2),shp_name1,lnres)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp;shp_name2 &nbsp; &nbsp;= &quot;/public/home/huanglei/map/shaanxi_city_l.shp&quot;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; prres=True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; prres@gsLineThicknessF = 1.0<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; prres@gsLineColor = &quot;black&quot;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp;id2 = new(3,graphic)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp;id2(0) = gsn_add_shapefile_polylines(wks,plot(0),shp_name2,prres)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp;id2(1) = gsn_add_shapefile_polylines(wks,plot(1),shp_name2,prres)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp;id2(2) = gsn_add_shapefile_polylines(wks,plot(2),shp_name2,prres)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;creat panel;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; resP &nbsp; &nbsp; = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; resP@gsnMaximize = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; resP@gsnFrame &nbsp;=False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; resP@gsnPanelLabelBar = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; resP@gsnPanelBottom = 0.05<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; resP@gsnPanelTop &nbsp; &nbsp; &nbsp; = 0.95 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; Make sure not too close to
<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; resP@gsnPanelBottom &nbsp; &nbsp;= 0.001 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; edge, so it maximizes better.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; resP@gsnPanelLabelBar &nbsp;= True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; resP@lbLabelFontHeightF = &nbsp;0.006<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; resP@gsnPanelFigureStrings= (/&quot;Total precipitation&quot;,&quot;RAINC&quot;,&quot;RAINNC&quot;/) ; add strings to panel<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; resP@amJust &nbsp; = &quot;TopRight&quot; <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; resP@gsnPanelFigureStringsFontHeightF = 0.005<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; resP@lbLabelAutoStride = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; resP@gsnPaperOrientation = &quot;landscape&quot;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; gsn_panel(wks,plot,(/1,3/),resP)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp;;psres = True &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;
<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; maximize_output(wks,psres) &nbsp;; calls draw and frame for you
<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">; draw(plot) &nbsp; &nbsp; &nbsp; ; This will draw the map and the shapefile outlines.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp; frame(wks)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">end<o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:&quot;Verdana&quot;,&quot;sans-serif&quot;">&nbsp;How can I slove the problem?</span><o:p></o:p></p>
</div>
</div>
</body>
</html>