Re: NCEP stageIV grib data plot

Dave,

Thanks for your reply.

Now I know lat2d@Lov and lon2d@Lov.
ncl_filedump tells everything in grib file.

Yongzuo

________________________________
From: David Brown [dbrown@ucar.edu]
Sent: Thursday, January 14, 2010 1:17 PM
To: Li, Yongzuo
Cc: Dennis Shea; ncl-talk@ucar.edu
Subject: Re: NCEP stageIV grib data plot

On Jan 14, 2010, at 10:09 AM, Li, Yongzuo wrote:

Dennis,

I made an almost perfect plot
ftp://ftp.caps.ou.edu/users/yli/precip-ncep-stageIV.pdf
using
ftp://ftp.caps.ou.edu/users/yli/ncep-st4.ncl
with
ftp://ftp.caps.ou.edu/users/yli/ST4.2008091307.01h.grib

Questions:

(1) Is my "ZOOM" method good?

(2) The map tilts to right because the center longitude
is -105
res@mpProjection = "Stereographic"
res@mpLimitMode = "Corners"
res@mpCenterLonF = lat2d@Lov
res@mpCenterLatF = 90.

(3) What does lat2d@Lov, the reference longitude, come?

Hi Yongzuo,
The Lov value comes from the grid description in the GRIB file. See the table for polar stererographic grids at:
http://www.nco.ncep.noaa.gov/pmb/docs/on388/tabled.html
 -dave

Yongzuo
________________________________________
From: Dennis Shea [shea@ucar.edu<mailto:shea@ucar.edu>]
Sent: Thursday, January 14, 2010 8:54 AM
To: Li, Yongzuo
Cc: ncl-talk@ucar.edu<mailto:ncl-talk@ucar.edu>
Subject: Re: NCEP stageIV grib data plot

If you did a

%> ncl_filedump ST4.2008091307.01h.grib

you would see that when NCL reads the grib file it
adds the long_name anf units attributes.

prc = a->A_PCP_GDS5_SFC_acc1h

delete(prc@long_name)
delete(prc@units)

Given your plot, maybe

"NCEP stageIV 1-h precip valid 20080913_"+HH
delete(prc@long_name)
delete(prc@units)
res@gsnCenterString = "1-h precip (mm)"

see gsnCenterString resources at:

http://www.ncl.ucar.edu/Document/Graphics/Resources/gsn.shtml

====
If you want to add a label to the label bar, see

http://www.ncl.ucar.edu/Applications/labelbar.shtml
Example 4

Good Luck

Li, Yongzuo wrote:
> Hi,
>
> I use
>
> ftp://ftp.caps.ou.edu/users/yli/ncep-st4.ncl
>
> to make
>
> ftp://ftp.caps.ou.edu/users/yli/precip-ncep-stageIV.pdf
>
> How can I remove "Total precipitation" and add description to color bar?
>
> I want to add subtitle.
>
> Thanks
>
> Yongzuo
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Jan 14 15:00:45 2010