20000306: gempak_pl15 bug?

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: address@hidden
>