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

[Sourav Chatterjee]

Hello NCL Users


I want to calculate the poleward energy transports along the latitude circles which is given by,


PAET([http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.phi.gif]) = [http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.int02pi.gif] [ CpVT([http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.lambda.gif], [http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.phi.gif] ) + gVZ([http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.lambda.gif], [http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.phi.gif] ) + LVQ([http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.lambda.gif], [http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.phi.gif] ) + VK([http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.lambda.gif], [http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.phi.gif] )] a cos([http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.phi.gif]) d[http://www.cgd.ucar.edu/cas/catalog/newbudgets/newbudgets.lambda.gif]

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





-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170905/a26847b6/attachment.html>