[ncl-talk] gsn_add_text, text not showing up?

Elizabeth Burakowski burakows at ucar.edu
Thu Sep 18 09:48:45 MDT 2014


Thanks for the prompt response, Adam.  It was the draw order and it is not
working properly when I move gsn_add_text above draw(plot_01).

Liz.


On Thu, Sep 18, 2014 at 9:16 AM, Adam Phillips <asphilli at ucar.edu> wrote:

> Hi Elizabeth,
> I think this is just a draw order problem. According to the gsn_add_text
> documentation, "The text strings will be drawn only when the plot is
> drawn.". So, I think all you need to do is move these five lines above the
> draw(plot_d01) call:
> tres = True
> tres at txFontHeightF = 0.15
> dum1 = gsn_add_text(wks,plot_d01,"d01, 36 km",-65,28,tres)
> dum2 = gsn_add_text(wks,plot_d01,"d02,12 km",-67,36,tres)
> dum3 = gsn_add_text(wks,plot_d01,"d03, 4 km",-69,41,tres)
> draw(plot_d01)
> ;---- set polygon line resources:
> lnres = True
>
> If that does not solve it please let ncl-talk know..
> Adam
>
>
>
> On Thu, Sep 18, 2014 at 9:12 AM, Elizabeth Burakowski <burakows at ucar.edu>
> wrote:
>
>> Hi,
>>
>> I am trying to use gsn_add_text to include some text on my overlaid
>> contour maps.  However, the text I enter does not show up in the map and I
>> do not receive any errors.  Any ideas why the text is not showing up?
>>
>> Thanks,
>>
>> Liz.
>>
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>>
>> ;==========================================================
>> ;                         The main code
>> ;==========================================================
>>
>> begin
>>
>> ;---WRF Domains directory
>>     dir       =
>> "/glade/scratch/burakows/summer_MP8LW1SW2CLMcurrent_200811-200911/run/"
>>
>>     wrfo_d01    = addfile(dir+"wrfout_d01_2008-11-01_00:00:00.nc","r")
>>     wrfo_d02    = addfile(dir+"wrfout_d02_2008-11-01_00:00:00.nc","r")
>>     wrfo_d03  = addfile(dir+"wrfout_d03_2008-11-01_00:00:00.nc","r")
>>
>>         d01_LC = wrfo_d01->LU_INDEX ; read present-day land use and
>> preserve metadata
>> d02_LC  = wrfo_d02->LU_INDEX
>> d03_LC  = wrfo_d03->LU_INDEX
>>
>> d01_LC at description = " "
>> d02_LC at description = " "
>> d03_LC at description = " "
>>
>> d01_LC at lat2d  = wrfo_d01->XLAT(0,:,:)             ; Assign lat to x
>>         d01_LC at lon2d  = wrfo_d01->XLONG(0,:,:)             ; Assign lon
>> to x
>> d01_lat2d     = wrfo_d01->XLAT(0,:,:)
>> d01_lon2d     = wrfo_d01->XLONG(0,:,:)
>> dims_d01      = dimsizes(d01_LC(0,:,:))
>> nlat_d01      = dims_d01(0)
>> nlon_d01      = dims_d01(1)
>>
>>         d02_LC at lat2d  = wrfo_d02->XLAT(0,:,:)
>>         d02_LC at lon2d  = wrfo_d02->XLONG(0,:,:)
>> d02_lat2d     = wrfo_d02->XLAT(0,:,:)
>> d02_lon2d     = wrfo_d02->XLONG(0,:,:)
>> dims_d02      = dimsizes(d02_LC(0,:,:))
>> nlat_d02      = dims_d02(0)
>> nlon_d02      = dims_d02(1)
>>         d03_LC at lat2d  = wrfo_d03->XLAT(0,:,:)
>>         d03_LC at lon2d  = wrfo_d03->XLONG(0,:,:)
>> d03_lat2d     = wrfo_d03->XLAT(0,:,:)
>> d03_lon2d     = wrfo_d03->XLONG(0,:,:)
>> dims_d03      = dimsizes(d03_LC(0,:,:))
>> nlat_d03      = dims_d03(0)
>> nlon_d03      = dims_d03(1)
>>
>> ;--- Define land use (LU_Index) categories for USGS 24 cat
>>
>>    info = (/                             \ ; USGS 24 cat
>>         " 1 Urban & Built-Up"       ,\ ; red
>>         " 2 Dryland Crop & Pasture"     ,\ ; lightgoldenrod1
>>         " 3 Irr. Crop & Pasture"       ,\ ; goldenrod1
>>         " 4 Mix Dry/Irr Crop & Pasture" ,\ ; orange
>>         " 5 Crop/Grass Mosaic"          ,\ ; tan2
>>         " 6 Crop/Wood Mosaic"           ,\ ; tan4 (too dark)
>>         " 7 Grassland"           ,\ ; khaki1
>>         " 8 Shrubland"           ,\ ; khaki3
>>         " 9 Mix Shrub/Grass"            ,\ ; bisque
>>         " 10 Savanna"        ,\ ; palegoldenrod
>>         " 11 Decid. Broadleaf"          ,\ ; palegreen
>>         " 12 Decid. Needleleaf"         ,\ ; olivedrab3
>>         " 13 Evergreen Broadleaf"       ,\ ; chartruese
>>         " 14 Evergreen Needleleaf"      ,\ ; darkgreen
>>         " 15 Mixed Forest"              ,\ ; seagreen4
>>         " 16 Water Bodies"              ,\ ; slategray
>>         " 17 Herbaceous Wetland"        ,\ ; mediumturquoise
>> " 18 Wooden Wetland"            ,\ ; midnightblue
>> " 19 Barren/Sparse"       ,\ ; moccasin
>> " 20 Herbaceous Tundra"         ,\ ; magenta
>> " 21 Wooded Tundra"          ,\ ; mediumorchid4
>> " 22 Mixed Tundra"            ,\ ; maroon4
>> " 23 Bare Ground Tundra"        ,\ ; mediumpurple
>> " 24 Snow or Ice"               /) ; white
>>
>>
>>    ninfo = dimsizes(info) ; # of categories
>>
>> ;----------------------------------------------------------------------
>> ; Plotting options section
>> ;----------------------------------------------------------------------
>>        pltType = "ps"                 ; plot type
>>        pltDir = "./"               ; plot directory
>>        pltName = "DomainMap" ; plot name (ps file)
>>        pltPath = pltDir+pltName         ; plot path
>>
>>        colorscheme =
>> (/"red","lightgoldenrod1","goldenrod1","orange","tan2","tan4",\
>>
>> "khaki1","khaki3","bisque","palegoldenrod","palegreen","olivedrab3",\
>>
>> "chartreuse","darkgreen","seagreen4","slategray1","mediumturquoise","midnightblue",\
>>      "moccasin","magenta","mediumorchid4","maroon4",\
>>      "mediumpurple","white"/)
>>         ncolors = dimsizes(colorscheme)
>>  if(ninfo.ne.ncolors) then
>> print("size mismatch lu cats="+ninfo+"  ,ncolors="+ncolors)
>> exit
>> end if
>>
>>    wks = gsn_open_wks(pltType,pltPath)              ; create workstation
>> for ps file
>>
>> cnres                  = True                ; plot mods desired
>>   cnres at gsnDraw          = False
>>   cnres at gsnFrame         = False
>>
>>     cnres at lbLabelBarOn     = False             ; turn on in panel plot
>> cnres at cnInfoLabelOn    = False
>>   cnres at cnFillOn         = True               ; color Fill
>>   cnres at cnFillMode       = "RasterFill"
>>   cnres at cnLinesOn        =  False             ; Turn off contour lines
>> cnres at cnLineLabelsOn   = False
>>   cnres at cnLevelSelectionMode = "ExplicitLevels" ; set explict contour
>> levels
>>   cnres at cnLevels         = (/2,3,4,5,6,7,8,\
>> 9,10,11,12,13,14,\
>> 15,16,17,18,19,20,\
>> 21,22,23,24/)
>>   cnres at cnFillPalette    = colorscheme        ; distinct colors for
>> categories
>>   cnres at gsnSpreadColors  = False              ; use each color
>> sequentially
>>
>> res = cnres ; copy contour resources from first plot
>>     res at mpFillOn              = False
>>     res at mpOutlineOn           = True
>>     res at mpOutlineBoundarySets = "AllBoundaries"
>>     res at mpProjection          = "LambertConformal"
>> res at mpOutlineDrawOrder  = "PostDraw" ; draw map outlines last
>> res at mpDataBaseVersion     = "LowRes"  ; State boundaries
>> res at mpGridAndLimbOn  = False
>> res at pmTickMarkDisplayMode = "Always"         ; Turn on map tickmarks
>>
>>     res at mpLambertParallel1F  = (max(d01_lat2d)+min(d01_lat2d))*0.3
>> res at mpLambertParallel2F  = (max(d01_lat2d)+min(d01_lat2d))*0.7
>> res at mpLambertMeridianF  = -72 ;
>> res at mpLimitMode           = "Corners"        ; required for Lambert
>> Conformal
>> res at mpLeftCornerLatF      = d01_lat2d(0,0)
>>     res at mpLeftCornerLonF      = d01_lon2d(0,0)
>>     res at mpRightCornerLatF     = d01_lat2d(nlat_d01-1,nlon_d01-1)
>>     res at mpRightCornerLonF     = d01_lon2d(nlat_d01-1,nlon_d01-1)
>> res at mpPerimOn  = True
>> res at mpPerimLineThicknessF = 2.0
>> res at mpPerimLineColor  = "black"
>>  res at lbTitleString   = " "
>>  res at tfDoNDCOverlay   = True  ; turn off transformation
>>
>> res at gsnLeftString         = ""                ; Turn off left subtitle
>>         res at gsnRightString        = ""                ; Turn off right
>> subtitle
>>
>> res at gsnAddCyclic = False ; regional grid
>>
>> res2 = cnres
>>   res2 at lbLabelBarOn = False  ; Labelbar already created in 1st plot
>>   res2 at gsnMaximize  = False  ; Use maximization from original plot
>>  res3 = cnres
>>   res3 at lbLabelBarOn = False  ; Labelbar already created in 1st plot
>>   res3 at gsnMaximize  = False  ; Use maximization from original plot
>>
>> ;---- Plot three domains and their land cover
>>
>>         plot_d01 = gsn_csm_contour_map(wks,d01_LC(0,:,:),res) ; create
>> plot of d01 present-day land cover
>> plot_d02 = gsn_csm_contour(wks,d02_LC(0,:,:),res2) ; create plot of d02
>> present-day land cover
>> plot_d03 = gsn_csm_contour(wks,d03_LC(0,:,:),res3) ; create plot of d03
>> present-day land cover
>>
>> overlay(plot_d01,plot_d02)
>> overlay(plot_d01,plot_d03)
>>
>> draw(plot_d01)
>>  ;---- set polygon line resources:
>> lnres = True
>> lnres at gsLineColor = "black"
>> lnres at gsLineThicknessF = 2.0
>>  gsn_polyline(wks,plot_d01,d02_lon2d(:,0),d02_lat2d(:,0),lnres)
>>
>> gsn_polyline(wks,plot_d01,d02_lon2d(:,nlon_d02-1),d02_lat2d(:,nlon_d02-1),lnres)
>> gsn_polyline(wks,plot_d01,d02_lon2d(0,:),d02_lat2d(0,:),lnres)
>>
>> gsn_polyline(wks,plot_d01,d02_lon2d(nlat_d02-1,:),d02_lat2d(nlat_d02-1,:),lnres)
>>  gsn_polyline(wks,plot_d01,d03_lon2d(:,0),d03_lat2d(:,0),lnres)
>>
>> gsn_polyline(wks,plot_d01,d03_lon2d(:,nlon_d03-1),d03_lat2d(:,nlon_d03-1),lnres)
>> gsn_polyline(wks,plot_d01,d03_lon2d(0,:),d03_lat2d(0,:),lnres)
>>
>> gsn_polyline(wks,plot_d01,d03_lon2d(nlat_d03-1,:),d03_lat2d(nlat_d03-1,:),lnres)
>>  tres = True
>> tres at txFontHeightF = 0.15
>> dum1 = gsn_add_text(wks,plot_d01,"d01, 36 km",-65,28,tres)
>>         dum2 = gsn_add_text(wks,plot_d01,"d02,12 km",-67,36,tres)
>>         dum3 = gsn_add_text(wks,plot_d01,"d03, 4 km",-69,41,tres)
>>
>>
>> frame(wks)
>>
>> end
>>
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
>
> --
> Adam Phillips
> Associate Scientist,  Climate and Global Dynamics Division, NCAR
> www.cgd.ucar.edu/staff/asphilli/   303-497-1726
>
> <http://www.cgd.ucar.edu/staff/asphilli>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140918/f0cda376/attachment.html 


More information about the ncl-talk mailing list