[ncl-talk] Text label on map projection

David Adams dave.k.adams at gmail.com
Tue Nov 11 14:31:59 MST 2014


Thank Alan, Rick
yes, that´s exactly what I needed to do.

Dave

On Tue, Nov 11, 2014 at 3:13 PM, Alan Brammer <abrammer at albany.edu> wrote:

> Think you just need to set gsnDraw and gsnFrame = False, then draw(plot)
> and frame(wks) after you add the text to the plot.
>
> https://www.ncl.ucar.edu/Document/Graphics/Resources/gsn.shtml#gsnDraw
> https://www.ncl.ucar.edu/Document/Graphics/Resources/gsn.shtml#gsnFrame
>
> https://www.ncl.ucar.edu/Document/Functions/Built-in/draw.shtml
> https://www.ncl.ucar.edu/Document/Functions/Built-in/frame.shtml
>
>
> Good Luck,
> Alan.
>
>
>
> --
> Alan Brammer
> PhD Candidate — African Easterly Waves and TC genesis
> University At Albany, NY
> www.atmos.albany.edu/student/abrammer
>
>
>
> From: David Adams <dave.k.adams at gmail.com> <dave.k.adams at gmail.com>
> Reply: David Adams <dave.k.adams at gmail.com>> <dave.k.adams at gmail.com>
> Date: November 11, 2014 at 15:51:33
> To: ncl-talk at ucar.edu <ncl-talk at ucar.edu>> <ncl-talk at ucar.edu>
> Subject:  [ncl-talk] Text label on map projection
>
>    Hi NCLers,
> I need to add simple text to GOES image, but the text doesn´t appear.
> It might be due to the funky 2d lat,lon coordinates of GOES images.  Any
> suggestions?
>
> thanks,
> Dave
>
>  ---------Code--------------------------------------------
>    diri = "./"
>   file_name = asciiread("file_name",-1,"string")
>   fili ="goes13_4_2013_196_2315_subset.nc"
>    f    = addfile(diri+fili,"r")
>    new_temp    = f->new_temp          ; float data(time, yc, xc) ;
>                            ; data:type = "IR" ;
> ;*********************************
> ; Fix the variable so it has recognizable missing data
> ;*********************************
>    lat  = f->new_lat
>    lon  = f->new_lon
>    lat at _FillValue = max(lat)
>    lon at _FillValue = lat at _FillValue
>    new_temp at _FillValue   = 255.
>                            ; associate coordinates with variable
>    new_temp at lat2d = lat           ; (yc, xc)
>    new_temp at lon2d = lon
>
>
> ;*********************************
> ; create plot
> ;*********************************
>    sfx     = get_file_suffix(fili,0)           ; use file name
>
> ;  write output file name using first time record of file list
>   out_file = str_get_cols(file_name(0),0,21)
>   output_file = out_file
>    pltName = output_file
>
>  pltType = "ps"                             ; "ps", "eps", "pdf", "png"
>    pltDir  = "./"
>    wks = gsn_open_wks(pltType, pltDir+pltName)
>   ;gsn_define_colormap(wks,"amwg")             ; choose a color map
>    gsn_define_colormap(wks,"BlAqGrYeOrReVi200"); choose a color map
>    res                     = True
>    res at cnFillOn            = True            ; turn on color
>    res at cnFillMode          = "RasterFill"    ; cell mode
>    res at cnLinesOn           = False           ; Turn off contour lines
>    res at gsnSpreadColors     = True            ; use full colormap
>    res at gsnAddCyclic        = False           ; data not cyclic
>    res at gsnMaximize         = True            ; ps, pdf, pdf
>    res at pmTickMarkDisplayMode = "Always"      ; use NCL default
>   ;res at lbOrientation       = "Vertical"      ; vertical label bar
>    res at lbLabelAutoStride   = True            ; let NCL decide spacing
>
>    res at mpMinLatF           = 22.50  ;    min(d at lat2d)    ; region to zoom
> in on
>    res at mpMaxLatF           = 33.80  ;    max(d at lat2d)
>    res at mpMinLonF           = -115.00   ;   min(d at lon2d)
>    res at mpMaxLonF           = -102.500 ; max(d at lon2d)
>    res at mpFillOn            = False
>   ;res at mpOutlineBoundarySets = "USStates"    ; turn on state boundaries
>   ;res at mpOutlineBoundarySets = "AllBoundaries"
>    res at mpOutlineBoundarySets = "National"    ; turn on country boundaries
>    res at trGridType          = "TriangularMesh"   ; Necessary b/c lat, lon
>
>    res at cnLevelSelectionMode = "ManualLevels"  ;"ExplicitLevels"
>    res at cnMinLevelValF       = 200.                            ; set the
> minimum contour level
>    res at cnMaxLevelValF       = 300.                          ; set the
> maximum contour level
>    res at cnLevelSpacingF      = 5.0                            ; set the
> contour interval
>    res at cnRasterSmoothingOn  = True
>  res at lbLabelStride        = 5.0             ; every other label bar label
>
>    res at gsnLeftString       = new_temp at type
>    res at tiMainString        = fili
>    plot = gsn_csm_contour_map_ce(wks,new_temp(:,:), res)
>
>    txid = new(1,graphic)                 ; one text string
>
>    txres               = True
>    txres at txFont        = "helvetica-bold"
>    txres at txFontHeightF = 0.10
>    txid = gsn_add_text (wks,plot,"BASC",-108.2,28.2 ,txres)
>
>
>
>
>
>
>
>   _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20141111/85ab898d/attachment.html 


More information about the ncl-talk mailing list