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

[Zilore Mumba]

Thanks a lot Ehsan, I will try it and let you know..

<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>

On Sat, Feb 19, 2022 at 8:01 PM Ehsan Taghizadeh <ehsantaghizadeh at yahoo.com>
wrote:

> Dear Zilore,
> I hope attached file could help you.
>
>
> *-----------------------------------------------------------------------------*
> *Sincerely,*
> *Ehsan Taghizadeh*
>
>
>
>
> On Saturday, February 19, 2022, 01:13:15 AM GMT+3:30, Zilore Mumba via
> ncl-talk <ncl-talk at mailman.ucar.edu> wrote:
>
>
> 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>
> <#m_4062106394686179226_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220219/4db94bd5/attachment.html>