[ncarg-talk] help
Mary Haley
haley at ucar.edu
Wed Mar 18 14:53:00 MDT 2015
res at cnLinesOn = False
Please post your NCL questions to ncl-talk at ucar.edu, and not
ncarg-talk at ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
--Mary
On Tue, Mar 17, 2015 at 7:06 AM, Gracielle Siqueira <
graciqueira17 at hotmail.com> wrote:
> Hi,
>
> Could you help me ?
> How delete of lines of contour map, call here "isolinhas" in portuguese.
>
> Thanks
>
> ------------------------------
> Date: Wed, 24 Dec 2014 08:30:49 -0700
> From: haley at ucar.edu
> To: liuy0813 at gmail.com
> CC: ncarg-talk at ucar.edu
> Subject: Re: [ncarg-talk] help for basemap overlay with contour_vector map
>
> In the future, please post your questions to ncl-talk and not
> ncarg-talk. This email list is for people using NCAR Graphics.
>
> You need to subscribe to ncl-talk before you can post:
>
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
> Regarding your problem: You cannot overlay two maps using NCL. You
> should have gotten an error when you tried to do this:
>
> overlay(map,plot1)
>
> because both "map" and "plot1" are maps.
>
> You can try using gsn_csm_vector_scalar instead of gsn_csm_vector_scalar_map,
> and then you will also need to remove all the map resources you are setting
> with vres1.
>
> Also, gsn_csm_contour should have worked just fine, and show the lat/lon
> coordinates, because you have lat2d and lon2d attached as attributes.
>
> Please see the two "try1.ncl" and "try2.ncl" scripts. Since you didn't
> provide the data file, I can't run them to guarantee they work. try1.ncl
> uses gsn_csm_vector_scalar, and try2.ncl uses both gsn_csm_contour and
> gsn_csm_vector.
>
> --Mary
>
> On Mon, Dec 22, 2014 at 6:41 PM, 刘勇 <liuy0813 at gmail.com> wrote:
>
>
> hi ,
>
> I want to draw a contour and vector plot on a self-define coastline,my
> code shows below,First I try to draw a basemap ,then use
> gsn_csm_vector_scalar_map to draw contour and vector map, finally try to
> overlay this two map,but the coastline of my first map were overlayed by
> the second map.
> I also try this way: First draw the basemap, then draw contour map
> through gsn_csm_contour(wks,out,res),but this doesn't show the lon&lat
> infomation of out,and can't overlay to my basemap.
> hope U can help me.
>
> thxs.
>
> liuy
>
>
> fig show below:
> http://oneputao.qiniudn.com/zonal.png
> http://oneputao.qiniudn.com/current.png
>
> code are there:
> ;*************************************************
>
> ; plot_xmbay.ncl
> ;
> ;*************************************************
> ;load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> ;load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "./ROMS_utils.ncl"
> undef("create_map")
> function create_map(wks,title,resolution)
> local a, res2
> begin
> res2 = True
> res2 at mpDataBaseVersion = resolution
> res2 at gsnMaximize = True
> res2 at gsnDraw = False
> res2 at gsnFrame = False
>
> res2 at mpOutlineOn = False
> res2 at mpFillOn = False
>
> ;---Turn on fancier tickmark labels.
> res2 at pmTickMarkDisplayMode = "Always"
>
> lon1=117.80
> lat1=24.20
> lon2=118.35+15.0/60.0
> lat2=24.4+25.0/60.0
> ;---Zoom in on area of interest
> res2 at mpLimitMode = "Corners"
> res2 at mpLeftCornerLatF = lat1
> res2 at mpLeftCornerLonF = lon1
> res2 at mpRightCornerLatF = lat2
> res2 at mpRightCornerLonF = lon2
> ;res2 at tiMainString = title
> ; res2 at mpGridAndLimbDrawOrder = "Predraw"
>
> ;---Create map.
> map = gsn_csm_map(wks,res2)
>
> return(map)
> end
>
> begin
> ;************************************************
> ; User settings
> ;************************************************
>
> fhis="roms_his_2010_nouse.nc"
> outfile = "roms_5"
> depth = -2.0
> rec = 1
>
> ;***********************************************
> ; Read data and interpolate
> ;***********************************************
> his = addfile (fhis,"r")
> lon2d = his->lon_rho
> lat2d = his->lat_rho
> out = roms_3d_interp(his,"salt",rec,depth)
> ur = roms_3d_interp(his, "u", rec, depth) ; they are automatically
> put to "rho"
> vr = roms_3d_interp(his, "v", rec, depth)
> angle = his->angle
> uvrot = uv_rot(ur,vr,angle)
>
> urr = uvrot(0,:,:)
> vrr = uvrot(1,:,:)
>
>
> delete(uvrot)
> delete(ur)
> delete(vr)
>
> urr at lat2d = lat2d
> urr at lon2d = lon2d
> vrr at lat2d = lat2d
> vrr at lon2d = lon2d
> out at lat2d = lat2d
> out at lon2d = lon2d
>
> minValue = 0.0
> maxValue = 35.0
> step = 1
> stride = 5
>
> ;************************************************
> ; create plot
> ;************************************************
> wks_type = "x11" ; or "eps"
> wks_type at wkOrientation = "Portrait"
> wks = gsn_open_wks (wks_type, outfile) ; open workstation
> gsn_define_colormap(wks, "BlAqGrYeOrRevi200")
> i = NhlNewColor(wks,0.8,0.8,0.8) ; add gray to
> colormap
>
> reso="HighRes"
> title="Coast draw test"
> map = create_map(wks,title,reso)
> data=asciiread("./ranges/bigzone_zhu.txt",-1,"float")
> temp=reshape(data,(/dimsizes(data)/2,2/))
> lat0=temp(:,1)
> lon0=temp(:,0)
> lnres = True
> lnres at gsLineThicknessF = 3.0
>
> lnres at gsLineColor = "brown"
>
> lnres at gsFillColor = "gray"
> lnres at gsLineColor = "black"
> ln0 = gsn_add_polygon(wks,map,lon0,lat0,lnres)
> ln2 = gsn_add_polyline(wks,map,lon0,lat0,lnres)
>
> ;
> ; Control appearance of map.
> ;
>
> vres1 = True ; plot mods desired
> vres1 at gsnDraw = True
> vres1 at gsnFrame = True
> vres1 at gsnMaximize = True ; Maximize plot in frame
> ; vres1 at gsnPaperOrientation = "Portrait"
> ; vres1 at cnFillDrawOrder = "PreDraw"
> vres1 at cnFillOn = True ; turn on color for
> contours
> vres1 at cnLinesOn = False ; turn off contour lines
> vres1 at cnLineLabelsOn = False ; turn off contour line
> labels
> vres1 at cnFillMode = "RasterFill"
> vres1 at gsnScalarContour = True ; contours desired
> vres1 at gsnSpreadColors = True ; use full color map
> vres1 at gsnSpreadColorEnd = -2
> ; vres1 at mpLandFillColor = "gray" ; set land to be gray
> vres1 at lbLabelBarOn = True
> vres1 at lbLabelStride = stride
> vres1 at cnLevelSelectionMode = "ManualLevels" ; set manual contour
> levels
> vres1 at cnMinLevelValF = minValue ; set min contour
> level
> vres1 at cnMaxLevelValF = maxValue ; set max contour
> level
> vres1 at cnLevelSpacingF = step ; set contour spacing
> vres1 at lbOrientation = "Vertical" ; /Vertical label bar
> vres1 at pmLabelBarOrthogonalPosF = -0.01 ; move label bar closer
> vres1 at pmLabelBarDisplayMode = "Always" ; Turn on a label bar.
> vres1 at lbPerimOn = False ; no box around it
> vres1 at lbBoxLinesOn = True ; Yes/No labelbar box
> lines.
> vres1 at vcRefMagnitudeF = 1.0 ; define vector ref
> mag
> vres1 at vcRefLengthF = 0.045 ; define length
> of vec ref
> vres1 at vcRefAnnoOrthogonalPosF = -1.0 ; move ref vector
> vres1 at vcRefAnnoParallelPosF = 0.108
> vres1 at vcRefAnnoArrowLineColor = "black" ; change ref vector
> color
> vres1 at vcRefAnnoArrowUseVecColor = False ; don't use vec
> color for ref
> vres1 at vcRefAnnoString1 = "1.0 m/s"
> vres1 at vcLabelsOn = False
> vres1 at vcLineArrowColor = "black"
> vres1 at vcRefAnnoOn = True
> vres1 at vcMonoLineArrowColor = True ; vec's colored by
> their mag
> vres1 at vcLineArrowHeadMaxSizeF = 0.008
> vres1 at vcLineArrowHeadMinSizeF = 0.0055
> vres1 at vcMinDistanceF = 0.0145 ; thin vectors ,
> xishu with big data
> vres1 at vcLineArrowThicknessF = 1.5 ; change vector
> thickness
> vres1 at vcGlyphStyle = "CurlyVector" ; turn on curly vectors
> lon_in1=117.80
> lat_in1=24.20
> lon_in2=118.35+15.0/60.0
> lat_in2=24.4+25.0/60.0
>
> ; MAP
> ; vres1 at mpProjection = "Mercator"
> vres1 at mpLimitMode = "Corners" ; choose range of
> map
> vres1 at mpLeftCornerLatF = lat_in1
> vres1 at mpLeftCornerLonF = lon_in1
> vres1 at mpRightCornerLatF = lat_in2
> vres1 at mpRightCornerLonF = lon_in2
> vres1 at mpDataBaseVersion = "HighRes" ; use high
> resolution coast
> vres1 at pmTickMarkDisplayMode = "Always" ; turn on tickmarks
>
> plot1= gsn_csm_vector_scalar_map(wks,urr,vrr,out,vres1)
> overlay(map,plot1)
> draw(map)
> frame(wks)
>
> end
> ------------------------------
>
> _______________________________________________
> ncarg-talk mailing list
> ncarg-talk at ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncarg-talk
>
>
>
> _______________________________________________ ncarg-talk mailing list
> ncarg-talk at ucar.edu http://mailman.ucar.edu/mailman/listinfo/ncarg-talk
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncarg-talk/attachments/20150318/af706a79/attachment-0001.html
More information about the ncarg-talk
mailing list