Overlay plot does not work

Dear NCL-talk,

I testest the suggestion below but if I do this way, I get the error
message that I have to define "line". I tried different versions but
they all did not work. Which dimension has "line" to be?

Thank you for your help!
Katja

Quoting Adam Phillips <asphilli@ucar.edu>:

[Hide Quoted Text]
Hi Katja,
Try putting gsn_add_polyline in your j do loop and create a new
graphical array to hold the gsn_add_polyline output:

line = new(ntime-2,graphic)
do j= 0,ntime-3
.....
   line(j) = gsn_add_polyline(wks,contour,(/x(j),x(j+1)/),(/y(j),y(j+1)/),pres)
end do
draw(contour)
frame(wks)
delete(wks)

If that does not help please respond to the ncl-talk email list and
let everyone know..
Adam

On 09/27/2013 05:38 AM, Katja Lohmüller wrote:
   Dear NCL-talk,

I'm trying to illustrate a figure with a contour plot in the
background and in the foreground, there should be a line drawn in
different colors. Now my problem is the function overlay(x,y) doesn't
work in my case because I'm using the procedure gsn_polyline to create
the line and an overlay needs a function as argument. If I'm using
line= gsn_add_polyline, the line is black and that's not what I want
it to be.

Any suggestions how to solve my problem?

Thanks for your help,
Katja.

The code is as follows:

;************************************************************
; Target Classification 2D
;************************************************************

;************************************************************
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"
;************************************************************

begin

;************************************************************
;creation of filenames for loop
;************************************************************

;do j=10,31

      day ="17"
      month ="02"
      year ="2013"
      rest ="_lindenberg_classification"
      directory
="/localdata/Lohmueller/BA/Daten_Lindenberg/Classification/"
      type ="nc"

      filename = directory+year+month+day+rest

;************************************************************
;read in data for contour plot
;************************************************************

      a = addfile (filename+"."+type , "r")

;a=addfile("/localdata/Lohmueller/BA/Daten_Lindenberg/Classification/20130101_lindenberg_classification.nc",
"r")

time=a->time
tar=a->target_classification
height=a->height

TAR = tar(height|:,time|:) ;switch atts

TAR@units=" "
TAR@long_name= " "
time@units=" "
time@long_name= " "
height@units=" "
height@long_name=" "
height@standard_name=" "

;************************************************************
;read in data for polyline
;************************************************************

    path = "/localdata/Lohmueller/BA/Amdardaten/Flutra01JanDLH_3_049.txt"
    data = asciiread(path,(/40,3/),"float")
    ntime=40

wks= gsn_open_wks("pdf","Notflug")

;************************************************************
;the data slice
;************************************************************

slice=TAR(:,:)

;**********************************************
;polyline parameters
;**********************************************

x= new(ntime,float)
y= new(ntime,float)
vi= new(ntime,float)

;**********************************************
;put values to variable for polyline
;**********************************************

do i=0,ntime-1
    x(i)=data(i,0) ;extract distance
    y(i)=data(i,1) ;extract height
    vi(i)=data(i,2) ;extract Icing
end do

;print(x)
;print(y)

    gsn_define_colormap(wks,"BkBlAqGrYeOrReViWh200")

;***********************************************
;resources contour plot
;***********************************************

    res = True
    res@gsnDraw = False
    res@gsnFrame = False

    res@gsnMaximize = True
    res@gsnPaperOrientation = "landscape"
    res@cnFillOn = True
    ;res@cnCellFillEdgeColor = "black"
    res@cnFillMode = "CellFill"
    res@cnLinesOn = False
    res@cnLineLabelsOn = False
    res@cnLevelSelectionMode = "ManualLevels"
                                 ; set manual contour levels

    res@lbLabelAlignment = "BoxCenters"
    res@cnMinLevelValF = 0.1
                                          ; set the minimum contour level
    res@cnMaxLevelValF = 9
                                        ; set the maximum contour level
    res@cnConstFLabelConstantSpacingF = 1
    res@cnFillColors =
(/201,20,38,54,71,91,102,119,134,156,170/) ; set the
colors to be used
    res@cnLevelSpacingF = 1

    res@tiYAxisString = "Hoehe (m)" ; y axis title
    res@tiXAxisString = "Strecke (m)" ;
x axis title
    res@tiYAxisFontHeightF = 0.01
    res@tiXAxisFontHeightF = 0.01

    res = True
    ;res@lbAutoManage = True
    res@lbLabelAutoStride = False
    res@lbPerimOn = False
    res@lbOrientation = "Vertical"
    res@lbLabelStride = 1
    res@lbLabelStrings = (/"Klarer
Himmel","Wolkentropfen","Niesel & Regen","Niesel/Regen &
Wolkentropfen", "Eis", "Eis & unterkuehlte Tropfen", "Schmelzendes
Eis", "Schmelzendes Eis & Wolkentropfen","Aerosole",
"Insekten","Aerosole & Insekten"/)
    res@lbLabelFontHeightF = 0.01

    res@tiMainString = "Generierter Notflug EDDV-EDDF
17.02.2013"
    res@trGridType = "TriangularMesh"

;**** y-Achsen Ausschnitt ****
    res@trYMaxF = 6500
    res@trYMinF = 0
;*****************************

;**************************************************
;resources Polyline
;**************************************************

pres = True
pres@gsLineThicknessF = 5.0
pres@gsLineColor = "black" ; color of lines

v=10
s=new(2879,float)
time2=new(2879,float)
time2=time*60*60

s=v*time2

   res@sfXArray = s
   res@sfYArray = height

contour = gsn_csm_contour(wks,slice,res)

;***************************************************
; Set up resources to attach legend to map.
;***************************************************

do j= 0,ntime-3
      if (vi(j).eq.4) then
      pres@gsLineColor= (/4/)
      else
      if(vi(j).eq.2)then
      pres@gsLineColor= (/2/)
      else
      if(vi(j).eq.5)then
      pres@gsLineColor= (/5/)
      else
      if(vi(j).eq.3)then
      pres@gsLineColor= (/3/)
      end if
      end if
      end if
      end if
      gsn_polyline(wks,(/x(j),x(j+1)/),(/y(j),y(j+1)/),pres)
     end do

;line= gsn_add_polyline(wks,contour,x,y,pres)

draw(contour)
frame(wks)

delete(wks)

;end do

end

-- 
Katja Lohmüller
Institut fuer Meteorologie und Klimatologie
Leibniz Universitaet Hannover
Herrenhaeuser Str. 2
30419 Hannover
Email: lohmueller@muk.uni-hannover.de
Mobil: 0049 176 61396963
Buero: 0049 511 762 2396
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk