<div dir="ltr">If you have a one-dimensional array (1D; eg: a 'time series') and a three-dimensional array (3D; eg: time series at multiple grid points), you *must* make the arrays 'conformant' for an array calculation<br><div><br>You have:<br></div><div>   CoD1 = ((data_single(:) - data_subset(:,:,:))/(data_<wbr>single(:) + data_subset(:,:,:)))^2<br><br></div><div>You did not include a 'printVarSummary' of these variables. Assuming data_single(ntim) and data_subset(ntim,nlat,mlon) then you must replicate (aka, propogate) the 1D array to the 3D array size. In NCL, this is called 'array conformance' . Please see and carefully read:<br><br><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/conform.shtml" target="_blank">https://www.ncl.ucar.edu/<wbr>Document/Functions/Built-in/<wbr>conform.shtml</a><br><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/conform_dims.shtml" target="_blank">https://www.ncl.ucar.edu/<wbr>Document/Functions/Built-in/<wbr>conform_dims.shtml</a><br><br>---<br></div><div>  DATA_SINGLE_3D = conform(data_subset, data_single, 0)    ; (ntim,nlat,mlon)<br></div><div>  printVarSummary(DATA_SINGLE_<wbr>3D)   ; Look at data<br></div><div><br>Do *not* use (:,:,:) in a computation... this is inefficient.<br><br>  CoD1 = ((DATA_SINGLE_3D - data_subset)/(DATA_SINGLE_3D+ data_subset))^2</div><div>  printVarSummary(CoD1)<br></div><div>  printMinMax(CoD1, 0)<br><br><br></div><div><br><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Nov 7, 2017 at 11:15 AM, Kunal Bali <span dir="ltr"><<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@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">Dear NCL user<div><br></div><div>Suppose, I have 2 datasets </div><div><span id="m_5515606383850329571gmail-24da4d6e-023a-4f1a-b118-d5e1a4c82fd0" class="m_5515606383850329571gmail-GINGER_SOFTWARE_mark">e.g.</span></div><div><br></div><div>data1 = of single point <span id="m_5515606383850329571gmail-6aca75ac-d298-40a8-985f-1cbbb96b5695" class="m_5515606383850329571gmail-GINGER_SOFTWARE_mark">lat</span> & <span id="m_5515606383850329571gmail-bff943b5-43c0-40fa-a663-b799819a52ff" class="m_5515606383850329571gmail-GINGER_SOFTWARE_mark">lon</span> (e.g 28N, 72E) with time steps 365 (having filling values)</div><div>data2 = of grid box (24-32N, 68-76E) with time steps 365  (having filling values)</div><div><br></div><div><span id="m_5515606383850329571gmail-66abf8dd-5002-4f45-affb-3d51dc239b49" class="m_5515606383850329571gmail-GINGER_SOFTWARE_mark">Now</span> I want to calculate the CoD from one single point to 400km grid box area and then lot the spatial map of the calculated product. </div><div><br></div><div><span class="m_5515606383850329571gmail-GINGER_SOFTWARE_mark">Formula</span> for CoD is GIVEN BELOW<img src="cid:ii_j9pxiu9x1_15f97ab70e5f7d73" width="467" height="98"><br>​</div><div><br></div><div>So I calculated in NCL as </div><div><br></div><div>CoD1 = ((data1 - data2<span id="m_5515606383850329571gmail-9c0766d1-b7ae-40ad-a6fd-7d16b19ee10a" class="m_5515606383850329571gmail-GINGER_SOFTWARE_mark">)</span>/(data1 + data2)<span id="m_5515606383850329571gmail-6dab5315-cc3c-41e8-b9d7-53268320e1c9" class="m_5515606383850329571gmail-GINGER_SOFTWARE_mark">)</span>^2</div><div>CoD_final = (sqrt<span id="m_5515606383850329571gmail-4532ee9b-4867-4186-806e-11e562eabf37" class="m_5515606383850329571gmail-GINGER_SOFTWARE_mark">(</span>sum<span id="m_5515606383850329571gmail-9c6d88cf-f4c7-419e-865b-49604afa5b43" class="m_5515606383850329571gmail-GINGER_SOFTWARE_mark">(</span>CoD1))<span id="m_5515606383850329571gmail-de678d04-615b-4668-963c-db988d6ae0e2" class="m_5515606383850329571gmail-GINGER_SOFTWARE_mark">)</span>/3</div><div><br></div><div><span id="m_5515606383850329571gmail-32f19815-56c6-4d4c-b70f-e51d0a0d030e" class="m_5515606383850329571gmail-GINGER_SOFTWARE_mark">please let me know is</span> it the <span id="m_5515606383850329571gmail-798c9baa-ab22-45e1-8ddb-6a1e673aa7fb" class="m_5515606383850329571gmail-GINGER_SOFTWARE_mark">right way</span> to do?</div><div><br></div><div>I tried to plot the final product CoD_final with the attached script, but error comes as</div><div><br></div><div><div>fatal:Minus: Number of dimensions do not match, can't continue</div><div>fatal:["Execute.c":8578]:<wbr>Execute: Error occurred at or near line 116 in file plot.ncl</div></div><div><br></div><div><br></div><div><br clear="all"><div><div class="m_5515606383850329571gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><span id="m_5515606383850329571gmail-f5df93db-c322-49af-98a9-077ea58741ad" class="m_5515606383850329571gmail-GINGER_SOFTWARE_mark">regards</span></div><span class="HOEnZb"><font color="#888888"><div>Kunal Bali<br></div><div><br></div></font></span></div></div></div></div></div>
</div></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></div>