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

swati shewale shewaleswati16 at gmail.com
Fri Sep 22 22:00:30 MDT 2017 

Thank you so much Mary....surely I will go through the link provided by
you....If come across any problem, will let you know

On Sat, Sep 23, 2017 at 12:36 AM, Mary Haley <haley at ucar.edu> wrote:

> 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/20170923/3f67a24a/attachment.html>

More information about the ncl-talk mailing list