<div dir="ltr">Hi Jonathan,<div>I think the issue is related to your last set of where calls:</div><div><div>cntrycodes_regionALL = where(cntrycodes_region01.eq.1,1,cntrycodes_regionALL)</div><div>cntrycodes_regionALL = where(cntrycodes_region02.eq.2,2,cntrycodes_regionALL)</div><div>cntrycodes_regionALL = where(cntrycodes_region03.eq.3,3,cntrycodes_regionALL)</div></div><div><br></div><div>Running small tests here, each call is likely overwriting the cntrycodes_regionALL array on each line. This is one of those situations where the coding looks right (to me at least), but where is behaving differently than I would expect. I think this has to do with how the where function interprets areas set to _FillValue. Try changing the above lines to this:</div><div><div>cntrycodes_regionALL = where(ismissing(cntrycodes_regionALL).and.cntrycodes_region01.eq.1,1,cntrycodes_regionALL)</div><div><div>cntrycodes_regionALL = where(ismissing(cntrycodes_regionALL).and.cntrycodes_region02.eq.2,2,cntrycodes_regionALL)</div></div><div><div>cntrycodes_regionALL = where(ismissing(cntrycodes_regionALL).and.cntrycodes_region03.eq.3,3,cntrycodes_regionALL)</div><div></div></div><div>(and make sure cntrcodes_regionALL has a _FillValue attribute set)..</div><div></div></div><div><br></div><div>Let ncl-talk know if that doesn't work.</div><div>Adam</div><div><br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Apr 16, 2021 at 9:05 AM Buzan, Jonathan via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">



<div>
Hi NCL-Talk,
<div><br>
</div>
<div>I have something conceptually wrong in my head with NCL.</div>
<div><br>
</div>
<div>I am attempting to take a .nc file with countries, and map them onto regions.</div>
<div><br>
</div>
<div>Country Map</div>
<div><br>
</div>
<div><img id="gmail-m_7432066726356615982gmail-m_1175420333317078302gmail-m_8195473246610769622gmail-m_8419946374867213108BC65C0BF-416D-453E-A515-B376427488BA" src="cid:178db6c25bc55e4f0bf1"></div>
<div><br>
</div>
<div><span style="white-space:pre-wrap"></span></div>
<div><br>
</div>
<div>I use:</div>
<div><br>
</div>
<div>
<div>cntrycodes!0 = "lat"</div>
<div>cntrycodes&lat = lat</div>
<div>cntrycodes!1 = "lon"</div>
<div>cntrycodes&lon = lon</div>
<div>cntrycodes@long_name = "Country Codes"</div>
<div>cntrycodes@missing_value = cntrycodes@_FillValue</div>
</div>
<div>…</div>
<div><br>
</div>
<div><br>
</div>
<div>
<div>cntrycodes_region04 = where(cntrycodes.eq.7,4,cntrycodes_region04)</div>
<div>cntrycodes_region04 = where(cntrycodes.eq.10,4,cntrycodes_region04)</div>
<div>cntrycodes_region04 = where(cntrycodes.eq.13,4,cntrycodes_region04)</div>
</div>
<div>…</div>
<div>
<div>cntrycodes_region04 = where(cntrycodes.eq.145,4,cntrycodes_region04)</div>
<div>cntrycodes_region04 = where(cntrycodes.eq.152,4,cntrycodes_region04)</div>
</div>
<div><br>
</div>
<div><br>
</div>
<div>I get:</div>
<div><img id="gmail-m_7432066726356615982gmail-m_1175420333317078302gmail-m_8195473246610769622gmail-m_841994637486721310886104144-C5F6-4BF1-9652-6FB7DE23FA2C" src="cid:178db6c25bcba22cfa12"></div>
<div><br>
</div>
<div>I do this for 16 regions.</div>
<div>Then I try to put it together:</div>
<div><br>
</div>
<div>
<div>cntrycodes_regionALL = where(cntrycodes_region01.eq.1,1,cntrycodes_regionALL)</div>
<div>cntrycodes_regionALL = where(cntrycodes_region02.eq.2,2,cntrycodes_regionALL)</div>
<div>cntrycodes_regionALL = where(cntrycodes_region03.eq.3,3,cntrycodes_regionALL)</div>
</div>
<div>…</div>
<div>
<div>cntrycodes_regionALL = where(cntrycodes_region15.eq.15,15,cntrycodes_regionALL)</div>
<div>cntrycodes_regionALL = where(cntrycodes_region16.eq.16,16,cntrycodes_regionALL)</div>
</div>
<div><br>
</div>
<div>And I get:</div>
<div><img id="gmail-m_7432066726356615982gmail-m_1175420333317078302gmail-m_8195473246610769622gmail-m_84199463748672131081624F096-F82F-4B18-A891-CE20D73D7192" src="cid:178db6c25bc23936be53"></div>
<div><br>
</div>
<div><br>
</div>
<div>I am not sure what I am doing wrong. But I have a feeling it is with how I am using the where function.</div>
<div><br>
</div>
<div>Cheers,</div>
<div>-Jonathan</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div></div>
</div>

_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist,  </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphilli/</a>   </font></span><span><font color="#888888">303-497-1726 </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>