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

Mary Haley haley at ucar.edu
Fri Sep 22 13:06:27 MDT 2017 

Please see this page:

http://www.ncl.ucar.edu/Applications/wrfgsn.shtml

and look at the first example (wrf_gsn_1.ncl and wrf_nogsn_1.ncl) and the
second example (wrf_gsn_2.ncl).

These show how to use gsn to plot data and it compares the plot with using
WRF functions.

In the wrf_gsn_1.ncl plot, you can see that the gsnLeft/RightString
resources are being set automatically, because the data has a description
and units attribute.

Try this script as a start on your own data, and also look at wrf_gsn_2.ncl
to see how to further customize your plot.

If you are unable to get the plot to look as nice as you want, email
ncl-talk at ucar.edu again with your new script and what you want it to look
like, and we can help further.

Good luck,

--Mary


On Fri, Sep 22, 2017 at 3:21 AM, swati shewale <shewaleswati16 at gmail.com>
wrote:

> Thank you Mary... its wrf file, I just renamed it for my convenience.
> Actually I am able to plot it using wrf function but unable to put
> left/right/centre strings, as those are related to gsn functions. so to
> overcome the difficulty of left/right/centre, I thought to make a plot
> using gsn functions injstead of wrf functions. but inbuilt wrf functions
> gives nicer plots than gsn one.
>
> On Fri, Sep 22, 2017 at 3:48 AM, Mary Haley <haley at ucar.edu> wrote:
>
>> It's hard for somebody to simply look at a script to see if it's correct
>> or not.  I do see some errors, but even if I correct the errors, I can't be
>> sure it's 100% correct without having access to your data.
>>
>> The problems are:
>>
>> [1]
>>
>> Your data file is called "patna_data.nc". You are using wrf_xxxx
>> functions to plot this data. Since WRF-ARW output files are normally called
>> something like "wrfout_d02_2005-08-29_00:00:00", I can't be sure that
>> you are actually plotting WRF-ARW data. You can only use wrf_xxx plotting
>> functions on WRF-ARW data.
>>
>> [2]
>>
>> You have:
>>
>>  t   = a->T2(0:7,:,:)
>>  ter = dim_avg_n( ter1, 0 )
>>  ter = ter - 273
>>
>> I don't know what your objective, so I don't know if the above is
>> correct.  However, you are calling "dim_avg_n" on "ter1" which doesn't
>> exist.  Did you mean:
>>
>>  ter = dim_avg_n( t, 0 )
>>
>> The above will take an average of the first 8 times teps of your data,
>> assuming that the leftmost dimension of T2 represents time.
>>
>> Further, you should use the _Wrap version of this function, so metadata
>> is maintained:
>>
>>  ter = dim_avg_n_Wrap( t, 0 )
>>
>> [3]
>>
>> You have:
>>
>> ter = ter - 273
>>
>> Since you are changing the units, you should update the "units" attribute
>> to something like this:
>>
>> ter at units = "degC"
>>
>>
>> The rest of the script looks okay, but again, the wrf_contour and
>> wrf_map_overlays procedures will only work if you are actually plotting
>> WRF-ARW data.
>>
>> --Mary
>>
>>
>>
>> On Thu, Sep 7, 2017 at 10:02 PM, swati shewale <shewaleswati16 at gmail.com>
>> wrote:
>>
>>> Thank you. I will try with your script. yesterday I tried it with below
>>> script. Please have a look and let me know, whether it is correct or not.
>>>
>>> On Fri, Sep 8, 2017 at 12:54 AM, Mary Haley <haley at ucar.edu> wrote:
>>>
>>>> If I understand you correctly, then you need to reshape your X array so
>>>> that it is NTIM x 8 x NLAT x NLON, and then you can average across the
>>>> second dimension from the left.
>>>>
>>>> I've attached a modified version of your script, which is UNTESTED. You
>>>> may need to modify it to do the right thing with regard to copying over the
>>>> required metadata.  I did this a little with copy_VarAtts, and copy_coords.
>>>>
>>>> --Mary
>>>>
>>>>
>>>>
>>>>
>>>> On Thu, Sep 7, 2017 at 3:42 AM, swati shewale <shewaleswati16 at gmail.com
>>>> > wrote:
>>>>
>>>>> Hiii,
>>>>> using ncrcat I have single file containing 24 time step (3hrly, 3
>>>>> days). Now I want to plot temperature variable (average for each day,
>>>>> spatial, panel plot). I am unable to take out time step for 1st day (0:7)
>>>>> and then do average of it...then second day (8:15) ...likewise for third
>>>>> day (16:23). I wrote a script, in that it takes average over all time
>>>>> steps. what changes should I make so that I can get average over specific
>>>>> time step? attached herewith script.
>>>>> summary of variable temperature is as follow.
>>>>>
>>>>> Variable: X
>>>>> Type: float
>>>>> Total Size: 5658552 bytes
>>>>>             1414638 values
>>>>> Number of Dimensions: 3
>>>>> Dimensions and sizes:    [Time | 23] x [south_north | 201] x
>>>>> [west_east | 306]
>>>>> Coordinates:
>>>>> Number Of Attributes: 6
>>>>>   FieldType :    104
>>>>>   MemoryOrder :    XY
>>>>>   description :    TEMP at 2 M
>>>>>   units :    K
>>>>>   stagger :
>>>>>   coordinates :    XLONG XLAT XTIME
>>>>>
>>>>>
>>>>>
>>>>> On Thu, Sep 7, 2017 at 3:04 PM, swati shewale <
>>>>> shewaleswati16 at gmail.com> wrote:
>>>>>
>>>>>> Thank you. I will try with it.
>>>>>>
>>>>>> On Wed, Sep 6, 2017 at 8:39 PM, Mary Haley <haley at ucar.edu> wrote:
>>>>>>
>>>>>>> The error message is telling you that there's something wrong with
>>>>>>> argument 0 passed to wrf_contour.  Argument 0 refers to the first argument
>>>>>>> in the list, which is "a(itime)" in this case.
>>>>>>>
>>>>>>> The issue is that when you use addfiles to open a series of files,
>>>>>>> you need to use square brackets, [], to access a single file, and not
>>>>>>> curved brackets, (). This is because addfiles returns a "List" object, and
>>>>>>> lists in NCL need to be access with square brackets.
>>>>>>>
>>>>>>> You have:
>>>>>>>
>>>>>>>         contour = wrf_contour(a(itime),wks,t,t_res)
>>>>>>>
>>>>>>> and this should be:
>>>>>>>
>>>>>>>         contour = wrf_contour(a[itime],wks,t,t_res)
>>>>>>>
>>>>>>> --Mary
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> On Mon, Sep 4, 2017 at 4:40 AM, 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
>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>
>>>>>
>>>>
>>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170922/ca9d60ad/attachment.html>

More information about the ncl-talk mailing list