[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