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"
Attachment:
vect_sst.ps
Description: PostScript document