Re: Overlay plot does not work

From: Katja Lohmüller <lohmueller_at_nyahnyahspammersnyahnyah>
Date: Sat Sep 28 2013 - 08:41:44 MDT

Oh, sorry, I muss have missed tha line while reading. But
unfortunately I have to say that it doesn't work either, if I do it
this way. The line is still black.

Thanks for your help!
Katja

Quoting Adam Phillips <asphilli@ucar.edu>:

> Hi Katja,
> Notice that in my original emaiI I put:
> line = new(ntime-2,graphic)
> above your do loop. That line is necessary to predefine the line array.
> Adam
>
>
>> On Sep 28, 2013, at 2:54, Katja Lohmüller
>> <lohmueller@muk.uni-hannover.de> wrote:
>>
>> 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
>

-- 
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 Sat Sep 28 08:41:55 2013

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