Re: how to calculate transport across the boundary and area integration for an area

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Mon Aug 05 2013 - 13:44:17 MDT

We appreciate it when people show what they have done 1st.
ncl-talk can not be expected to write code.

===

If variables are type "short" use short2flt(...)
===
    f = addfile("...","r")

    latS = -10
    lonN = 20
    lonL = 150
    lonR = 190

    u = f->U(:,:,{latS:latN},{lonL:lonR})
    q = f->Q(:,:,{latS:latN},{lonL:lonR})

    uq = u*q
    copy_VarCoords(u, uq)
    uq@long_name = "zonal moisture transport"
    uq@units = "(m/s)(kg/kg)" or "(m/s)(g/kg)"

same for "v"

The 'box' boundaries are the edges.

===
 From the NCL Wkshop class notes: See Processing: conform

See:
http://www.ncl.ucar.edu/Document/Functions/Contributed/dpres_plevel_Wrap.shtml

T(ntim, klev, nlat,mlon), dp(klev) ; T(0 , 1 , 2 , 3)
dpT = conform (T, dp, 1) ; dpT(ntim,klev,nlat,mlon)
T_wgtAve = dim_sum_n (T*dpT, 1)/dim_sum_n(dpT, 1) ; wgted vertical

On 8/5/13 11:46 AM, Bithi De wrote:
> Hi,
> I was calculating moisture transport using NCEP Reanalysis2 data. I
> need to calculate transport across each of the boundaries of a box
> containing my area of interest.
>
> I also need to calculate area integration for the same box.
> How is it possible to do using NCL?
>
> Thanks,
> Bithi
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Aug 5 13:44:20 2013

This archive was generated by hypermail 2.1.8 : Wed Aug 07 2013 - 13:59:00 MDT