[ncl-talk] converting pressure levels to hybrid height coordiantes ( output of HadGEM2-ES).

[Dennis Shea]

>From the file and dump you sent offline:

double lev(lev) ;
lev:bounds = "lev_bnds" ;
lev:units = "m" ;
lev:axis = "Z" ;
lev:positive = "up" ;
lev:long_name = "hybrid height coordinate" ;
lev:standard_name = "atmosphere_hybrid_height_coordinate" ;
lev:formula = "z(k,j,i) = a(k) + b(k)*orog(j,i)" ;
lev:formula_terms = "a: lev b: b orog: orog" ;

The "b" coefficient variable  is present

double b(lev) ;
b:long_name = "vertical coordinate formula term: b(k)" ;

b = 0.99771648645401, 0.990881502628326, 0.979542553424835,
0.9637770652771,
    0.943695485591888, 0.919438362121582, 0.891178011894226,
    0.859118342399597, 0.823493480682373, 0.784570515155792,
    0.742646217346191, 0.698050200939178, 0.651142716407776,
    0.602314412593842, 0.55198872089386, 0.500619947910309,
0.44869339466095,
    0.39672577381134, 0.34526526927948, 0.294891387224197,
0.24621507525444,
    0.199878215789795, 0.156554222106934, 0.116947874426842,
    0.0817952379584312, 0.0518637150526047, 0.0279368180781603,
    0.0107164792716503, 0.00130179093685001, 0, 0, 0, 0, 0, 0, 0, 0, 0 ;


However, there is no "a" coefficient variable present (nor is there an
"a_bnds" coefficient variable).

If you send another file, please only send one time-step containing the 'a'
and 'be' coefficients.
... that is all that is needed

Also, what might a sample pressure level variable look like"

Is it isobaric levels p(lev) or a 4D pressure variable: p(time,lev,lat,lon)

If the former [ p(lev] ] just include in an email; if 4D please send as a
netCDF file.

THX


On Mon, Jun 15, 2015 at 6:16 AM, chao.tang <chao.tang at univ-reunion.fr>
wrote:

> Hello everyone,
>
>
> I was searching for the solution to converting pressure levels to hybrid
> height coordinates (unit: m, output of HadGEM2-ES).
>
> Any idea would be really appreciated.
>
> thanks in advance, best wishes.
>
>
>
> *Best regards,*
> *Chao TANG*
> *--------------------------------------*
> *--------------------------------------**---------------------*
> *Laboratoire d'Energétique, d'Electronique et Procédés*
> *Unité de Recherche EA 4079*
> Université de La Réunion, LE2P
> 15 avenue René Cassin CS 92003
> 97744 ST DENIS CEDEX 9
> Mobile: 0693 91 57 38 <secretariat.le2p at univ-reunion.fr>
> Courriel: chao.tang at univ-reunion.fr  <secretariat.le2p at univ-reunion.fr>
> *--------------------------------------*
> *--------------------------------------**---------------------*
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150615/74a64a10/attachment.html