Psbasemap boundaries problems

Q&A
#1

I get the problem in plooting boundaries in psbasemap. It alway shows 4 boundaries (W,N,S,E) no matter the change the value. Please helps me to fix it.
This is my code to plot my maps

HYCOM_Feb_2014-2018.D26
HYCOM_Feb_2014-2018.D263234×2308 111 KB
:
gmtset FONT_ANNOT_PRIMARY 20p,Helvetica,black FONT_LABEL 20p
gmtset FONT_LABEL 20p MAP_TICK_LENGTH 5p MAP_TICK_PEN 2p MAP_FRAME_PEN 0.08p
gmtset PROJ_LENGTH_UNIT INCH PS_MEDIA a0

makecpt -CGMT_EASNFS.cpt -T0/120/10 -D >UOHC.cpt

psbasemap -R130/165/0/25 -JQ180/9 -BWS -Bxa -Bya -X1 -Y2 -K -P> %1.ps

psxy %1 -R -CUOHC.cpt -J -Ss0.05 -K -O >>%1.ps

pscoast -J -R -W1 -G200 -Di -O -K >>%1.ps

psxy bwp201902_C5_WUTIP.bt1km -R -J -O -W0.04i,0/0/0 -K >>%1.ps
psxy bwp201902_C5_WUTIP.bt1km -R -J -O -Ctc_scale1.cpt -Sc0.15 -W0.01i -K >>%1.ps
pstext bwp201902_C5_WUTIP.bt1km.date -R -J -F+f10p,Helvetica,black,-=0.5p,white -B5 -K -O >>%1.ps
gmtset FONT_LABEL 20p MAP_TICK_LENGTH 5p MAP_TICK_PEN 2p MAP_FRAME_PEN 2p

psscale -CUOHC.cpt -R -J -DJCR+o0.5/0i+w7.2i/0.15i+e -Bxa20f10 -By+l(m) -O -K >>%1.ps

rem ----------- plot legend ---------------
echo 149 25.2 0 > temp2
echo 152 25.2 63 >> temp2
echo 155 25.2 83 >> temp2
echo 158 25.2 96 >> temp2
echo 161 25.2 114 >> temp2
echo 164 25.2 136 >> temp2

psxy temp2 -R -J -O -Ctc_scale1.cpt -Sc0.2 -Y0.2 -K -N >> %1.ps
psxy temp2 -R -J -O -W0.03i,0/0/0 -Sc0.2 -K -N >> %1.ps

echo 149 26.5 18 0 0 2 TD/TS> temp2
echo 152 26.5 18 0 0 2 Cat.1>> temp2
echo 155 26.5 18 0 0 2 Cat.2>> temp2
echo 158 26.5 18 0 0 2 Cat.3>> temp2
echo 161 26.5 18 0 0 2 Cat.4>> temp2
echo 164 26.5 18 0 0 2 Cat.5>> temp2
pstext temp2 -R -J -Ctc_scale1.cpt -O -K -N >>%1.ps

echo 130 25.8 20,1,black LT %1 > name
pstext name -R -J -W -F+f+j -O -N >>%1.ps

psconvert %1.ps -Tg -A -P
del %1.ps name temp2

Best regards,

#2

Removing -B5 from the above command may solve your problem.

1 Like
#3

thank you very much. Appreciate that!

#4

Hello, in this map I also got the problem with the appearance of W, S. It should be shown 0N or 100W. But I didn’t. How I fix it?

#5

Can you please show me your new figure so that I can understand the problem?

#6

HYCOM_Wutip_2019021812.D20
HYCOM_Wutip_2019021812.D203258×2269 147 KB

Here also the new code:
gmtset FONT_ANNOT_PRIMARY 20p,Helvetica FONT_LABEL 20p
gmtset FONT_LABEL 20p MAP_TICK_LENGTH 5p MAP_TICK_PEN 2p MAP_FRAME_PEN 1p
gmtset PROJ_LENGTH_UNIT INCH PS_MEDIA a0

makecpt -CGMT_EASNFS.cpt -T0/240/20 -D >UOHC.cpt

psbasemap -R130/165/0/25 -JQ180/9 -Bxa5f5 -Bya5f5 -BWeSn -X1 -Y2 -K -P> %1.ps

psxy %1 -R -CUOHC.cpt -J -Ss0.05 -K -O >>%1.ps

pscoast -J -R -W1 -G200 -Di -O -K >>%1.ps

psxy bwp201902_C5_WUTIP.bt1km -R -J -O -W0.04i,0/0/0 -K >>%1.ps
psxy bwp201902_C5_WUTIP.bt1km -R -J -O -Ctc_scale1.cpt -Sc0.15 -W0.01i -K >>%1.ps
pstext bwp201902_C5_WUTIP.bt1km.date -R -J -F+f18p,Helvetica,white,-=~2p,black,+a30+j -K -O >>%1.ps

psscale -CUOHC.cpt -R -J -DJCR+o0.58/0i+w6i/0.15i+e -Bxa40f20 -By+lm -O -K >>%1.ps

rem ----------- plot legend ---------------
echo 149 25.2 0 > temp2
echo 152 25.2 63 >> temp2
echo 155 25.2 83 >> temp2
echo 158 25.2 96 >> temp2
echo 161 25.2 114 >> temp2
echo 164 25.2 136 >> temp2

psxy temp2 -R -J -O -Ctc_scale1.cpt -Sc0.2 -Y0.2 -K -N >> %1.ps
psxy temp2 -R -J -O -W0.03i,0/0/0 -Sc0.2 -K -N >> %1.ps

echo 149 26.5 18 0 0 2 TD/TS> temp2
echo 152 26.5 18 0 0 2 Cat.1>> temp2
echo 155 26.5 18 0 0 2 Cat.2>> temp2
echo 158 26.5 18 0 0 2 Cat.3>> temp2
echo 161 26.5 18 0 0 2 Cat.4>> temp2
echo 164 26.5 18 0 0 2 Cat.5>> temp2
pstext temp2 -R -J -Ctc_scale1.cpt -O -K -N >>%1.ps

echo 130 25.8 20,1,black LT %1 > name
pstext name -R -J -F+f+j -O -N >>%1.ps

psconvert %1.ps -Tg -A -P
del %1.ps name temp2

#7

Do you want to have suffix like N, S, W and E? If so, you should set the parameter FORMAT_GEO_MAP.

#8

I also fixed that. Thanks a lot!