<div dir="ltr">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=

<b style="color:rgb(128,0,128);font-family:"Times New Roman";font-size:large">T(850 mb) + Td(850 mb) - 2[T(500 mb)].</b><div><font color="#000000" face="arial, sans-serif"><b style="">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.</b></font></div><div><b style=""><font color="#000000" style="" face="arial, sans-serif">Alternatively if someone has a script which calculates TTI, I will appreciate it if they can share.</font></b></div><div><font color="#800080" face="Times New Roman" size="4"><b><br></b></font><div>begin<br>;---Open WRF output file.<br>  dir      = "/home/zmumba/DA/OUTPUT/2021081900/noda/"<br>  filename = "wrfout_d01_2021-08-19_00:00:00"<br>  f = addfile(dir + filename + ".nc","r")<br>  QVAPOR = f->QVAPOR(0,6,:,:)<br>    P      = f->P(0,6,:,:)<br>    PB     = f->PB(0,6,:,:)<br>    p      = P + PB<br>    td    = wrf_td( p , QVAPOR )<br>  printVarSummary(td)<br>  exit<br>  print(td)<br>end<br></div></div></div><div id="DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2"><br>
<table style="border-top:1px solid #d3d4de">
        <tr>
        <td style="width:55px;padding-top:13px"><a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail" target="_blank"><img src="https://ipmcdn.avast.com/images/icons/icon-envelope-tick-round-orange-animated-no-repeat-v1.gif" alt="" width="46" height="29" style="width: 46px; height: 29px;"></a></td>
                <td style="width:470px;padding-top:12px;color:#41424e;font-size:13px;font-family:Arial,Helvetica,sans-serif;line-height:18px">Virus-free. <a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail" target="_blank" style="color:#4453ea">www.avast.com</a>
                </td>
        </tr>
</table><a href="#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2" width="1" height="1"></a></div>