Wanhong,
As I said, the following method will perform OA on the surface file and you
should be able to plot a contour and/or fill plot, but I have not been able to
plot point data (CTYPE = P in GDCNTR) yet, and I'm not really sure why.
in OAGRID
GDFILE = oagrid.gem
DELTAN = 0.40
DELTAX = 0.20
DELTAY = 0.20
GRDAREA = 36;-90;39.25;-81.75
EXTEND = 2;2;2;2
DTAAREA =
SOURCE = SF
SNFILE =
SFFILE = dev.gem
SNPARM =
SFPARM = ptot
DATTIM = /0800
LEVELS = 10
MAXGRD = 10
r
you should see a file in your present working directory called oagrid.gem (mine
is 4.6k)
in OABSFC:
SFFILE = dev.gem
GDFILE = oagrid.gem
SFPARM = ptot
DATTIM = /0800
DTAAREA =
GUESS =
GUESFUN =
GAMMA = 0.3
SEARCH = 20
NPASS = 2
QCNTL =
OABND =
GDATTIM = last
GFUNC = ptot
GLEVEL = 0
GVCORD = none
r
oagrid.gem should now be about 6.1k.
to confirm with GDINFO:
GDFILE = oagrid.gem
LSTALL = YES
OUTPUT = T
GDATTIM = last
GLEVEL = 0
GVCORD = none
GFUNC = ptot
r
you should see:
Number of grids in file: 1
Maximum number of grids in file: 10
NUM TIME1 TIME2 LEVL1 LEVL2 VCORD PARM
1 110208/0800 0 NONE PTOT
in GDCNTR:
GDATTIM = last
GLEVEL = 0
GVCORD = none
GFUNC = ptot
GDFILE = oagrid.gem
CINT =
LINE = 3
MAP = 10/1/2
MSCALE = 0
TITLE =
DEVICE = xw
SATFIL = 0
RADFIL = 0
IMCBAR =
PROJ = mer
GAREA = 36;-90;39.25;-81.75
IJSKIP = 0
CLEAR = y
PANEL = 0
TEXT = 1/21//hw
SCALE = -1
LATLON =
HILO =
HLSYM =
CLRBAR = 24/H/LL//.80;.02
CONTUR = 0
SKIP =
FINT = 0
FLINE = 10-20
CTYPE = c
LUTFIL = 0
STNPLT =
r
I've attached a gif image of what I see (with CTYPE = C)
Again, when I change to CTYPE = P it's not plotting grid point values as it
should. I'm still trying to figure out why.
Michael James
Unidata
> Hi Michael,
> I really appreciate your help. We will remove numbers on maps if the
> floating data
> can not be ploted on the map using sfmap. Please email me the scripts
> using OAGRID and
> OABSFC if you could(even you don't finish that script), I will test that
> and we will
> use that later.
>
> Thanks.
>
> Wanhong Wang
> System Admin/Prog.
> UK Weather Center
>
>
>
>
> On 2011-2-9 21:39, Unidata GEMPAK Support wrote:
> > Unfortunately floating points are not able to be plotted in SFMAP.
> >
> > Two solutions I have are 1) use SFPARM = ptot*100 to display precip to
> > hundredths of inches but without the decimal point, or 2) use OAGRID and
> > OABSFC to grid the surface precip obs to plot on top of contoured model
> > forecast values.
> >
> > I've done option 2 here on my end but can't seem to display point values
> > (TYPE = P) for the PTOT grid and I'm not sure why. If option 1 works for
> > you then no bother, otherwise I'll have to look at this again tomorrow.
> >
> > best,
> >
> > Michael
> >
> >
> >
> >> Hi Michael James,
> >>
> >> I got error message when try to upload the files. I put all file
> >> on following link, you may get all you want.
> >>
> >> http://weather.uky.edu/ldm-images/shared/
> >>
> >> dev.gem 09-Feb-2011 18:59 22K
> >> dev.grd 09-Feb-2011 18:59 38K
> >> gemglb.nts 09-Feb-2011 18:59 22K
> >> kypm30d_prec.fin 09-Feb-2011 18:59 2.2K
> >> kypm30d_prec.gif 09-Feb-2011 19:00 29K
> >> kypm_prec.csh 09-Feb-2011 19:00 3.5K
> >> kypmsites.tbl 09-Feb-2011 19:00 5.4K
> >>
> >>
> >> We got map, but we want print the rainfall real data on the map,
> >> not as integer...
> >>
> >> http://weather.uky.edu/ldm-images/shared/kypm30d_prec.gif
> >>
> >>
> >>
> >>
> >>
> >>
> >> -------- Original Message --------
> >> Subject: Fwd: [GEMPAK #ZGN-126493]: Fwd: Fwd: KY PM 7 Days data and map
> >> Date: Wed, 9 Feb 2011 18:26:57 -0500
> >> From: Priddy, Kenny T<address@hidden>
> >> To: Wang, Wanhong<address@hidden>
> >>
> >>
> >>
> >>
> >>
> >> Sent from my iPhone
> >>
> >> Begin forwarded message:
> >>
> >>> *From:* Unidata GEMPAK Support<address@hidden
> >>> <mailto:address@hidden>>
> >>> *Date:* February 9, 2011 2:46:55 PM EST
> >>> *To:* "Priddy, Kenny T"<address@hidden<mailto:address@hidden>>
> >>> *Cc:* "address@hidden
> >>> <mailto:address@hidden>"
> >>> <address@hidden<mailto:address@hidden>>
> >>> *Subject:* *[GEMPAK #ZGN-126493]: Fwd: Fwd: KY PM 7 Days data and map*
> >>> *Reply-To:* "address@hidden
> >>> <mailto:address@hidden>"
> >>> <address@hidden<mailto:address@hidden>>
> >>>
> >>> I'm not sure if you attached the mesonet obs grid (dev.grd?) or if it
> >>> was stripped by our mail server, which happens occasionally.
> >>>
> >>> Can you upload the file for me at the GEMPAK support repository?
> >>>
> >>> http://motherlode.ucar.edu/repository/alias/gempakuploads/
> >>>
> >>> Thanks,
> >>>
> >>> Michael James
> >>> Unidata
> >>>
> >>>
> >>>> Can y'all help? I tried our local authority....Gilbert....but no luck.
> >>>> We want to print rainfall values of our Mesonet sites using a gempak
> >>>> script...
> >>>> but all we can get to print is interger values of rainfall. Is that one
> >>>> of the parameters
> >>>> in the script? Any help is appreciated.
> >>>>
> >>>> -------- Original Message --------
> >>>> Subject: Fwd: KY PM 7 Days data and map
> >>>> Date: Wed, 09 Feb 2011 12:52:03 -0500
> >>>> From: Tom Priddy<address@hidden<mailto:address@hidden>>
> >>>> To: Gilbert Sebenste<address@hidden
> >>>> <mailto:address@hidden>>, wwang
> >>>> <address@hidden<mailto:address@hidden>>
> >>>>
> >>>>
> >>>>
> >>>> From: Tom Priddy
> >>>>
> >>>> Hey, Gilbert. Here's a question we've been scratching our heads on....
> >>>>
> >>>> We want to print rainfall values of our Mesonet sites using a gempak
> >>>> script...
> >>>> but all we can get to print is interger values. Is that one of the
> >>>> parameters
> >>>> in the script?
> >>>>
> >>>>
> >>>> -------- Original Message --------
> >>>> Subject: KY PM 7 Days data and map
> >>>> Date: Wed, 9 Feb 2011 12:44:22 -0500
> >>>> From: Wang, Wanhong<address@hidden<mailto:address@hidden>>
> >>>> To: Priddy, Kenny T<address@hidden<mailto:address@hidden>>
> >>>>
> >>>>
> >>>>
> >>>>
> >>>>
> >>>>
> >>>> how to change parameters to print a number(1.22) as decimal
> >>>> fraction(1.22 or 1.2) not as integer(1)?
> >>>>
> >>>>
> >>>> The maps:
> >>>>
> >>>>
> >>>>
> >>>>
> >>>>
> >>>> the programs used for creating maps:
> >>>>
> >>>> gdcntr<< EOF6
> >>>> \$RESPOND = no
> >>>> \$MAPFIL = hipowo.cia + hicnus.gsf
> >>>> GDATTIM = /0800
> >>>> GLEVEL = 0
> >>>> GVCORD = none
> >>>> GFUNC = ptot
> >>>> GDFILE = dev.grd
> >>>> CINT = ${int}
> >>>> LINE = 2/1
> >>>> MAP = 10/1/2 + 10/10/1
> >>>> TITLE = 24/-2/${ddate2} - ${titleinfo} (inches)
> >>>> DEVICE = ${device}
> >>>> SATFIL =
> >>>> RADFIL =
> >>>> PROJ = mer
> >>>> GAREA = 36;-90;39.25;-81.75
> >>>> CLEAR = yes
> >>>> PANEL =
> >>>> TEXT = 1/21//hw
> >>>> SCALE = 0
> >>>> LATLON = 0
> >>>> HILO = 0
> >>>> HLSYM = 0
> >>>> CLRBAR = 24/H/LL//.80;.02
> >>>> CONTUR = 0
> >>>> SKIP =
> >>>> FINT = ${int}
> >>>> FLINE = 101;21-30;4
> >>>> CTYPE = c/f
> >>>> LUTFIL = 0
> >>>> STNPLT = 0
> >>>> run
> >>>>
> >>>> exit
> >>>> EOF6
> >>>>
> >>>> sfmap<< EOF5a
> >>>> \$RESPOND = no
> >>>> \$MAPFIL = hipowo.cia + hicnus.gsf
> >>>> AREA = 36;-90;39.25;-81.75
> >>>> GAREA = 36;-90;39.25;-81.75
> >>>> SATFIL = 0
> >>>> RADFIL = 0
> >>>> SFPARM = ptot
> >>>> DATTIM = /0800
> >>>> SFFILE = dev.gem
> >>>> COLORS =
> >>>> MAP = 10/1/2 + 10/10/1
> >>>> LATLON =
> >>>> MARKER = 0
> >>>> TITLE = 24/-1/${titleinfo2}
> >>>> CLEAR = no
> >>>> PANEL = 0
> >>>> DEVICE = ${device}
> >>>> PROJ = mer
> >>>> FILTER = no
> >>>> TEXT = 1/21//hw
> >>>> SKPMIS = 0
> >>>> LUTFIL = 0
> >>>> STNPLT = 0|31/5/1.25/1.25|${tbl}
> >>>> run
> >>>>
> >>>> exit
> >>>> EOF5a
> >>>>
> >>>>
> >>>> attached data file
> >>>>
> >>>>
> >>>>
> >>>>
> >>>>
> >>>
> >>> Ticket Details
> >>> ===================
> >>> Ticket ID: ZGN-126493
> >>> Department: Support GEMPAK
> >>> Priority: Normal
> >>> Status: Open
> >>>
> >>
> >
> > Ticket Details
> > ===================
> > Ticket ID: ZGN-126493
> > Department: Support GEMPAK
> > Priority: Normal
> > Status: Open
> >
>
>
Ticket Details
===================
Ticket ID: ZGN-126493
Department: Support GEMPAK
Priority: High
Status: OpenAttachment:
meso.gif
Description: GIF image
