[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
19990127: time range missing in gdplot
- Subject: 19990127: time range missing in gdplot
- Date: Thu, 28 Jan 1999 13:30:23 -0700
Steve,
The root of the problem is the number of grid times in the file (341).
GEMPAK has 2 array sizes defined in $GEMPAKHOME/include/gemprm.xxx,
LLMXTM (maximum number of times in a file, default 200) and
LLMXGT (maximum number of grid times, default 1000). My interpretation
of those variables says you should be allowd to exceed 200 times in
your file, however, the routine gr_tlst has LLMXTM as the array
size rather than LLMXGT. So, the problem is really that the number of
times exceeds the array sizes and they are getting stomped on.
The complete fix will be to find all usage of the number of times and make
them agree.
Here is a quick fix for the immediate problem at hand.
edit $GEMPAKHOME/src/gemlib/gr/grtlst.f
--> change timlst(LLMXTM) to timlst(LLMXGT)
CHARACTER*20 timlst (LLMXGT), timout*40
--> for consitency, set LLMXGT in GDGTIM call...its not used in that
--> subroutine however...go figure!
C* Get list of times in data file.
C
CALL GD_GTIM ( igdfln, LLMXGT, timlst, ntimin, ier )
Now rebuild the subroutine:
cd $GEMPAKHOME/src/gemlib/gr/
make clean
make
make clean
....Now here is why gdplot is a little different in its screw up than
the other programs you were using...it makes the tgfind call directly.
Edit $GEMPAKHOME/src/programs/gd/gdplot/gdbtim.f
--> change LLMXTM to LLMXGT in these instances
CHARACTER*36 time, timlst (LLMXGT), tstrt, tstop, ttt, timout
and
CALL GD_GTIM ( iflno, LLMXGT, timlst, ntimin, ier )
and rebuild all the gd programs
cd $GEMPAKHOME/src/programs/gd
make clean
make all
make install
make clean
This quick fix will work for the internal gr_tlst and tg_find calls since
even though you have 341 timies in your file, when you plot
gdattim=f00031-f12031-01200, you are only asking for 11 times in the
range, so you are no where near exceeding the 200 time maximum in
other locations at this time. The bombing out right now is simply because
if can't get all the times from the file to search for your range.
As I said, the complete fix will be to ensure all sizes are in agreement
so that if you ask for a range of greater than 200 times, then
things will behave properly (of course you will hit the frame
limit in looping in gdplot at 32 or so frames).
Let me know if this works for you, or other problems still exist.
Steve Chiswell
Unidata User SUpport
On Thu, 28 Jan 1999, Steven L. Mullen wrote:
> Steve,
>
> Just login in to my account on cyclone.
>
> Steve
>
>
> Steve Chiswell wrote:
>
> > Steve,
> >
> > If you can post the file, I'd be willing to take a look at it.
> > Seems like the time range is causing a problem. It looks like the
> > first example should stop at f06031, but I see output
> > for f10831 and f12031. Thats not what I would expect.
> > Its possible to have problems in a range or list is there
> > are any times in the file with Time1:time2 defined.
> > I also notice a level of -9999 for one of the parameters,
> > and not sure how things like that might affect the lookup.
> >
> > Steve
> >
> > On Wed, 27 Jan 1999, Steven L. Mullen wrote:
> >
> > > Steve,
> > >
> > > Thanks for the excellent tip on ps directories. I would like
> > > to tap your genius status again.
> > >
> > > We have another particularity with our installation of gempak.
> > > What might cause gdplot to not be able to plot certain times
> > > for fields that definitely exist (I can gdlist them and gddiag
> > > them) in a gempak grid management file? I have a couple of
> > > gdplots outputs, a gdlist output, and a portion of a gdlist output.
> > >
> > > Your advice is appreciated.
> > >
> > > Steve
> > >
> > > GDFILE Grid file
> > > /cuser/mullen/for/sanders/plots/gem/semi1.gem
> > > GDATTIM Grid date/time f00031-f06031-01200
> > > GLEVEL Grid level 500
> > > GVCORD Grid vertical coordinate pres
> > > PANEL Panel loc/color/dash/width/regn 0
> > > SKIP Skip_cntr/skip_plt_x;skip_plt_y 0
> > > SCALE Scalar scale / vector scale 15!-1
> > > GFUNC Scalar grid
> > > cros(grad(sm9s(zbar)),grad(mag(grad(sm9s(zbar)))!sm9s(zbar)
> > > CTYPE Contour type: C/F f!c
> > > CONTUR Subbox/smooth 0
> > > CINT Contour interval/min/max !6
> > > LINE Line color/type/width/label !1/1/4
> > > FINT Fill interval/min/max
> > > -10;-9;-8;-7;-6;-5;-4;-3;-2;-1;0;1;2;3;4;5;6;7;8;9;10
> > > FLINE Fill colors 30-21;0;0;20-10
> > > HILO Color/symbol/rng/rad/cnt/intp
> > > HLSYM HILO txt size/posn/font/wdth/hw
> > > CLRBAR Color/ornt/anch/x;y/ln;wd/freq 1/h/ll/.01;.01/.8;.02
> > > GVECT Vector grid
> > > WIND Wind symbol/siz/wdth/typ/hdsz am1
> > > REFVEC Mag;x;y;txtsiz/font/wdth/HW;labl
> > > TITLE Title color/line/title 1
> > > TEXT Text size/font/width/hw flag 1
> > > CLEAR Clear screen flag YES
> > > GAREA Graphics area 5;130;45;-75
> > > PROJ Map projection/angles/margins|dr lcc/60;-180;30
> > > MAP Map color/dash/width 8/1/2
> > > LATLON Line color/dash/width/label/inc 8/30/3
> > > DEVICE Device|name|x size;y size|color XW
> > > STNPLT Txc|mkc/mkt/mks/mkwd/hwfl|stnfil
> > > SATFIL Satellite image filename(s)
> > > RADFIL Radar image filename(s)
> > > LUTFIL Enhancement lookup table filenam
> > > GEMPAK-GDPLOT>gdattim=f00031-f12031-01200
> > > GEMPAK-GDPLOT>r
> > >
> > > PARAMETERS FOR GDPLOT:
> > >
> > > Grid file: /cuser/mullen/for/sanders/plots/gem/semi1.gem
> > > TIMES: 760106/1200F00031
> > > 760106/1200F01231
> > > 760106/1200F02431
> > > 760106/1200F03631
> > > 760106/1200F04831
> > > 760106/1200F06031
> > > 760106/1200F10831
> > > 760106/1200F12031
> > >
> > > GRID IDENTIFIER:
> > > PLOT LEVEL VCORD PARM WPARM
> > > 1 500 PRES CROS(GRAD(SM9S(ZBAR)),GRAD(
> > > 2 500 PRES SM9S(ZBAR)
> > >
> > > GAREA: 5;130;45;-75
> > > Parameters requested:
> > > GDFILE,GDATTIM,GLEVEL,GVCORD,PANEL,SKIP,SCALE,GFUNC,
> > >
> > > CTYPE,CONTUR,CINT,LINE,FINT,FLINE,HILO,HLSYM,CLRBAR,GVECT,WIND,REFVEC,TITLE,
> > >
> > > TEXT,CLEAR,GAREA,PROJ,MAP,LATLON,DEVICE,STNPLT,SATFIL,RADFIL,LUTFIL.
> > > GEMPAK-GDPLOT>
> > >
> > >
> > > GDFILE Grid file
> > > /cuser/mullen/for/sanders/plots/gem/semi1.gem
> > > GDATTIM Grid date/time f00031-f12031-01200
> > > GLEVEL Grid level 500
> > > GVCORD Grid vertical coordinate pres
> > > PANEL Panel loc/color/dash/width/regn 0
> > > SKIP Skip_cntr/skip_plt_x;skip_plt_y 0
> > > SCALE Scalar scale / vector scale 15!-1
> > > GFUNC Scalar grid
> > > cros(grad(sm9s(zbar)),grad(mag(grad(sm9s(zbar)))!sm9s(zbar)
> > > CTYPE Contour type: C/F f!c
> > > CONTUR Subbox/smooth 0
> > > CINT Contour interval/min/max !6
> > > LINE Line color/type/width/label !1/1/4
> > > FINT Fill interval/min/max
> > > -10;-9;-8;-7;-6;-5;-4;-3;-2;-1;0;1;2;3;4;5;6;7;8;9;10
> > > FLINE Fill colors 30-21;0;0;20-10
> > > HILO Color/symbol/rng/rad/cnt/intp
> > > HLSYM HILO txt size/posn/font/wdth/hw
> > > CLRBAR Color/ornt/anch/x;y/ln;wd/freq 1/h/ll/.01;.01/.8;.02
> > > GVECT Vector grid
> > > WIND Wind symbol/siz/wdth/typ/hdsz am1
> > > REFVEC Mag;x;y;txtsiz/font/wdth/HW;labl
> > > TITLE Title color/line/title 1
> > > TEXT Text size/font/width/hw flag 1
> > > CLEAR Clear screen flag YES
> > > GAREA Graphics area 5;130;45;-75
> > > PROJ Map projection/angles/margins|dr lcc/60;-180;30
> > > MAP Map color/dash/width 8/1/2
> > > LATLON Line color/dash/width/label/inc 8/30/3
> > > DEVICE Device|name|x size;y size|color XW
> > > STNPLT Txc|mkc/mkt/mks/mkwd/hwfl|stnfil
> > > SATFIL Satellite image filename(s)
> > > RADFIL Radar image filename(s)
> > > LUTFIL Enhancement lookup table filenam
> > > GEMPAK-GDPLOT>gdattim=f08431
> > > GEMPAK-GDPLOT>r
> > > [TG 2] The time 760106/1200F08431 is not found in the data set.
> > > [TG -5] Invalid time range.
> > > [TG 0]
> > > Parameters requested:
> > > GDFILE,GDATTIM,GLEVEL,GVCORD,PANEL,SKIP,SCALE,GFUNC,
> > >
> > > CTYPE,CONTUR,CINT,LINE,FINT,FLINE,HILO,HLSYM,CLRBAR,GVECT,WIND,REFVEC,TITLE,
> > >
> > > TEXT,CLEAR,GAREA,PROJ,MAP,LATLON,DEVICE,STNPLT,SATFIL,RADFIL,LUTFIL.
> > > GEMPAK-GDPLOT>
> > >
> > >
> > >
> > > GDATTIM Grid date/time f08431
> > > GLEVEL Grid level 500
> > > GVCORD Grid vertical coordinate pres
> > > GFUNC Scalar grid
> > > cros(grad(sm9s(zbar)),grad(mag(grad(sm9s(zbar)))
> > > GDFILE Grid file
> > > /cuser/mullen/for/sanders/plots/gem/semi1.gem
> > > GAREA Graphics area 25;160;55;-150
> > > PROJ Map projection/angles/margins|dr lcc/60;-180;30
> > > SCALE Scalar scale / vector scale 15
> > > OUTPUT Output device/filename t
> > > GEMPAK-GDLIST>r
> > >
> > >
> > > GDLIST PARAMETERS:
> > >
> > > Grid file: /cuser/mullen/for/sanders/plots/gem/semi1.gem
> > > GRID IDENTIFIER:
> > > TIME1 TIME2 LEVL1 LEVL2 VCORD PARM
> > > 760106/1200F08431 500 PRES CROSGRADUGRA
> > >
> > > GAREA: 25;160;55;-150
> > > SCALE FACTOR : 10**15
> > > OUTPUT: TERM/
> > >
> > > MINIMUM AND MAXIMUM VALUES -15.78 10.79
> > > Enter <cr> to accept parameters or type EXIT:
> > >
> > >
> > > Grid file: /cuser/mullen/for/sanders/plots/gem/semi1.gem
> > > GRID IDENTIFIER:
> > > TIME1 TIME2 LEVL1 LEVL2 VCORD PARM
> > > 760106/1200F08431 500 PRES CROSGRADUGRA
> > > AREA: 25;160;55;-150 GRID SIZE:
> > > 24 24
> > > COLUMNS: 9 17 ROWS: 17 23
> > >
> > > Scale factor: 10**15
> > >
> > >
> > > COLUMN: 9 10 11 12 13 14
> > > 15 16
> > > 17
> > > ROW 23 .14 .06 -.54 -1.70 -2.41 -1.98
> > > -1.30 -.91
> > > -.53
> > > ROW 22 .21 1.10 1.48 .75 -.82 -2.01
> > > -2.67 -2.47
> > > -1.52
> > > ROW 21 .22 3.15 5.66 6.13 3.02 -2.00
> > > -5.81 -6.15
> > > -4.00
> > > ROW 20 .27 5.07 9.59 10.79 5.27 -3.94
> > > -11.38 -12.54
> > > -7.87
> > > ROW 19 -.91 3.78 8.45 9.73 4.31 -5.37
> > > -14.51 -15.78
> > > -7.78
> > > ROW 18 -3.38 -.53 3.16 5.42 3.12 -3.45
> > > -10.14 -9.05
> > > .41
> > > ROW 17 -5.00 -3.12 -.91 2.26 3.10 .55
> > > -2.27 .12
> > > 4.24
> > > Parameters requested:
> > > GDATTIM,GLEVEL,GVCORD,GFUNC,GDFILE,GAREA,PROJ,SCALE,
> > > OUTPUT.
> > > GEMPAK-GDLIST>
> > >
> > >
> > > Portion of gdlist output of
> > > /cuser/mullen/for/sanders/plots/gem/semi1.gem
> > >
> > > 511 760106/1200F08427 0 PRES PMSL
> > > 512 760106/1200F08427 500 PRES HGHT
> > > 513 760106/1200F08428 0 PRES PMSL
> > > 514 760106/1200F08428 500 PRES HGHT
> > > 515 760106/1200F08429 0 PRES PMSL
> > > 516 760106/1200F08429 500 PRES HGHT
> > > 517 760106/1200F08430 0 PRES PMSL
> > > 518 760106/1200F08430 500 PRES HGHT
> > > 519 760106/1200F08431 0 PRES PMSL
> > > 520 760106/1200F08431 500 PRES HGHT
> > > 521 760106/1200F08431 500 PRES VORADV
> > >
> > > 522 760106/1200F08431 500 PRES ZBAR
> > > 523 760106/1200F08431 -9999 PRES
> > > THCKADV
> > > 524 760106/1200F09601 0 PRES PMSL
> > > 525 760106/1200F09601 500 PRES HGHT
> > > --
> > > -------------------------------------------------------
> > > {} Steven L. Mullen
> > > {} {} Dept. of Atmospheric Sciences
> > > {} {} {} University of Arizona
> > > {} {} {} 1118 E. 4th St.
> > > {}{}{} PO Box 210081
> > > {} Tucson, Arizona 85721-0081
> > > {} Tel:(520) 621-6842
> > > {} Fax:(520) 621-6833
> > > {} Email: address@hidden
> > > -------------------------------------------------------
> > >
> > >
> > >
>
> --
> -------------------------------------------------------
> {} Steven L. Mullen
> {} {} Dept. of Atmospheric Sciences
> {} {} {} University of Arizona
> {} {} {} 1118 E. 4th St.
> {}{}{} PO Box 210081
> {} Tucson, Arizona 85721-0081
> {} Tel:(520) 621-6842
> {} Fax:(520) 621-6833
> {} Email: address@hidden
> -------------------------------------------------------
>
>
>