Re: WRF: Max & Min Temperature calculation

From: Adam Phillips <asphilli_at_nyahnyahspammersnyahnyah>
Date: Mon Mar 08 2010 - 10:01:53 MST

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-talk
Received 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