<html><head></head><body><div style="font-family:Helvetica Neue, Helvetica, Arial, sans-serif;font-size:16px;"><div style="font-family:Helvetica Neue, Helvetica, Arial, sans-serif;font-size:16px;"><div></div>
<div>Dears,</div><div>Pardon me to forward my questions twice. May I ask how using attached script for multiple files?</div><div><br></div><div>Sincerely</div><div>Ehsan</div><div><br></div>
</div><div id="ydpb4fa9243yahoo_quoted_5055561477" class="ydpb4fa9243yahoo_quoted">
<div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;">
<div>----- Forwarded Message -----</div>
<div><b>From:</b> Ehsan Taghizadeh <ehsantaghizadeh@yahoo.com></div><div><b>To:</b> Ncl-talk <ncl-talk@ucar.edu>; alesall@iri.columbia.edu <alesall@iri.columbia.edu></div><div><b>Sent:</b> Tuesday, August 14, 2018, 7:03:20 PM GMT+4:30</div><div><b>Subject:</b> Re: [ncl-talk] Fw: area average</div><div><br></div>
<div><div id="ydpb4fa9243yiv0406576422"><div><div style="font-family:bookman old style, new york, times, serif;font-size:16px;"><div>Thank you for your helps.</div><div>Dear Alessandra</div><div>Thank you for you comments. A reference of my data could be this page (<a shape="rect" href="https://nsidc.org/data/SPL2SMP_E/versions/2" rel="nofollow" target="_blank">https://nsidc.org/data/SPL2SMP_E/versions/2</a>). However Dennis's comments<span><span style="color:rgb(0, 0, 0);"> about my data</span></span> are true. I've attached a picture of these data for date 1 April 2018.</div><div>Also the result of PrintVarSummary(Var):</div><div><br clear="none"></div><div><span></span><div><font face="times new roman, new york, times, serif" size="2"> Variable: var</font></div><div><font face="times new roman, new york, times, serif" size="2"> Type: float</font></div><div><font face="times new roman, new york, times, serif" size="2"> Total Size: 1124096 bytes</font></div><div><font face="times new roman, new york, times, serif" size="2"> 281024 values</font></div><div><font face="times new roman, new york, times, serif" size="2"> Number of Dimensions: 1</font></div><div><font face="times new roman, new york, times, serif" size="2"> Dimensions and sizes: [DIM_000 | 281024]</font></div><div><font face="times new roman, new york, times, serif" size="2"> Coordinates:</font></div><div><font face="times new roman, new york, times, serif" size="2"> Number Of Attributes: 6</font></div><div><font face="times new roman, new york, times, serif" size="2"> _FillValue : -9999</font></div><div><font face="times new roman, new york, times, serif" size="2"> coordinates : /Soil_Moisture_Retrieval_Data/latitude /Soil_Moisture_Retrieval_Data/longitude</font></div><div><font face="times new roman, new york, times, serif" size="2"> long_name : Representative soil moisture measurement for the Earth based grid cell.</font></div><div><font face="times new roman, new york, times, serif" size="2"> units : cm**3/cm**3</font></div><div><font face="times new roman, new york, times, serif" size="2"> valid_max : 0.5</font></div><div><font face="times new roman, new york, times, serif" size="2"> valid_min : 0.02</font></div><div><br clear="none"></div></div><div>Please let me know if I should explain more.</div><div><br clear="none"></div><div>I've used Dennis's suggestion, and it worked well. However may I know 2 extra points:</div><div>1- I've tried to use attached script to get average from multiple files, but an error as below has occurred:</div><div><br clear="none"></div><div><span style="background-color:rgb(253, 248, 105);"><font face="courier new, courier, monaco, monospace, sans-serif" size="2"><span>fatal:Dimension sizes of left hand side and right hand side of assignment do not match</span><br clear="none"></font></span></div><div><span><span style="background-color:rgb(253, 248, 105);"><font face="courier new, courier, monaco, monospace, sans-serif" size="2">fatal:["Execute.c":8637]:Execute: Error occurred at or near line 61 in file SPL2SMP_E_areaave_mod1.ncl</font></span><br clear="none"></span></div><div><br clear="none"></div><div>This error could be expected because of the below lines in loop:<br clear="none"><font face="times new roman, new york, times, serif" size="2"><span> var<span style="white-space:pre-wrap;"> </span>= SPL2_E_f->$var_path$</span><br clear="none"></font></div><div><span><font face="times new roman, new york, times, serif" size="2"><span> lat<span style="white-space:pre-wrap;"> </span>= SPL2_E_f->$lat_path$</span><br clear="none"></font></span></div><div><span><span><font face="times new roman, new york, times, serif" size="2"><span> lon<span style="white-space:pre-wrap;"> </span>= SPL2_E_f->$lon_path$</span><br clear="none"></font></span></span></div><div><span><span><font face="times new roman, new york, times, serif" size="2"><span><span> sm_avg(i) = avg(var(ism))</span><br clear="none"></span></font></span></span></div><div><span><span><font face="times new roman, new york, times, serif" size="2"><span><span><br clear="none"></span></span></font></span></span></div><div>However using "<font face="times new roman, new york, times, serif" size="2">:=</font>" couldn't help so.</div><div>May I ask help again about using above lines in loop.</div><div><br clear="none"></div><div>2- May I ask how compute areal average of these data? I mean which function of ncl could help me about that (just for curiously!).</div><div><br clear="none"></div><div>Sincerely</div><div>Ehsan</div><div><br clear="none"></div><div class="ydpb4fa9243yiv0406576422yahoo_quoted" id="ydpb4fa9243yiv0406576422yahoo_quoted_5180311432">
<div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;">
<div class="ydpb4fa9243yiv0406576422yqt6806569319" id="ydpb4fa9243yiv0406576422yqtfd10056"><div>
On Monday, August 13, 2018, 8:23:47 AM GMT+4:30, Dennis Shea <shea@ucar.edu> wrote:
</div>
<div><br clear="none"></div>
<div><br clear="none"></div>
<div><div id="ydpb4fa9243yiv0406576422"><div><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 clear="none"></div><div><br clear="none"></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 clear="none"><br clear="none"></div><div>Your code comments the following:<br clear="none"></div><div><br clear="none"><b>/;</b><br clear="none"> --37.80-----<br clear="none"> | |<br clear="none"> 45.8 46.03<br clear="none"> | |<br clear="none"> ---37.45----<br clear="none"><b>;/</b><br clear="none"><br clear="none"></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 clear="none"><br clear="none"></div><div> latS = 37.45<br clear="none"></div><div> latN = 37.80<br clear="none"></div><div> lonL = 45.8<br clear="none"></div><div> lonR = 4603<br clear="none"><br clear="none"></div><div> sm = f->soil_moisture ; [*]<br clear="none"></div><div> lat = f->latitude ; [*]<br clear="none"></div><div> lon = f->longitude ; [*]<br clear="none"></div><div><br clear="none"></div><div> ism = <b><a shape="rect" href="http://www.ncl.ucar.edu/Document/Functions/Built-in/ind.shtml" rel="nofollow" target="_blank">ind</a></b>(<b>.not</b><a shape="rect" href="http://.ismissing" rel="nofollow" target="_blank"><b>.ismissing</b></a>(sm) .and. lat.ge.latS .and. lat.le.latN \<br clear="none"></div><div> .and. lon.ge.lonL .and. lon.le.lonR)<br clear="none"></div><div> npts = <a shape="rect" href="http://www.ncl.ucar.edu/Document/Functions/Built-in/dimsizes.shtml" rel="nofollow" target="_blank"><b>dimsizes</b></a>(ism)<br clear="none"></div><div><br clear="none"></div><div> sm_avg = <a shape="rect" href="http://www.ncl.ucar.edu/Document/Functions/Built-in/avg.shtml" rel="nofollow" target="_blank"><b>avg</b></a>(sm(ism))<br clear="none"></div><div>============================== ====================<br clear="none">Side points:<br clear="none"><br clear="none"></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 clear="none"><br clear="none"></div> xGrid = <b><a shape="rect" href="http://www.ncl.ucar.edu/Document/Functions/Built-in/conform_dims.shtml" rel="nofollow" target="_blank">conform_dims</a></b>( (/nlat,mlon/), x1d) ; (nlat,mlon); similar to <a shape="rect" href="http://www.ncl.ucar.edu/Document/Functions/Built-in/onedtond.shtml" rel="nofollow" target="_blank"><b>onedtond</b></a><br clear="none"><div><div><div><br clear="none">===<br clear="none">The documentation for <b><a shape="rect" href="http://www.ncl.ucar.edu/Document/Functions/Built-in/wgt_areaave.shtml" rel="nofollow" target="_blank">wgt_areaave</a></b> should be clarified. This function expects a <b>rectilinear</b> grid.<br clear="none"><br clear="none">======<br clear="none"><div class="ydpb4fa9243yiv0406576422gmail_extra"><div class="ydpb4fa9243yiv0406576422gmail_quote">On Sun, Aug 12, 2018 at 4:13 PM, alesall <span dir="ltr"><<a shape="rect" href="mailto:alesall@iri.columbia.edu" rel="nofollow" target="_blank">alesall@iri.columbia.edu</a>></span> wrote:<br clear="none"><blockquote class="ydpb4fa9243yiv0406576422gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex;"><div class="ydpb4fa9243yiv0406576422yqt8774489905" id="ydpb4fa9243yiv0406576422yqt11037"><div><div><br clear="none"></div>Hi Ehsan, <div><br clear="none"></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 clear="none"></div><div>warm regards</div><div><br clear="none"></div><div><br clear="none"></div><div><br clear="none"></div><div><br clear="none"><div>
<div>-- <br clear="none">Alessandra Giannini<br clear="none">IRI for Climate and Society - The Earth Institute at Columbia University<br clear="none">P.O. Box 1000, Palisades NY 10964-8000<br clear="none">phone/fax: +1 845 680-4473/4864 - email: <a shape="rect" href="mailto:alesall@iri.columbia.edu" rel="nofollow" target="_blank">alesall@iri.columbia.ed u</a></div>
</div>
<div><br clear="none"><blockquote type="cite"><div><div class="ydpb4fa9243yiv0406576422m_-5502896329409997538gmail-h5"><div>On Aug 12, 2018, at 3:16 PM, Ehsan Taghizadeh <<a shape="rect" href="mailto:ehsantaghizadeh@yahoo.com" rel="nofollow" target="_blank">ehsantaghizadeh@yahoo.com</a>> wrote:</div><br clear="none" class="ydpb4fa9243yiv0406576422m_-5502896329409997538gmail-m_5826833313638763522Apple-interchange-newline"></div></div><div><div><div class="ydpb4fa9243yiv0406576422m_-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 clear="none"></div><div>Sincerely</div><div>Ehsan</div><div><br clear="none"></div><div><br clear="none"></div>
</div><div class="ydpb4fa9243yiv0406576422m_-5502896329409997538gmail-m_5826833313638763522ydpa8de3a6fyahoo_quoted" id="ydpb4fa9243yiv0406576422m_-5502896329409997538gmail-m_5826833313638763522ydpa8de3a6fyahoo_quoted_4594534684">
<div>
<div>----- Forwarded Message -----</div>
<div><b>From:</b> Ehsan Taghizadeh <<a shape="rect" href="mailto:ehsantaghizadeh@yahoo.com" rel="nofollow" target="_blank">ehsantaghizadeh@yahoo.com</a>></div><div><b>To:</b> Ncl-talk <<a shape="rect" href="mailto:ncl-talk@ucar.edu" rel="nofollow" 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 clear="none"></div>
<div><div id="ydpb4fa9243yiv0406576422m_-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 clear="none"></div><div><br clear="none"></div><div>Could I ask any help to have are average from 1 dimension data?!</div><div><br clear="none"></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></span></span></span></span></span></font><div><span style="background-color:rgb(173,215,115);"> Copyright (C) 2015 Free Software Foundation, Inc.</span></div><div><span><span><br clear="none"></span></span></div><div><span><span>I'll be thankful for any help.</span></span><br clear="none"></div><div><span><span><br clear="none"></span></span></div><div><span><span>Sincerely</span></span></div><div><span><span>Ehsan</span></span></div></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><SPL2SMP_E_areaave.ncl></span>_______ ______________________________ __________<br clear="none">ncl-talk mailing list<br clear="none"><a shape="rect" href="mailto:ncl-talk@ucar.edu" rel="nofollow" target="_blank">ncl-talk@ucar.edu</a><br clear="none">List instructions, subscriber options, unsubscribe:<br clear="none"><a shape="rect" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="nofollow" target="_blank">http://mailman.ucar.edu/mailma n/listinfo/ncl-talk</a><br clear="none"></div></blockquote></div><br clear="none"></div></div></div><br clear="none">______________________________ _________________<br clear="none">
ncl-talk mailing list<br clear="none">
<a shape="rect" href="mailto:ncl-talk@ucar.edu" rel="nofollow" target="_blank">ncl-talk@ucar.edu</a><br clear="none">
List instructions, subscriber options, unsubscribe:<br clear="none">
<a shape="rect" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="nofollow" target="_blank">http://mailman.ucar.edu/mailma n/listinfo/ncl-talk</a><br clear="none">
<br clear="none"></blockquote></div><br clear="none"></div></div></div></div></div></div></div></div>
</div></div><div class="ydpb4fa9243yiv0406576422yqt6806569319" id="ydpb4fa9243yiv0406576422yqtfd48578">
</div></div></div></div></div></div>
</div>
</div></div></body></html>