Re: about polar graph?

[Jonathan Callahan <callahan@xxxxxxxxxxxxx>]

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
> =====================================================