[Thread Prev][Thread Next][Index]

Re: mp_stereographic



Hi John,
Thanks for the answer.  I'd just like to add that if you want to change
the sign you  can run the mp_stereographic.jnl script as is, and then use
the opposite sign when you draw the plot

shade  var, -1* x_page, -1* y_page

Ansley Manke


John Donners wrote:

> 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

--
Ansley Manke  Pacific Marine Environmental Laboratory  Seattle WA




[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement