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