Thanks for your patience. I got the right image after adding =
data(::-1,: ).
Can I define a specific domain ? The following sentence do not work for =
me.
res@mpMinLonF = 120
res@mpMaxLonF = 145
res@mpMinLatF = 30
res@mpMaxLatF = 60
On Jun 11, 2014, at 1:07, Wei Huang <huangwei@ucar.edu> wrote:
> Jinqing,
>
> This is a satellite image, we have to tune to make things look right.
> Attached are a script, and a plot generated.
> You probably need to consult with a satellite image expert to see
> what those parameters mean.
>
> Regards,
>
> Wei
>
>
>
> On Tue, Jun 10, 2014 at 8:50 AM, 刘金卿 =
<liujq_0912@126.com> wrote:
> 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>
>
>
>
> =
邮件带有附件预览=
接，若您转发或=
复此邮件时不希望=
对方预览附件，建=
您手动删除链接=
> 共有 2 个附件
> T.png(450K)
> 极速下载 在线预览=
> fy2c.ncl(4K)
> 极速下载
> <T.png><fy2c.ncl>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jun 10 10:54:37 2014
This archive was generated by hypermail 2.1.8 : Wed Jul 23 2014 - 15:33:46 MDT