[Thread Prev][Thread Next][Index]
Re: about polar graph?
Kun,
I'd recommend using the new map projection functionality which is
available in ferret version 5 and higher. Here's a simple script to do
what you want. Information on map projections is at:
http://ferret.wrc.noaa.gov/Ferret/Documentation/Users_Guide/v50/chapter6.html#V5.0_projection
set win/size=.3
use coads_climatology
set region/x=-180:180/y=45:90/l=1
go mp_stereographic_north
set grid sst
go mp_aspect
! create a blank underlay which is slightly larger
! (y=45:90) so that there is room to draw the
! labels afterwards
shade/noaxes/nokey/pal=white sst, x_page, y_page
! now for the region we're interested in
set region/x=-180:180/y=50:90/l=1
shade/over/noaxes/key/pal=no_red sst, x_page, y_page
go mp_fland
go mp_graticule 0 360 30 50 90 5 1
! label positions tweaked to look good
go mp_label 180 49 0 0 .12 180^o
go mp_label -120 45 0 0 .12 "120^o W"
go mp_label -60 45 0 0 .12 "60^o W"
go mp_label 0 47 0 0 .12 0^o
go mp_label 60 46 0 0 .12 "60^o E"
go mp_label 120 46 0 0 .12 "120^o E"
go mp_label 45 50.5 0 45 .12 "50^o N"
go mp_label 45 60.5 0 45 .12 "60^o N"
-- Jonathan Callahan
Yang Kun wrote:
>
> Dear Ferreters,
>
> Just a simple question. But I didn't know how to modify
> the script.
> I want to plot a polar graph for some variables in North
> Pole & South Pole. Below is my script. and it works well.
> Question is: i just want to see 50s:90s & 50n:90n in polar
> graphs, but not the whole earth. I tried as below, but it
> doesn't works. What should i do? I think it is only a minor
> problem, but i can not figure it out.
>
> Thanks!
> --Kun Yang
> +++++++++++++++++++++++++++++++++++++++++++++
> set data test_new
> cancel region/x
> set region/l=1
> set window 1
> set region/y=50n:90n
> set mode metafile sicthe1.plt
> set wind/aspect=1:ax/size=.5
> go polar_2d fill/pal=no_red sicthe[y=50n:90n,l=1:12@ave] 100
> go polar_fland 60 gray
> go polar_land " " 1.0 ! limit lines to a radius of 1
> go polar_grid_fancy
> cancle mode metafile
> set window 2
> set region/y=50s:90s
> set mode metafile sicthe2.plt
> set wind/aspect=1:ax/size=.5
> go polar_2d fill/pal=no_red sicthe[y=50s:90s,l=1:12@ave] 100,,,,,-1
> go polar_fland 60 gray
> go polar_land green,1.0
> go polar_grid_fancy
> cancle mode metafile
>
> =====================================================
> Kun Yang, Researching Assistant
> Dept. of Oceanography, Texas A&M Uni. College Station
> Mailing Addr: 806F, Oran Circle, Bryan, TX, 77801
> Tel: (409)-268-0963 (h), (409)-845-3238 (o)
> Email addr: yangk@ocean.tamu.edu; yangk@netease.com
> =====================================================
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement