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

20010314: plotting data values with MARK in sfmap



Dave,

MARKER is not a parameter in SFMAP in 5.6+.
Instead, you can specify MARK in the sfparm value.
(this allows you to specify marker based on value
as well as color coding the markers using the COLOR=() syntax).

I see your script plotting SGNL uses MARKER. THat may be the simple problem.

Chiz

On Wed, 14 Mar 2001, David Knight wrote:

> ----------
> X-Sun-Data-Type: text
> X-Sun-Data-Description: text
> X-Sun-Data-Name: text
> X-Sun-Charset: us-ascii
> X-Sun-Content-Lines: 39
> 
> 
> Chiz,
>     The program seems to run, but data is not plotted...
> sflist sees the data and prints it fine. I'll attach a sample
> script. When I use gempak54upc_pl15 and program=ltgmap
> it works fine. When I use gempak56a_upc and sfmap
> I get no lightning data plotted. Oh I called the
> script with sfcity dal 14 18. I made the marker bigger
> in case that was it (since I noticed that 56 seems to use
> smaller fonts by default that 54)...
> 
> Thanks for your help
> David
> > 
> > 
> > Dave,
> > 
> > I increased the number of stations in all of GEMPAK to 30,000
> > so as not to need ltglist and ltgmap (which were just sfmap and sflist
> > with the larger station numbers defined).
> > 
> > I can run sfmap and sflist here. What problems are you seeing?
> > 
> > Chiz
> > 
> > On Wed, 14 Mar 2001, David Knight wrote:
> > 
> > > 
> > > Hi Steve,
> > >      There appears to be no ltgmap in gempak5.6. Did I miss something,
> > > or should it no longer be needed. Near as I can see sfmap will not read
> > > our lightning files. Before I adjusted sfmap to be able to read them
> > > I thought I'd check with you.
> > > Thanks
> > > David
> > > 
> > > 
> > 
> > 
> ----------
> X-Sun-Data-Type: cshell-script
> X-Sun-Data-Description: cshell-script
> X-Sun-Data-Name: sfcity
> X-Sun-Charset: us-ascii
> X-Sun-Content-Lines: 129
> 
> #! /bin/csh -f
> #
> #  D. Knight
> #
> # 
> # simple script to display surface observations in x-window using gempak
> #
> # csh usage: sfcity [city] [dd] [hh]
> #     where: city is three digit city ID - default is ALB
> #            dd is two digit day - default is today
> #            hh is two digit hour - default is most recent hour
> #
> #     note:  city can actually be any "legal" gempak geographic area
> #            i.e. a state name, ny co etc. or a country us uk
> #            or a zoomed in view i.e. alb+ or zoomed out view alb-
> #
> # if only one paramter is given is assumed to be: city
> # if two parameters are given they are assumed to be: city hh
> #
> #  first lets set up the path
> #  set path for normal systems stuff
> set path=( . ~/bin /bin /usr/sbin /etc )
> #  include in path most locally installed software
> set path=($path /usr/local/bin)
> #  define path and environmental variables for gempak
> source /unidata/gempak54upc_pl15/Gemenviron
> #source /unidata/gempak56a_upc/Gemenviron
> set program=ltgmap
> #set program=sfmap
> ###################################################
> #   set year and month
> 
>         set yy=`date -u +%y`
>         set YY=`date -u +%Y`
>         set mm=`date -u +%m`
> 
> ############ PARSE INPUT  ##########################
> 
> switch ($#argv)
>  case 0:
>       set city=alb
>         set dd=`date -u +%d`
>         set hh=`date -u +%H`
>         set mi=`date -u +%M`
>       if ($mi < 8) set hh = `expr $hh - 1`
>         breaksw
>  case 1:
>         set city = $1
>         set dd=`date -u +%d`
>         set hh=`date -u +%H`
>         set mi=`date -u +%M`
>       if ($mi < 8) set hh = `expr $hh - 1`
>         breaksw
>  case 2:
>       set city = $1
>       set hh = $2
>         set dd=`date -u +%d`
>         breaksw
>  case 3:
>       set city = $1
>       set dd = $2
>       set hh = $3
>         breaksw
>  default:
>       echo " usage: sfcity [city] [dd] [hh] "
>       echo " where:  city is three digit city ID - default is ALB"
>       echo "         dd is two digit day - default is today "
>       echo "         hh is two digin hour - default is most recent hour "
>         exit(1)
>  endsw
> 
>       set h1 = `expr $hh - 1`
> 
>     echo  running
> set nldnfile=`ls $nldn/$YY$mm$dd$hh*.gem`
> 
> 
> $program << endsfmap 
> \$mapfil=hipowo.gsf
> area  = @$city/$city
> garea = $city
> sfparm= skyc:.7;tmpf;vsby;smsl;ptnd:.7;dwpf;p24i*100;wsym:.7;BRBK:.7:1:13;snow
> dattim= $yy$mm$dd/$hh
> sffile= $gemdat/$yy$mm$dd.sf
> colors= 31;2;7;5;6;3;1;17;1;1
> map   = 4
> latlon= 0
> marker= 0
> title = 31/-3
> clear = y
> panel = 0
> dev   = xw|sfcity
> proj  = mer
> filter= 0.5 
> filter= 1.5 
> text  = 0.7
> skpmis= no
> r
> 
> sffile=$gemdat/$yy$mm$dd$hh.sb
> sfparm= skyc:.7;tmpf;vsby;smsl;ptnd:.7;dwpf;wsym:.7;BRBK:.7:1:13
> colors= 31;2;7;5;6;3;17;1
> clear = n
> map   = 0
> title = 0
> r
> 
> area  = $city
> SFPARM= sgnl<0
> SFFILE= $nldnfile
> COLORS= 0
> MARKER= 29/15/1/5
> MARKER= 29/15/4/5
> TITLE = 29/-2/ $h1 - $hh NLDN Neg Flash
> r
> 
> SFPARM= sgnl>0
> MARKER= 10/1/.6/2
> MARKER= 10/1/6/2
> TITLE = 10/-1/ $h1 - $hh NLDN Pos Flash
> r
> 
> endsfmap
> #
> 
> echo " Don't forget to clean up by running gpend "
> 
> exit (0)
> 
>