Hi,
sometimes you miss some stupid mistakes if you work too long on the
same thing.. Thank you very much, your hind to change the colors
solved my problem.
Have a great day!
Katja
Zitat von Adam Phillips <asphilli@ucar.edu>:
> Hi Katja,
> You are specifying the BkBlAqGrYeOrReWh200 colormap:
> http://www.ncl.ucar.edu/Document/Graphics/ColorTables/BkBlAqGrYeOrReViWh200.shtml
> and then within your if else statements you asking for color indices
> 2,3,4 or 5. Looking at the colormap those 4 indices are black or
> close to black. Try changing those 4 numbers to say 50, 100, 150, or
> 200 and see if the line colors switch from black. If they do not,
> then I would put in print statements in your if/else statements to
> see if they are getting executed...
> Hope that helps.
> Adam
>
>> On Sep 29, 2013, at 11:32, Katja Lohmüller
>> <lohmueller@muk.uni-hannover.de> wrote:
>>
>> I did it the same way like it is done in the snippet you sent me.
>> My script is now 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=" "
>>
>> ;print(TAR(:,2))
>>
>> ;************************************************************
>> ;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 (normally:extract
>> latitude(Breite(N)))
>> y(i)=data(i,1) ;extract height (normally:extract
>> longitude(Länge(E)))
>> 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@lbOrientation = "Default"
>> 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.
>> ;***************************************************
>> line= new(ntime-2,graphic)
>>
>> 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
>> line(j) =
>> gsn_add_polyline(wks,contour,(/x(j),x(j+1)/),(/y(j),y(j+1)/),pres)
>> ;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
>>
>> Quoting Adam Phillips <asphilli@ucar.edu>:
>>
>>> Hi Katja,
>>> Did you place your numerous if else statements that define
>>> gsLineColor within the do loop? When I put "....." in the snippet
>>> of code I sent you yesterday I assumed you would leave your if
>>> else statements in there..
>>>
>>> If the if else statements were in your script and you are still
>>> getting black lines please respond to the ncl-talk email list with
>>> your complete script attached along with a copy of the resulting
>>> plot.
>>> Adam
>>>
>>>> On Sep 28, 2013, at 8:41, Katja Lohmüller
>>>> <lohmueller@muk.uni-hannover.de> wrote:
>>>>
>>>> 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
>>
>>
>>
>> --
>> 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
>> <Notflug.png>
>> _______________________________________________
>> 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-talkReceived on Mon Sep 30 09:30:07 2013
This archive was generated by hypermail 2.1.8 : Tue Oct 01 2013 - 14:41:43 MDT