<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=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<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:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
pre
{mso-style-priority:99;
mso-style-link:"HTML Preformatted Char";
margin:0in;
margin-bottom:.0001pt;
font-size:10.0pt;
font-family:"Courier New";}
span.HTMLPreformattedChar
{mso-style-name:"HTML Preformatted Char";
mso-style-priority:99;
mso-style-link:"HTML Preformatted";
font-family:Consolas;}
span.EmailStyle21
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
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-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal">Hi Dennis,<o:p></o:p></p>
<p class="MsoNormal">Thanks for your response. Nice to hear from you and hope you are doing well and safe from recent wildfire in north boulder.
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I’m not sure that I understand the bin_sum function specially the second part (‘user must perform averaging’).<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">For my case, dlon, dlat and d are from csv file and lat/lon from netcdf file. See blow their dimension. I would like create new variable at lat/lon grid by summing all d based on number of dlon, dlat within each of lat/lon grid.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">For example – if we have 50 dlat/dlon in one latXlon grid, then new variable in lat/lon grid = sum of d from all 50 dlat/dlon.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">With regards<o:p></o:p></p>
<p class="MsoNormal">-Dave<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p></o:p></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Variable: dlat<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Type: float<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Total Size: 1392092 bytes<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"> 348023 values<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Number of Dimensions: 1<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Dimensions and sizes: [348023]<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Coordinates:<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Number Of Attributes: 1<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"> _FillValue : 9.96921e+36<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Variable: dlon<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Type: float<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Total Size: 1392092 bytes<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"> 348023 values<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Number of Dimensions: 1<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Dimensions and sizes: [348023]<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Coordinates:<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Number Of Attributes: 1<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"> _FillValue : 9.96921e+36<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Variable: d<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Type: float<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Total Size: 1392092 bytes<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"> 348023 values<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Number of Dimensions: 1<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Dimensions and sizes: [348023]<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Coordinates:<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Number Of Attributes: 1<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"> _FillValue : 9.96921e+36<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><b>Output grid –<o:p></o:p></b></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Variable: lat<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Type: float<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Total Size: 744 bytes<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"> 186 values<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Number of Dimensions: 1<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Dimensions and sizes: [lat | 186]<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Coordinates:<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"> lat: [-46.375..-0.125]<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Number Of Attributes: 4<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"> units : degree_north<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"> standard_name : latitude<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"> long_name : latitude<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"> bounds : lat_bnds<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Variable: lon<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Type: float<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Total Size: 964 bytes<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"> 241 values<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Number of Dimensions: 1<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Dimensions and sizes: [lon | 241]<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Coordinates:<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"> lon: [100.125..160.125]<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%">Number Of Attributes: 4<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"> units : degree_east<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"> standard_name : longitude<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"> long_name : longitude<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#203864;mso-style-textfill-fill-color:#203864;mso-style-textfill-fill-alpha:100.0%"> bounds : lon_bnds<o:p></o:p></span></p>
<p class="MsoNormal"><b>From:</b> Dennis Shea <shea@ucar.edu> <br>
<b>Sent:</b> Sunday, October 18, 2020 2:25 PM<br>
<b>To:</b> Paimazumder, Debasish <Debasish.Paimazumder@aig.com><br>
<b>Cc:</b> Barry Lynn <barry.h.lynn@gmail.com>; ncl-talk <ncl-talk@ucar.edu><br>
<b>Subject:</b> [EXTERNAL] Re: [ncl-talk] accumulate data on grid<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><b><span style="color:gray">This message is from an external sender; be cautious with links and attachments.</span></b><o:p></o:p></p>
</div>
<div>
<div>
<p class="MsoNormal">Maybe the <b>"bin" <a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.ncl.ucar.edu_Applications_binning.shtml&d=DwMFaQ&c=kn4_INW_mBCDHV_xJEVJkg&r=fpudUXThXLQwAKQPDCnmVaSzm2dwE9OxcWi7MuJ0jbo&m=6F7M48QNHan0cjDH_KWLNBKvEM8Z4MLy0AVrfi5jQo8&s=esBoEZ_hYSHMOH9ETaF6M17fqzbclWD-fQGCOTkzhL0&e=">
examples</a> /</b> <b>functions</b>: <b><a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.ncl.ucar.edu_Document_Functions_Built-2Din_bin-5Fsum.shtml&d=DwMFaQ&c=kn4_INW_mBCDHV_xJEVJkg&r=fpudUXThXLQwAKQPDCnmVaSzm2dwE9OxcWi7MuJ0jbo&m=6F7M48QNHan0cjDH_KWLNBKvEM8Z4MLy0AVrfi5jQo8&s=xnFvaC2e9cElHG3_pafK084oquAu_2dJsejynvv5C-k&e=">
bin_sum</a></b> ,<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.ncl.ucar.edu_Document_Functions_Built-2Din_bin-5Favg.shtml&d=DwMFaQ&c=kn4_INW_mBCDHV_xJEVJkg&r=fpudUXThXLQwAKQPDCnmVaSzm2dwE9OxcWi7MuJ0jbo&m=6F7M48QNHan0cjDH_KWLNBKvEM8Z4MLy0AVrfi5jQo8&s=9NeB2QchXaY0pnTSSuyn5Qmj-7oUsC6EfwD4EXGQNYA&e="><b>
bin_avg</b></a><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Assuming the data are randomly located: d(npts), dlat(npts), dlon(npts)
<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><b>**AND**</b> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">the target grid is <b>RECTILINEAR</b> with grid locations/sizes: glon[mlon], glat[nlat]<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<pre> gbin = <a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.ncl.ucar.edu_Document_Functions_Built-2Din_new.shtml&d=DwMFaQ&c=kn4_INW_mBCDHV_xJEVJkg&r=fpudUXThXLQwAKQPDCnmVaSzm2dwE9OxcWi7MuJ0jbo&m=6F7M48QNHan0cjDH_KWLNBKvEM8Z4MLy0AVrfi5jQo8&s=R3wqMSXgnUywYdf9AP9u9NrETGphvJvVf7GItjprpxI&e="><strong><span style="font-family:"Courier New"">new</span></strong></a> ( (/nlat,mlon/), float ) <o:p></o:p></pre>
<pre> gknt = <a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.ncl.ucar.edu_Document_Functions_Built-2Din_new.shtml&d=DwMFaQ&c=kn4_INW_mBCDHV_xJEVJkg&r=fpudUXThXLQwAKQPDCnmVaSzm2dwE9OxcWi7MuJ0jbo&m=6F7M48QNHan0cjDH_KWLNBKvEM8Z4MLy0AVrfi5jQo8&s=R3wqMSXgnUywYdf9AP9u9NrETGphvJvVf7GItjprpxI&e="><strong><span style="font-family:"Courier New"">new</span></strong></a> ( (/nlat,mlon/), integer) <o:p></o:p></pre>
<pre><o:p> </o:p></pre>
<pre> gbin = 0.0 ; initialization<o:p></o:p></pre>
<pre> gknt = 0 <o:p></o:p></pre>
</div>
<div>
<pre> <strong><span style="font-family:"Courier New"">bin_sum</span></strong>(gbin,gknt,glon,glat,dlon,dlat,d)<br><br>;*****************************************************************<o:p></o:p></pre>
<pre>; User must perform averaging<o:p></o:p></pre>
<pre>;*****************************************************************<o:p></o:p></pre>
<pre> ; avoid division by 0<o:p></o:p></pre>
<pre> gknt = <a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.ncl.ucar.edu_Document_Functions_Built-2Din_where.shtml&d=DwMFaQ&c=kn4_INW_mBCDHV_xJEVJkg&r=fpudUXThXLQwAKQPDCnmVaSzm2dwE9OxcWi7MuJ0jbo&m=6F7M48QNHan0cjDH_KWLNBKvEM8Z4MLy0AVrfi5jQo8&s=4SbR5DjHbZc8obQVy_1-w-LH_brR9Bnv0Ixv8K46mA0&e="><strong><span style="font-family:"Courier New"">where</span></strong></a>(gknt.eq.0 , gknt@_FillValue, gknt)<o:p></o:p></pre>
<pre><o:p> </o:p></pre>
<pre> gbin = gbin/gknt<o:p></o:p></pre>
<pre> <o:p></o:p></pre>
<pre> gbin!0 = "lat"<o:p></o:p></pre>
<pre> gbin!1 = "lon"<o:p></o:p></pre>
<pre> gbin&lat = glat<o:p></o:p></pre>
<pre> gbin&lon = lon<o:p></o:p></pre>
<pre><o:p> </o:p></pre>
<pre> <a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.ncl.ucar.edu_Document_Functions_Contributed_copy-5FVarCoords.shtml&d=DwMFaQ&c=kn4_INW_mBCDHV_xJEVJkg&r=fpudUXThXLQwAKQPDCnmVaSzm2dwE9OxcWi7MuJ0jbo&m=6F7M48QNHan0cjDH_KWLNBKvEM8Z4MLy0AVrfi5jQo8&s=H8H1S0QA3naQXQKhJSuJAT2QNQyKzwCTxbukEdlFtI8&e="><strong><span style="font-family:"Courier New"">copy_VarCoords</span></strong></a>(gbin, gknt) ; copy coords<o:p></o:p></pre>
<pre><o:p> </o:p></pre>
<pre> if (<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.ncl.ucar.edu_Document_Functions_Built-2Din_isfilevaratt.shtml&d=DwMFaQ&c=kn4_INW_mBCDHV_xJEVJkg&r=fpudUXThXLQwAKQPDCnmVaSzm2dwE9OxcWi7MuJ0jbo&m=6F7M48QNHan0cjDH_KWLNBKvEM8Z4MLy0AVrfi5jQo8&s=pbnNxst9G_K0_R0bUV0jHQfBbWTAjBmkwShzIGuJjwE&e="><strong><span style="font-family:"Courier New"">isfilevaratt</span></strong></a>(f, vNam, "long_name")) then<o:p></o:p></pre>
<pre> gbin@long_name = "BINNED: "+vNam<o:p></o:p></pre>
<pre> gknt@long_name = "BINNED COUNT: "+vNam<o:p></o:p></pre>
<pre> end if<o:p></o:p></pre>
<pre><o:p> </o:p></pre>
<pre> if (<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.ncl.ucar.edu_Document_Functions_Built-2Din_isfilevaratt.shtml&d=DwMFaQ&c=kn4_INW_mBCDHV_xJEVJkg&r=fpudUXThXLQwAKQPDCnmVaSzm2dwE9OxcWi7MuJ0jbo&m=6F7M48QNHan0cjDH_KWLNBKvEM8Z4MLy0AVrfi5jQo8&s=pbnNxst9G_K0_R0bUV0jHQfBbWTAjBmkwShzIGuJjwE&e="><strong><span style="font-family:"Courier New"">isfilevaratt</span></strong></a>(f, vNam, "units")) then<o:p></o:p></pre>
<pre> gbin@units = f->$vNam$@units<o:p></o:p></pre>
<pre> end if<o:p></o:p></pre>
<pre style="margin-bottom:12.0pt"><o:p> </o:p></pre>
<pre> printVarSummary(gbin)<br> printMinMax(gbin,1)<br> print("=========")<br> print("=========")<br> printVarSummary(gknt) <br> printMinMax(gbin,1)<br><br>====<o:p></o:p></pre>
<pre>Note: If the target grid is <b>CURVILINEAR</b>: [ glat(nlat,mlon), glon(nlat,mlon], there is no NCL function to accomplish the task.<o:p></o:p></pre>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">On Sun, Oct 18, 2020 at 11:32 AM Paimazumder, Debasish via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.ucar.edu</a>> wrote:<o:p></o:p></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in">
<div>
<p class="MsoNormal">Dave- <span style="color:black">getind_latlon2d needs two-dimensional lat/lon while my output grid 1-d lat/lon. </span>
<o:p></o:p></p>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="color:black">Interpolation from location to 1-d grid could be an option but we won’t get the accumulation (sum) at the output grid.</span><o:p></o:p></p>
<div>
<p class="MsoNormal">Best regards, <o:p></o:p></p>
<div>
<p class="MsoNormal">Dave<o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"><br>
<br>
<o:p></o:p></p>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<p class="MsoNormal" style="margin-bottom:12.0pt">On Oct 18, 2020, at 12:54 PM, Barry Lynn <<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>> wrote:<o:p></o:p></p>
</blockquote>
</div>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<p class="MsoNormal"> <o:p></o:p></p>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><b><span style="color:gray">This message is from an external sender; be cautious with links and attachments.</span></b><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Right, that can be done. Or, you can interpolate the data from the fine to coarse grid.
<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">I am not sure if the latter will conserve the averages like doing the averaging in a loop.<o:p></o:p></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">On Sun, Oct 18, 2020 at 6:58 PM Dave Allured - NOAA Affiliate via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<o:p></o:p></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in">
<div>
<p class="MsoNormal">Use the function <b>getind_latlon2d</b> to get the (i,j) grid point indices for the point data coordinates. Then do a single loop over the point data to accumulate sums and counts on the output grid.
<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">On Sun, Oct 18, 2020 at 9:09 AM Paimazumder, Debasish via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<o:p></o:p></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in">
<div>
<p class="MsoNormal">Hi Barry, <o:p></o:p></p>
<div>
<p class="MsoNormal">Thanks for your response. Actually I don’t have time dimension. I have location level data (at high resolution).<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">For example, x(m) with lat(m) and lon (m) where m = 300k location <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">My output grid is 32x32 km and I would like to do following <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">1. Estimate the number of locations are in each 32x32 km grid<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">2. Accumulate (sum) the x at each output grid.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">For example- if I have 100 locations at one of the 32x32km grid, I will sum all x from 100 locations and assign the sum to the output grid.<o:p></o:p></p>
<div>
<p class="MsoNormal">Best regards, <o:p></o:p></p>
<div>
<p class="MsoNormal">Dave<o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"><br>
<br>
<o:p></o:p></p>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<p class="MsoNormal">On Oct 18, 2020, at 3:50 AM, Barry Lynn <<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>> wrote:<o:p></o:p></p>
</blockquote>
</div>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<div>
<p class="MsoNormal">Hi: <o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">If you have two grids of the same dimensions, but different times, you can simply add them together.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">If you have grids of a different size, you can interpolate the data to this grid, one grid at a time.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">See for example:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.ncl.ucar.edu_Applications_regrid.shtml&d=DwMFaQ&c=kn4_INW_mBCDHV_xJEVJkg&r=fpudUXThXLQwAKQPDCnmVaSzm2dwE9OxcWi7MuJ0jbo&m=IaFm8AuXStzIPZXMbIRmp6fjt5zy74eAUUQzoQb5pBc&s=tUIMDu-yehcnE2-mHTaX2NhZiAg5eBHJsqf3e1Bmbl8&e=" target="_blank">https://www.ncl.ucar.edu/Applications/regrid.shtml</a><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">This may also help.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.ncl.ucar.edu_Document_Functions_Built-2Din_dim-5Fsum-5Fn.shtml&d=DwMFaQ&c=kn4_INW_mBCDHV_xJEVJkg&r=fpudUXThXLQwAKQPDCnmVaSzm2dwE9OxcWi7MuJ0jbo&m=IaFm8AuXStzIPZXMbIRmp6fjt5zy74eAUUQzoQb5pBc&s=4tdKFBxFBOUpMYFj95WyO8MUdUcYRniTWsJIvCCbhAE&e=" target="_blank">https://www.ncl.ucar.edu/Document/Functions/Built-in/dim_sum_n.shtml</a><o:p></o:p></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">On Sun, Oct 18, 2020 at 7:59 AM Paimazumder, Debasish via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<o:p></o:p></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in">
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Hi NCL users<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">I have location level data and 2d lat, lon grid and I’m trying to accumulate the data on 2d lat, lon grid.<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">For example
<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Input - x (n) with lat(n), lon(n)<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Output grid - lat(i)xlon(j)<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">I’m trying to estimate number of location within each grid and accumulate x at each grid.
<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Is there any function in NCL that will solve this problem?<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">With regards<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">-Debasish
<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><b><span style="font-size:12.0pt;color:#00B0F0">Debasish PaiMazumder, PhD |</span></b><span style="font-size:12.0pt;color:#666666">Tel
<a href="tel:%2B1%20215%20255%206371" target="_blank">+1 215 255 6052</a> | Cell <a href="tel:%2B1%20718%20913%206946" target="_blank">
+1 267 449 5226</a> |</span><u><span style="font-size:12.0pt;color:blue"><a href="mailto:Debasish.Paimazumder@aig.com" target="_blank">Debasish.Paimazumder@aig.com</a></span></u><o:p></o:p></p>
</div>
</div>
</blockquote>
</div>
</div>
</blockquote>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
<p class="MsoNormal">_______________________________________________<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://urldefense.proofpoint.com/v2/url?u=https-3A__mailman.ucar.edu_mailman_listinfo_ncl-2Dtalk&d=DwMFaQ&c=kn4_INW_mBCDHV_xJEVJkg&r=fpudUXThXLQwAKQPDCnmVaSzm2dwE9OxcWi7MuJ0jbo&m=MrHhJNEEUYaIu_m4tHYe3BagQB4y_F9zUKivzqyU6Qc&s=OAPVqMfiOsSTA-yHW1I4powc9uQpzwob8irpoxiBBy0&e=" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><o:p></o:p></p>
</blockquote>
</div>
<p class="MsoNormal"><br clear="all">
<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<p class="MsoNormal">-- <o:p></o:p></p>
<div>
<div>
<div>
<div>
<p class="MsoNormal">Barry H. Lynn, Ph.D <o:p></o:p></p>
<div>
<div>
<p class="MsoNormal">Senior Associate Scientist, Lecturer,<o:p></o:p></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="color:#888888">The Institute of the Earth Science, <br>
The Hebrew University of Jerusalem, <br>
Givat Ram, Jerusalem 91904, Israel </span><o:p></o:p></p>
</div>
<p class="MsoNormal"><span style="color:#888888">Tel: 972 547 231 170<br>
Fax: (972)-25662581</span><o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">C.E.O, Weather It Is, LTD<br>
Weather and Climate Focus<br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__weather-2Dit-2Dis.com&d=DwMFaQ&c=kn4_INW_mBCDHV_xJEVJkg&r=fpudUXThXLQwAKQPDCnmVaSzm2dwE9OxcWi7MuJ0jbo&m=MrHhJNEEUYaIu_m4tHYe3BagQB4y_F9zUKivzqyU6Qc&s=wTDx59WcqkIpQWlxXtecOdqmP_2sM2GEo2mXT7ZAaPw&e=" target="_blank">http://weather-it-is.com</a><br>
Jerusalem, Israel<br>
Local: 02 930 9525<br>
Cell: 054 7 231 170<br>
Int-IS: x972 2 930 9525<o:p></o:p></p>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
<p class="MsoNormal">_______________________________________________<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://urldefense.proofpoint.com/v2/url?u=https-3A__mailman.ucar.edu_mailman_listinfo_ncl-2Dtalk&d=DwMFaQ&c=kn4_INW_mBCDHV_xJEVJkg&r=fpudUXThXLQwAKQPDCnmVaSzm2dwE9OxcWi7MuJ0jbo&m=6F7M48QNHan0cjDH_KWLNBKvEM8Z4MLy0AVrfi5jQo8&s=JHsKhKBBboyoxNbA9b8hgUBDyt3F_iBYSmaDbFnDfy0&e=" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><o:p></o:p></p>
</blockquote>
</div>
</div>
</body>
</html>