[ncl-talk] Calculating Dew point temperature at given level in a WRF file.

[Zilore Mumba]

Can someone help me with a simple script to calculate dew point temperature
(td) at say 500hPa on a WRF file. Specifically I am trying to calculate
Total total index which is given as TTI= *T(850 mb) + Td(850 mb) - 2[T(500
mb)].*
*To start with I have the script below which is supposed to give me td.
Somehow I am not able to go beyond this. print(td) prints what looks like
surface dew points at all grid points. Whatever number I put instead of 6
in my script, the temperatures printed seem to be the same.*
*Alternatively if someone has a script which calculates TTI, I will
appreciate it if they can share.*

begin
;---Open WRF output file.
  dir      = "/home/zmumba/DA/OUTPUT/2021081900/noda/"
  filename = "wrfout_d01_2021-08-19_00:00:00"
  f = addfile(dir + filename + ".nc","r")
  QVAPOR = f->QVAPOR(0,6,:,:)
    P      = f->P(0,6,:,:)
    PB     = f->PB(0,6,:,:)
    p      = P + PB
    td    = wrf_td( p , QVAPOR )
  printVarSummary(td)
  exit
  print(td)
end

<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>
Virus-free.
www.avast.com
<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>
<#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220218/5ec1698e/attachment.html>