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