hi Basit,
you are almost there. After you convert those strings to integers (year,
month, day etc individually for every time of output), use those in
cd_inv_calendar() to get the coordinate value of 'time' axis. Now, to do
this, you'll require to chose a 'units'. Examples are there inside that
routine. Use that same units as 'units' attribute of the variable (1-d
array) you've just calculated at every time of output. This variable can
now be the input to time_axis_labels().
regards,
-- *Dr Arindam Chakraborty* Assistant Professor Centre for Atmospheric and Oceanic Sciences Indian Institute of Science Bangalore - 560 012, INDIA. Tel: +91-80-22933074 http://caos.iisc.ernet.in/ <http://caos.iisc.ernet.in/faculty/arch> On Mon, May 28, 2012 at 4:04 PM, BasitAli Khan <BasitAli.Khan@kaust.edu.sa>wrote: > Dear Chakraborty, > > Thanks for the prompt reply. > Well, .. When you say numeric value (not a string) it means the entire > date and time value should be an integer …. right? The date and time > value is one string (can be obtained as one string) for every history time > in wrfout. So what I do is pick the numeric part of date and time from that > string and convert it into integer. But then I would have 6 different > variables I.e year, month, day, hour, minutes and seconds. I can't > obviously add them to get them in time_label_axis, I can only put them in > an array however, time_axis_label won't accept any 2 dimensional array for > its 0th argument. > > Sorry, I am still far from reaching any solution. If you have some ncl > script written for wrfoutput, that probably could help me in getting this > done. > > Cheers, > ---- > > Basit A. Khan, Ph.D.**** > > Postdoctoral Research Fellow**** > > Division of Physical Sciences & Engineering**** > > Office# 3204, Level 3, Building 1,**** > > King Abdullah University of Science & Technology**** > > 4700 King Abdullah Blvd, Box 2753, Thuwal 23955 –6900,**** > > Kingdom of Saudi Arabia.**** > > ** ** > > Office: +966(0)2 808 0276, Mobile: +966(0)5 0860 3617**** > > E-mail: basitali.khan@kaust.edu.sa**** > > Skype name: basit.a.khan > > From: Arindam Chakraborty <arch@caos.iisc.ernet.in> > Reply-To: "arch13@gmail.com" <arch13@gmail.com> > Date: Mon, 28 May 2012 12:35:31 +0300 > To: KAUST <basitali.khan@kaust.edu.sa> > Cc: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu> > Subject: Re: how to use "time_axis_labels"? > > hi Basit, > make sure that the 0th argument to time_axis_label() is an 1-d array > containing the 'time' values as numbers. The numbers must have been > calculated from a origin with a proper units, specified by the "units" > attribute to the 0th argument. In your case, please make sure that > 'newdate' is a numerical variable and not a string, and has a units > attribute something of this form: "days since yyyy-mm-dd hh:mm:ss". Check > cd_inv_calendar() for more details about this "units" attribute. > > regards, > -- > *Dr Arindam Chakraborty > * > > > On Mon, May 28, 2012 at 2:32 PM, BasitAli Khan <BasitAli.Khan@kaust.edu.sa > > wrote: > >> Hi NCL users, >> I have been trying to use time_axis_labels to plot an xy-plot from my >> wrfoutput however, I am getting the following error. >> >> fatal:Argument type mismatch on argument (0) of (time_axis_labels) can >> not coerce >> >> I looked at NCL mail archive and tried Saji Hameed's solution but that >> also didn't work for me. I would greatly appreciate if some could help to >> fix this. My script for the xy plot is copied below. >> >> >> 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/wrf/WRFUserARW.ncl" >> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" >> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl" >> load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/time_axis_labels.ncl" >> >> begin >> >> f1 = addfile("/project/k14/basit/wrf-sims/wrf_153/times.nc","r") >> f2 = addfile("/project/k14/basit/wrf-sims/wrf_153/tauaer2.nc","r") >> wks = gsn_open_wks ("x11","outFname") >> >> nTimes = dimsizes(times_f1) >> >> date_f1 = times_f1 >> year_f1 = str_get_cols(date_f1,0,3) >> month_f1 = str_get_cols(date_f1,5,6) >> day_f1 = str_get_cols(date_f1,8,9) >> hour_f1 = str_get_cols(date_f1,11,12) >> minutes_f1 = str_get_cols(date_f1,14,15) >> >> aoT2 = wrf_user_getvar(f2,"TAUAER2", -1) >> aoT2_new = aoT2(:,0,300,300) >> >> newdate=year_f1+month_f1+day_f1+hour_f1+minutes_f1 >> resplot =True >> restick = True >> restick@ttmFormat="%Y%N%D%H%M" >> newdate@units="200606010000" >> restick@ttmAxis ="XB" >> print(newdate) >> >> time_axis_labels(newdate,resplot,restick) >> >> plot = gsn_csm_xy(wks, aoT2_new,newdate,resplot) >> >> end >> >> >> >> Thanks and best regards, >> ---- >> >> Basit A. Khan, Ph.D.**** >> >> Postdoctoral Research Fellow**** >> >> Division of Physical Sciences & Engineering**** >> >> Office# 3204, Level 3, Building 1,**** >> >> King Abdullah University of Science & Technology**** >> >> 4700 King Abdullah Blvd, Box 2753, Thuwal 23955 –6900,**** >> >> Kingdom of Saudi Arabia.**** >> >> ** ** >> >> Office: +966(0)2 808 0276, Mobile: +966(0)5 0860 3617**** >> >> E-mail: basitali.khan@kaust.edu.sa**** >> >> Skype name: basit.a.khan >> >> -- >> This message has been scanned for viruses and >> dangerous content by *MailScanner* <http://www.mailscanner.info/>, and >> is >> believed to be clean. >> _______________________________________________ >> 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
Received on Mon May 28 05:43:16 2012
This archive was generated by hypermail 2.1.8 : Wed May 30 2012 - 09:25:29 MDT