Re: can't see the global map

On Fri, 14 Aug 2009, [gb2312] Ëﳬ wrote:

> Dear friends
> I have finished the simulation.But when I use the NCL to plot fields from the
> global WRF run .I still could only got a thin long bar which appear to be the
> meridian and message as followed:
> "
> fatal:VectorPlotDraw: VVECTR - VECTOR NDC LENGTH TOO GREAT
> fatal:VectorPlotDraw: error drawing vectors
> fatal:VectorPlotDraw: draw error
> fatal:PlotManagerDraw: error in plot draw
> fatal:_NhlPlotManagerDraw: Draw error
> fatal:VectorPlotDraw: VVECTR - VECTOR NDC LENGTH TOO GREAT
> fatal:VectorPlotDraw: error drawing vectors
> fatal:VectorPlotDraw: draw error
> fatal:PlotManagerDraw: error in plot draw
> fatal:_NhlPlotManagerDraw: Draw error
> "
>
> I have already reduced the time to 1s and set the fft_filter_lat to be 90.I
> checked the wrfoutput data ,they seemed to be right.
> And I tried another script(The script is attached to the letter.) which doesn't
> have the vectorplotdraw , then I didn't recieve error message.
> But I just can only see a thin line as before.
> Dose it mean that there is something else caused I can't see the whole map? I
> compiled the ncl by myself.Whether the compiling caused the error? What should I
> do ?

This kind of problem is hard to debug without being able to look at
your data. There's something fundamentally wrong here, and it could be
any number of things, like a problem with the wrfout file, a problem
with the wrf_xxx scripts, or a simple user error somewhere.

For one, if you are not getting a map plot, global map parameters on
the file may be wrong.

Second, the errors you are getting about the vector data makes me
believe something is wrong with the vector data.

For debugging purposes, I would create a *very simple* program to make
sure it is not your map parameters or vector data that's the problem.

This program should draw three frames:

   1. A map plot that uses wrf_map.
   2. A vector plot that uses gsn_vector to draw the vectors.
   3. A vector plot that uses wrf_vector to draw the vectors.

The third frame may look odd, because "wrf_vector" is not really
meant to be used in this fashion.

Here's the program:

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"

begin
   a = addfile("./wrfout.nc","r")
   u10 = wrf_user_getvar(a,"U10",1)
   v10 = wrf_user_getvar(a,"V10",1)

   wks = gsn_open_wks("x11","test")

   res = True

; First frame
   map = wrf_map(wks,a,res)
   draw(map)
   frame(wks)

; Second frame
   plot = gsn_vector(wks,u10,v10,res)

; Third frame
   vector = wrf_vector(a,wks,u10,v10,res)
   draw(vector)
   frame(wks)

end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Aug 14 2009 - 10:23:11 MDT