20000306: gempak_pl15 bug?

>From: "David J. Knight" <address@hidden>
>Organization: .
>Keywords: 200003061841.LAA01106

>Chiz,
>     yep, gfunc=res
>works fine.
>
>now try gfunc=mul(res,3600)
>This fails for me.
>
>Thanks
>David


David,

the mul(res,3600) works for me with gdattim = f000:f024 and glevel=700:850:

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

 GAREA:     uslcc                                                               
     SCALE: -3

 MINIMUM AND MAXIMUM VALUES  -188.68   205.52

 LINE CONTOURS: 

 LEVELS:      -160.00  -120.00   -80.00   -40.00     0.00    40.00    80.00
 COLORS:            5        5        5        5        5        5        5
 LINTYP:            1        1        1        1        1        1        1
 LINWID:            1        1        1        1        1        1        1
 LABEL:             1        1        1        1        1        1        1

 LEVELS:       120.00   160.00   200.00
 COLORS:            5        5        5
 LINTYP:            1        1        1
 LINWID:            1        1        1
 LABEL:             1        1        1
etc...

Chiz

>
>> 
>> 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 routi
> nes
>> > ,
>> >> 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.ge
> m
>> >> > 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.alban
> y.e
>> > du
>> >> >
>> >> 
>> >> *************************************************************************
> ***
>> >> Unidata User Support                                    UCAR Unidata Prog
> ram
>> >> (303)497-8644                                                  P.O. Box 3
> 000
>> >> address@hidden                                   Boulder, CO 80
> 307
>> >> -------------------------------------------------------------------------
> ---
>> >> Unidata WWW Service                        http://www.unidata.ucar.edu/  
>    
>> >> *************************************************************************
> ***
>> >> 
>> >
>> 
>> ****************************************************************************
>> 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/     
>> ****************************************************************************
>