vertical integration

From: Jeanne Thibeault <jeanne.thibeault_at_nyahnyahspammersnyahnyah>
Date: Sat Jan 28 2012 - 06:47:42 MST

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
Received on Sat Jan 28 18:47:54 2012

This archive was generated by hypermail 2.1.8 : Thu Feb 02 2012 - 03:10:31 MST