[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