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