[ncl-talk] how to do vertical interpolation from hybrid to sigma levels

[Dennis Shea]

There is no explicit function to do this. You must do multiple steps.

See: *Model Vertical Coordinate*
<https://www.cesm.ucar.edu/models/atm-cam/docs/usersguide/node25.html>
==> p(k) = a(k)*p0 + b(k)*ps

You can see the values of the CAM hybrid coefficients [top-to-bottom]  via

%> *ncdump* -v hyam,hybm CAM.nc

The *hybm *are the 'sigma' level portion of the hybrid (pressure-sigma)
coefficients.
The hybm are not the classic standalone:  sigma_classic=p/psfc   or
p/p0

[1] calculate pressures at CAM5 levels:  *pres_hybrid_ccm*
<https://www.ncl.ucar.edu/Document/Functions/Built-in/pres_hybrid_ccm.shtml>
: p_cam = *pres_hybrid_ccm*(....)  ; (time,klvl,nlat,mlon)
[2] read a variable you wish to interpolate from hybrid to sigma:   eg
t_cam = fcam->T             ; (time,klvl,nlat,mlon)
[3] sigma_cam = p_cam/p_sfc    or   p_cam/p0
[4] Use *int2p*
<https://www.ncl.ucar.edu/Document/Functions/Built-in/int2p.shtml> or
*int2p_Wrap*
<https://www.ncl.ucar.edu/Document/Functions/Contributed/int2p_Wrap.shtml>
to interpolate to user desired sigma levels ( SIGMA(KSIGMA) ).
The documentation and function were originally developed using pressure.
Hence, the use of pressure in the documentation. However, as I recall,
*int2p* actually is generic.

     SIGMA = (/....../)     ; (KSIGMA)    <== user specified SIGMA levels
     SIGMA!0 = "SIGMA"
     SIGMA at long_name = "sigma level: sigma=p/psfc"    ; "sigma level:
sigma=p/p0"
     SIGMA at units = ""
;  unitless

     t_SIGMA = *int2p_Wrap*
<https://www.ncl.ucar.edu/Document/Functions/Contributed/int2p_Wrap.shtml>
(sigma_cam ,t_cam, SIGMA,  *1* )     ; (ntim,KSIGMA,nlat,mlon)
     *printVarSummary*(t_SIGMA)


On Sun, Jan 24, 2021 at 6:25 PM Hui Ding via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Dear Sir or Madam,
> I am trying to interpolate monthly CAM5 diabatic data from hybrid level to
> sigma level? Do you know how to do it? Thank you!
> Best regards,
> Hui Ding
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210125/8220f625/attachment.html>