Thank you, you have found a bug. poly/coord_ax=Z/key/over/nolab/lev=30/pal=default jlon+PLMxpolymark*sym_size, jlat+PLMypolymark*sym_size, tt This command is set in the Ferret script server/jnls/insitu/insitu_poly_xy.jnl -------------------- poly/coord_ax=Z/key/over/nolab/lev=$3"50"/pal=$4"rnb2" jlon+PLMxpolymark*sym_size, jlat+PLMypolymark*sym_size, $2 IF `($xstart) LT 0` THEN poly/coord_ax=Z/key/over/nolab/lev=$3"50"/pal=$4"rnb2" jlon-360+PLMxpolymark*sym_size, jlat+PLMypolymark*sym_size, $2 ENDIF IF `($xend) GT 360` THEN poly/coord_ax=Z/key/over/nolab/lev=$3"50"/pal=$4"rnb2" jlon+360+PLMxpolymark*sym_size, jlat+PLMypolymark*sym_size, $2 ENDIF -------------------------- $2 refers to the variable name 'tt', but in insitu_setup.jnl, the variable has been subsampled due to it's size and renamed 'jtt'. that is why the number of values(18547) does not match the number of polygons(4637). Notice the other arguments to the polygon command include jlon and jlat instead of xax and yax - these are our subsampled axis variables at the correct length of 4637. SOLUTION: change the Ferret script arg reference $2 to j$2 (and the correct variable will be used in the polygon command. poly/coord_ax=Z/key/over/nolab/lev=$3"50"/pal=$4"rnb2" jlon+PLMxpolymark*sym_size, jlat+PLMypolymark*sym_size, j$2 I will submit a bug report and correct the scripts. Thank you, Joe ------------------------------------------- cch wrote:
|