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

[Dennis Shea]

Oops!

    Integral ===> Summation  (not an average)

It should be:

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

The multiply by the appropriate latitudinal area: cos(lat)*dlon

    paet_zonal = paet_zonal*a*cos(rad*lat)*dlon
    paet_zonal at long_name = "....."
    paet at units = "..."

----------------

Also, many of the energy quantities are available at:

ftp ftp.cgd.ucar.edu
anonymous
your_email
cd archive/BUDGETS
ls
cd ERAI
ls
get ERAI.TEDIV.1979-2015.nc
get ...
...
quit

------------------


-----


On Mon, Sep 18, 2017 at 1:59 PM, Dennis Shea <shea at ucar.edu> wrote:

> 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/nc
>> lscripts/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/fb56077c/attachment.html>