Re: confusion about the Argument of cz2ccm

Dennis,

Thanks!
Yes, the size of "bnds" is 2.
The a_bnds are
  0.000 0.000
  0.000 0.5038E-02
  0.5038E-02 0.1114E-01
...
  0.5104E-04 0.2662E-04
  0.2662E-04 0.9869E-05
a_bnds(k,0) and a_bnds(k,1) are the lower and upper bounds of level k.

On Tue, Feb 4, 2014 at 3:25 PM, Dennis Shea <shea@ucar.edu> wrote:

> [1] The documentation for hyai and hybi has been changed.
> Thank you for calling this our attention.
>
> [2] You have
>
>
> > double lev ( lev )
> > bounds : lev_bnds
> [snip]
>
> Then you have
>
>
> double a_bnds ( lev, bnds )
> long_name : vertical coordinate formula term: a(k+1/2)
>
> You do not indicate the size 'lev' and there is no information
> on the named dimension 'bnds'.
>
> **IF** the size of 'bnds' is 2, then maybe
>
> a_bnds(:,0) for hyai and b_bnds(:,1) for hybi?
>
> This assume the the 2 refers to "upper" and "lower"
>
> a_bnds(k,0)
> ai(k)
> a_bnds(k,1)
>
> I would print these out to verify
>
>
>
>
> On 2/4/14, 12:50 PM, Andy Show wrote:
>
>> Dear NCL users,
>>
>> http://www.ncl.ucar.edu/Document/Functions/Built-in/cz2ccm.shtml
>> function cz2ccm (
>> ps : numeric,
>> phis : numeric,
>> tv : numeric,
>> p0 [1] : numeric,
>> hyam [*] : numeric,
>> hybm [*] : numeric,
>> hyai [*] : numeric,
>> hybi [*] : numeric
>> )
>> ...
>> hyam
>> A one-dimensional array of hybrid coordinate coefficients for base
>> pressures (layer midpoints) ordered bottom-to-top.
>> hybm
>> A one-dimensional array of hybrid coordinate coefficients for surface
>> pressures (layer midpoints) ordered bottom-to-top.
>> hyai
>> A one-dimensional array of hybrid A coefficients (layer midpoints) ordered
>> bottom-to-top.
>> hybi
>> A one-dimensional array of hybrid B coefficients (layer midpoints) ordered
>> bottom-to-top.
>> ...
>> 1) hyai and hybi should be the coefficients on the layer interface or
>> midpoint?
>>
>> 2) In CMIP5 6hr outputs (eg. GFDL-CM3), there are two pairs of a&b and
>> a_bnds&b_bnds coefficients in the hybrid coordinate. Then how to use
>> cz2ccm
>> to calculate the geopotential height, a_bnds(:,0) for hyai and b_bnds(:,0)
>> for hybi?
>>
>> double lev ( lev )
>> bounds : lev_bnds
>> units : 1
>> axis : Z
>> positive : down
>> long_name : hybrid sigma pressure coordinate
>> standard_name : atmosphere_hybrid_sigma_
>> pressure_coordinate
>> formula : p = a*p0 + b*ps
>> formula_terms : p0: p0 a: a b: b ps: ps
>> double a ( lev )
>> long_name : vertical coordinate formula term: a(k)
>> double b ( lev )
>> long_name : vertical coordinate formula term: b(k)
>> double a_bnds ( lev, bnds )
>> long_name : vertical coordinate formula term: a(k+1/2)
>> double b_bnds ( lev, bnds )
>> long_name : vertical coordinate formula term: b(k+1/2)
>>
>> Thanks!
>>
>>
>>
>> _______________________________________________
>> 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 Tue Feb 4 13:38:41 2014