[ncl-talk] Text label on map projection

[David Adams]

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)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20141111/ab35c6ad/attachment.html