HI Sahidul,
Using dim_max/dim_min (or dim_max_n/dim_min_n if you have version
5.1.1), this should not be that bad. For example:
-----------------------------
fils = systemfunc("ls wrfout*")
a = addfiles(fils+".nc","r")
temp = a[:]->T ; T = temperature or theta? If theta use wrf_tk to
convert if desired
dailymax = temp(0::24,:,:,:) ; copy metadata
dailymax = dailymax@_FillValue
dailymin = dailymax
cntr = 0
do gg = 0,dimsizes(temp&Time)-1
dailymax(cntr,:,:,:) = dim_max_n(temp(cntr:cntr+23,:,:,:),0)
dailymin(cntr,:,:,:) = dim_min_n(temp(cntr:cntr+23,:,:,:),0)
cntr = cntr+24
end do
-----------------------------
That will calculate the max/min temperature for every grid point at
every level. I am not sure whether for your model data if temperature is
output or whether theta is output. If theta is, and you want to convert
it to temperature use wrf_tk:
http://www.ncl.ucar.edu/Document/Functions/Built-in/wrf_tk.shtml
Good luck,
Adam
On 03/08/2010 09:25 AM, Sahidul wrote:
> Dear Sir/NCL Users,
>
> I would like to calculate maximum and Minimum temperature in a day(24 hours).
>
> I am having WRF output files saved in 1 hour interval. Like below ( 1
> day forecast files = 24)
>
> wrfout_d01_2010-03-05-00:00:00 (Times = 2010-03-05-00:00:00)
> (temperature -> T (Time, lev,lat,lon) )
> wrfout_d01_2010-03-05-01:00:00 (Times = 2010-03-05-01:00:00)
> wrfout_d01_2010-03-05-02:00:00 (Times = 2010-03-05-02:00:00)
> ......
> .....
> wrfout_d01_2010-03-05-23:00:00 (Times = 2010-03-05-23:00:00)
>
>
> please help me to calculate it. If anyone has any scripts for this,
> kindly provide me.
>
> Thanking you,
>
> Sahidul
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
-- __________________________________________________ Adam Phillips asphilli@ucar.edu National Center for Atmospheric Research tel: (303) 497-1726 Climate and Global Dynamics Division fax: (303) 497-1333 P.O. Box 3000 Boulder, CO 80307-3000 http://www.cgd.ucar.edu/cas/asphilli _______________________________________________ ncl-talk mailing list List instructions, subscriber options, unsubscribe: http://mailman.ucar.edu/mailman/listinfo/ncl-talkReceived on Mon Mar 8 10:01:57 2010
This archive was generated by hypermail 2.1.8 : Mon Mar 08 2010 - 12:07:42 MST