Hi, Wei,
Thanks for your reply, but the picture still seems wrong compared with =
the image from website quick view (attachment).
And I add
res@mpDataSetName = "Earth..4"
res@mpDataBaseVersion = "MediumRes"
res@mpOutlineSpecifiers = (/"China:states","India","Taiwan”/=
)
in the script, but it dose not work. Is there anyway to turn on the =
boundaries for this data?
On Jun 10, 2014, at 0:41, Wei Huang <huangwei@ucar.edu> wrote:
> Jinqing,
>
> My may try this script:
>
> ;************************************************
> fn = "FY2C_TBB_IR1_NOM_20070816_1200.hdf.h5"
> f = addfile(fn, "r")
> vNam = getfilevarnames(f) ; all variables names =
on file
> print(vNam)
> ;************************************************
> ;pltType = "x11"
> pltType = "png"
> pltName = "T"
> ;************************************************
> ; plotting parameters
> ;************************************************
> wks = gsn_open_wks (pltType,pltName) ; open workstation
> ;gsn_define_colormap(wks,"amwg") ; choose colormap
> ;gsn_define_colormap(wks,"gsdtol")
> ;****************************************************
> res = True ; plot mods desired
> res@gsnMaximize = True ; make ps/eps/pdf large
> res@gsnSpreadColors = True
>
> res@cnFillOn = True
> ;res@cnLinesOn = True
> ;res@cnLineLabelsOn = True
> res@cnFillMode = "RasterFill" ; faster
> ;res@cnRasterSmoothingOn = True
> ;****************************************************
> x=f->FY2C_TBB_Hourly_Product(:,:)
>
> x = where(x.lt.150.0, -999.0, x)
> x@_FillValue = -999.0
>
> printVarSummary(x)
> plot = gsn_csm_contour(wks, x, res)
>
>
> Wei
>
>
>
> On Mon, Jun 9, 2014 at 4:56 AM, 刘金卿 =
<liujq_0912@126.com> wrote:
> Hi, I try to plot a hdf file, but got the following warnings:
> warning:_NhlCreateSplineCoordApprox: Attempt to create spline =
approximation for X axis failed: consider adjusting trXTensionF value
> warning:IrTransInitialize: error creating spline approximation for =
trXCoordPoints; defaulting to linear
>
> The picture seems wrong too. The data has been transferred on the ftp, =
and the script is in attachment. Any suggestions will be appreciate.
> ftp> put FY2C_TBB_IR1_NOM_20070816_1200.hdf
>
> -Liu
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> =
邮件带有附件预览=
接，若您转发或=
复此邮件时不希望=
对方预览附件，建=
您手动删除链接=
> 共有 1 个附件
> T.png(458K)
> 极速下载 在线预览=
> <T.png>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jun 10 02:51:01 2014
This archive was generated by hypermail 2.1.8 : Wed Jul 23 2014 - 15:33:46 MDT