Hi,
How can the axis tick marker be defined as integer for gsn_xy plot?
I got 980.0 982.0 ... on y-axis. I just want to have 980 , 982 on y-axis.
Thanks.
Yongzuo
________________________________________
From: ncl-talk-bounces@ucar.edu [ncl-talk-bounces@ucar.edu] on behalf of ncl-talk-request@ucar.edu [ncl-talk-request@ucar.edu]
Sent: Sunday, January 10, 2010 7:54 PM
To: ncl-talk@ucar.edu
Subject: ncl-talk Digest, Vol 74, Issue 25
Send ncl-talk mailing list submissions to
ncl-talk@ucar.edu
To subscribe or unsubscribe via the World Wide Web, visit
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
or, via email, send a message with subject or body 'help' to
ncl-talk-request@ucar.edu
You can reach the person managing the list at
ncl-talk-owner@ucar.edu
When replying, please edit your Subject line so it is more specific
than "Re: Contents of ncl-talk digest..."
Today's Topics:
1. Re: subtraction function (jwsmith@ucar.edu)
2. Re: Y-axis label and values (kalim ullah)
----------------------------------------------------------------------
Message: 1
Date: Sun, 10 Jan 2010 14:25:47 -0700
From: jwsmith@ucar.edu
Subject: Re: subtraction function
To: Dennis Shea <shea@ucar.edu>
Cc: ncl-talk <ncl-talk@ucar.edu>
Message-ID: <20100110142547.2kd3llso9w48k44c@web3.acd.ucar.edu>
Content-Type: text/plain; charset=UTF-8; DelSp="Yes"; format="flowed"
co(Time, bottom_top, south_north, west_east)
I want to take a ten - day average of co at 12z. I then want to
subtract this average from the individual 12 z time in that period.
This will be the anamolous field. I want to make a vertical cross
section of the anamolies.
When I did the subtraction in nco using ncbo and ncdiff. The
difference of co was taken, but it also seems to have taken the
difference of the XLAT and XLON. I think this happened because co is a
variable with lat/lon (bottom_top) coordinates, (I could be
misinterpreting the use of ncdiff and ncbo).
Hencem I was asking if the subtraction could be done in ncl. I hope I
have provided more clarity this time.
Thanks,
Jonathan
Quoting Dennis Shea <shea@ucar.edu>:
> I agree with Huang. It is not clear what you want to do.
> It does sound like conform is what ou need.
>
> If you have:
>
> x3(Times,lat,lon) and xBase(lat,lon)
>
> xDiff = x3 - conform(x3, xBase, (/1,2/) )
> copy_VarMeta(x3, xDiff)
>
> http://www.ncl.ucar.edu/Document/Functions/Built-in/conform.shtml
>
> Good Luck
>
> ph0007@ustc.edu wrote:
>> Hi, Jonathan
>>
>> I'm not really sure what you want to do. I guess the function "conform"
>> could help you.
>>
>> Please see the document about details.
>>
>> Huang Ping
>>
>>
>> 2010-01-10
>> ------------------------------------------------------------------------
>> ------------------------------------------------------------------------
>> *????* jwsmith
>> *?????* 2010-01-10 13:16:56
>> *????* ncl-talk
>> *???*
>> *???* [ncl-talk] subtraction function
>> Hello ncl
>> I see that there is a sum function in ncl. Is there a subtraction
>> fuction? I need to subtract a wrfout that contains the average of
>> fields from an single wrf model time. The ncbo and ncdiff functions
>> in nco did not work because the files contain lat/lon coordinates. Is
>> it possible to something like this in ncl?
>> Jonathan
>> Jonathan Wynn Smith
>> Doctoral Student in Howard University
>> Program in Atmospheric Sciences (HUPAS)
>> 2008 ASP Graduate Student Visitor
>> 2006 ASP Summer Colloquium Participant
>> 408 Thirkield Building
>> 2355 6th St. NW
>> Washington, DC 20059
>> E-mail: jw_smith@howard.edu or jwsmith9@gmail.com
>> Cell Phone: 336-601-4563
>> ----------------------------------------------------------------
>> This message was sent using IMP, the Internet Messaging Program.
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>> ------------------------------------------------------------------------
>>
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
Jonathan Wynn Smith
Doctoral Student in Howard University
Program in Atmospheric Sciences (HUPAS)
2008 ASP Graduate Student Visitor
2006 ASP Summer Colloquium Participant
408 Thirkield Building
2355 6th St. NW
Washington, DC 20059
E-mail: jw_smith@howard.edu or jwsmith9@gmail.com
Cell Phone: 336-601-4563
----------------------------------------------------------------
This message was sent using IMP, the Internet Messaging Program.
------------------------------
Message: 2
Date: Sun, 10 Jan 2010 17:54:09 -0800 (PST)
From: kalim ullah <kalimrid@yahoo.com>
Subject: Re: Y-axis label and values
To: Mary Haley <haley@ucar.edu>
Cc: ncl-talk@ucar.edu
Message-ID: <275943.50743.qm@web53708.mail.re2.yahoo.com>
Content-Type: text/plain; charset="utf-8"
Dear Mary, thanks alot for your kind help.
regrads,
kalim
--- On Fri, 1/8/10, Mary Haley <haley@ucar.edu> wrote:
From: Mary Haley <haley@ucar.edu>
Subject: Re: Y-axis label and values
To: "kalim ullah" <kalimrid@yahoo.com>
Cc: ncl-talk@ucar.edu
Date: Friday, January 8, 2010, 7:21 AM
Dear?
Kalim,
You can use ut_calendar to convert your time coordinate values to something more manageable, like years, months, and days, and then you can use these arrays to determine which part of the Y axis you want to label.
I changed your script to show how you might label the first day of every month. You can modify this to label the midpoint of every month, using the days_in_month calendar as you did before. You'll still want to use ut_calendar so you can retrieve the months and days.
You can use this same method to determine where you want minor tick marks.
--Mary
On Jan 6, 2010, at 8:19 PM, kalim ullah wrote:
Dear Mary,
Thanks a lot for spending your precious time on my problem.?I think you are right my time coordinate array may be wrong. When I ignore the (res@tmYLMode = "Explicit?) it works and the label appear on Y-axis like (1779600- 1780200) but I don?t know how to convert it into days and month label? ?So, for your conveyance my plot, NCEP data file and zAvgTime print report is attached given below. Your guidance in this regard will be appreciated.
thanks
kalim
*****************************
Variable: z
Type: float
Total Size: 10828800 bytes
2707200 values
Number of Dimensions: 3
Dimensions and sizes: [time | 150] x [lat | 94] x [lon | 192]
Coordinates:
time: [1779456..1780350]
lat: [88.542..-88.542]
lon: [ 0..358.125]
Number Of Attributes: 18
long_name : 6-hourly Mean of Latent Heat Net Flux at Surface
unpacked_valid_range : ( -2420, 2420 )
actual_range : ( -667.2, 1884.8 )
units : W/m^2
_FillValue : -32767
precision : 1
least_significant_digit : 0
GRIB_id : 121
GRIB_name : LHTFL
var_desc : Latent Heat Net Flux
dataset : NCEP/DOE AMIP-II Reanalysis (Reanalysis-2)
level_desc : Surface
statistic : Mean
parent_stat : Individual Obs
standard_name : surface_upward_latent_heat_flux
cell_methods : time: mean
valid_range : ( -2420, 2419.9 )
missing_value : -32767
?
Variable: zAvgTime
Type: float
Total Size: 115200 bytes
28800 values
Number of Dimensions: 2
Dimensions and sizes: [time | 150] x [lon | 192]
Coordinates:
time: [1779456..1780350]
lon: [ 0..358.125]
Number Of Attributes: 19
missing_value : -32767
valid_range : ( -2420, 2419.9 )
cell_methods : time: mean
standard_name : surface_upward_latent_heat_flux
parent_stat : Individual Obs
statistic : Mean
level_desc : Surface
dataset : NCEP/DOE AMIP-II Reanalysis (Reanalysis-2)
var_desc : Latent Heat Net Flux
GRIB_name : LHTFL
GRIB_id : 121
least_significant_digit : 0
precision : 1
_FillValue : -32767
units : W/m^2
actual_range : ( -667.2, 1884.8 )
unpacked_valid_range : ( -2420, 2420 )
long_name : 6-hourly Mean of Latent Heat Net Flux at Surface
average_op_ncl : dim_avg over dimension: lat
--- On Wed, 1/6/10, Mary Haley <haley@ucar.edu> wrote:
?
If these values are not in the range of 0.0 to 144.8, then this is what your problem is.
--Mary
On Jan 5, 2010, at 7:50 PM, kalim ullah wrote:
Dear Mary, thanks for your response. Actually, I want to plot as example (hov_7.ncl:) but after some mudification I got the?figure?without label(i.e.month names) on y-axis?,so my plot is attached herewith.
thanks
kalim
--- On Tue, 1/5/10, Mary Haley <haley@ucar.edu> wrote:
From: Mary Haley <haley@ucar.edu>
Subject: Re: Y-axis label and values
To: "kalim ullah" <kalimrid@yahoo.com>
Cc: ncl-talk@ucar.edu
Date: Tuesday, January 5, 2010, 8:59 AM
Kalim,
If your tickmark labels are not showing up in "explicit" mode, then usually the problem is that the tickmark values you're using are not within the range of your axis values.
Since I can't run your script, I'm not sure what your Y axis values are. ?You need to look at the arrays zAvgTime&time and major_ticks and make sure they are in the same data space.
major_ticks should be set to values that are within the range of your time values.
--Mary
On Jan 4, 2010, at 7:09 PM, kalim ullah wrote:
Dear NCL Users,
I?m trying to plot time vs lon figure but after plotting there is no Y-axis labels and values appear and also no error. Please could you help me in this regard? My script is attached given below:
Thanks in advance
kalim
******************************
?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
;=============================================
; data processing
;=============================================
?? f = addfile("lhtfl.sfc.gauss.2007.nc","r")
?? z = short2flt(f->lhtfl(0:149,:,:))???? ;AUG(1-5)
?? printVarSummary(z)
?? zAvgTime = dim_avg_Wrap(z(time|:,lon|:,{lat|0:40}))
?? printVarSummary(zAvgTime)
?
?? wks? = gsn_open_wks ("x11", "hov" )
?? gsn_define_colormap(wks,"BlAqGrYeOrRe")
?? ;*******************************
?? res??????????? = True????????????????????????
?? res@cnFillOn????????????? = True?
?? months????? = (/1,? 2,??? 3,?? 4,?? 5/)?? ???
?? years?????? = (/2007,2007,2007,2007,2007/)
?? tk_per_mon? = 5 ; Tickmarks per month.
?? nmonths???? = dimsizes(months)
?? ndays?????? = nmonths * tk_per_mon
?
?? day_in_mon? = new(nmonths,float)?????????? ; Number of days in each month
?? do i=0,nmonths-1
?? day_in_mon(i) = days_in_month(years(i),months(i))*1.
?? end do
?
?? major_ticks = new(ndays,float)??? ; Array to hold major tickmark values.
?? imonth????? = 0
?? days_accum? = 0.
?? ii4???????? = ispan(1,4,1)*1.
?? do i=0,ndays-1,5
?? major_ticks(i) = days_accum
?? major_ticks(i+1:i+4) = major_ticks(i) + ii4*(day_in_mon(imonth)/tk_per_mon)
?? days_accum???? = days_accum + day_in_mon(imonth)
?? imonth = imonth + 1
?? end do
?? ;*********************************
?? res@cnLinesOn??????? ?????= False
?? res@gsnSpreadColors?????? = True
?? res@cnLevelSelectionMode = "ManualLevels"
?? res@cnMinLevelValF?????? = -10.?????????????? ; min level
?? res@cnMaxLevelValF?????? =? 100.?????????????? ; max leve
?? res@lbOrientation???????? = "Vertical"
?? res@gsnLeftString???????? = "Latent Heat Net Flux(W/m^2)"
?? res@gsnRightString?????? = ""
?? res@tmYLMode?????????? = "Explicit"???????? ; Define own tick mark labels.
?? res@tmYLMinorValues??? = major_ticks(0:ndays-1:5)
?? res@tmYLValues???????? = major_ticks
?? res@tmYLLabels???????? = new(ndays,string)
?? res@tmYLLabels(:)???? = ""
?
; Label mid tickmark with the month label.
?? res@tmYLLabels(3::5)? = (/"JAN", "FEB", "MAR", "APR", "MAY"/)
?
?? res@tmYLMinorOutwardLengthF = 0.02??? ; Increase length and size of
?? res@tmYLMinorLengthF??????? = 0.02??? ; minor tickmarks
?? res@tmYLMinorThicknessF???? = 2.0
?? res@tmYLMajorOutwardLengthF = 0.01
?? res@tmYLMajorLengthF??????? = 0.01
?? plot = gsn_csm_hov(wks, zAvgTime(:,{30:120}) , res )
?? end
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
<hov.pdf>
<Lhov3.pdf><lhtfl.sfc.gauss.2007.nc>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100110/849ae607/attachment.html
------------------------------
_______________________________________________
ncl-talk mailing list
ncl-talk@ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
End of ncl-talk Digest, Vol 74, Issue 25
****************************************
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jan 11 07:37:54 2010
This archive was generated by hypermail 2.1.8 : Tue Jan 12 2010 - 15:38:20 MST