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

20021116: 20021115: GEMPAK scripts

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: 20021116: 20021115: GEMPAK scripts
  • Date: Sat, 16 Nov 2002 11:20:19 -0700
Werner,

Check to make sure that you have correctly set the panel for the first plot and 
are not
just overwriting the first display. Your problem sounds like a scripting error 
since
it is not the case that you cannot do multiple overlays within a single 
invocation.

Steve Chiswell

>From: Werner Wintels <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200211160540.gAG5ekL23504

>
>This is only a segment of the script.  I did run gpend at the end of the 
>script, as I do at the end of all my scripts.
> 
>The point is that I have to exit and reenter gempak or only the first run 
>will plot.  I have this problem with gif files as well. 
>
>Sorry for not including the entire script, but it was rather long and I 
>wanted to get to the point. 
>
>Thanks again.
>
>On Fri, 15 Nov 2002, Unidata Support wrote:
>
>> 
>> 
>> Werner,
>> 
>> You show below that you are using the PS device.
>> I don't see that you have run "gpend" following
>> you scripts. Without running gpend, you have not
>> flushed the output buffer and appended the 
>> closing sequence with "showpage" in the postscript file.
>> 
>> Remember, that until you have run gpend, you can still
>> append more plots to the file, and only after you have
>> run gpend will your drawing commands be flushed to the file
>> from the remaining buffer and have the complete trailing postscript
>> in the file which tells the printer that the page is complete.
>> 
>> A postscript plot will have the vector drawing commands buffered and
>> written a chunk at a time. A gif file, on the otherhand will
>> have nothing until the gpend command is issued since the gif file
>> can only be created from the final raster after all drawing has finished.
>> 
>> Steve Chiswell
>> 
>> 
>> 
>> >The following only plots the first panel:
>> >
>> >gdcntr << endgdcntr
>> > GDFILE = 02102612_eta211.gem
>> > GDATTIM  = F000
>> >
>> ><snip>
>> >
>> > clear = n
>> > title = 1/-2/1000 -500 hPa Thickness (dashed) (6 dam)
>> >r
>> >
>> > GDFILE = 02102512_eta211.gem
>> > GDATTIM = F024
>> > panel = 2
>> >r
>> >
>> >exit
>> >endgdcntr
>> 
>> 
>> 
>> 
>> 
>> 
>> 
>> 
>> 
>> >From: Werner Wintels <address@hidden>
>> >Organization: UCAR/Unidata
>> >Keywords: 200211152040.gAFKeQL11677
>> 
>> >
>> >There seems to be some inconsistency in the performance of gempak 
>> >scripts. On some platforms (and some versions of GEMPAK) you need to exit 
>> >and re-enter the GEMPAK programs each time in order to produce 
>> >multi-panel or multi contour plots. 
>> >
>> >For example, this script segment will produce plots in panel 1 and panel 2:
>  
>> >
>> >
>> >gdcntr << endgdcntr
>> > GDATTIM = F000
>> > GLEVEL = 850
>> > panel  = 1
>> > GFUNC = mixr
>> > GDFILE = 02102612_eta211.gem
>> > CINT     = 1
>> >LINE     = 1/1/1/1
>> > MAP      = 1
>> > DEVICE   = ps
>> > PROJ     = STR/90;-80;0
>> >GAREA    = 35;-90;50;-60
>> > CLEAR    = n
>> > PANEL    = 3
>> > TEXT     = .7
>> > SCALE    = 0
>> > LATLON   = 1/10
>> > HILO     =
>> > HLSYM    = 3.2;1/3/3/6
>> > CLRBAR   =
>> > CONTUR   =
>> > FINT     = 2;3;4;5;6;8;10
>> >  FLINE    = 0;19;16;13;11;8;19;5
>> > CTYPE    = C/F
>> > LUTFIL   = 
>> > STNPLT   =
>> >r
>> >
>> >exit 
>> >endgdcntr 
>> >
>> >gdcntr << endgdcntr
>> > GDFILE   = 02102512_eta211.gem
>> > GDATTIM  = F024
>> > panel    = 2
>> >r
>> >
>> >exit
>> >endgdcntr
>> >
>> >The following only plots the first panel:
>> >
>> >gdcntr << endgdcntr 
>> > GDFILE = 02102612_eta211.gem
>> > GDATTIM  = F000
>> >
>> ><snip>
>> >
>> > clear = n
>> > title = 1/-2/1000 -500 hPa Thickness (dashed) (6 dam)
>> >r 
>> >
> > GDFILE = 02102512_eta211.gem
>> > GDATTIM = F024
>> > panel = 2
>> >r
>> >
>> >exit 
>> >endgdcntr
>> >
>> >
>> >I'm not sure why the second script only works some of the time.  
>> >
>> >Thanks
>> >
>> >Werner Wintels
>> >
>> 
>> ****************************************************************************
>> Unidata User Support                                    UCAR Unidata Program
>> (303)497-8643                                                  P.O. Box 3000
>> address@hidden                                   Boulder, CO 80307
>> ----------------------------------------------------------------------------
>> Unidata WWW Service                        http://www.unidata.ucar.edu/     
>> ****************************************************************************
>> 
>
>********************************
>Werner Wintels 
>Atmospheric and Oceanic Sciences
>McGill University    
>Montreal, Canada
>********************************
>