<div dir="ltr"><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div>Oops!  <br><br>    Integral ===> Summation  (not an average)<br><br></div>It should be:<br><br>    paet_zonal = dim_sum_n_Wrap(paet, 2)   ; (time,lat)<br><br>The multiply by the appropriate latitudinal area: cos(lat)*dlon<br><br></div>    paet_zonal = paet_zonal*a*cos(rad*lat)*dlon<br></div>    paet_zonal@long_name = "....."<br></div>    paet@units = "..."<br><br>----------------<br><br></div>Also, many of the energy quantities are available at:<br><br></div>ftp <a href="http://ftp.cgd.ucar.edu">ftp.cgd.ucar.edu</a><br></div>anonymous<br></div>your_email<br></div>cd archive/BUDGETS<br></div>ls<br></div>cd ERAI<br></div>ls<br></div>get <a href="http://ERAI.TEDIV.1979-2015.nc">ERAI.TEDIV.1979-2015.nc</a><br></div>get ...<br>...<br></div>quit<br><br>------------------<br></div><div><div><div><div><div><br><div><div><div><div><div><div><div><div><div><div><div><br>-----<br><br></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">On Mon, Sep 18, 2017 at 1:59 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div>If you have:<br><br></div>paet(time,lat,lon) ; dimension numbers (0,1,2)<br><br></div>and, you want "to calculate the poleward <span class="m_5565865308374687230gmail-il">energy</span> transports along the latitude circle"<br><br></div>Not clear what you mean but this sounds like the zonal means<br></div><div><br></div>   paet_zonal = dim_avg_n_Wrap(paet, 2)   ; (time,lat)<br><br></div>The multiply by the appropriate latitudinal area.<br><div><div><br><br> </div></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Mon, Sep 4, 2017 at 11:33 PM, Sourav Chatterjee <span dir="ltr"><<a href="mailto:srvsxc@outlook.com" target="_blank">srvsxc@outlook.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="h5">




<div dir="ltr">
<div id="m_5565865308374687230m_-5893429167856747811divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<p>Hello NCL Users</p>
<p><br>
</p>
<p>I want to calculate the poleward energy transports along the latitude circles which is given by,</p>
<p><br>
</p>
<p><a name="m_5565865308374687230_m_-5893429167856747811_Sec11"></a></p>
<p><i><span style="color:rgb(0,0,0)"></span>PAET</i>(<img src="http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.phi.gif">) =
<img src="http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.int02pi.gif">[
<i>C_p</i><font color="red"><span style="color:rgb(0,0,0)">VT</span></font><span style="color:rgb(0,0,0)">(</span><img src="http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.lambda.gif">,
<img src="http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.phi.gif">) + <i>
g</i><font color="red"><span style="color:rgb(0,0,0)">VZ</span></font><span style="color:rgb(0,0,0)">(</span><img src="http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.lambda.gif"><span style="color:rgb(0,0,0)">,
</span><img src="http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.phi.gif"><span style="color:rgb(0,0,0)">) +
</span><i><span style="color:rgb(0,0,0)">L</span></i><font color="red"><span style="color:rgb(0,0,0)">VQ</span></font><span style="color:rgb(0,0,0)">(</span><img src="http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.lambda.gif"><span style="color:rgb(0,0,0)">,
</span><img src="http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.phi.gif"><span style="color:rgb(0,0,0)">) +
</span><font color="red"><span style="color:rgb(0,0,0)">VK</span></font><span style="color:rgb(0,0,0)">(</span><img src="http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.lambda.gif"><span style="color:rgb(0,0,0)">,
</span><img src="http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.phi.gif"><span style="color:rgb(0,0,0)">)</span>]
<i>a</i> cos(<img src="http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.phi.gif">) d<img src="http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.lambda.gif">
<br>
</p>
<p><a href="http://www.cgd.ucar.edu/cas/catalog/newbudgets/" class="m_5565865308374687230m_-5893429167856747811OWAAutoLink" id="m_5565865308374687230m_-5893429167856747811LPlnk976476" target="_blank">http://www.cgd.ucar.edu/cas/ca<wbr>talog/newbudgets/</a></p>
<div id="m_5565865308374687230m_-5893429167856747811LPBorder_GT_15045877138990.4027121073145895" style="margin-bottom:20px;overflow:auto;width:100%;text-indent:0px">
<table id="m_5565865308374687230m_-5893429167856747811LPContainer_15045877138920.9649448193477005" style="width:90%;background-color:rgb(255,255,255);overflow:auto;padding-top:20px;padding-bottom:20px;margin-top:20px;border-top:1px dotted rgb(200,200,200);border-bottom:1px dotted rgb(200,200,200)" cellspacing="0">
<tbody>
<tr style="border-spacing:0px" valign="top">
<td id="m_5565865308374687230m_-5893429167856747811TextCell_15045877138940.6346520655902038" style="vertical-align:top;padding:0px;display:table-cell" colspan="2">
<div id="m_5565865308374687230m_-5893429167856747811LPRemovePreviewContainer_15045877138940.5495064427456213"></div>
<div id="m_5565865308374687230m_-5893429167856747811LPTitle_15045877138940.4038262243975589" style="color:rgb(0,120,215);font-weight:400;font-size:21px;font-family:"wf_segoe-ui_light","Segoe UI Light","Segoe WP Light","Segoe UI","Segoe WP",Tahoma,Arial,sans-serif;line-height:21px">
<a id="m_5565865308374687230m_-5893429167856747811LPUrlAnchor_15045877138960.6416793234186516" style="text-decoration:none" href="http://www.cgd.ucar.edu/cas/catalog/newbudgets/" target="_blank">Vertically Integrated Mass, Moisture, Heat, and Energy ...</a></div>
<div id="m_5565865308374687230m_-5893429167856747811LPMetadata_15045877138960.3301664613330533" style="margin:10px 0px 16px;color:rgb(102,102,102);font-weight:400;font-family:"wf_segoe-ui_normal","Segoe UI","Segoe WP",Tahoma,Arial,sans-serif;font-size:14px;line-height:14px">
<a href="http://www.cgd.ucar.edu" target="_blank">www.cgd.ucar.edu</a></div>
<div id="m_5565865308374687230m_-5893429167856747811LPDescription_15045877138980.7495977679656846" style="display:block;color:rgb(102,102,102);font-weight:400;font-family:"wf_segoe-ui_normal","Segoe UI","Segoe WP",Tahoma,Arial,sans-serif;font-size:14px;line-height:20px;max-height:100px;overflow:hidden">
Vertically Integrated Mass, Moisture, Heat, and Energy Budget Products Derived from the NCEP/NCAR Reanalysis March 2003 Climate Analysis Section, CGD, NCAR</div>
</td>
</tr>
</tbody>
</table>
</div>
I have the mass corrected values for VT,VZ,VQ,VK (time*lat*lon*) on gaussian lat lon grid.
<p></p>
<p></p>
<p><br>
</p>
<p>I want to retain the time dimension.</p>
<p><br>
</p>
<p><br>
</p>
<p></p>
<div>l = 2.5104*10^6<br>
a = 6.37122*10^6<br>
cp = 1004.64<br>
g = 9.81<br>
<br>
 pi = 3.14159265<br>
 re = 6.371*10^6            ; radius of earth<br>
 coef = re/1.e15       ; in PW<br>
<br>
<br>
ntime = dimsizes(vk(:,0,0)) ; 444<br>
nlat = dimsizes(vk(0,:,0)) ; 256<br>
nlon = dimsizes(vk(0,0,:)) ; 512<br>
dlon = 2.*pi/nlon   ; dlon in radians<br>
lat = vt&lat<br>
time = vt&time<br>
<br>
lat&lat = lat<br>
<br>
gw =  latGauWgt(nlat, "lat", "gaussian weights", "")<br>
<br>
<br>
energy = cp*vt+g*vz+l*vq+vk<br>
copy_VarCoords(vt,energy)<br>
printVarSummary(netflux)<br>
<br>
</div>
<p></p>
<p>I have tried the example for ocean given in <a href="http://www.mit.edu/course/13/13.715/ncarg-5.5.1/lib/ncarg/nclscripts/csm/ccm_func.ncl" class="m_5565865308374687230m_-5893429167856747811OWAAutoLink" id="m_5565865308374687230m_-5893429167856747811LPlnk18119" target="_blank">
http://www.mit.edu/course/13/1<wbr>3.715/ncarg-5.5.1/lib/ncarg/nc<wbr>lscripts/csm/ccm_func.ncl</a>.</p>
<p><br>
</p>
<p>But I am not getting correct values with that. Can anyone help on this!<br>
</p>
<p><br>
</p>
<p><br>
</p>
<p>Thanks</p><span class="m_5565865308374687230HOEnZb"><font color="#888888">
<p>Sourav<br>
</p>
<p><br>
</p>
<p><br>
</p>
<p><br>
</p>
<p><br>
</p>
<p><br>
</p>
</font></span></div>
</div>

<br></div></div><span class="">______________________________<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></span></blockquote></div><br></div>
</blockquote></div><br></div>