Dear Sir or madam:
I want to draw the cold front symbol in a map .I use the below script and it turns out the cold symbol with the orientations
northeast
but I want the blue triangle directe towards the reverse side . how can I write the command ,hope for your help!
make this symbol from to
thank you !
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
begin
cmap = (/ \
(/ 1., 1., 1. /), \ ; 0 - White background.
(/ 0., 0., 0. /), \ ; 1 - Black foreground.
(/ 1., 0., 0. /), \ ; 2 - Red.
(/ 0., 0., 1. /), \ ; 3 - Blue.
(/ 0., 1., 1. /), \ ; 4 - Cyan.
(/ 1., 0., 1. /) \ ; 5 - Magenta.
/)
wks = gsn_open_wks("ps","weather_sym")
gsn_define_colormap(wks,cmap)
res = True
res@gsnFrame = False
;
; Define size of map in frame.
;
res@vpXF = 0.06
res@vpYF = 0.90
res@vpWidthF = 0.88
res@vpHeightF = 0.88
res@mpFillOn = False
res@mpOutlineBoundarySets = "National"
res@mpNationalLineColor = 1
res@mpGridAndLimbOn = True
res@mpCenterLatF = 40.
res@mpCenterLonF = 110.
res@mpLimitMode = "LatLon" ; limit map via lat/lon
res@mpMinLatF = 24.0 ; map area
res@mpMaxLatF = 50.0 ; latitudes
res@mpMinLonF = 90 ; and
res@mpMaxLonF = 140 ; longitudes
map = gsn_map(wks,"LambertEqualArea",res)
wmsetp("ezf",1) ; Tell wmap we are using an existing map projection.
wmsetp("lin",1.0) ; Line width of front curve.
wmsetp("fro","cold"); Specify stationary front.
wmsetp("cfc",3) ; Use blue for the triangles.
; wmsetp("wfc",2) ; Use red for the bumps.
wmsetp("swi",0.02) ; Increase the size of the bumps and triangles.
xlat = (/ 40.65, 39.2, 38.18, 37.67, 37.85 , 38,37/) ; Latitudes.
xlon = (/127.1, 124.06, 121, 117.3, 113.5,110,109/) ; Longitudes.
wmdrft(wks, xlat, xlon)
frame(wks)
cmd="convert -geometry 800x800 -density 300 -trim "+fon+".eps "+fon+".png"
print(cmd)
system(cmd)
end
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2.1.8 : Mon Oct 08 2012 - 15:54:16 MDT