Re: Re: [ncl-talk] trYMinF,trYMaxF can not work

From: 朱学明 <557xiaozhu_at_nyahnyahspammersnyahnyah>
Date: Mon Dec 16 2013 - 09:15:59 MST

Hi Mary,
   Thank you for your replying.
   I have figured out the issue. It was my fault. So just discard this problem please.
   Regards.

2013-12-17



Xueming Zhu 朱学明
Key Laboratory of Research on Marine Hazards Forecasting (LoMF), SOA
National Marine Environmental Forecasting Center (NMEFC)
No.8, Dahuisi Road, Haidian District, Beijing, 100081
People's Republic of China
Tel:+86-10-82481923



发件人:Mary Haley <haley@ucar.edu>
发送时间:2013-12-17 00:07
主题:Re: [ncl-talk] trYMinF,trYMaxF can not work
收件人:"朱学明"<557xiaozhu@163.com>
抄送:"ncl-talk@ucar.edu USERS"<ncl-talk@ucar.edu>

Hello,

I couldn't reproduce the problem using the data and script you provided. For starters, you only gave me one dimension of the "dataL" array, and it looks like it is a 2D array in reality.
Also, I don't have the "spec@phase" values that are used to generate "phase".

  
--Mary

On Dec 14, 2013, at 9:22 PM, 朱学明 <557xiaozhu@163.com> wrote:

> Hello,
> I had got some trouble about the function gsn_csm_xy2. The figure (please see attachment) got from the code as follow looks like wired very much. I do not know the reason for that. It looks like that theresL@trYMinF and resL@trYMaxF do not work, or any other reason?
> I had also attached the output files, freq1.dat, freq.dat, coher.dat. The variables's value look like reasonable. But why the figure...? I am confused.
>
> Anyone can help me to take a look at the code? Thank you so much.
>
>
> ;************************************************
> ; set function arguments
> ;************************************************
> iopt = 0 ; detrending iopt: 0=>remove mean 1=>remove mean and detrend
> jave = 5 ; smoothing periodogram: (0 <= sm <= ??.) should be at least 3 and odd
> pct = 0.10 ; percent tapered: (0.0 <= pct <= 1.0) 0.10 common.
> ;************************************************
> ; calculate spectrum
> ;************************************************
> spec = specxy_anal(xn,yn,iopt,jave,pct)
> ;; spec@frq=1/spec@frq
> ;************************************************
> ; plotting parameters that remain constant
> ;************************************************
> wks = gsn_open_wks("png",picname) ; Opens a ps file
> ; res = True ; no plot mods desired
> ; res@gsnFrame = False ; required for panel plots
> ; res@gsnDraw = False ; required for panel plots
> ; res@tiXAxisString = "Period (day)" ; xaxis ; Frequency (cycles/day)
>
> ; res@trYLog = True
> ; res@trXLog = True
> ;***********************************************
> ; create plot of coherance
> ;************************************************
> dataL=new((/2,dimsizes(spec@coher)/),typeof(spec@coher))
> dataL(0,:)=spec@coher
> dataL(1,:)=spec@coher_probability(1)
>
> resL = True
> resL@gsnMaximize = True
> ; resL@vpWidthF = 1.0
> ; resL@trXMinF = 0.001
> ; resL@trXMaxF = 0.01
> ; resL@trYMinF = min(dataL(0,:))
> ; resL@trYMaxF = max(dataL(0,:))
> resL@tiYAxisString = "Coherance" ; yaxis
> resL@xyLineThicknesses = (/1.0,2.0/)
> resL@xyLineColors = (/"black","red"/)
> print(min(dataL(0,:)))
> print(max(dataL(0,:)))
> ;***********************************************
> ; create plot of phase
> ;************************************************
> asciiwrite("probability.dat",spec@coher_probability)
> phase=spec@phase
> phase@_FillValue= 1e20
> phase=where(spec@coher.le.spec@coher_probability(1),phase@_FillValue,phase)
>
> resR = True
> ; resR@trXMinF = 0.001
> ; resR@trXMaxF = 0.01
> resR@trYMinF = 0. ;
> resR@trYMaxF = 60.
> resR@tiYAxisString = "Phase" ; yaxis
> resR@xyMarkLineModes = "Markers" ; choose which have markers
> resR@xyMarkers = 16 ; choose type of marker
> resR@xyMarkerColor = "red" ; Marker color
> resR@xyMarkerSizeF = 0.005 ; Marker size (default 0.01)
> ; asciiwrite("coher.dat",dataL(0,:))
> ; asciiwrite("freq.dat",spec@frq)
>
> plot=gsn_csm_xy2(wks,spec@frq,dataL,phase,resL,resR); create plot
> delete(wks)
> ; resP = True ; panel mods desired
> ; resP@gsnMaximize = True ; blow up plot
> ; resP@gsnPanelRowSpec = True ; tell panel whar order to plot
> ; gsn_panel(wks,plot,(/1,2/),resP) ; create panel plot
> ; asciiwrite("freq.dat",spec@frq)
> system("convert -alpha off -background white -density 300 -trim "+picname+".eps"+" "+picname+wks_type)
>
> 2013-12-15
> Xueming Zhu 朱学明
> Key Laboratory of Research on Marine Hazards Forecasting (LoMF), SOA
> National Marine Environmental Forecasting Center (NMEFC)
> No.8, Dahuisi Road, Haidian District, Beijing, 100081
> People's Republic of China
> Tel:+86-10-82481923
> <cross_spec_TK-PCM_seasonal.png><probability.dat><coher.dat><freq.dat>_______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Dec 16 09:31:22 2013

This archive was generated by hypermail 2.1.8 : Wed Dec 18 2013 - 14:29:51 MST