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

20010314: 20010314: plotting data values with MARK in sfmap

David,

The info on mark is spread between the online help fir "sfmap" and
"colors".

You can still run the program twice as before, once for dots
and once for + signs like:

sfparm = mark:15:4:1;sgnl>0
colors = 30;0
clear = y
r

sfparm = mark:1:1:1;sgnl<0
colors = 16;0
clear = n
r

In this case, there really isn't much new in 5.6.

However, you could run the program just once for color coding the dots
(orange negative, purple positive) like:
sfparm = mark:15:4:1
colors = (0/16;30/sgnl)
clear = y
r

The above will color code the marker based on the value of SGNL.
The SGNL values less than or equal to 0 will set the MARK to color 16.
The SGNL values geater than 0 will set the MARK to color 30. The <= or >=
can be configured in COLORS with /L or /U.
Of course, that is only useful if you don't mind using the same marker type.

Steve Chiswell
Unidata User Support



>From: David Knight <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200103142143.f2ELh0L15194

>
>Hi Chiz,
>     Oh. OK. I see MARK as a new parameter in sfparm.
>But all the help for sfparm says is
> MARK - Markers
>Not very helpful ;-(  Can you tell me how to use this parameter?
>Say for example that I want to make a purple dot for pos flashes
>and a orange plus sign for negative flashes, centered where the
>strike occured.
>
>You were right, sfmap can plot the lightning data.
>Now I just have to figure out how to get it to just
>show the strikes.
>
>How are you doing this (or are you)
>
>Thanks for your help. Looks like there are going to be a few
>adjustments in transitioning to 5.6
>
>David
>> 
>> 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 rea
> d
>> > > > 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)
>> > 
>> > 
>> 
>> 
>