overlay plot doesn't work

From: Katja Lohmüller <lohmueller_at_nyahnyahspammersnyahnyah>
Date: Fri Sep 27 2013 - 05:38:21 MDT

  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
Received on Fri Sep 27 05:38:29 2013

This archive was generated by hypermail 2.1.8 : Tue Oct 01 2013 - 14:41:43 MDT