20000306: gempak_pl15 bug?

David,

I created a RES variable here in one of my grid files:

 GRID IDENTIFIER: 
    TIME1             TIME2         LEVL1 LEVL2   VCORD PARM
000306/1200F000    000306/1200F024    700  850     PRES RES

Then ran gdcntr:
 GDATTIM  = f000:f024
 GLEVEL   = 700:850
 GVCORD   = pres
 GFUNC    = res

This all worked here.

Since you mentioned that you are running pl15 and compared to pl6,
I assume this means that you built pl15 in a new clean separate directory
from the source distribution. Is that correct?

I might need to get you to upload a sample of your grid so that I can
see if I can duplicate your results here. It doen't appear that
I can duplicate your results by creating a grid from gddiag as I did above.

Steve Chiswell
Unidata User Support


>From: "David J. Knight" <address@hidden>
>Organization: .
>Keywords: 200003061655.JAA27660

>
>Nope,
>     That doesn't work either. I tried spcifying
>gdattim, glevel and gvcord individually,
>then as part of the function, then both.
>All three gave the same error
>David
>> 
>> 
>> David,
>> 
>> With your specification of gdattim, glevel and gvcord,
>> you should be able to use:
>> GFUNC=mul(RES,3600)
>> Since the ^, @ and % you have are the same as what
>> your definitions above are. Does that work? if not, there may
>> be something more fundemental than the mul() routine.
>> 
>> The primary things that were upgraded to pl15 were the 4 digit year routines
> ,
>> so there may be something there that is causing problems.
>> 
>> Steve Chiswell
>> Unidata User Support
>> 
>> 
>> 
>> 
>> >From: "David J. Knight" <address@hidden>
>> >Organization: .
>> >Keywords: 200003061412.HAA21958
>> 
>> >Hi,
>> >     we seem to be having a problem with gempak54upc_pl15
>> >
>> >when I run gdcntr with
>> >GDATTIM  = 940706/00F21:940706/00F24
>> > GLEVEL   = 925:850
>> > GVCORD   = pres
>> > GFUNC    = mul(RES^940706/00F21:940706/00F24@925:850%pres,3600)
>> > GDFILE   = /net/acacia/acacia1/data/nocera/mm5/mmoutp2.type2b2thermo.gem
>> > CINT     = 1/-1/-1
>> > LINE     = 6/2/4/2
>> > MAP      = 3//5
>> > TITLE    = 1/-2/940706/00F24/925mb Thermo Res per h
>> > DEVICE   = ps
>> > SATFIL   = 
>> > RADFIL   = 
>> > PROJ     = ced
>> > GAREA    = -35;-75;-5;-53
>> > CLEAR    = n
>> > PANEL    = 4
>> > TEXT     = 1/22//HW
>> > SCALE    = 1
>> > LATLON   = 1/10/1/2;2/10;10
>> > HILO     = 
>> > HLSYM    = 
>> > CLRBAR   = 
>> > CONTUR   = 0
>> > SKIP     = 0
>> > FINT     = 1500/1500/6000
>> > FLINE    = 0;8;8;8;8
>> > CTYPE    = C
>> > LUTFIL   = 
>> > STNPLT   = 
>> >
>> >we get the error:
>> >[DG -23]  LEVEL @ in MUL is invalid.
>> >
>> >
>> >with 
>> >GFUNC    = RES^940706/00F21:940706/00F24@925:850%pres
>> >
>> >it runs fine.
>> >
>> >using gempak5.4 pl6 with
>> > GFUNC    = mul(RES^940706/00F21:940706/00F24@925:850%pres,3600)
>> >also works fine
>> >
>> >Any suggestions? I guess the obvious is to include
>> >the 3600 when the grid is created, but that doesn't solve
>> >what appears to be the underlying problem
>> >
>> >Thanks
>> >David
>> >David Knight
>> >Department of Earth and Atmospheric Sciences   Tel: (518)-442-4204
>> >SUNYA   ES-228                                 Fax: (518)-442-4494
>> >Albany, NY  12222                              Email: knight@atmos.albany.e
> du
>> >
>> 
>> ****************************************************************************
>> Unidata User Support                                    UCAR Unidata Program
>> (303)497-8644                                                  P.O. Box 3000
>> address@hidden                                   Boulder, CO 80307
>> ----------------------------------------------------------------------------
>> Unidata WWW Service                        http://www.unidata.ucar.edu/     
>> ****************************************************************************
>> 
>