[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

20000306: gempak_pl15 bug?

This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.

  • Subject: 20000306: gempak_pl15 bug?
  • Date: Mon, 06 Mar 2000 09:45:42 -0700
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
>