<div dir="ltr"><div><div><div>The data values on '<b>SMAP_L2_</b>SM_P_E_*h5' are from satellite <b>swaths</b>. "<b>L2</b>"<br></div><div><br></div>The data variable ("soil_moisture") and the location variables: "latitude" and "longitude" are all <b>one-dimensional</b>. It looks like the swath values and the associated latitude & longitude locations were just made into one-dimensional arrays. <br><br></div><div>Your code comments the following:<br></div><div><br><b>/;</b><br> --37.80-----<br> | |<br> 45.8 46.03<br> | |<br> ---37.45----<br><b>;/</b><br><br></div><div>There is <b>no need </b>to perform and areal average when using such a small area [specifically such a small latitudinal range]. .The following should suffice:<br><br></div><div> latS = 37.45<br></div><div> latN = 37.80<br></div><div> lonL = 45.8<br></div><div> lonR = 4603<br><br></div><div> sm = f->soil_moisture ; [*]<br></div><div> lat = f->latitude ; [*]<br></div><div> lon = f->longitude ; [*]<br></div><div><br></div><div> ism = <b><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/ind.shtml" target="_blank">ind</a></b>(<b>.not</b><a href="http://.ismissing" target="_blank"><b>.ismissing</b></a>(sm) .and. lat.ge.latS .and. lat.le.latN \<br></div><div> .and. lon.ge.lonL .and. lon.le.lonR)<br></div><div> npts = <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/dimsizes.shtml"><b>dimsizes</b></a>(ism)<br></div><div><br></div><div> sm_avg = <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/avg.shtml"><b>avg</b></a>(sm(ism))<br></div><div>==============================<wbr>====================<br>Side points:<br><br></div>Alessandra's suggestions would be appropriate if the original data were on a rectilinear grid and then had been placed onto 1D arrays. Then, if you knew the appropriate nlat and mlon sizes you could have used<br><br></div> xGrid = <b><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/conform_dims.shtml" target="_blank">conform_dims</a></b>( (/nlat,mlon/), x1d) ; (nlat,mlon); similar to <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/onedtond.shtml" target="_blank"><b>onedtond</b></a><br><div><div><div><br>===<br>The documentation for <b><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/wgt_areaave.shtml" target="_blank">wgt_areaave</a></b> should be clarified. This function expects a <b>rectilinear</b> grid.<br><br>======<br><div class="gmail_extra"><div class="gmail_quote">On Sun, Aug 12, 2018 at 4:13 PM, alesall <span dir="ltr"><<a href="mailto:alesall@iri.columbia.edu" target="_blank">alesall@iri.columbia.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div><br></div>Hi Ehsan, <div><br></div><div>can you reorganize your data so that it is 2-dimensional, in latitude and longitude?</div><div>If your “var” measures something which varies in latitude and longitude, you need to figure out how the data in it is ordered/organized. </div><div>Is the fastest dimension longitude or latitude?</div><div>You can use functions like “onedtond” to break down a vector of length nlon x nlat into a matrix of dimensions nlon, nlat.</div><div>And then assign names/values of latitude and longitude to the dimensions accordingly.</div><div><br></div><div>warm regards</div><div><br></div><div><br></div><div><br></div><div><br><div>
<div>-- <br>Alessandra Giannini<br>IRI for Climate and Society - The Earth Institute at Columbia University<br>P.O. Box 1000, Palisades NY 10964-8000<br>phone/fax: +1 845 680-4473/4864 - email: <a href="mailto:alesall@iri.columbia.edu" target="_blank">alesall@iri.columbia.ed<wbr>u</a></div>
</div>
<div><br><blockquote type="cite"><div><div class="m_-5502896329409997538gmail-h5"><div>On Aug 12, 2018, at 3:16 PM, Ehsan Taghizadeh <<a href="mailto:ehsantaghizadeh@yahoo.com" target="_blank">ehsantaghizadeh@yahoo.com</a>> wrote:</div><br class="m_-5502896329409997538gmail-m_5826833313638763522Apple-interchange-newline"></div></div><div><div><div class="m_-5502896329409997538gmail-h5"><div><div style="font-family:bookman old style,new york,times,serif;font-size:16px"><div style="font-family:bookman old style,new york,times,serif;font-size:16px"><div></div>
<div>Hi,</div><div>I'm sorry to forward my question again! Do you think it is possible doing this or I should leave it?</div><div>I'll be thankful for any suggestion.</div><div><br></div><div>Sincerely</div><div>Ehsan</div><div><br></div><div><br></div>
</div><div id="m_-5502896329409997538gmail-m_5826833313638763522ydpa8de3a6fyahoo_quoted_4594534684" class="m_-5502896329409997538gmail-m_5826833313638763522ydpa8de3a6fyahoo_quoted">
<div style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px;color:rgb(38,40,42)">
<div>----- Forwarded Message -----</div>
<div><b>From:</b> Ehsan Taghizadeh <<a href="mailto:ehsantaghizadeh@yahoo.com" target="_blank">ehsantaghizadeh@yahoo.com</a>></div><div><b>To:</b> Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>></div><div><b>Sent:</b> Thursday, August 2, 2018, 8:00:22 PM GMT+4:30</div><div><b>Subject:</b> area average</div><div><br></div>
<div><div id="m_-5502896329409997538gmail-m_5826833313638763522ydpa8de3a6fyiv7902342373"><div><div style="font-family:bookman old style,new york,times,serif;font-size:16px"><div>Hi,</div><div>I've tried to compute area average using "<span>wgt_areaave</span>". However it seems my data are 1 dimensional. So the following error occurred, as it is expected:</div><div><span style="background-color:rgb(253,248,105)"> fatal:wgt_areaave: the input array must have at least 2 dimensions</span><br></div><div><br></div><div>Could I ask any help to have are average from 1 dimension data?!</div><div><br></div><div>I attached my script and put input files (<span>SPL2SMP_E.tar</span>) in ftp. Other information are:</div><div> <span style="background-color:rgb(173,215,115)"><span><span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px">$</span></span> <span><span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px">ncl -V</span></span></span></div><div><span><span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px;background-color:rgb(173,215,115)"> <span>6.5.0</span></span></span></div><div><span><span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px"><span style="background-color:rgb(173,215,115)">$ <span><span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px">uname -a</span></span></span></span></span></div><div><font face="verdana, sans-serif" color="#333333"><span style="font-size:13.3333px;background-color:rgb(173,215,115)"> <span>Linux localhost.localdomain 3.10.0-514.2.2.el7.x86_64 #1 SMP Tue Dec 6 23:06:41 UTC 2016 x86_64 x86_64 x86_64 GNU/Linux</span></span></font></div><div><font face="verdana, sans-serif" color="#333333"><span style="font-size:13.3333px"><span style="background-color:rgb(173,215,115)">$ <span><span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px">gcc --version</span></span></span></span></font></div><div><span style="background-color:rgb(173,215,115)"><font face="verdana, sans-serif" color="#333333"><span style="font-size:13.3333px"><span><span><span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px"> <span></span></span></span></span></span></font><span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px">gcc (GCC) 4.8.5 20150623 (Red Hat 4.8.5-28)</span></span><font face="verdana, sans-serif" color="#333333"><span style="font-size:13.3333px"><span><span><span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px"><span><div><span style="background-color:rgb(173,215,115)"> Copyright (C) 2015 Free Software Foundation, Inc.</span></div><div><span><span style="font-family:"old times",serif;font-size:16px"><br></span></span></div><div><span><span style="font-family:"old times",serif;font-size:16px">I'll be thankful for any help.</span></span><br></div><div><span><span style="font-family:"old times",serif;font-size:16px"><br></span></span></div><div><span><span style="font-family:"old times",serif;font-size:16px">Sincerely</span></span></div><div><span><span style="font-family:"old times",serif;font-size:16px">Ehsan</span></span></div></span></span></span></span></span></font></div><div style="font-family:bookman old style,new york,times,serif;font-size:16px"></div></div></div></div></div>
</div>
</div></div></div></div></div><span id="m_-5502896329409997538gmail-m_5826833313638763522cid:7576075c-30ea-f3cd-0666-3efb5f1294cc@yahoo.com"><SPL2SMP_E_areaave.ncl></span>_______<wbr>______________________________<wbr>__________<br>ncl-talk mailing list<br><a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>List instructions, subscriber options, unsubscribe:<br><a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br></div></blockquote></div><br></div></div><br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div></div></div></div></div>