<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Rick, it worked! Thank you so much for your time and investigation into this! I deleted the islands in QGIS, ran the script again, and now it is working beautifully!
<div class=""><br class="">
</div>
<div class="">Best regards,</div>
<div class="">Toni</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class=""><b style="color: rgb(51, 51, 153); background-color: rgb(255, 255, 255);" class="">Toni Klemm</b><br style="color: rgb(51, 51, 153); background-color: rgb(255, 255, 255);" class="">
<span style="color: rgb(51, 51, 153); background-color: rgb(255, 255, 255);" class="">Ph.D. Candidate</span><br style="color: rgb(51, 51, 153); background-color: rgb(255, 255, 255);" class="">
<span style="color: rgb(51, 51, 153); background-color: rgb(255, 255, 255);" class="">The University of Oklahoma</span><br style="color: rgb(51, 51, 153); background-color: rgb(255, 255, 255);" class="">
<span style="color: rgb(51, 51, 153); background-color: rgb(255, 255, 255);" class="">South Central Climate Science Center</span><br style="color: rgb(51, 51, 153); background-color: rgb(255, 255, 255);" class="">
</div>
<div class=""><span style="color: rgb(51, 51, 153); background-color: rgb(255, 255, 255);" class="">phone: 405 325 1272</span></div>
<div class=""><span style="color: rgb(51, 51, 153); background-color: rgb(255, 255, 255);" class=""><a href="http://www.toni-klemm.de" class="">www.toni-klemm.de</a></span></div>
<div class=""><a href="https://www.eccforum.org" class="">Early Career Climate Forum</a></div>
<div class="">Twitter:&nbsp;<a href="https://twitter.com/ToniKlemm" class="">@toniklemm</a></div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br class="Apple-interchange-newline">
</div>
<br class="">
<div>
<blockquote type="cite" class="">
<div class="">On May 25, 2017, at 3:15 PM, Rick Brownrigg &lt;<a href="mailto:brownrig@ucar.edu" class="">brownrig@ucar.edu</a>&gt; wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div dir="ltr" class="">
<div class="">
<div class="">
<div class="">
<div class="">
<div class="">
<div class="">
<div class="">
<div class="">
<div class="">Hi Toni,<br class="">
<br class="">
</div>
Sorry for the delay...this took a bit to figure out.<br class="">
<br class="">
</div>
This *may* be a bug in gc_inout, or perhaps the documentation is misleading. The issue is that although the shapefile says it contains 1 feature, that feature is comprised of 6 &quot;segments&quot;:&nbsp; 1 for the large 4 state area and 5 for the islands off the coast of
 Texas (see the attached png).&nbsp; This apparently is confusing the in/out tests.<br class="">
<br class="">
</div>
I don't know of a good way of programmatically working around this.&nbsp; As a hack (an extreme hack!), I know from looking at the shapefile that the large 4-state area's coordinates begin at index 56, and extend to the end of the x/y arrays. So if I change your
 script from:<br class="">
<br class="">
&nbsp; shp = addfile(&quot;./<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__4-5Fstates-5Fmerged.nc&amp;d=DwMFaQ&amp;c=qKdtBuuu6dQK9MsRUVJ2DPXW6oayO8fu4TfEHS8sGNk&amp;r=trcSnQNnUT2ddBXcpd358A&amp;m=yJKWezSUje6we6fxpkzLjE7jtiGQ1kGaXuaMupFyefY&amp;s=7NPENOYYTc5LXWJHVC9CzI7gs6KoCRKqxMV-orMtQfo&amp;e=" class="">4_states_merged.nc</a>&quot;,&quot;r&quot;)<br class="">
&nbsp; shp_lat = shp-&gt;y<br class="">
&nbsp; shp_lon = shp-&gt;x&#43;360&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; adding 360 to get shapefile longitudes in same value range as data longitudes<br class="">
<br class="">
</div>
to<br class="">
<br class="">
&nbsp; shp = addfile(&quot;./<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__4-5Fstates-5Fmerged.nc&amp;d=DwMFaQ&amp;c=qKdtBuuu6dQK9MsRUVJ2DPXW6oayO8fu4TfEHS8sGNk&amp;r=trcSnQNnUT2ddBXcpd358A&amp;m=yJKWezSUje6we6fxpkzLjE7jtiGQ1kGaXuaMupFyefY&amp;s=7NPENOYYTc5LXWJHVC9CzI7gs6KoCRKqxMV-orMtQfo&amp;e=" class="">4_states_merged.nc</a>&quot;,&quot;r&quot;)<br class="">
&nbsp; shp_lat = shp-&gt;y(56:)<br class="">
&nbsp; shp_lon = shp-&gt;x(56:)&#43;360&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; adding 360 to get shapefile longitudes in same value range as data longitudes<br class="">
<br class="">
</div>
Then it works, as seen in the second png. <br class="">
<br class="">
</div>
The docs for gc_inout() suggest one should be able to pass in multiple polygons to test against, but as practical matter, I don't see how they would be delineated in the lat/lon arguments, so I suspect the docs are wrong, or maybe the implementation incomplete.<br class="">
<br class="">
</div>
If you have access to something like ArcGIS that can edit shapefiles, the better solution would be to create a simpler mask excluding the island polygons.<br class="">
<br class="">
</div>
Wish I had a better answer, but I hope that helps.<br class="">
</div>
Rick<br class="">
<br class="">
<br class="">
<div class="">
<div class="">
<div class="">
<div class="">
<div class=""><br class="">
</div>
</div>
</div>
</div>
</div>
</div>
<div class="gmail_extra"><br class="">
<div class="gmail_quote">On Wed, May 24, 2017 at 3:19 PM, Klemm, Toni <span dir="ltr" class="">
&lt;<a href="mailto:toni@ou.edu" target="_blank" class="">toni@ou.edu</a>&gt;</span> wrote:<br class="">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="">
<div style="word-wrap:break-word" class="">Hi Rick,
<div class=""><br class="">
</div>
<div class="">Thanks for getting back to me!</div>
<div class=""><br class="">
</div>
<div class="">I got it to work in the meantime using masking example 9 (<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.ncl.ucar.edu_Applications_Scripts_mask-5F9.ncl&amp;d=DwMFaQ&amp;c=qKdtBuuu6dQK9MsRUVJ2DPXW6oayO8fu4TfEHS8sGNk&amp;r=trcSnQNnUT2ddBXcpd358A&amp;m=yJKWezSUje6we6fxpkzLjE7jtiGQ1kGaXuaMupFyefY&amp;s=ylM244pXeKWXrgrLrfc52xrVmOM31ChELZkdcBRprSw&amp;e=" target="_blank" class="">https://www.ncl.ucar.edu/<wbr class="">Applications/Scripts/mask_9.<wbr class="">ncl</a>),
 which looks similar to example 4 on the page you linked (<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__ncl.ucar.edu_Applications_Scripts_shapefiles-5F4.ncl&amp;d=DwMFaQ&amp;c=qKdtBuuu6dQK9MsRUVJ2DPXW6oayO8fu4TfEHS8sGNk&amp;r=trcSnQNnUT2ddBXcpd358A&amp;m=yJKWezSUje6we6fxpkzLjE7jtiGQ1kGaXuaMupFyefY&amp;s=6JyGs2nQ3Pgnkmm4jcgRn8kB0Ks0IPIHtcNyKNEPwQw&amp;e=" target="_blank" class="">http://ncl.ucar.edu/<wbr class="">Applications/Scripts/<wbr class="">shapefiles_4.ncl</a>).
 The masking worked, but for some reason there is large wedge with missing data right in the middle of my plot now. I imagine it has somewhere to do with the “<i class="">if (MASK_INSIDE) then</i>” loop in my script below, but I can’t figure it out. Do you
 have any suggestions? I also attached my original data file, the resulting datafile (_<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__shp.nc&amp;d=DwMFaQ&amp;c=qKdtBuuu6dQK9MsRUVJ2DPXW6oayO8fu4TfEHS8sGNk&amp;r=trcSnQNnUT2ddBXcpd358A&amp;m=yJKWezSUje6we6fxpkzLjE7jtiGQ1kGaXuaMupFyefY&amp;s=9z-ZdsoHo55bNj8sxbMN37AfvDWtsfEEJBhFNHPs7qw&amp;e=" target="_blank" class="">shp.nc</a>)
 and plot, and the .nc file of my study area. This shapefile was created in ArcMap 10.3 from the states.* files in the&nbsp;<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.ncl.ucar.edu_Applications_Data_-23shp&amp;d=DwMFaQ&amp;c=qKdtBuuu6dQK9MsRUVJ2DPXW6oayO8fu4TfEHS8sGNk&amp;r=trcSnQNnUT2ddBXcpd358A&amp;m=yJKWezSUje6we6fxpkzLjE7jtiGQ1kGaXuaMupFyefY&amp;s=VUSaP80OlWKgR2W9u6MrYascvWy6DMd8D4Jr1ocakfI&amp;e=" target="_blank" class="">NCL
 shapefile set</a>&nbsp;and converted to .nc using ncl_convert2nc.&nbsp;</div>
<div class=""><br class="">
</div>
<div class="">Again, thank you very much!</div>
<div class=""><br class="">
</div>
<div class="">Best,</div>
<div class="">Toni</div>
<div class=""><br class="">
</div>
<div class="">Script:</div>
<div class=""><br class="">
</div>
<div class="">
<div class=""><i class="">; ******************************<wbr class="">******************************<wbr class="">*****</i></div>
<div class=""><i class="">; Create map of GFDL-OBS bias averages by FORECAST MONTH 1980-2014</i></div>
<div class=""><i class="">;&nbsp;</i></div>
<div class=""><i class="">; Use SHAPEFILE OF STUDY REGION to define subset and plot</i></div>
<div class=""><i class="">; ******************************<wbr class="">******************************<wbr class="">*****</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">load &quot;$NCARG_ROOT/lib/ncarg/<wbr class="">nclscripts/csm/gsn_code.ncl&quot;</i></div>
<div class=""><i class="">load &quot;$NCARG_ROOT/lib/ncarg/<wbr class="">nclscripts/csm/gsn_csm.ncl&quot;</i></div>
<div class=""><i class="">load &quot;$NCARG_ROOT/lib/ncarg/<wbr class="">nclscripts/csm/contributed.<wbr class="">ncl&quot;</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">; ******************************<wbr class="">******************************</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">begin</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">do target_month = 1,12<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>&nbsp;; do-loop through the 12 calendar months</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; <span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>if (target_month .lt. 10) then</i></div>
<div class=""><i class="">&nbsp;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>&nbsp; target_month_zero = &quot;0&quot; &#43; target_month<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; turns &quot;1&quot;...&quot;9&quot; into &quot;01&quot;...&quot;09&quot; to match filename structure (months have 2 digits&quot;</i></div>
<div class=""><i class="">&nbsp;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>end if</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp;if (target_month .gt. 9) then</i></div>
<div class=""><i class="">&nbsp; &nbsp; target_month_zero = target_month<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; leaves &quot;10...&quot;12&quot; unchanged</i></div>
<div class=""><i class="">&nbsp;end if</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp;print(&quot;**** &nbsp;START: target_month &quot; &#43; target_month &#43; &quot; &nbsp;****&quot;)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; ; ******************************<wbr class="">***********************</i></div>
<div class=""><i class="">&nbsp; ; READ IN FILE INFORMATION BY TARGET MONTH</i></div>
<div class=""><i class="">&nbsp; ; ******************************<wbr class="">***********************</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">; *****</i></div>
<div class=""><i class="">print(&quot;1 - Read in file information&quot;)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; bias_f = addfile(&quot;~/GFDL/TAS/DATAFILES/<wbr class="">monthly/5b-4alt_GFDL_TAS_bias/<wbr class="">avg_by_target_month/stack/<wbr class="">GFDL_TAS_bias_1980-2014_<wbr class="">target_&quot; &#43; target_month_zero &#43; &quot;.nc&quot;,&quot;r&quot;)</i></div>
<div class=""><i class="">&nbsp; print(&quot;Finished: 1/8 Read in multiple files&quot;)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; ; ******************************<wbr class="">*****************</i></div>
<div class=""><i class="">&nbsp; ; READ IN THE VARIABLE (precipitation)</i></div>
<div class=""><i class="">&nbsp; ; ******************************<wbr class="">*****************</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; print(&quot;2 - Read in variable and dimensions and average all years&quot;)</i></div>
<div class=""><i class="">&nbsp; lat = (/bias_f-&gt;lat/)</i></div>
<div class=""><i class="">&nbsp; lon = (/bias_f-&gt;lon/)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; printMinMax(lat, False)</i></div>
<div class=""><i class="">&nbsp; printMinMax(lon, False)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; nlat = dimsizes(lat)</i></div>
<div class=""><i class="">&nbsp; lat!0 = &quot;lat&quot;</i></div>
<div class=""><i class="">&nbsp; lat&amp;lat = lat</i></div>
<div class=""><i class="">&nbsp; lat@long_name = &quot;latitude&quot;</i></div>
<div class=""><i class="">&nbsp; lat@units = &quot;degrees_north&quot;</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; mlon = dimsizes(lon)</i></div>
<div class=""><i class="">&nbsp; lon!0 = &quot;lon&quot;</i></div>
<div class=""><i class="">&nbsp; lon&amp;lon = lon</i></div>
<div class=""><i class="">&nbsp; lon@long_name = &quot;longitude&quot;</i></div>
<div class=""><i class="">&nbsp; lon@units = &quot;degrees_east&quot;</i></div>
<div class=""><i class="">&nbsp;&nbsp;</i></div>
<div class=""><i class="">&nbsp; GFDL_bias_3D = rm_single_dims(bias_f-&gt;GFDL_<wbr class="">OBS_bias)<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span><span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap"></span>; extract precipitation into new variable, output is mm/day</i></div>
<div class=""><i class="">&nbsp; GFDL_bias_2D = dim_avg_n(GFDL_bias_3D,0) &nbsp; ; averaging all bias for GFDL and LTM made in that month for 1985-2010</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; min_bias = decimalPlaces(min(GFDL_bias_<wbr class="">2D),2,True)</i></div>
<div class=""><i class="">&nbsp; max_bias = decimalPlaces(max(GFDL_bias_<wbr class="">2D),2,True)</i></div>
<div class=""><i class="">&nbsp; avg_bias = decimalPlaces(avg(GFDL_bias_<wbr class="">2D),2,True)</i></div>
<div class=""><i class="">&nbsp;&nbsp;</i></div>
<div class=""><i class="">; &nbsp;printVarSummary(GFDL_bias_2D)</i></div>
<div class=""><i class="">&nbsp;&nbsp;</i></div>
<div class=""><i class="">; &nbsp;GFDL_bias_2D@_FillValue = 1.0E20</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">; &nbsp;printVarSummary(GFDL_bias_2D)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">; &nbsp;printVarSummary(GFDL_bias_3D)</i></div>
<div class=""><i class="">; &nbsp;printVarSummary(GFDL_bias_2D)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">; ******************************<wbr class="">***********</i></div>
<div class=""><i class="">; ***** Read in shapefile information ***** (source: <a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.ncl.ucar.edu_Applications_Scripts_mask-5F9.ncl&amp;d=DwMFaQ&amp;c=qKdtBuuu6dQK9MsRUVJ2DPXW6oayO8fu4TfEHS8sGNk&amp;r=trcSnQNnUT2ddBXcpd358A&amp;m=yJKWezSUje6we6fxpkzLjE7jtiGQ1kGaXuaMupFyefY&amp;s=ylM244pXeKWXrgrLrfc52xrVmOM31ChELZkdcBRprSw&amp;e=" target="_blank" class="">
https://www.ncl.ucar.edu/<wbr class="">Applications/Scripts/mask_9.<wbr class="">ncl</a>)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; MASK_INSIDE = True &nbsp; &nbsp;; Whether to mask data inside or outside the given geographical area.</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; shp = addfile(&quot;~/States_shapefile_<wbr class="">data/NCL/4_states_merged.shp&quot;,<wbr class="">&quot;r&quot;)</i></div>
<div class=""><i class="">&nbsp; shp_lat = shp-&gt;y</i></div>
<div class=""><i class="">&nbsp; shp_lon = shp-&gt;x&#43;360 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; adding 360 to get shapefile longitudes in same value range as data longitudes</i></div>
<div class=""><i class="">; &nbsp;printVarSummary(shp_lat)</i></div>
<div class=""><i class="">; &nbsp;printVarSummary(shp_lon)</i></div>
<div class=""><i class="">&nbsp; printMinMax(shp_lat, False)</i></div>
<div class=""><i class="">&nbsp; printMinMax(shp_lon, False)</i></div>
<div class=""><i class="">&nbsp;&nbsp;</i></div>
<div class=""><i class="">&nbsp; nshp = dimsizes(shp_lon)</i></div>
<div class=""><i class="">&nbsp;&nbsp;</i></div>
<div class=""><i class="">&nbsp; min_shp_lat = min(shp_lat)</i></div>
<div class=""><i class="">&nbsp; max_shp_lat = max(shp_lat)</i></div>
<div class=""><i class="">&nbsp; min_shp_lon = min(shp_lon)</i></div>
<div class=""><i class="">&nbsp; max_shp_lon = max(shp_lon)</i></div>
<div class=""><i class="">&nbsp;&nbsp;</i></div>
<div class=""><i class="">&nbsp;&nbsp;</i></div>
<div class=""><i class="">; ***** Create new subset based on shapefile ***** (source:
<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.ncl.ucar.edu_Applications_Scripts_mask-5F9.ncl&amp;d=DwMFaQ&amp;c=qKdtBuuu6dQK9MsRUVJ2DPXW6oayO8fu4TfEHS8sGNk&amp;r=trcSnQNnUT2ddBXcpd358A&amp;m=yJKWezSUje6we6fxpkzLjE7jtiGQ1kGaXuaMupFyefY&amp;s=ylM244pXeKWXrgrLrfc52xrVmOM31ChELZkdcBRprSw&amp;e=" target="_blank" class="">
https://www.ncl.ucar.edu/<wbr class="">Applications/Scripts/mask_9.<wbr class="">ncl</a>)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; if(MASK_INSIDE) then</i></div>
<div class=""><i class="">;---Start with data filled in.</i></div>
<div class=""><i class="">&nbsp; &nbsp; data_mask = GFDL_bias_2D</i></div>
<div class=""><i class="">&nbsp; else</i></div>
<div class=""><i class="">;---Start with data all missing</i></div>
<div class=""><i class="">&nbsp; data_mask = new(dimsizes(GFDL_bias_2D),<wbr class="">typeof(GFDL_bias_2D),GFDL_<wbr class="">bias_2D@_FillValue)</i></div>
<div class=""><i class="">&nbsp; copy_VarCoords(GFDL_bias_2D,<wbr class="">data_mask)</i></div>
<div class=""><i class="">&nbsp; end if</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">; ***** Determine the boundary lat/lons of the subset region within the lat/lon box&nbsp;</i></div>
<div class=""><i class="">&nbsp; ilt_mn = ind(min_shp_lat.gt.lat)</i></div>
<div class=""><i class="">&nbsp; ilt_mx = ind(max_shp_lat.lt.lat)</i></div>
<div class=""><i class="">&nbsp; iln_mn = ind(min_shp_lon.gt.lon)</i></div>
<div class=""><i class="">&nbsp; iln_mx = ind(max_shp_lon.lt.lon)</i></div>
<div class=""><i class="">&nbsp; ilt1 &nbsp; = ilt_mn(dimsizes(ilt_mn)-1) &nbsp; &nbsp;; Start of lat box</i></div>
<div class=""><i class="">&nbsp; iln1 &nbsp; = iln_mn(dimsizes(iln_mn)-1) &nbsp; &nbsp;; Start of lon box</i></div>
<div class=""><i class="">&nbsp; ilt2 &nbsp; = ilt_mx(0) &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; End of lat box</i></div>
<div class=""><i class="">&nbsp; iln2 &nbsp; = iln_mx(0) &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; End of lon box</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; if (MASK_INSIDE) then</i></div>
<div class=""><i class="">; ***** Put missing values in the areas that we want masked. (source:
<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.ncl.ucar.edu_Applications_Scripts_mask-5F9.ncl&amp;d=DwMFaQ&amp;c=qKdtBuuu6dQK9MsRUVJ2DPXW6oayO8fu4TfEHS8sGNk&amp;r=trcSnQNnUT2ddBXcpd358A&amp;m=yJKWezSUje6we6fxpkzLjE7jtiGQ1kGaXuaMupFyefY&amp;s=ylM244pXeKWXrgrLrfc52xrVmOM31ChELZkdcBRprSw&amp;e=" target="_blank" class="">
https://www.ncl.ucar.edu/<wbr class="">Applications/Scripts/mask_9.<wbr class="">ncl</a>)</i></div>
<div class=""><i class="">&nbsp; &nbsp; do ilt=ilt1,ilt2</i></div>
<div class=""><i class="">&nbsp; &nbsp; &nbsp; do iln=iln1,iln2</i></div>
<div class=""><i class="">&nbsp; &nbsp; &nbsp; &nbsp; if(gc_inout(lat(ilt),lon(iln),<wbr class="">shp_lat,shp_lon)) then</i></div>
<div class=""><i class="">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; data_mask(ilt,iln) = data_mask@_FillValue &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; replace data inside and outside study region with missing values</i></div>
<div class=""><i class="">&nbsp; &nbsp; &nbsp; &nbsp; end if</i></div>
<div class=""><i class="">&nbsp; &nbsp; &nbsp; end do</i></div>
<div class=""><i class="">&nbsp; &nbsp; end do</i></div>
<div class=""><i class="">&nbsp; else</i></div>
<div class=""><i class="">; ***** Put data back in the areas that we don't want masked. (source:
<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.ncl.ucar.edu_Applications_Scripts_mask-5F9.ncl&amp;d=DwMFaQ&amp;c=qKdtBuuu6dQK9MsRUVJ2DPXW6oayO8fu4TfEHS8sGNk&amp;r=trcSnQNnUT2ddBXcpd358A&amp;m=yJKWezSUje6we6fxpkzLjE7jtiGQ1kGaXuaMupFyefY&amp;s=ylM244pXeKWXrgrLrfc52xrVmOM31ChELZkdcBRprSw&amp;e=" target="_blank" class="">
https://www.ncl.ucar.edu/<wbr class="">Applications/Scripts/mask_9.<wbr class="">ncl</a>)</i></div>
<div class=""><i class="">&nbsp; &nbsp; do ilt=ilt1,ilt2</i></div>
<div class=""><i class="">&nbsp; &nbsp; &nbsp; do iln=iln1,iln2</i></div>
<div class=""><i class="">&nbsp; &nbsp; &nbsp; &nbsp; if(gc_inout(lat(ilt),lon(iln),<wbr class="">shp_lat,shp_lon)) then</i></div>
<div class=""><i class="">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; data_mask(ilt,iln) = GFDL_bias_2D(ilt,iln) &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; replace missing values *inside* study region with data again</i></div>
<div class=""><i class="">&nbsp; &nbsp; &nbsp; &nbsp; end if</i></div>
<div class=""><i class="">&nbsp; &nbsp; &nbsp; end do</i></div>
<div class=""><i class="">&nbsp; &nbsp; end do</i></div>
<div class=""><i class="">&nbsp; end if</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; ; ******************************<wbr class="">*****************</i></div>
<div class=""><i class="">&nbsp; ; WRITE THE AVERAGES INTO NEW netCDF FILE</i></div>
<div class=""><i class="">&nbsp; ; ******************************<wbr class="">*****************</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; print(&quot;3 - Write the averages into new netCDF file&quot;)</i></div>
<div class=""><i class="">&nbsp;&nbsp;</i></div>
<div class=""><i class="">&nbsp; &nbsp; netCDF1 = False<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; Output format is NetCDF</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; if (netCDF1) then</i></div>
<div class=""><i class="">&nbsp; &nbsp; diro = &quot;~/GFDL/TAS/DATAFILES/monthly/<wbr class="">5b-4alt_GFDL_TAS_bias/avg_by_<wbr class="">target_month/nc_shp/&quot;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; output in subdirectory &quot;subset&quot; within input directory</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; &nbsp; filename = &quot;GFDL_TAS_bias_1980-2014_<wbr class="">target_month_&quot; &#43; target_month_zero &#43; &quot;_shp&quot;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; output filename for subset using initialization year and month to differentiate</i></div>
<div class=""><i class="">&nbsp; &nbsp; filo = filename &#43;&quot;.nc&quot;</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; end if</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; setfileoption(&quot;nc&quot;,&quot;Format&quot;,&quot;<wbr class="">LargeFile&quot;)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; if (netCDF1) then</i></div>
<div class=""><i class="">&nbsp; &nbsp; system(&quot;/bin/rm -f &quot;&#43; diro &#43; filo)<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; removes old files with the same name if they are present</i></div>
<div class=""><i class="">&nbsp; &nbsp; ncdf = addfile(diro &#43; filo,&quot;c&quot;)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; &nbsp; setfileoption(ncdf,&quot;<wbr class="">DefineMode&quot;,True)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; &nbsp; ; create attributes</i></div>
<div class=""><i class="">&nbsp; &nbsp; fAtt = True</i></div>
<div class=""><i class="">&nbsp; &nbsp; fAtt@title = &quot;Data: GFDL FLOR-B01 for the Contiguous U.S., Forecast Month &quot; &#43; target_month_zero</i></div>
<div class=""><i class="">&nbsp; &nbsp; fAtt@source = &quot;Vecchi et al. 2014&quot;</i></div>
<div class=""><i class="">&nbsp; &nbsp; fAtt@Conventions = &quot;None&quot;</i></div>
<div class=""><i class="">&nbsp; &nbsp; fAtt@creation_date = systemfunc(&quot;date&quot;)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; &nbsp; fileattdef(ncdf,fAtt)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; &nbsp; dimNames = (/&quot;target_month&quot;,&quot;lat&quot;,&quot;lon&quot;/)</i></div>
<div class=""><i class="">&nbsp; &nbsp; dimSizes = (/dimsizes(target_month),<wbr class="">dimsizes(lat),dimsizes(lon)/)</i></div>
<div class=""><i class="">&nbsp; &nbsp; dimUnlim = (/True,False,False/)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; &nbsp; filedimdef(ncdf,dimNames,<wbr class="">dimSizes,dimUnlim)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; &nbsp; filevardef(ncdf,&quot;target_month&quot;<wbr class="">,typeof(target_month),(/&quot;<wbr class="">target_month&quot;/))</i></div>
<div class=""><i class="">&nbsp; &nbsp; filevardef(ncdf,&quot;lat&quot;,typeof(<wbr class="">lat),(/&quot;lat&quot;/))</i></div>
<div class=""><i class="">&nbsp; &nbsp; filevardef(ncdf,&quot;lon&quot;,typeof(<wbr class="">lon),(/&quot;lon&quot;/))</i></div>
<div class=""><i class="">&nbsp; &nbsp; filevardef(ncdf,&quot;data_mask&quot;,<wbr class="">typeof(data_mask),(/&quot;lat&quot;,&quot;<wbr class="">lon&quot;/))</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; &nbsp; ; now write all the variables to the file</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; &nbsp; ncdf-&gt;target_month = (/target_month/)</i></div>
<div class=""><i class="">&nbsp; &nbsp; ncdf-&gt;lat = (/lat/)</i></div>
<div class=""><i class="">&nbsp; &nbsp; ncdf-&gt;lon = (/lon/)</i></div>
<div class=""><i class="">&nbsp; &nbsp; ncdf-&gt;data_mask = (/data_mask/)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; end if</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; ; ******************************<wbr class="">*****************</i></div>
<div class=""><i class="">&nbsp; ; CREATE MAP</i></div>
<div class=""><i class="">&nbsp; ; ******************************<wbr class="">*****************</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; ; Define boundaries for map</i></div>
<div class=""><i class="">&nbsp; data_mask!0 = &quot;lat&quot;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; latitude information come from variable lat_subdomain</i></div>
<div class=""><i class="">&nbsp; data_mask&amp;lat = lat</i></div>
<div class=""><i class="">&nbsp; data_mask!1 = &quot;lon&quot;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; longitude information come from variable lon_subdomain</i></div>
<div class=""><i class="">&nbsp; data_mask&amp;lon = lon</i></div>
<div class=""><i class="">&nbsp; data_mask@long_name = &quot;Range: &quot; &#43; min_bias &#43; &quot; - &quot; &#43; max_bias &#43; &quot; mm/day&quot; &nbsp; &nbsp;; map subtitle, printed on the top left</i></div>
<div class=""><i class="">&nbsp; data_mask@units = &quot;Avg.: &quot; &#43; avg_bias &#43; &quot; mm/day&quot;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; map units, printed on the top right</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; diro = &quot;~/GFDL/TAS/DATAFILES/monthly/<wbr class="">5b-4alt_GFDL_TAS_bias/avg_by_<wbr class="">target_month/png_shp/&quot;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; output directory</i></div>
<div class=""><i class="">&nbsp; wks = gsn_open_wks(&quot;png&quot;, diro &#43; &quot;GFDL_TAS_bias_1980-2014_<wbr class="">target_month_&quot; &#43; target_month_zero &#43; &quot;_shp&quot;)<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; map name and file type</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; cmap = read_colormap_file(&quot;<wbr class="">BlueGreen14&quot;)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; res &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= True<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; plot mods desired</i></div>
<div class=""><i class="">&nbsp; res@tiMainString &nbsp; &nbsp; &nbsp; &nbsp; = &quot;Avg. GFDL bias, Forecast Month &quot; &#43; target_month_zero &#43; &quot;, 1980-2014, no bias correction&quot;
<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; main title</i></div>
<div class=""><i class="">&nbsp; res@cnFillOn &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = True<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; turn on color fill, works with gsn_csm_map() function 3 lines down</i></div>
<div class=""><i class="">&nbsp; res@gsnAddCyclic &nbsp; &nbsp; &nbsp; &nbsp; = False</i></div>
<div class=""><i class="">&nbsp; res@cnLinesOn &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= False<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; contour lines of the map fill on or off &nbsp;</i></div>
<div class=""><i class="">&nbsp; res@mpProjection &nbsp; &nbsp; &nbsp; &nbsp; = &quot;Mercator&quot;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; Mercator map projection</i></div>
<div class=""><i class="">&nbsp; res@cnFillPalette &nbsp; &nbsp; &nbsp; &nbsp;= cmap<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; reverse color map</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; res@mpLimitMode &nbsp; &nbsp;= &quot;LatLon&quot;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; define plotted area by lats and lons, namely those of the selection (20-50N, 230-300E)</i></div>
<div class=""><i class="">&nbsp; res@mpMaxLatF &nbsp; &nbsp; &nbsp;= 41.5</i></div>
<div class=""><i class="">&nbsp; res@mpMinLatF &nbsp; &nbsp; &nbsp;= 25</i></div>
<div class=""><i class="">&nbsp; res@mpMaxLonF &nbsp; &nbsp; &nbsp;= 267</i></div>
<div class=""><i class="">&nbsp; res@mpMinLonF &nbsp; &nbsp; &nbsp;= 250</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; res@mpOutlineBoundarySets = &quot;GeophysicalAndUSStates&quot;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; lines and boundaries of USA and US States</i></div>
<div class=""><i class="">&nbsp; res@<wbr class="">mpGeophysicalLineThicknessF = 2 &nbsp; &nbsp;
<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>&nbsp; &nbsp; &nbsp; &nbsp;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; thickness of the USA lines</i></div>
<div class=""><i class="">&nbsp; res@mpUSStateLineThicknessF &nbsp; &nbsp; = 2<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; thickness of the state lines</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; res@mpSpecifiedFillColors &nbsp;= (/0,100/) &nbsp; &nbsp; &nbsp;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>&nbsp; &nbsp; &nbsp;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; fill with background color, works with gsn_map() 2 lines down</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">; &nbsp;res@cnLevelSelectionMode = &quot;ManualLevels&quot;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; manually set the contour levels with the following 3 resources</i></div>
<div class=""><i class="">; &nbsp;res@cnMinLevelValF &nbsp;= -2 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>&nbsp; &nbsp; &nbsp;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; set the minimum contour level</i></div>
<div class=""><i class="">; &nbsp;res@cnMaxLevelValF &nbsp;= 4 &nbsp; &nbsp; &nbsp; <span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; set the maximum contour level</i></div>
<div class=""><i class="">; &nbsp;res@cnLevelSpacingF = 0.25 &nbsp; &nbsp; &nbsp; &nbsp;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; set the interval between contours</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; ; ******************************<wbr class="">*****************</i></div>
<div class=""><i class="">&nbsp; ; PLOT MAP</i></div>
<div class=""><i class="">&nbsp; ; ******************************<wbr class="">*****************</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; print (&quot;5 - Plot map&quot;)</i></div>
<div class=""><i class="">&nbsp;&nbsp;</i></div>
<div class=""><i class="">&nbsp; plot = gsn_csm_contour_map(wks,data_<wbr class="">mask, res)<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; universal map projection, define in res@mpProjection</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; print(&quot;**** ALL DONE: target_month &quot; &#43; target_month &#43; &quot; ****&quot;)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">&nbsp; delete(bias_f)<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>&nbsp;<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>&nbsp;; delete before the next loop or program will fail</i></div>
<div class=""><i class="">&nbsp; delete(GFDL_bias_3D)<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; delete before the next loop or program will fail</i></div>
<div class=""><i class="">&nbsp; delete(GFDL_bias_2D)<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; delete before the next loop or program will fail</i></div>
<div class=""><i class="">&nbsp; delete(data_mask) &nbsp; <span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; delete before the next loop or program will fail</i></div>
<div class=""><i class="">&nbsp;&nbsp;</i></div>
<div class=""><i class="">end do<span class="m_3164076619572964805x_Apple-tab-span" style="white-space:pre-wrap">
</span>; end the loop through the 12 calendar months (1-12)</i></div>
<div class=""><i class=""><br class="">
</i></div>
<div class=""><i class="">end</i></div>
</div>
<span class="">
<div class=""><br class="">
</div>
<div class="">
<div class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div class=""><b style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class="">Toni Klemm</b><br style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class="">
<span style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class="">Ph.D. Candidate</span><br style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class="">
<span style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class="">The University of Oklahoma</span><br style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class="">
<span style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class="">South Central Climate Science Center</span><br style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class="">
</div>
<div class=""><span style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class="">phone:
<a href="tel:(405)%20325-1272" value="&#43;14053251272" target="_blank" class="">405 325 1272</a></span></div>
<div class=""><span style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class=""><a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.toni-2Dklemm.de&amp;d=DwMFaQ&amp;c=qKdtBuuu6dQK9MsRUVJ2DPXW6oayO8fu4TfEHS8sGNk&amp;r=trcSnQNnUT2ddBXcpd358A&amp;m=yJKWezSUje6we6fxpkzLjE7jtiGQ1kGaXuaMupFyefY&amp;s=WQc3GB-ur2qv3crEfh486hHrTNM9_W2wW_lRnDzS6vw&amp;e=" target="_blank" class="">www.toni-klemm.de</a></span></div>
<div class=""><a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.eccforum.org&amp;d=DwMFaQ&amp;c=qKdtBuuu6dQK9MsRUVJ2DPXW6oayO8fu4TfEHS8sGNk&amp;r=trcSnQNnUT2ddBXcpd358A&amp;m=yJKWezSUje6we6fxpkzLjE7jtiGQ1kGaXuaMupFyefY&amp;s=JfnBVqnY8ot6dhKlHP7_rdIQpRnK7RTyvk5rkvu2PQc&amp;e=" target="_blank" class="">Early
 Career Climate Forum</a></div>
<div class="">Twitter:&nbsp;<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__twitter.com_ToniKlemm&amp;d=DwMFaQ&amp;c=qKdtBuuu6dQK9MsRUVJ2DPXW6oayO8fu4TfEHS8sGNk&amp;r=trcSnQNnUT2ddBXcpd358A&amp;m=yJKWezSUje6we6fxpkzLjE7jtiGQ1kGaXuaMupFyefY&amp;s=hNxdFfUOc4W-DuBFvw3cfz-9NKnq4b6JMIW-8fnrJ9A&amp;e=" target="_blank" class="">@toniklemm</a></div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</span></div>
<div style="word-wrap:break-word" class="">
<div class="">
<div id="m_3164076619572964805x_AppleMailSignature" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div class=""></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<div style="word-wrap:break-word" class="">
<div class="">
<div id="m_3164076619572964805x_AppleMailSignature" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div class=""></div>
<div class=""></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<div class=""><span id="cid:848bb57a-8d3f-4131-8b70-e31ed6753761@namprd03.prod.outlook.com">&lt;GFDL_TAS_bias_1980-2014_target_month_01_shp.png&gt;</span>
</div>
<div class="">
<div class="h5">
<div style="word-wrap:break-word" class="">
<div class="">
<div id="m_3164076619572964805x_AppleMailSignature" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div class=""></div>
<div class=""></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<div style="word-wrap:break-word" class="">
<div class="">
<div class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class="">
<div class=""></div>
<div class=""><br class="">
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br class="m_3164076619572964805x_Apple-interchange-newline">
</div>
<br class="">
<div class="">
<blockquote type="cite" class="">
<div class="">On May 24, 2017, at 12:35 PM, Rick Brownrigg &lt;<a href="mailto:brownrig@ucar.edu" target="_blank" class="">brownrig@ucar.edu</a>&gt; wrote:</div>
<br class="m_3164076619572964805x_Apple-interchange-newline">
<div class="">
<div dir="ltr" class="">
<div class="">
<div class="">
<div class="">Hi Toni,<br class="">
<br class="">
</div>
Sorry for the delay in response. There are several examples of using shapefiles to mask out data, if that's what are wanting to do. In particular, look at examples 4, 5, and 21 at:<br class="">
<br class="">
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__ncl.ucar.edu_Applications_shapefiles.shtml&amp;d=DwMFaQ&amp;c=qKdtBuuu6dQK9MsRUVJ2DPXW6oayO8fu4TfEHS8sGNk&amp;r=trcSnQNnUT2ddBXcpd358A&amp;m=_2unS1RvDS7egGpkWknJW0V9KMHhv2GWDckbCG7ZOcQ&amp;s=VqnvWQCKomt1QM_SAC2MNal5QUSSlnZmBwXdI3_KjO0&amp;e=" target="_blank" class="">http://ncl.ucar.edu/<wbr class="">Applications/shapefiles.shtml</a><br class="">
<br class="">
</div>
II think there's a utility script they make use that you'll need to download.&nbsp; Hopefully this gets you at least started -- please post back the list if you have other questions.<br class="">
<br class="">
</div>
Rick<br class="">
</div>
<div class="m_3164076619572964805x_gmail_extra"><br class="">
<div class="m_3164076619572964805x_gmail_quote">On Mon, May 22, 2017 at 10:10 PM, Klemm, Toni
<span dir="ltr" class="">&lt;<a href="mailto:toni@ou.edu" target="_blank" class="">toni@ou.edu</a>&gt;</span> wrote:<br class="">
<blockquote class="m_3164076619572964805x_gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word" class="">Good evening,
<div class=""><br class="">
</div>
<div class="">I am analyzing seasonal climate forecasts skill for the U.S. Great Plains, specifically Texas, Oklahoma, Kansas, and Colorado. So far I’ve been using lat/lon to define my study region. However, this rectangular box includes parts of Mexico and
 the Gulf of Mexico that are outside my study region and which have extreme outliers that distort my results.</div>
<div class=""><br class="">
</div>
<div class="">To get more accurate results, I’d like to use a shapefile of my 4 states (<font face="Courier" class="">TX_OK_KS_CO.shp</font>, created in QGIS 2.18 from a U.S. states shapefile) to define my study region. However, I’m not sure how to interpret
 and use the results I’m getting. I used &nbsp; “<font face="Courier" class="">ncl_convert2nc TX_OK_KS_CO.shp -nc2c</font>“ &nbsp;to create &nbsp; “<font face="Courier" class="">TX_OK_KS_CO.nc</font>“ .</div>
<div class=""><br class="">
</div>
<div class="">“<font face="Courier" class="">printVarSummary</font>(<span style="font-family:Courier" class="">TX_OK_KS_CO.n<wbr class="">c)</span>” and ”<font face="Courier" class="">PrintMinMax</font>“ of the x-coordinate (longitude) as well as y-coordinate
 (latitude) look like this:</div>
<div class=""><br class="">
</div>
<div class="">
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class="">Variable: shp</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class="">Type: file</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class="">File path:<span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap">
</span>/home/States_shapefile_data/GA<wbr class="">DM/TX_OK_KS_CO.nc</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class="">Number of global attributes:<span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap">
</span>11</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class="">Number of dimensions:<span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap">
</span>5</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class="">Number of variables:<span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap">
</span>4</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class="">(0)<span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap">
</span>min=0 &nbsp; max=0</span></div>
</div>
<div class=""><br class="">
</div>
<div class=""><font face="Courier" class="">ncdump -h TX_OK_KS_CO.nc</font>&nbsp;returns this:</div>
<div class=""><br class="">
</div>
<div class="">
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class="">netcdf TX_OK_KS_CO {</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class="">dimensions:</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>geometry = 2 ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>segments = 2 ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>num_features = 4 ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>num_segments = 262 ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>num_points = 102749 ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class="">variables:</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>int geometry(num_features, geometry) ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>int segments(num_segments, segments) ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>double x(num_points) ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>double y(num_points) ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201);min-height:16px" class="">
<span class=""></span><br class="">
</div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class="">// global attributes:</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>:segs_numPnts = 1 ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>:segs_xyzIndex = 0 ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>:geom_numSegs = 1 ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>:geom_segIndex = 0 ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>:geometry_type = &quot;polygon&quot; ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>:layer_name = &quot;TX_OK_KS_CO&quot; ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>:creation_date = &quot;Mon May 22 17:20:46 EDT 2017&quot; ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>:NCL_Version = &quot;6.4.0&quot; ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>:system = &quot;Linux
<a href="http://ou.edu/" target="_blank" class="">ou.edu</a> x86_64 #1 SMP Fri Mar 3 00:04:05 UTC 2017 x86_64 x86_64 x86_64 GNU/Linux&quot; ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>:Conventions = &quot;None&quot; ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class=""><span class="m_3164076619572964805x_m_-6640439190133930147Apple-tab-span" style="white-space:pre-wrap"></span>:title = &quot;NCL: convert-OGR-to-netCDF&quot; ;</span></div>
<div style="margin:0px;font-size:13px;line-height:normal;font-family:Courier;color:rgb(255,255,255);background-color:rgb(43,102,201)" class="">
<span class="">}</span></div>
</div>
<div class=""><br class="">
</div>
<div class="">I am not sure how to use this information, or if this is even the information I should be getting, to define the study region when extracting data from my global forecast dataset. I am fairly sure the shapefile is correct. Opening it again in
 QGIS showed the four states with their outlines and filling, and nothing else.</div>
<div class=""><br class="">
</div>
<div class="">I appreciate any help!</div>
<div class=""><br class="">
</div>
<div class="">Thank you very much,</div>
<div class="">Toni</div>
<span class="m_3164076619572964805x_HOEnZb"><font color="#888888" class="">
<div class=""><br class="">
</div>
<div class=""><br class="">
<div class="">
<div class="">
<div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word" class="">
<div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word" class="">
<div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word" class="">
<div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word" class="">
<div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word" class="">
<div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word" class="">
<div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word" class="">
<div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word" class="">
<div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word" class="">
<div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word" class="">
<div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word" class="">
<div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word" class="">
<div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word" class="">
<div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word" class="">
<div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word" class="">
<div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word" class="">
<div class=""><b style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class="">Toni Klemm</b><br style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class="">
<span style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class="">Ph.D. Candidate</span><br style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class="">
<span style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class="">The University of Oklahoma</span><br style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class="">
<span style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class="">South Central Climate Science Center</span><br style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class="">
</div>
<div class=""><span style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class="">phone:
<a href="tel:(405)%20325-1272" value="&#43;14053251272" target="_blank" class="">405 325 1272</a></span></div>
<div class=""><span style="color:rgb(51,51,153);background-color:rgb(255,255,255)" class=""><a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.toni-2Dklemm.de&amp;d=DwMFaQ&amp;c=qKdtBuuu6dQK9MsRUVJ2DPXW6oayO8fu4TfEHS8sGNk&amp;r=trcSnQNnUT2ddBXcpd358A&amp;m=_2unS1RvDS7egGpkWknJW0V9KMHhv2GWDckbCG7ZOcQ&amp;s=QAW2_xJ9gp71I8wBvMBw6JI7Kb7j7d8KF0tlyNuhEYE&amp;e=" target="_blank" class="">www.toni-klemm.de</a></span></div>
<div class=""><a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.eccforum.org&amp;d=DwMFaQ&amp;c=qKdtBuuu6dQK9MsRUVJ2DPXW6oayO8fu4TfEHS8sGNk&amp;r=trcSnQNnUT2ddBXcpd358A&amp;m=_2unS1RvDS7egGpkWknJW0V9KMHhv2GWDckbCG7ZOcQ&amp;s=BsFxGtcxiOomgvIRjvEfgPA7YppuvrO0n-TW0SLSiok&amp;e=" target="_blank" class="">Early
 Career Climate Forum</a></div>
<div class="">Twitter:&nbsp;<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__twitter.com_ToniKlemm&amp;d=DwMFaQ&amp;c=qKdtBuuu6dQK9MsRUVJ2DPXW6oayO8fu4TfEHS8sGNk&amp;r=trcSnQNnUT2ddBXcpd358A&amp;m=_2unS1RvDS7egGpkWknJW0V9KMHhv2GWDckbCG7ZOcQ&amp;s=B5YLwKe0WRmXO4D2nzE_HgiPLcFnIFvPX6Z-kf3KSrw&amp;e=" target="_blank" class="">@toniklemm</a></div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br class="m_3164076619572964805x_m_-6640439190133930147Apple-interchange-newline">
</div>
<br class="">
</div>
</div>
</font></span></div>
<br class="">
______________________________<wbr class="">_________________<br class="">
ncl-talk mailing list<br class="">
<a href="mailto:ncl-talk@ucar.edu" target="_blank" class="">ncl-talk@ucar.edu</a><br class="">
List instructions, subscriber options, unsubscribe:<br class="">
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.ucar.edu_mailman_listinfo_ncl-2Dtalk&amp;d=DwMFaQ&amp;c=qKdtBuuu6dQK9MsRUVJ2DPXW6oayO8fu4TfEHS8sGNk&amp;r=trcSnQNnUT2ddBXcpd358A&amp;m=_2unS1RvDS7egGpkWknJW0V9KMHhv2GWDckbCG7ZOcQ&amp;s=g_J7lY2OPhdO6e1uZkyCAbj2G1XmDbXoRk01EanWTdI&amp;e=" rel="noreferrer" target="_blank" class="">http://mailman.ucar.edu/mailma<wbr class="">n/listinfo/ncl-talk</a><br class="">
<br class="">
</blockquote>
</div>
<br class="">
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
<span id="cid:9AB0C65A-79D4-4769-9115-1C5C9384234B@ou.edu">&lt;texasCoast.png&gt;</span><span id="cid:CC4E8DD6-3F00-4473-8DCF-18F6D52663ED@ou.edu">&lt;GFDL.png.png&gt;</span></div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>