[Thread Prev][Thread Next][Index]
Vectors on polar stereographic plots
I am currently trying to produce some polar stereographic plots with
vectors. I have produced the same plots in a 'normal' projection and
had no problems. I have also been able to produce 'fill' plots on the
stereographic plots with no problems. However, when producing vector
plots I only get arrows from regions where my measured field has its
highest values. At some point in the code it seems that the small
values are being filtered out (and they're not that much smaller):
set mode metafile
use u10_2000_01_djf.cdf
set grid u10
let pi=4*atan(1)
let wind=((u10^2)+(v10^2))^0.5
let ang=atan(v10/u10)
let dir=if u10 gt 0 then ((pi/2)-ang) else (((3*pi)/2)-ang)
let Ts=(5.1e-4)*(wind^3)
let Tsx=Ts*sin(dir)
let Tsy=Ts*cos(dir)
set window/aspect=2
set viewport lower
set region/x=0:360/y=-90:0/l=1:90
go mp_stereographic_south
!fill/levels=(0,1.5,0.1)/nolab/title="Seasonally-Averaged Integrated
Stokes Transport for DJF 2000/01 based on the math of McWilliams and
Restrepo, Southern Hemisphere" Ts[t=0:1@ave]*mp_mask, x_page, y_page
vector/length=0.5 Tsx[t=0:1@ave]*mp_mask,
Tsy[t=0:1@ave]*mp_mask,x_page,y_page
go mp_fland
label/nouser `($PPL$XLEN)*0.5`,-0.4,0,0,.05 Seasonally-Averaged
Integrated Stokes Transport for DJF 2000/01 based on the math of
McWilliams and Restrepo, Southern Hemisphere
set viewport upper
set region/x=0:360/y=0:90/l=1:90
go mp_stereographic_north
!fill/levels=(0,1.5,0.1)/nolab/title="Seasonally-Averaged Integrated
Stokes Transport for DJF 2000/01 based on the math of McWilliams and
Restrepo, Northern Hemisphere" Ts[t=0:1@ave]*mp_mask, x_page, y_page
vector/length=0.1 Tsx[t=0:1@ave]*mp_mask, Tsy[t=0:1@ave]*mp_mask,
x_page, y_page
go mp_fland
label/nouser `($PPL$XLEN)*0.5`,-0.4,0,0,.05 Seasonally-Averaged
Integrated Stokes Transport for DJF 2000/01 based on the math of
McWilliams and Restrepo, Northern Hemisphere
cancel mode metafile
Hope someone can help!!
Thank you,
Stephen Nash
--
__________________________________
Stephen Nash
Department of Meteorology
The University of Reading
Early Gate
PO BOX 243
Reading, RG6 6BB
s.d.nash@reading.ac.uk
07766075885
__________________________________
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement