<div dir="ltr">Hi Melanie,<div><span style="font-size:12.8px">clmMonTLL and </span><span style="font-size:12.8px">calcMonAnomTLL expect the input arrays to be 3-dimensional and ordered time x latitude x longitude. Your x2 array is a 4-dimensional array. Try doing this to eliminate the depth dimension:</span></div><div><span style="font-size:12.8px"> x2Clm = clmMonTLL(x2(:,0,:,:)) ; (12,lat,lon;)</span></div><div style="font-size:12.8px">x2Anom = calcMonAnomTLL(x2(:,0,:,:), xClm) ; (time, lat,lon)</div><div><br></div><div>alternatively, you can use rmMonAnnCycTLL to do the same thing:</div><div>x2Anom = rmMonAnnCycTLL(x2(:,0,:,:))</div><div><br></div><div><a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/rmMonAnnCycTLL.shtml">https://www.ncl.ucar.edu/Document/Functions/Contributed/rmMonAnnCycTLL.shtml</a><br><br></div><div>Hope that helps!</div><div>Adam</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Sep 1, 2017 at 12:58 PM, Melanie O' hanoly <span dir="ltr"><<a href="mailto:mel.hanoly@gmail.com" target="_blank">mel.hanoly@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Adam..<div><br></div><div>Thank you so much for help me..</div><div>But now have some problems in calculating anomalies</div><div><br></div><div><div>Variable: x</div><div>Type: double</div><div>Total Size: 29760000 bytes</div><div> 3720000 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[time | 372] x [lat | 100] x [lon | 100]</div><div>Coordinates: </div><div> time: [ 0..16260480]</div><div> lat: [-89.09999999999999..89.<wbr>10000000000002]</div><div> lon: [ 0..356.4]</div><div>Number Of Attributes: 2</div><div> _FillValue :<span style="white-space:pre-wrap"> </span>-999000000</div><div> missing_value :<span style="white-space:pre-wrap"> </span>-999000000</div><div><br></div><div>Variable: x2</div><div>Type: float</div><div>Total Size: 14880000 bytes</div><div> 3720000 values</div><div>Number of Dimensions: 4</div><div>Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[time | 372] x [depth | 1] x [latitude | 100] x [longitude | 100]</div></div><div><br></div><div><div>fatal:Number of dimensions in parameter (0) of (clmMonTLL) is (4), (3) dimensions were expected </div><div>fatal:["Execute.c":8640]:<wbr>Execute: Error occurred at or near line 86 in file eof.ncl</div><div><br></div><div>fatal:Number of dimensions in parameter (0) of (calcMonAnomTLL) is (4), (3) dimensions were expected </div><div>fatal:["Execute.c":8640]:<wbr>Execute: Error occurred at or near line 88 in file eof.ncl</div><div><br></div><div>This is my script</div><div><br></div><div>"; ==============================<wbr>==============================<wbr>==</div><div>; compute climatology and Anomalies 2 file</div><div>; ==============================<wbr>==============================<wbr>=</div><div> (86) x2Clm = clmMonTLL(x2) ; (12,lat,lon;)</div><div> (87) </div><div> (88) x2Anom = calcMonAnomTLL(x2, xClm) ; (time, lat,lon)</div><div><br></div><div> printVarSummary(x2Anom) </div><div> printMinMax(x2Anom, True) </div><div>; ==============================<wbr>==============================<wbr>=====" </div><div><br></div><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Aug 22, 2017 at 5:42 PM, Adam Phillips <span dir="ltr"><<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Melanie,<div>regCoef operates on the rightmost dimensions of the two input arrays, and expects those two rightmost dimensions to have the same size. See the documentation here:</div><div><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/regCoef-1.shtml" target="_blank">https://www.ncl.ucar.edu/Docum<wbr>ent/Functions/Built-in/regCoef<wbr>-1.shtml</a><br></div><div>referring specifically to the Arguments - y section near the top.</div><div><br></div><div>Thus, assuming your EOF/principal component timeseries was computed over the TIME1 dimension, you should change this:</div><span><div>eof_regres(ne,:,:) = (/ regCoef(eof_ts(ne,:), x2(TIME1|:,DEPTH1_1|:,LAT|:,LO<wbr>N|:))/)<br></div></span><div>to this:</div><div>eof_regres(ne,:,:) = (/ regCoef(eof_ts(ne,:), x2(DEPTH1_1|0,LAT|:,LON|:,TIME<wbr>1|:))/)</div><div>(Note that I also changed "DEPTH1_1|:" to "DEPTH1_1|0" so that the first and only index is chosen, that way the array returned to the lefthand side of the equal sign will be dimensioned lat (100) x lon (100).)</div><div><br></div><div>If you were using regCoef_n, you could do this:</div><div>eof_regres(ne,:,:) = (/ regCoef_n(eof_ts(ne,:), x2(TIME1|:,DEPTH1_1|0,LAT|:,LO<wbr>N|:),0,0) /)<br></div><div><br></div><div>The last two arguments represent the dimension number to operate on for the eof_ts and x2 arrays. In this specific case, you are index subscripting eof_ts from a 2D array to a 1D array, and thus regCoef_n will only see a 1D array and the 1st dimension should be operated on. Finally, you are dimension reordering x2, and you don't really need to do this if you were to use regCoef_n. You just need to indicate the correct dimension to operate on with the last provided argument.</div><div><br></div><div>Hope that all helps. If you have any further questions please respond to the ncl-talk email list.</div><div>Adam </div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="m_5235427741404856551h5">On Tue, Aug 22, 2017 at 2:19 PM, Melanie O' hanoly <span dir="ltr"><<a href="mailto:mel.hanoly@gmail.com" target="_blank">mel.hanoly@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="m_5235427741404856551h5"><div dir="ltr"><div style="font-size:12.8px"><pre style="white-space:pre-wrap;color:rgb(0,0,0)"><i>Dear NCL users, I just have a little knowledge on NCL. So, could anyone help me?</i></pre><pre style="white-space:pre-wrap;color:rgb(0,0,0)"><i>I have two data with 372 timestep, one has 3D and the other 4D, someone please help me with this error?</i></pre><pre style="white-space:pre-wrap;color:rgb(0,0,0)"><i style="font-family:arial,sans-serif">"</i><font color="#000000" style="font-family:arial,sans-serif"><i>fatal:regCoef: The rightmost dimension of x must be equal to the rightmost dimension of y"</i></font><br></pre><pre style="white-space:pre-wrap;color:rgb(0,0,0)"><font color="#000000" style="font-family:arial,sans-serif"><i>The summary is bellow</i></font></pre><font color="#000000" face="monospace"><i><div style="white-space:pre-wrap"><font color="#000000" face="monospace"><i>This is the scrip line error </i></font></div><div style="white-space:pre-wrap"><font color="#000000" face="monospace"><i><br></i></font></div><div><font color="#000000" face="monospace" style="white-space:pre-wrap"><i>"</i></font><span style="white-space:pre-wrap">; Regress</span></div><span style="white-space:pre-wrap">; ==============================<wbr>==============================<wbr>=====
eof_regres = eof ; create an array w meta data
do ne=0,neof-1
eof_regres(ne,:,:) = (/ regCoef(eof_ts(ne,:), x2(TIME1|:,DEPTH1_1|:,LAT|:,LO<wbr>N|:))/)
end do
printVarSummary(eof_regres)
printMinMax(eof_regres, True)</span></i></font></div><div style="font-size:12.8px"><font color="#000000" face="monospace"><span style="white-space:pre-wrap"><i>"</i></span></font></div><div style="font-size:12.8px"><font color="#000000" face="monospace"><span style="white-space:pre-wrap"><i><br></i></span></font></div><div style="font-size:12.8px"><font color="#000000" face="monospace"><i><span style="white-space:pre-wrap">If i use regCoef_n, shows this error "fatal:syntax error: function regCoef_n expects 4 arguments, got 2"<br></span></i><i><span style="white-space:pre-wrap"> </span></i></font></div><span style="font-size:12.8px">Summary</span><div style="font-size:12.8px"><div>Copyright (C) 1995-2017 - All Rights Reserved</div><div> University Corporation for Atmospheric Research</div><div> NCAR Command Language Version 6.4.0</div><div> The use of this software is governed by a License Agreement.</div><div> See <a href="http://www.ncl.ucar.edu/" target="_blank">http://www.ncl.ucar.edu/</a> <wbr>for more details.</div><div><br></div><div>Variable: x</div><div>Type: double</div><div>Total Size: 29760000 bytes</div><div> 3720000 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes:<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>[time | 372] x [lat | 100] x [lon | 100]</div><div>Coordinates: </div><div> time: [ 0..16260480]</div><div> lat: [-89.09999999999999..89.100000<wbr>00000002]</div><div> lon: [ 0..356.4]</div><div>Number Of Attributes: 2</div><div> _FillValue :<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>-999000000</div><div> missing_value :<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>-999000000</div><div><br></div><div>Variable: x2</div><div>Type: double</div><div>Total Size: 29760000 bytes</div><div> 3720000 values</div><div>Number of Dimensions: 4</div><div>Dimensions and sizes:<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>[TIME1 | 372] x [DEPTH1_1 | 1] x [LAT | 100] x [LON | 100]</div><div>Coordinates: </div><div> TIME1: [1928.06849316134..1958.561643<wbr>829649]</div><div> DEPTH1_1: [17.5..17.5]</div><div> LAT: [-89.1..89.1]</div><div> LON: [ 0..356.4]</div><div>Number Of Attributes: 4</div><div> long_name :<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>DIC[L=1:372]</div><div> _FillValue :<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>-9.999999999999999e+33</div><div> missing_value :<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>-9.999999999999999e+33</div><div> history :<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>From dic_uvic</div><div><br></div><div>Variable: eof_ts</div><div>Type: double</div><div>Total Size: 2976 bytes</div><div> 372 values</div><div>Number of Dimensions: 2</div><div>Dimensions and sizes:<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>[evn | 1] x [time | 372]</div><div>Coordinates: </div><div> evn: [1..1]</div><div> time: [ 0..16260480]</div><div>Number Of Attributes: 3</div><div> ts_mean :<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>1.196080271862835e-13</div><div> matrix :<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>covariance</div><div> _FillValue :<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>-999000000</div><div><br></div><div>fatal:regCoef: The rightmost dimension of x must be equal to the rightmost dimension of y</div><div>fatal:["Execute.c":8640]:Execu<wbr>te: Error occurred at or near line 128 in file eof_cfsr.ncl</div><div><br></div><div>Variable: eof_regres</div><div>Type: double</div><div>Total Size: 80000 bytes</div><div> 10000 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes:<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>[evn | 1] x [lat | 100] x [lon | 100]</div><div>Coordinates: </div><div> evn: [1..1]</div><div> lat: [-89.09999999999999..89.100000<wbr>00000002]</div><div> lon: [ 0..356.4]</div><div>Number Of Attributes: 6</div><div> _FillValue :<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>-999000000</div><div> method :<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>transpose</div><div> matrix :<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>covariance</div><div> pcvar :<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>13.29658</div><div> eval :<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>634047.31020702</div><div> eval_transpose :<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>24386.43500796231</div><div>(0)<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span> </div><div>(0)<span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-m_1757016587382215491gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>min=-0.04332826535433597 max=0.02833814180143831</div></div><div class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-yj6qo m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-ajU" style="margin:2px 0px 0px;font-size:12.8px"><div id="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-:10x" class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-ajR"><img class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-ajT" src="https://ssl.gstatic.com/ui/v1/icons/mail/images/cleardot.gif"></div></div><span class="m_5235427741404856551m_-6559124948943822219HOEnZb"><font color="#888888"><span class="m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-HOEnZb m_5235427741404856551m_-6559124948943822219m_7784903099633192123gmail-adL" style="font-size:12.8px"><font color="#888888"><div><br></div><div><br></div><div><b><font color="#ff00ff">Melanie O'Hanoly</font></b></div><div><b><font color="#ff00ff">Oceangrapher</font></b></div><div><b><font color="#ff00ff">PhD canditate at <span style="line-height:18.2px">Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research</span></font></b></div></font></span></font></span></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><span class="m_5235427741404856551HOEnZb"><font color="#888888"><br><br clear="all"><span class="HOEnZb"><font color="#888888"><div><br></div>-- <br><div class="m_5235427741404856551m_-6559124948943822219gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist, </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphill<wbr>i/</a> </font></span><span><font color="#888888"><a href="tel:(303)%20497-1726" value="+13034971726" target="_blank">303-497-1726</a> </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>
</font></span></font></span></div>
</blockquote></div><br></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><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist, </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphilli/</a> </font></span><span><font color="#888888">303-497-1726 </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>
</div>