[ncl-talk] Wrfoutput plots empty workstation

afwande juliet afwandej965 at gmail.com
Mon Nov 2 18:54:52 MST 2015


ok mary
thanks a lot
this problem invalid data  and time appears to occur when I average my
files using CDO first before ingesting them in NCL or another graphical tool
thanks once again

On Mon, Nov 2, 2015 at 12:18 PM, Mary Haley <haley at ucar.edu> wrote:

> Thanks for providing the data and script.
>
> First, I was unable to even run your script because of this line:
>
>   times = wrf_user_getvar(a,"Time",-1)  ; get all times in the file
>
>
> There is no variable "Time" on the file.
>
> I fixed this with:
>
>   times  = wrf_user_list_times(a)
>   ntimes = dimsizes(times)         ; number of times in the file
>   printVarSummary(times)
>
> Second, I was unable to read "slp" off the file, because apparently you
> have some invalid data around time step index 94. I got an error:
>
> Troubles finding level  100 above ground.
> Problems first occur at (   1,   1)
> Surface pressure =    0.0 hPa.
> STOP Error_in_finding_100_hPa_up
>
> In order to debug this, instead of trying to read "slp" for all time steps
> at once, I moved the "slp" code to inside the do loop, so I could retrieve
> slp for one timestep at a time.  This is the timestep where it failed:
>
>      Working on time: 2000-04-12_21:00:00
>      STOP Error_in_finding_100_hPa_up
>      Troubles finding level  100 above ground.
>      Problems first occur at (   1,   1)
>      Surface pressure =    0.0 hPa.
>
> This is something that wrfhelp at ucar.edu (cc'ed here) may be able to help
> you with.  I've attached the modified version of your code, so you can at
> least see that your script appears to be working for the earlier time steps.
>
> --Mary
>
>
>
> On Mon, Nov 2, 2015 at 10:06 AM, afwande juliet <afwandej965 at gmail.com>
> wrote:
>
>> Hi mary and all
>> I have tried to work on it but don't know why;
>>
>> I have sent data and script via ftp
>>
>> On Thu, Oct 29, 2015 at 10:31 AM, Mary Haley <haley at ucar.edu> wrote:
>>
>>> Afwande,
>>>
>>> I forgot to mention that you don't want to *xxx*@gsnFrame and *xxx*@gsnDraw
>>> to False.  This is not necessary in your case.
>>>
>>> If this still doesn't fix your problem, then in addition to the data
>>> file, which I already requested, please email me a new copy of your script,
>>> so I can try it here.
>>>
>>> --Mary
>>>
>>>
>>> On Thu, Oct 29, 2015 at 9:55 AM, afwande juliet <afwandej965 at gmail.com>
>>> wrote:
>>>
>>>> thanks
>>>> Mary and Dennis,
>>>> I have tried both suggestions from Dennis and Mary but it doesn't plot.
>>>> The workstation is opened with no plot and error message. I have checked
>>>> the data and realized it has only upto 18hours run but whatever the case
>>>> there should be a plots
>>>>
>>>> thanks
>>>>
>>>> On Thu, Oct 29, 2015 at 7:54 AM, Mary Haley <haley at ucar.edu> wrote:
>>>>
>>>>> Afwande,
>>>>>
>>>>> My guess is that your script never got past this line, because NCL
>>>>> doesn't know what type of file your wrf output file is:
>>>>>
>>>>>   a = addfile("./wrfout_d01_2000-04-01_06:00:00","r")
>>>>>
>>>>> Didn't your script give you an error when you tried to run it?
>>>>>
>>>>> You need to add a ".nc" to the end of the filename (note: is NOT
>>>>> necessary to rename the file itself), so that NCL knows to open it as a
>>>>> NetCDF file:
>>>>>
>>>>>   a = addfile("./wrfout_d01_2000-04-01_06:00:00.nc","r")
>>>>>
>>>>> --Mary
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> On Thu, Oct 29, 2015 at 5:27 AM, afwande juliet <afwandej965 at gmail.com
>>>>> > wrote:
>>>>>
>>>>>> Dear NCL users,
>>>>>>  I am trying to plot wrf output for temp, u and v wind  dewpoint temp
>>>>>> in one plot but i only get empty workstation, could someone help me know
>>>>>> why . I am using script from WRF-NCL page and the runs is for about 13days
>>>>>> of 3hourly. I have attached the script
>>>>>> thanks
>>>>>> Afwande
>>>>>>
>>>>>> _______________________________________________
>>>>>> 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/20151102/5616ba78/attachment.html 


More information about the ncl-talk mailing list