[ncl-talk] Interpolate WRF eta levels to pressure levels
Dennis Shea
shea at ucar.edu
Tue Jan 7 21:13:57 MST 2020
***If*** the WRF file has a variable
P(Time,bottom_top,south_north,west_east)
then those are the pressure levels at each 'eta' coordinate at each grid
point and time.
If say you want to interpolate the variable
T(Time,bottom_top,south_north,west_east)
then something like:
T = f->T
P = f->P
printVarSummary(T)
printMinMax(T,0)
print("====")
printVarSummary(P)
printMinMax(P,0)
print("====")
po = (/1000,900,850,700,600,500,400,300,200,100,50,10,5/)*100.0
po!0 = "po"
po at long_name = "Pressure Level"
po at units = "Pa"
To =* int2p_n_Wrap*
<http://www.ncl.ucar.edu/Document/Functions/Contributed/int2p_n_Wrap.shtml>
(P,T,po,1,1)
printVarSummary(To)
printMinMax(To,0)
print("====")
====
If P is not on the WRF file, then you will have to get help from WRF
support.
On Tue, Jan 7, 2020 at 9:12 AM Anahita Amiri Farahani via ncl-talk <
ncl-talk at ucar.edu> wrote:
> Hi all,
>
> Is there any function/command in NCL to convert WRF eta levels to pressure
> levels?
>
> Thanks,
> Ana
> _______________________________________________
> 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/20200107/c2ea0646/attachment.html>
More information about the ncl-talk
mailing list