<div dir="ltr">Thanks a lot Ehsan, I will try it and let you know.. <br></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><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Feb 19, 2022 at 8:01 PM Ehsan Taghizadeh <<a href="mailto:ehsantaghizadeh@yahoo.com">ehsantaghizadeh@yahoo.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div style="font-family:"bookman old style","new york",times,serif;font-size:16px"><div><div dir="ltr">Dear Zilore,</div><div dir="ltr">I hope attached file could help you.</div><div><br></div><div><div dir="ltr" style="font-family:"old times",serif;font-size:16px"><div><div dir="ltr" style="color:rgb(0,0,0);font-family:"old times",serif;font-size:16px"><i><font size="2">-----------------------------------------------------------------------------</font></i></div><div dir="ltr" style="color:rgb(0,0,0);font-family:"old times",serif;font-size:16px"><i><font size="2">Sincerely,</font></i></div><div dir="ltr" style="color:rgb(0,0,0);font-family:"old times",serif;font-size:16px"><i><font size="2">Ehsan Taghizadeh</font></i></div><div dir="ltr" style="color:rgb(0,0,0);font-family:"old times",serif;font-size:16px"><br></div></div><br></div></div></div>
        <div><br></div><div><br></div>
        
        </div><div id="gmail-m_4062106394686179226yahoo_quoted_5436635089">
            <div style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px;color:rgb(38,40,42)">
                
                <div>
                    On Saturday, February 19, 2022, 01:13:15 AM GMT+3:30, Zilore Mumba via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:
                </div>
                <div><br></div>
                <div><br></div>
                <div><div id="gmail-m_4062106394686179226yiv1044055519"><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:New;font-size:large">T(850 mb) + Td(850 mb) - 2[T(500 mb)].</b><div><font color="#000000" face="arial, sans-serif"><b>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><font color="#000000" 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="gmail-m_4062106394686179226yiv1044055519DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2"><br>
<table style="border-top:1px solid rgb(211,212,222)">
        <tbody><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" rel="noreferrer noopener" 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; min-height: 29px;"></a></td>
                <td style="width:470px;padding-top:12px;color:rgb(65,66,78);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" style="color:rgb(68,83,234)" rel="noreferrer noopener" target="_blank">www.avast.com</a>
                </td>
        </tr>
</tbody></table><a href="#m_4062106394686179226_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2" rel="noreferrer noopener"></a></div>
</div>_______________________________________________<br>ncl-talk mailing list<br><a href="mailto:ncl-talk@mailman.ucar.edu" rel="noreferrer noopener" target="_blank">ncl-talk@mailman.ucar.edu</a><br>List instructions, subscriber options, unsubscribe:<br><a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer noopener" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></div>
            </div>
        </div></div></blockquote></div>