[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