[Thread Prev][Thread Next][Index]

Fw: poly_vectors.jnl



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


[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement