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 fi=
le
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, an=
d
> 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
>
>
Received on Mon Jun 09 04:42:18 2014
This archive was generated by hypermail 2.1.8 : Wed Jul 23 2014 - 15:33:46 MDT