Re: Vibeta function

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Tue Nov 12 2013 - 10:58:12 MST

  Your pressure levels are:

100000,92500,85000,70000,60000,50000,40000,30000

I speculate that over the Himalaya region, certain
grid points have missing (_FillValue) data at 50000
and below (100000).

The messages are 'warning' messages. I don't understand the 'fatal'
message.

  fatal:vibeta: ier = -999
> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 44
in file

The underlying code has:

c make sure there are at least several levels.

       if (nlmax.le.3) then
           ier = -999
           vint = xmsg
           return
       endif

This is the only place the -999 occurs.

On 11/12/13 7:00 AM, Guilherme Martins wrote:
> Hi Dennis,
>
> Now, when I run the script appear:
>
> warning:vibeta: there must be at least three levels with data above the
> surface
> warning:vibeta: there must be at least three levels with data above the
> surface
> warning:vibeta: there must be at least three levels with data above the
> surface
> warning:vibeta: there must be at least three levels with data above the
> surface
> fatal:vibeta: ier = -999
> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 44 in file
> fluxo.q.ncl
>
> Guilherme.
>
>
> 2013/11/12 Dennis Shea <shea@ucar.edu>
>
>>
>> Assuming pc(time,lat,lon) and divc(time,lev,lat,lon)
>>
>> then ....
>>
>> fint = vibeta(p,divc(time|:,lat|:,lon|:,lev|:),linlog,pc,pbot,ptop)
>>
>>
>> On 11/11/13 10:28 PM, Guilherme Martins wrote:
>>
>>> Hi all,
>>>
>>> I'm using the vibeta function, but when I run my script I get the error:
>>>
>>> fatal:vibeta: If psfc is not a constant, then its dimension sizes must be
>>> equal to all but the last (rightmost) dimension sizes of x
>>> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 43 in
>>> file
>>> fluxo.q.ncl
>>>
>>> The NCL told me the error, but I can't understand it. Can anyone help me?
>>>
>>> My script is in attach.
>>>
>>> My ncl is 6.1.2.
>>>
>>> Thanks,
>>>
>>> Guilherme.
>>>
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Nov 12 10:58:22 2013

This archive was generated by hypermail 2.1.8 : Fri Nov 22 2013 - 09:36:32 MST