Re: How to converts from pressure levels to sigma coordinates

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Wed, 03 Sep 2008 07:57:30 -0600

The p0 = 1

should be
     p0 = 1000 ; hPa

Dennis Shea wrote:
> 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
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Sep 03 2008 - 07:57:30 MDT

This archive was generated by hypermail 2.2.0 : Wed Sep 03 2008 - 11:24:29 MDT