<div dir="ltr"><div><div><div><div>Dear Dennis Sir,<br><br></div>Thanks for this explanation. <br><br></div>After changing/or reorder the dimension name then script ran successfully but with a warming<br><span style="color:rgb(255,0,0)"><b>warning:esccr: Non-fatal conditions encountered: all missing or constant values<br><br></b></span></div><span style="color:rgb(255,0,0)"><font color="#000000">After searching this query most of the users said<br><br>The warning massage is generally occurred if one or several time, <br>series have all the same values. The correlation coefficient is <br>calculated by dividing standard deviation. If all values are same, <br>the standard deviation is zero.<br><br><br></font></span></div><span style="color:rgb(255,0,0)"><font color="#000000">One user also said it could be ignored (<a href="https://www.ncl.ucar.edu/Support/talk_archives/2010/0175.html">https://www.ncl.ucar.edu/Support/talk_archives/2010/0175.html</a>)</font><span style="color:rgb(0,0,0)">. <br>So can we really ignore this warning?</span><b> </b><span style="color:rgb(0,0,0)">if yes, then is there any possibility to remove this warning.</span><b><br></b></span><div><div><div><div><br><br><br></div><div><b>Query regarding :</b><br><b>The last dimension of x must be equal to the last dimension of y</b><br></div><div>In the previous mail, this query removed by changing the order of dimensions name. But if the only dimension size is given not dimension name then how can we set the dimension. <br><br></div><div>Such as<br><br>Variable: CoD_final<br>Type: float<br>Total Size: 5076540 bytes<br>            1269135 values<br>Number of Dimensions: 2<br><b>Dimensions and sizes:    [1185] x [1071]</b><br>Coordinates: <br>Number Of Attributes: 1<br>  _FillValue :    -28672<br>(0)    min=0   max=1<br><br>Variable: ccr<br>Type: float<br>Total Size: 5076540 bytes<br>            1269135 values<br>Number of Dimensions: 3<br><b>Dimensions and sizes:    [1185] x [1071] x [1]</b><br>Coordinates: <br>Number Of Attributes: 1<br>  _FillValue :    -28672<br><span style="color:rgb(0,0,255)"><b>(0)    min=-2.13598   max=2.74001  ===> also, I am getting correaltion values greater than 1</b></span><br>    <br><br></div><div><span style="color:rgb(0,0,255)"><b>I tried to calculate covarience using<br>ccv = escovc(CoD_final,ccr(time|:, lat|:, lon|:))<br></b></span></div><div><span style="color:rgb(0,0,255)"><b>then fatal:(lon) is not a dimension name in variable (ccr), could not determine dimension number<br></b></span></div><div><span style="color:rgb(0,0,255)"><b><br>the I tried<br>ccv = escovc(CoD_final,ccr(:,:,0))</b></span><br></div><div><br></div><div>But after trying this I got several warning<br></div><div><br>warning:escovc: Non-fatal conditions encountered in series or xstd = 0.0<br>(0)    is_valid_latlon2d_attr: Warning: The 'lat2d' attribute must either be<br>(0)    the same dimension sizes as the data, or one element larger in both directions.<br>(0)    Your data will most likely not be overlaid on the map correctly.<br>(0)    check_for_y_lat_coord: Warning: Data either does not contain<br>(0)    a valid latitude coordinate array or doesn't contain one at all.<br>(0)    A valid latitude coordinate array should have a 'units'<br>(0)    attribute equal to one of the following values: <br>(0)        'degrees_north' 'degrees-north' 'degree_north' 'degrees north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'<br>(0)    is_valid_latlon2d_attr: Warning: The 'lon2d' attribute must either be<br>(0)    the same dimension sizes as the data, or one element larger in both directions.<br>(0)    Your data will most likely not be overlaid on the map correctly.<br>(0)    check_for_lon_coord: Warning: Data either does not contain<br>(0)    a valid longitude coordinate array or doesn't contain one at all.<br>(0)    A valid longitude coordinate array should have a 'units'<br>(0)    attribute equal to one of the following values: <br>(0)        'degrees_east' 'degrees-east' 'degree_east' 'degrees east' 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'<br>(0)    Error: scalar_field: If the input data is 1-dimensional, you must set sfXArray and sfYArray to 1-dimensional arrays of the same length.<br>warning:create: Bad HLU id passed to create, ignoring it<br><br><br> <br></div></div></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>---<br></div><div><font size="2"><span>Kunal Bali<br></span></font></div><br><div><br></div><div><p style="margin:0px;border-collapse:collapse;font-family:Tahoma,Verdana;font-size:12px"><font color="#1F497D"><br></font></p></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
<br><div class="gmail_quote">On Fri, Feb 23, 2018 at 8:51 PM, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.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 dir="ltr"><div>As we have said before, you must read the error message and then look at the function being called and the day being used. Further, looking at the Examples helps a lot.<br><br>===<br><span style="color:rgb(255,0,0)"><span class="gmail-">fatal:esccr: the rightmost dimension of x and y must be the same<br>fatal:["Execute.c":8640]:Execu<wbr>te: Error occurred at or near line 141 in file ccr.ncl<br></span>===<br><br></span></div><span style="color:rgb(255,0,0)">look at line 141<br><br>ccr = esccr(ts1,ts2,maxlag)         <wbr>          ; calc cross correlations<br><br clear="all"></span><div><div><div><div class="gmail-m_8407371779927078199gmail-m_1846718283174335763gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>at lines 133,134, you have<br><br>  ts1 = data_single(:)<br>  ts2 = data_subset(:,:,:)<br><br><b><span style="color:rgb(255,0,0)">Variable: data_single</span>      ==> ts1<span class="gmail-"><br>Dimensions and sizes:    [time | 96]</span></b><br><br><b><span style="color:rgb(255,0,0)">Variable: data_subset</span>      ==> ts2<span class="gmail-"><br>Dimensions and sizes:    [time | 96] x [lat | 1185] x [lon | 1071]<br><br></span></b></div><div>As noted in the function documentation,<br><i>x</i>: An array of any numeric type or size. The rightmost dimension is usually time.
<i></i>
<i><br><br>y: </i>An array of any numeric type or size. The rightmost dimension is usually time. 
<span style="color:rgb(0,0,255)"><b>The size of the rightmost dimension must be the same as <i>x</i>.
</b></span>
<p class="gmail-m_8407371779927078199gmail-indent">====</p><p class="gmail-m_8407371779927078199gmail-indent">Clearly. the "rightmost" dimension of 'ts1' [size=1071'] does not match the size of ts1' [size=96]<br></p><p class="gmail-m_8407371779927078199gmail-indent"> Example 2 in the documentation, show how the address this situation. Use, NCL's dimension reordering:<br></p><p class="gmail-m_8407371779927078199gmail-indent">ccr = esccr(ts1,ts2(lat|:,lon|:,<wbr>time|:), maxlag)<br></p><br></div><div><b></b></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="gmail-h5">On Fri, Feb 23, 2018 at 4:07 AM, Kunal Bali <span dir="ltr"><<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div class="gmail-h5"><div dir="ltr"><div><div>Dear NCL user,<br><br></div>I am using the attached script to calculate the cross-correlation between single grid cell to neighbor grid cells with time. <br></div><br><div><div><div>I am getting an error<br><span style="color:rgb(255,0,0)">fatal:esccr: the rightmost dimension of x and y must be the same<br>fatal:["Execute.c":8640]:Execu<wbr>te: Error occurred at or near line 141 in file ccr.ncl<br><br clear="all"></span><div><div><div class="gmail-m_8407371779927078199m_1846718283174335763gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>---</div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>The summary is here:<br><br><br></div><div> <b>Copyright (C) 1995-2017 - All Rights Reserved<br> University Corporation for Atmospheric Research<br> NCAR Command Language Version 6.4.0<br> The use of this software is governed by a License Agreement.<br> See <a href="http://www.ncl.ucar.edu/" target="_blank">http://www.ncl.ucar.edu/</a> for more details.<br><br>Variable: totc<br>Type: float<br>Total Size: <a href="tel:(202)%20508-6976" value="+12025086976" target="_blank">2025086976</a> bytes<br>            506271744 values<br>Number of Dimensions: 3<br>Dimensions and sizes:    [time | 96] x [lat | 2112] x [lon | 2497]<br>Coordinates: <br>            time: [35430.125..104813.<a href="tel:079%201666%206667" value="+917916666667" target="_blank">7916666667</a>]<br>            lat: [5.0078125..37.9921875]<br>            lon: [  59..  98]<br>Number Of Attributes: 8<br>  long_name :    AOT at 0.55 micron<br>  hdfeos_name :    Optical_Depth_055<br>  projection :    Albers Conical Equal_Area<br>  unit :    None<br>  _FillValue_original :    -28672<br>  _FillValue :    -28672<br>  missing_value_original :    -28672<br>  missing_value :    -28672<br><br><span style="color:rgb(255,0,0)">Variable: data_subset</span><br>Type: float<br>Total Size: 487347840 bytes<br>            121836960 values<br>Number of Dimensions: 3<br>Dimensions and sizes:    [time | 96] x [lat | 1185] x [lon | 1071]<br>Coordinates: <br>            time: [35430.125..104813.<a href="tel:079%201666%206667" value="+917916666667" target="_blank">7916666667</a>]<br>            lat: [23.5078125..5.0078125]<br>            lon: [72.78125..89.5]<br>Number Of Attributes: 8<br>  missing_value :    -28672<br>  missing_value_original :    -28672<br>  _FillValue :    -28672<br>  _FillValue_original :    -28672<br>  unit :    None<br>  projection :    Albers Conical Equal_Area<br>  hdfeos_name :    Optical_Depth_055<br>  long_name :    AOT at 0.55 micron<br>(0)    nearest location: k=0  n=541 m=1291  13.4609375  79.17187499999999<br>(0)    -----<br><br><span style="color:rgb(255,0,0)">Variable: data_single</span><br>Type: float<br>Total Size: 384 bytes<br>            96 values<br>Number of Dimensions: 1<br>Dimensions and sizes:    [time | 96]<br>Coordinates: <br>            time: [35430.125..104813.<a href="tel:079%201666%206667" value="+917916666667" target="_blank">7916666667</a>]<br>Number Of Attributes: 10<br>  lon :    79.171875<br>  lat :    13.4609375<br>  missing_value :    -28672<br>  missing_value_original :    -28672<br>  _FillValue :    -28672<br>  _FillValue_original :    -28672<br>  unit :    None<br>  projection :    Albers Conical Equal_Area<br>  hdfeos_name :    Optical_Depth_055<br>  long_name :    AOT at 0.55 micron<br>(0)    -----</b><br><span style="color:rgb(255,0,0)">fatal:esccr: the rightmost dimension of x and y must be the same<br>fatal:["Execute.c":8640]:Execu<wbr>te: Error occurred at or near line 141 in file ccr.ncl</span><br clear="all"><div><div><div class="gmail-m_8407371779927078199m_1846718283174335763gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><br><br>---<br></div><div><font size="2"><span>Kunal Bali<br></span></font></div><br><div><br></div><div><p style="margin:0px;border-collapse:collapse;font-family:Tahoma,Verdana;font-size:12px"><font color="#1F497D"><br></font></p></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</div></div></div></div></div>
<br></div></div>______________________________<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>
</blockquote></div><br></div></div>