[ncl-talk] Poleward Atmospheric Energy Transport-Integration along longitudes

Dennis Shea shea at ucar.edu
Mon Sep 18 13:59:01 MDT 2017

```If you have:

paet(time,lat,lon) ; dimension numbers (0,1,2)

and, you want "to calculate the poleward energy transports along the
latitude circle"

Not clear what you mean but this sounds like the zonal means

paet_zonal = dim_avg_n_Wrap(paet, 2)   ; (time,lat)

The multiply by the appropriate latitudinal area.

On Mon, Sep 4, 2017 at 11:33 PM, Sourav Chatterjee <srvsxc at outlook.com>
wrote:

> Hello NCL Users
>
>
> I want to calculate the poleward energy transports along the latitude
> circles which is given by,
>
>
> *PAET*() = [ *Cp*VT(, ) + * g*VZ(, ) + *L*VQ(, ) + VK(, )] *a* cos() d
>
> http://www.cgd.ucar.edu/cas/catalog/newbudgets/
> Vertically Integrated Mass, Moisture, Heat, and Energy ...
> <http://www.cgd.ucar.edu/cas/catalog/newbudgets/>
> www.cgd.ucar.edu
> Vertically Integrated Mass, Moisture, Heat, and Energy Budget Products
> Derived from the NCEP/NCAR Reanalysis March 2003 Climate Analysis Section,
> CGD, NCAR
> I have the mass corrected values for VT,VZ,VQ,VK (time*lat*lon*) on
> gaussian lat lon grid.
>
>
> I want to retain the time dimension.
>
>
>
> l = 2.5104*10^6
> a = 6.37122*10^6
> cp = 1004.64
> g = 9.81
>
>  pi = 3.14159265
>  re = 6.371*10^6            ; radius of earth
>  coef = re/1.e15       ; in PW
>
>
> ntime = dimsizes(vk(:,0,0)) ; 444
> nlat = dimsizes(vk(0,:,0)) ; 256
> nlon = dimsizes(vk(0,0,:)) ; 512
> dlon = 2.*pi/nlon   ; dlon in radians
> lat = vt&lat
> time = vt&time
>
> lat&lat = lat
>
> gw =  latGauWgt(nlat, "lat", "gaussian weights", "")
>
>
> energy = cp*vt+g*vz+l*vq+vk
> copy_VarCoords(vt,energy)
> printVarSummary(netflux)
>
> I have tried the example for ocean given in http://www.mit.edu/course/13/
> 13.715/ncarg-5.5.1/lib/ncarg/nclscripts/csm/ccm_func.ncl.
>
>
> But I am not getting correct values with that. Can anyone help on this!
>
>
>
> Thanks
>
> Sourav
>
>
>
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170918/1cde9c93/attachment.html>
```