[ncl-talk] function vinth2p

[Dennis Shea]

I would not use 'vinth2p' with 'fake' vertical vertical hybrid 'A'
coefficients.
'vinth2p' was developed for the CAM hybrid model.

  *p(k) = a(k)*p0 + b(k)*ps*

Since, 'ps' (surface pressure) is variable the 'top' calculated pressure
would also be variable.

Attached is a script that, I think, does what you want.
===

NOTE: Automatic loading of NCL library scripts:

>From 6.2.0 onward, you no longer need to include the following load
commands at the top of your NCL script:

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

>From 6.4.0 onwards, the following scripts are automaticall loaded:

    "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
    "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
    "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
    "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
    "$NCARG_ROOT/lib/ncarg/nclscripts/csm/bootstrap.ncl"
    "$NCARG_ROOT/lib/ncarg/nclscripts/csm/extval.ncl"
    "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
'


On Tue, Apr 25, 2017 at 7:04 PM, Raj K. Rai <rajkrai at gmail.com> wrote:

> Hi,
> I am using the function vinth2p(data,fake,sigma,pres,psfc,1,1000,1,True)
> to convert the WRF data that are in sigma level to pressure level (from
> 1000 mbar to 100 mbar). I guess this function will take default value for
> the pressure at the top (if so what is its value?) and bottom pressure from
> psfc. How can I use 20 mbar for the top pressure that corresponds to the
> top sigma level?
>
> Thanks for your help.
> Raj
>
> _______________________________________________
> 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/20170426/5ffd82d4/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: wrf_vint.ncl
Type: text/x-ncl
Size: 2491 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170426/5ffd82d4/attachment.bin