Re: question about dpres_hybrid_ccm

     lev = 30 ; mid-layers
     ilev = 31 ; interface layers
     double hyai(ilev) ;
     double hybi(ilev) ;
     double hyam(lev) ;
     double hybm(lev) ;

     float PS(time, lat, lon)
     float U(time, lev, lat, lon)
     float Q(time, lev, lat, lon) ;

Are you using hyai and hybi (the interface layers)
or hyam, hybm? You should be using hyai and hybi

  f = addfile("","r")
  p0 = f->P0
  ps = f->PS
  hyai = f->hyai
  hybi = f->hybi
  dp = dpres_hybrid_ccm (ps,p0,hyai,hybi)
  printVarSummary(dp)

On 2/27/12 2:41 PM, Jeanne Thibeault wrote:
> Hello NCL users,
>
> I am trying to use dpres_hybrid_ccm in a vertical integration
> calculation (following Example 3 at
> http://www.ncl.ucar.edu/Document/Functions/Built-in/dpres_hybrid_ccm.shtml)
> and am running into a problem with the dimensions when I try to multiply
> qu by dp. I see in the example that the number of levels in dp is lev-1,
> but it's not clear how to deal with the mismatch in levels between the
> original variable and dp. Any help is appreciated.
>
> Thank you,
> Jeanne
>
> I have monthly qu that I calculated from 6hrly data:
>
> Variable: qu
> Type: float
> Total Size: 15667200 bytes
> 3916800 values
> Number of Dimensions: 4
> Dimensions and sizes: [time | 60] x [lev | 24] x [lat | 40] x [lon | 68]
> Coordinates:
> time: [40150.25..41945.25]
> lev: [0.9961083927371855..0.003650292822203918]
> lat: [1.011235955056179..79.88764044943821]
> lon: [181.25..348.75]
>
>
> I also have monthly ps that I calculated using the 6 hrly data:
> Variable: ps
> Type: float
> Total Size: 652800 bytes
> 163200 values
> Number of Dimensions: 3
> Dimensions and sizes: [time | 60] x [lat | 40] x [lon | 68]
> Coordinates:
> time: [40150.25..41945.25]
> lat: [1.011235955056179..79.88764044943821]
> lon: [181.25..348.75]
> Number Of Attributes: 5
> average_op_ncl : dim_avg_n over dimension(s): time
> standard_name : surface_air_pressure
> cell_methods : time: point
> units : Pa
> long_name : Surface Air Pressure
>
> I use the following line to calculate dp, but there is a mismatch
> between the number of levels in qu and dp.
> dp = dpres_hybrid_ccm (ps,p0,hyai,hybi) ; Pa [kg/(m s2)]
>
> If I try to multiply qu*dp (as in the example) I get the following error
> fatal:Mul: Dimension size, for dimension number 1, of operands does not
> match, can't continue
>
>
>
> Variable: dp
> Type: float
> Total Size: 15014400 bytes
> 3753600 values
> Number of Dimensions: 4
> Dimensions and sizes: [60] x [23] x [40] x [68]
> Coordinates:
> (0) ERROR: copy_VarCoords: dimension sizes do not match
>
>
>
>
> _______________________________________________
> 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 Mon Feb 27 19:54:07 2012