Hi Tony, The /graticule qualifier is probably the best answer for you, unless you really want to use a map projection. When you use /GRATICULE on the plot commands, it puts the lines at all of the large tics on the axes. Use the PPL XAXIS and YAXIS commands to set the tic interval. fill/graticule/set my_var ppl xaxis 200,360,10 ppl yaxis -30,30,10 ppl fill ! contour/levels , with no arguments to /levels, uses the same ! levels as used in the previous FILL command. contour/nolab/overlay/levels my_var go fland 5 The mp_graticule script works only if you have called the plot commands as curvilinear plots, with three arguments. The script mp_graticule.jnl is expecting the plot to be in the page coordinates defined by the projection defined by mp_plate_caree.jnl, not latitude-longitude coordinates. ... go mp_plate_caree fill/nolab/levels=(($ppl$zmin))($8,$9,$11)(($ppl$zmin)) $12, x_page,y_page ! plot memoire des contours contour/nolab/overlay/levels=($8,$9,$13) $12, x_page, y_page ! masque de Terre go mp_fland 5 ! grille lon-lat go mp_graticule $2 $3 10 $4 $5 10 In this plotting mode, you lose the latitude-longitude labels on the axes, and you would need to use label commands to put them on individually (as in the script mp_stereo_demo.jnl, which is explained in the FAQ on creating publication- quality map projection plots: http://www.ferret.noaa.gov/Ferret/FAQ/custom_plots/polar_stereographic.html Tony Jolibois wrote:
|