[ncl-talk] Error while handling multiple WRF output files

Tue Sep 5 22:19:51 MDT 2017

Hi Swati:

This program plots radar reflectivity from a single wrf output file at
different times, over the eastern US, in a panel format.  You need to
change the variable to the one (e.g, REFL10) that you have in your output.

To change the location of the plot, change the values within lats/lons.

Of course, you can change the variable you wish to plot.

I have included the program I use to plot gfs temperatures. It requires
that you run write_times.F, first (by the way, through interactions with
NCL-Talk I modified a similar program to obtain the time of each gefs file
directly from the file itself -- you can find that in the e-mail chain). I
include this because it has a nice way of setting the temperature intervals
and color bar.

After you obtain the result you want from a single wrf output file, please
let me know if you have trouble reading more than output file at a time.
The plot_gfs program shows you how to do this.

Barry

On Wed, Sep 6, 2017 at 7:06 AM, swati shewale <shewaleswati16 at gmail.com>
wrote:

> Hiii Bary,
> Please send me the script, I will look into it sideby. Thanks
>
> On Wed, Sep 6, 2017 at 9:35 AM, swati shewale <shewaleswati16 at gmail.com>
> wrote:
>
>> Hiii Rick,
>> Keeping your suggestion (of changing "a(itime)" to "a[itime]" in the call
>> to wrf_contour()) in the script gives me error as follow
>>
>> fatal:Number of dimensions in parameter (2) of (wrf_contour) is (3), (2)
>> dimensions were expected
>> fatal:["Execute.c":8567]:Execute: Error occurred at or near line 44 in
>> file temp.ncl
>>
>> (0)    Error: gsn_panel: all of the plots passed to gsn_panel appear to
>> be invalid
>>
>>
>> On Wed, Sep 6, 2017 at 5:55 AM, Barry Lynn <barry.h.lynn at gmail.com>
>> wrote:
>>
>>> Hi: did Rick's suggestion solve your problem? If not, I can send you my
>>> program that uses the same approach - but for radar from a single wrf
>>> output file -- to map panels from different times.
>>>
>>>
>>> On Tuesday, September 5, 2017, Rick Brownrigg <brownrig at ucar.edu> wrote:
>>>
>>>> I *think* the issue is that addfiles() returns a list of
>>>> file-variables, not an array. You might try changing "a(itime)" to
>>>> "a[itime]" in the call to wrf_contour().
>>>>
>>>> Rick
>>>>
>>>> On Tue, Sep 5, 2017 at 12:29 AM, Barry Lynn <barry.h.lynn at gmail.com>
>>>> wrote:
>>>>
>>>>> Are you able to plot the temperature variable for a single WRF run?
>>>>>
>>>>> On Tue, Sep 5, 2017 at 9:16 AM, swati shewale <
>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>
>>>>>> Thank you for reciprocation so quickly. I want to plot temperature
>>>>>> variable.
>>>>>>
>>>>>> On Mon, Sep 4, 2017 at 5:19 PM, Barry Lynn <barry.h.lynn at gmail.com>
>>>>>> wrote:
>>>>>>
>>>>>>> Hi:
>>>>>>> Can you tell me what you are trying to plot (without me needing to
>>>>>>> look...).
>>>>>>>
>>>>>>> Which variable?
>>>>>>>
>>>>>>> I might have a program you can use with just a small modification.
>>>>>>>
>>>>>>> Barry
>>>>>>>
>>>>>>> On Mon, Sep 4, 2017 at 1:40 PM, swati shewale <
>>>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>>>
>>>>>>>> Hiiii
>>>>>>>>
>>>>>>>> I want to do panel plotting using WRF output. I have total 8 files,
>>>>>>>> each contains single time step (one day, 3 hrly datafiles so total 8
>>>>>>>> files). out of eight times I want to plot 6 time step (panel plot) for
>>>>>>>> variable slp and temperature.
>>>>>>>>
>>>>>>>> It shows below error.
>>>>>>>>
>>>>>>>> fatal:Argument type mismatch on argument (0) of (wrf_contour) can
>>>>>>>> not coerce
>>>>>>>> fatal:["Execute.c":8567]:Execute: Error occurred at or near line
>>>>>>>> 44 in file temp.ncl
>>>>>>>> (0)    Error: gsn_panel: all of the plots passed to gsn_panel
>>>>>>>> appear to be invalid
>>>>>>>>
>>>>>>>>  Attached herewith script.
>>>>>>>>
>>>>>>>> Thank you in advance. Any help will be appreciated
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> _______________________________________________
>>>>>>>> ncl-talk mailing list
>>>>>>>> ncl-talk at ucar.edu
>>>>>>>> List instructions, subscriber options, unsubscribe:
>>>>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> --
>>>>>>> Barry H. Lynn, Ph.D
>>>>>>> Senior Lecturer,
>>>>>>> The Institute of the Earth Science,
>>>>>>> The Hebrew University of Jerusalem,
>>>>>>> Givat Ram, Jerusalem 91904, Israel
>>>>>>> Tel: 972 547 231 170
>>>>>>> Fax: (972)-25662581
>>>>>>>
>>>>>>> C.E.O, Weather It Is, LTD
>>>>>>> Weather and Climate Focus
>>>>>>> http://weather-it-is.com
>>>>>>> Jerusalem, Israel
>>>>>>> Local: 02 930 9525
>>>>>>> Cell: 054 7 231 170
>>>>>>> Int-IS: x972 2 930 9525
>>>>>>> US 914 432 3108 <(914)%20432-3108>
>>>>>>>
>>>>>>
>>>>>>
>>>>>
>>>>>
>>>>> --
>>>>> Barry H. Lynn, Ph.D
>>>>> Senior Lecturer,
>>>>> The Institute of the Earth Science,
>>>>> The Hebrew University of Jerusalem,
>>>>> Givat Ram, Jerusalem 91904, Israel
>>>>> Tel: 972 547 231 170
>>>>> Fax: (972)-25662581
>>>>>
>>>>> C.E.O, Weather It Is, LTD
>>>>> Weather and Climate Focus
>>>>> http://weather-it-is.com
>>>>> Jerusalem, Israel
>>>>> Local: 02 930 9525
>>>>> Cell: 054 7 231 170
>>>>> Int-IS: x972 2 930 9525
>>>>> US 914 432 3108 <(914)%20432-3108>
>>>>>
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> ncl-talk at ucar.edu
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>
>>>>>
>>>>
>>>
>>> --
>>> Barry H. Lynn, Ph.D
>>> Senior Lecturer,
>>> The Institute of the Earth Science,
>>> The Hebrew University of Jerusalem,
>>> Givat Ram, Jerusalem 91904, Israel
>>> Tel: 972 547 231 170
>>> Fax: (972)-25662581
>>>
>>> C.E.O, Weather It Is, LTD
>>> Weather and Climate Focus
>>> http://weather-it-is.com
>>> Jerusalem, Israel
>>> Local: 02 930 9525
>>> Cell: 054 7 231 170
>>> Int-IS: x972 2 930 9525
>>> US 914 432 3108 <(914)%20432-3108>
>>>
>>>
>>
>

-- 
Barry H. Lynn, Ph.D
Senior Lecturer,
The Institute of the Earth Science,
The Hebrew University of Jerusalem,
Givat Ram, Jerusalem 91904, Israel
Tel: 972 547 231 170
Fax: (972)-25662581

C.E.O, Weather It Is, LTD
Weather and Climate Focus
http://weather-it-is.com
Jerusalem, Israel
Local: 02 930 9525
Cell: 054 7 231 170
Int-IS: x972 2 930 9525
US 914 432 3108
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170906/3ad2f31d/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: wrf_gsn_radar_Xkm.ncl
Type: application/octet-stream
Size: 9567 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170906/3ad2f31d/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: plot_gfs_t2.ncl
Type: application/octet-stream
Size: 9491 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170906/3ad2f31d/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: write_times.F
Type: application/octet-stream
Size: 1650 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170906/3ad2f31d/attachment-0002.obj>