Re: mp_stereographic

[John Donners <donners@xxxxxxx>]

Hi Gabriel,

I had a quick look at the script mp_stereographic_south.jnl, and if 
you change the two following lines, I think you get the representation
as 
you would like it to be:

let/quiet x_page = (-1.0) * mp_rho * sin(mp_lambda-mp_lambda0) ! eqn
(108)
let/quiet y_page = (-1.0) * mp_rho * cos(mp_lambda-mp_lambda0) ! eqn
(109)

into

let/quiet x_page = mp_rho * sin(mp_lambda-mp_lambda0) ! eqn (108)
let/quiet y_page = mp_rho * cos(mp_lambda-mp_lambda0) ! eqn (109)

With regards,

John

On Thu, 24 Jul 2003, Gabriel Clauzet wrote:

> Dear Jon and Ferreters,
> 
> I condenced my script using the ETOPO dataset, the script run well but
I can't invert the projection, to have the Antarctica at the bottom of
the figure and the equator at the top. If any one have an idea to help
me please tell me.
> 
> Best Regards 
> Gabriel
> 
> !!!!!!!!!!!!! PLOT SOUTH ATLANTIC OCEAN !!!!!!!!!!!!!!!!!!!!!!!!!!!!
> 
> SET WINDOW/SIZE=0.6
> USE etopo40
> SET GRID rose
> 
> ! Define a region in the southern hemisphere.(80w:0e/10s:80s)
> LET xmin = 280; LET xmax = 360
> LET ymin = -80; LET ymax = -10
> LET xdel = 10; LET ydel = 10
> 
> LET mid_lon = `xmin+(xmax-xmin)/2`
> LET mid_lat = `ymin+(ymax-ymin)/2`
> SET REGION/X=`xmin`:`xmax`/Y=`ymin`:`ymax`
> 
> GO mp_stereographic_south `mid_lon` `mid_lat`
> LET masked_rose = rose * mp_mask
> GO mp_aspect
> 
> FILL/NOAXES/NOLAB/PALETTE=land_sea/SET_UP masked_rose, x_page, y_page
> PPL SHAKEY 1, 0, 0.1, 2, , , 2.5, 8.5, 7.5, 8.2
> PPL FILL
> 
> GO mp_land
> 
> GO mp_graticule `xmin` `xmax` `xdel` `ymin` `ymax` `ydel` 
> 
> REPEAT/K=`xmin`:`xmax`:`xdel` (GO mp_label `k` `ymax+2` 0 `mid_lon-k`
0.12 "@CR`k`")
> REPEAT/K=`-1*(ymin+ydel)`:`-1*(ymax-ydel)`:`2*ydel` (GO mp_label
`xmin-2` `-1*k` -1 `mid_lon-xmin` 0.12 "@CR`-1*k`")
> REPEAT/K=`-1*(ymin+ydel)`:`-1*(ymax-ydel)`:`2*ydel` (GO mp_label
`xmax+2` `-1*k`  1 `mid_lon-xmax` 0.12 "@CR`-1*k`")
> 
> !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!    END
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
>  
>   ----- Original Message ----- 
>   From: Jonathan Callahan 
>   To: Gabriel Clauzet 
>   Cc: ferret_users@ferret.wrc.noaa.gov 
>   Sent: Wednesday, July 23, 2003 7:32 PM
>   Subject: Re: mp_stereographic
> 
> 
>   Gabriel,
> 
>   If you could condense your script to as few lines as possible and
have it use the default 'coads_climatology' dataset it'll be much easier
for us to see what's going wrong.
> 
> 
>   -- Jon
> 
> 
>   Gabriel Clauzet wrote:
> 
>     Dear Ferreters: 
>     I am trying to do contour plot on polar stereographic projection
using go mp_stereographic_south.jnl. I can put the Antarctica in the
bottom of the figure and the equator in the upper part of the figure.
Ferret always plot the figure in a North Hemispheric view.
> 
>     Does any one have an idea to help me ?
> 
>     Thanks
> 
>     Gabriel
> 

------------------------------------------
John Donners
Royal Netherlands Meteorological Institute
Dept. Oceanographic Research

P.O. Box 201
Wilheminalaan 10
3730 AE  De Bilt
tel. nr. : ++31 30 2206701
fax. nr. : ++31 30 2202570
e-mail   : donners@knmi.nl
--------------------------------------------------------------
Zie ook/see also: http://www.knmi.nl/maildisclaimer.html