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

[swati shewale]

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/8c373139/attachment.html>