[ncl-talk] wrf_dbz

[Dennis Shea]

What does 'press1' look like?

*printVarSummary*
<https://www.ncl.ucar.edu/Document/Functions/Built-in/printVarSummary.shtml>
(press1)
-------------
Note: you can have the following WRF function calculate 'dbz' directly.

 dbz = *wrf_user_getvar*
<https://www.ncl.ucar.edu/Document/Functions/Built-in/wrf_dbz.shtml>(file1,"dbz",-1)

 printVarSummary(dbz)


On Thu, Dec 24, 2020 at 9:18 AM David W Reynolds via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> When attempting to run my ncl script to execute wrf_dbz using this
> function I get the following error.
>
>
>
> fatal:wrf_dbz: The prs array must have at least 3 dimensions
>
> fatal:["Execute.c":8575]:Execute: Error occurred at or near line 266 in
> file MAIN_plot_Reflectivity_Txt.ncl
>
>
>
> Line 266 is the line trying to run the function.  Can anyone help with
> this error.  I am very new to ncl so apologize for my ignorance.
>
>
>
> 263     ; Calculate Reflectivity
>
>     264     ivarint = 0  ; or 1
>
>     265     iliqskin = 1  ; or 1
>
>     266     dbz = wrf_dbz( press1, tk1, qvapor1, qrain, qsnow, qgraup,
> ivarint, iliqskin)
>
>
>
> Here is how I read the data.  Trying to output a vertical profile of
> reflectivity over a given lat long.
>
>
>
> ; Get fields from file
>
>     pbtemp1 = file1->$PBTEMPNAME1$(TIME1,:,i1,j1)
>
>     qvapor1 = file1->$QVAPORNAME1$(TIME1,:,i1,j1)
>
>     bapress1 = file1->$BAPRESNAME1$(TIME1,:,i1,j1)
>
>     pbpress1 = file1->$PBPRESNAME1$(TIME1,:,i1,j1)
>
>     bageopot1 = file1->$BAGEOPNAME1$(TIME1,:,i1,j1)
>
>     pbgeopot1 = file1->$PBGEOPNAME1$(TIME1,:,i1,j1)
>
>     uwrf1 = file1->$UCOMPNAME1$(TIME1,:,:,:)
>
>     vwrf1 = file1->$VCOMPNAME1$(TIME1,:,:,:)
>
>     qrain = file1->$QRAINNAME1$(TIME1,:,i1,j1)
>
>     qsnow = file1->$QSNOWNAME1$(TIME1,:,i1,j1)
>
>     qgraup = file1->$QGRAUPELNAME1$(TIME1,:,i1,j1)
>
>
>
>
>
> Thanks and Happy Holidays!!
>
>
>
> Dave
>
>
>
>
>
> Sent from Mail <https://go.microsoft.com/fwlink/?LinkId=550986> for
> Windows 10
>
>
> _______________________________________________
> 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/20201224/117ee8ec/attachment.html>