Re: Fw: poly_vectors.jnl

[Ansley Manke <ansley.b.manke@xxxxxxxx>]

Hi Gabriel,
Are there missing values in the temperature variable? If the velocity 
field exists but sst is missing then you would see arrows but  they 
would be unfilled.You can check this by comparing what the plot looks 
like if you use a command like

  missing(sst,0)

instead of sst in the script: At the end of your script:

LET press = YSEQUENCE( MISSING( sst[X=-70:20:2,Y=-40:10:2], 0) )


Gabriel Clauzet wrote:

> > Hi all,
> > I´m using the poly_vectors.jnl  function to plot the surface current 
> > fields
> > colored according the SST for the South Atlantic region.
> > My problem is that the function does not color some arrows. I tried to
> > change the number of arrows, the resolution of my data but some 
> > arrows are
> > still blank.
> > Does anyone have some idea ??
> >
> > here is my script:
> >
> > use "/data/dados/UVELcli350-399.nc"
> > use "/data/dados/VVELcli350-399.nc"
> > use "/data/dados/TEMPcli350-399.nc"
> >
> > let u=uvel[k=1,l=@ave,d=1]
> > let v=vvel[k=1,l=@ave,d=2]
> > let sst=temp[k=1,l=@ave,d=3]
> >
> > SET WINDOW/siz=0.6
> >
> > SET REGION/X=-70:20/Y=-40:10
> >
> > GO basemap x=-70:20 y=-40:10 20
> > LABEL/NOUSER `($ppl$xlen)/2`, `($ppl$yorigin)-1`, 0, 0, 0.16, "U,V"
> >
> > LET x_vec = X[GX=u,X=-70:20:2] + 0*Y[GY=u,Y=-40:10:2]
> > LET y_vec = 0*X[GX=u,X=-70:20:2] + Y[GY=u,Y=-40:10:2]
> > LET u_vec = uu1[X=-70:20:2,Y=-40:10:2]
> > LET v_vec = vv[X=-70:20:2,Y=-40:10:2]
> >
> > GO poly_vectors.jnl x_vec y_vec u_vec v_vec 8 "arrow"
> >
> > LET press = YSEQUENCE(sst[X=-70:20:2,Y=-40:10:2])
> > POLYGON/OVER/NOAX/NOLAB/KEY/LINE/PAL=light_rainbow plt_x_arrow,
> > plt_y_arrow, press
> >
> > GO poly_arrow_key "Surface Currents", "8 m/s"