Re: How to converts from pressure levels to sigma coordinates

I think you should be able to use

http://www.ncl.ucar.edu/Document/Functions/Built-in/pres2hybrid.shtml

    xp = ; variable to be interpolated
    psfc = ; surface pressure (Pa) or (hPa)

    hyb = sigma ; sigma coefficients
    hya = h ; create array same size/type as hya
    hya = 0.0 ; set to zero

    p0 = 100000 ; reference pressure (Pa)
or
    p0 = 1 ; hPa

    plev = ; match units of p0 and psfc
    plev_at_units = "Pa" or "hPa"

    x_sigma = pres2hybrid_Wrap(plev,psfc,p0,xp,hya,hyb,0)
    printVarSummary( x_sigma)

If retention of meta data is desired, use

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

x_sigma = pres2hybrid_Wrap(plev,psfc,p0,xp,hya,hyb,0)
printVarSummary(x_sigma)

Good luck
Lawrence wrote:
> Hi all,
>
> Is there any fuction in NCL that converts from pressure levels to sigma
> coordinates? Thank you !
>
> Regards,
>
> Lawrence
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Sep 03 2008 - 07:43:32 MDT