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

Barry Lynn barry.h.lynn at gmail.com
Wed Sep 6 01:37:02 MDT 2017


Hi:

I think it looks fine.

(Be sure to include ncl-talk in correspondence, so there is record of the
e-mail chain).

Barry

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

> Hiii
> I did club all the 3 hrly output datafile from wrf output (total 8 files)
> using ncrcat so that i get a single file, which now include 8 time step. to
> handle this processed file, script is available. I hope the way i did might
> be correct one. I am getting the plot attached here along with script.
> please let me know whether it is correct one or not. Thanks to all for
> participation and cooperation. As I am new to NCL, I may post queries
> several times.
>
>
>
> On Wed, Sep 6, 2017 at 10:15 AM, swati shewale <shewaleswati16 at gmail.com>
> wrote:
>
>> Thanks a lot Barry. I will go through it, if come across any problem,
>> will let you know
>>
>> On Wed, Sep 6, 2017 at 9:49 AM, Barry Lynn <barry.h.lynn at gmail.com>
>> wrote:
>>
>>> 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 <%28914%29%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 <%28914%29%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 <%28914%29%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>
>>>
>>
>>
>


-- 
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/04b28ae9/attachment.html>


More information about the ncl-talk mailing list