Re: vertical integration

http://www.ncl.ucar.edu/Document/Functions/Built-in/conform.shtml

qu [t|1] x [p|17] x [latitude|72] x [longitude | 96]
    0 1 2 f

qu_dp = qu*conform(qu, dp, 1) ; 'dp' matches the 'p' dimension

Good Luck

> Hello ncl-talk,
>
> I hope this is an easy question. I am trying to do a vertical
> integration following the examples listed under dpres_plevel and
> dpres_hybrid_ccm (I have datasets with both types of pressure units).
>
> One set has pre-calculated monthly qu and qv for jja. I can't get
> past multiplying the variables by dp. I was trying to follow Example
> 3; http://www.ncl.ucar.edu/Document/Functions/Contributed/
> dpres_plevel_Wrap.shtml.
>
> This is the code I am trying to use:
> qu = a1->uq
> qv = a2->vq
> lev = qu&p
>
> printVarSummary(qu)
>
> psfc = 1018
> ptop = 10
>
> ; dp(ntim,klev,nlat,mlon)
> dp = dpres_plevel(lev, psfc, ptop, 0)
> printVarSummary(dp)
> g = 9.81
>
> qu_dp = qu*dp **This is where I get the error message.
> qv_dp = qv*dp
>
>
> Here is the ncl printout for qu.
> Variable: qu
> Type: float
> Total Size: 470016 bytes
> 117504 values
> Number of Dimensions: 4
> Dimensions and sizes: [t | 1] x [p | 17] x [latitude | 72] x
> [longitude | 96]
> Coordinates:
> t: [39825..39825]
> p: [10..1000]
> latitude: [88.75..-88.75]
> longitude: [1.875..358.125]
> Number Of Attributes: 11
> source : Unified Model Output:
> name : uq
> title : u*q
> date : 01/12/59
> time : 00:00
> long_name : u*q
> units : m kg s-1 kg-1
> missing_value : 2e+20
> _FillValue : 2e+20
> valid_min : -0.1634946
> valid_max : 0.273821
>
>
> Variable: dp
> Type: float
> Total Size: 68 bytes
> 17 values
> Number of Dimensions: 1
> Dimensions and sizes: [17]
> Coordinates:
> fatal:Mul: Number of dimensions do not match, can't continue
>
>
> I see that the problem with multiplying by dp has to do with the
> number of dimensions, so is there something missing from my code (and
> maybe the examples) that I can try?
>
> Thank you for your help,
>
> Jeanne
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sat Jan 28 23:50:51 2012