<div dir="ltr"><div class="gmail_default" style="font-size:small">Lin,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I'm confused by the comment on this line:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><div class="gmail_default"> Opt@SrcGridMask = where(.not.ismissing(thetao(0,:,:)),1,0) ; thetao is the source grids, with data zero in land, no missing </div><div><br></div><div>When you say "data zero in land, no missing", do you mean that the "thetao" variable is equal to 0 whenever it is over land, and that it contains no missing values? If so, then using "ismissing" on a variable that doesn't contain missing values is never going to return True.</div><div><br></div><div>Instead, you want something like this:</div><div><br></div><div><div style="font-size:12.8px"> Opt@SrcGridMask = where(thetao.ne.0,1,0)</div></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">This creates a mask array that is equal to 1 at locations where thetao is NOT equal to 0, and equal to 0 at locations where thetao is equal to 0. </div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">Just as an example, this line should produce the exact same mask array:</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><div style="font-size:small"><div style="font-size:12.8px"> Opt@SrcGridMask = where(thetao.eq.0,0,1)</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">--Mary</div><div style="font-size:12.8px"><br></div><div><br></div></div></div><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Dec 19, 2016 at 11:53 AM, 林祥 <span dir="ltr"><<a href="mailto:xianglin72@icloud.com" target="_blank">xianglin72@icloud.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Hi, all<div><br></div><div> I have many CMIP5 sea ice concentration data and want to regrid them into the same grid coordinate as HadISST.</div><div><br></div><div> Since the land area own missing_value (_FillValue) in HadiSST, I also want the regridded data show missing values in land .</div><div><br></div><div>I use the following scripts:</div><div> </div><div> <span style="font-family:Menlo;font-size:14px;font-variant-ligatures:no-common-ligatures"> fi = </span><span style="font-family:Menlo;font-size:14px;font-variant-ligatures:no-common-ligatures;color:rgb(52,187,199)">addfile</span><span style="font-family:Menlo;font-size:14px;font-variant-ligatures:no-common-ligatures">(</span><span style="color:rgb(195,55,32);font-family:Menlo;font-size:14px;font-variant-ligatures:no-common-ligatures">"HadISST_ice.nc"</span><span style="font-family:Menlo;font-size:14px;font-variant-ligatures:no-common-ligatures">,</span><span style="color:rgb(195,55,32);font-family:Menlo;font-size:14px;font-variant-ligatures:no-common-ligatures">"r"</span><span style="font-family:Menlo;font-size:14px;font-variant-ligatures:no-common-ligatures">)</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures"> hsic = fi->sic</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures"> hlon = fi->longitude</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures"> hlat = fi->latitude(</span><span style="font-variant-ligatures:no-common-ligatures;color:#ce7924">{</span><span style="font-variant-ligatures:no-common-ligatures">-</span><span style="font-variant-ligatures:no-common-ligatures;color:#c33720">90</span><span style="font-variant-ligatures:no-common-ligatures">.:-</span><span style="font-variant-ligatures:no-common-ligatures;color:#c33720">50</span><span style="font-variant-ligatures:no-common-ligatures">.</span><span style="font-variant-ligatures:no-common-ligatures;color:#ce7924">}</span><span style="font-variant-ligatures:no-common-ligatures">)</span></div><div> …….</div><div> </div><div> Opt@SrcGridMask = where(.not.ismissing(thetao(0,<wbr>:,:)),1,0) ; thetao is the source grids, with data zero in land, no missing </div><div> Opt@DstGridMask = where(.not.ismissing(hsic(0,:,<wbr>:)),1,0) ; hsic is the HadISST SIC data, with missing values in land</div><div> </div><div> <span style="font-family:Menlo;font-size:14px">Opt@DstGridLat </span><span style="font-family:Menlo;font-size:14px"> </span><span style="font-family:Menlo;font-size:14px">= hlat</span></div><div><span style="font-family:Menlo;font-size:14px"> </span><span style="font-family:Menlo;font-size:14px">Opt@DstGridLon </span><span style="font-family:Menlo;font-size:14px"> </span><span style="font-family:Menlo;font-size:14px">= hlon</span> </div><div>……..</div><div><br></div><div> <span style="font-family:Menlo;font-size:14px"> </span><span style="font-family:Menlo;font-size:14px">thetao_regrid = ESMF_regrid(thetao,Opt)</span></div><div><br></div><div>However, I found the regridded data “thetao_regrid" show zero values, not missing values in land.</div><div><br></div><div><br></div><div> How to make the regridded data have missing values for the same grids as HadISST?</div><div><br></div><div><br></div><div>Thank</div><div><br></div><div><br></div><div>Lin </div><div><br></div><div> </div><div><br></div></div><br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>