Re: wmsetp

Insert a call:

 wmsetp("rev",1)

after:

   wmsetp("fro","cold"); Specify stationary front.

Fred Clare

On Sep 28, 2012, at 8:10 PM, nlcheng <cnl88@163.com> wrote:

> 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 <截图1.png> to <截图2.png>
>
> 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
>
>
>
>
> <weather_sym_4.ncl>_______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sat Sep 29 16:30:36 2012