Re: Why wind vectors could not be plotted in a sub-area with the same ncl codes?

Hi Dennis:

Thanks for your help.

I have ftped the file to the specified directory. It is a compressed file: xadala_pi00002011050215_utc.nc.gz

Would you please let me know if you cannot find the file? Thanks!

Ed

--- On Fri, 7/6/12, Dennis Shea <shea@ucar.edu> wrote:

> From: Dennis Shea <shea@ucar.edu>
> Subject: Re: Why wind vectors could not be plotted in a sub-area with the same ncl codes?
> To: "Yang Yang" <yang123yang@yahoo.com>
> Cc: ncl-talk@ucar.edu
> Date: Friday, July 6, 2012, 7:22 AM
> Nothing obviously wrong see 'diff'
> below'.
>
>
> It might be best to send the file to:
>
> ftp ftp.cgd.ucar.edu
> anonymous
> <use your email address for the password>
> cd incoming
> put xadala_pi00002011050215_utc.nc
> quit
>
> Send email to ncl-talk
> after you have successfully completed the transfer.
>
> ================================
> Rather than a double do loop
>
>     do i=100, 500
>     do j=100, 500
>    
>    xx=(v10m(i,j+1)-v10m(i,j-1))/200.
>    
>    yy=(u10m(i+1,j)-u10m(i-1,j))/200.
>     vor(i,j)=xx-yy
>     end do
>     end do
>
> Why not use the builtin func
>     dd  = 200.0
>     xx  = center_finite_diff_n
> (v10m,200,False,0, 1)
>     yy  = center_finite_diff_n
> (u10m,200,False,0, 0)
>     vor = xx-yy
> ================================
> %> diff plot_1.ncl plot_2.ncl
>
> 91c91
> <     res1@gsnLeftString = "Plot_1 "
> ---
> >     res1@gsnLeftString = "Plot_2 "
> 129c129
> <     wks = gsn_open_wks
> ("eps","plot_1")        ; open
> workstation
> ---
> >     wks = gsn_open_wks
> ("eps","plot_2")        ; open
> workstation
> 139,142c139,142
> <   res1@mpLeftCornerLatF   
>    = -36.85
> <   res1@mpLeftCornerLonF   
>    = 174.54
> <   res1@mpRightCornerLatF   
>   = -36.6
> <   res1@mpRightCornerLonF   
>   = 174.85
> ---
> >   res1@mpLeftCornerLatF   
>    = -36.75
> >   res1@mpLeftCornerLonF   
>    = 174.64
> >   res1@mpRightCornerLatF   
>   = -36.65
> >   res1@mpRightCornerLonF   
>   = 174.75
> 149c149
> <     inc = 10     
>                 ;
> plot at every "inc" grid point
> ---
> >     inc = 2     
>                 ;
> plot at every "inc" grid point
> =================================
>
>
> On 7/5/12 3:44 PM, Yang Yang wrote:
> > Hi:
> >
> > I posted the issue two days ago. Perhaps I didn't
> describe it well. I have not received any response. So I
> repost it here. Sorry!
> >
> > I used the same ncl codes to plot wind vectors in an
> area (~30x30km) and a sub-area (~10x10km) of it. The only
> differences are the parameters controlling the area to plot.
> However, the wind vectors can be plotted for the large area
> (this means my data are OK) but not on the sub-area. Can
> anyone help me fix the problem?
> >
> > plot_1.png was plotted by plot_1.ncl for the large
> area. plot_2.png was plotted for the sub-area by
> plot_2.ncl.The resolution of the data is 100m.
> >
> > Thank you very much!
> >
> > Ed
> >
> > --- On Wed, 7/4/12, Yang Yang<yang123yang@yahoo.com> 
> wrote:
> >
> >> From: Yang Yang<yang123yang@yahoo.com>
> >> Subject: Odd plot of wind vectors
> >> To: ncl-talk@ucar.edu
> >> Date: Wednesday, July 4, 2012, 8:14 AM
> >> Hi:
> >>
> >> I have a regional model output (100m resolution).
> When I
> >> ploted an area about 30x30km (plot_1.png attached)
> using
> >> plot_1.ncl, the plot seems OK. However, when I
> tried to plot
> >> an area about 10x10km within the 30x30km area, the
> plotted
> >> wind vectors are quite odd (plot_2.png by
> plot_2.ncl), and
> >> wind vectors were not plotted in most areas. Can
> anyone help
> >> me fix the problem? Thanks!
> >>
> >> Ed
> >> -----Inline Attachment Follows-----
> >>
> >> _______________________________________________
> >> 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 Jul 5 22:57:17 2012