[ncl-talk] Help Required in plotting average 2m temperature for 24 hours for 30 days in NCL
Muhammad Omer Mughal
m.mughal1 at postgrad.curtin.edu.au
Mon Oct 16 07:53:46 MDT 2017
Hi Mary
How can we overlay the average of wind speed vectors in the following code
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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "/data/muhdomer/NSCC/scratch/WRF_myfunc.ncl"
begin
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; input from WRF
; The WRF ARW input file.
wrf_file1 = systemfunc("ls /data/muhdomer/Build_WRF/Martilli/glade/u/home/martilli/wrf_3.8/wrfout_d05_2016-04-01_00:00:00")
wrf_file2 = systemfunc("ls /data/muhdomer/Build_WRF/Martilli/glade/u/home/martilli/wrf_3.8/WRF_Forest/run/wrfout_d05_2016-04-01_00:00:00")
; This needs to have a ".nc" appended, so just do it.
a = addfiles(wrf_file1+".nc","r")
b = addfiles(wrf_file2+".nc","r")
ListSetType(b,"cat")
LH1 = addfiles_GetVar(a,wrf_file1,"T2")
LH2 = addfiles_GetVar(b,wrf_file2,"T2")
;LH1 = wrf_user_getvar(a,"rh2",-1)
;LH2 = wrf_user_getvar(b,"rh2",-1)
;printVarSummary(LH1)
;printVarSummary(LH2)
LH = LH1 - LH2
;LH=LH-273.16
; What times and how many time steps are in the data set?
; times = wrf_user_list_times(a) ; get times in the file
times = a[:]->Times
ntimes = dimsizes(times) ; number of times in the file
ntimes_plot = ntimes(0)
lh3d = new((/29,129,210/), float)
lh2d = new((/129,210/), float)
c = 0
do i=time_serie, ntimes_plot-2,24
lh2d(:,:) = LH(i,:,:)
lh3d(c,:,:) = lh2d(:,:)
c = c + 1
end do
lh_avg = dim_avg_n(lh3d,0)
; We generate plots, but what kind do we prefer?
type = "png"
wks = gsn_open_wks(type,"LH_d05_diff_3_1_"+hour)
; gsn_define_colormap(wks,"spread_15lev")
; Set some basic resources
res = True
;res at ContourParameters = (/ 273., 310., 5./)
res at cnFillOn = True
res at cnLevelSelectionMode = "ExplicitLevels"
; res at cnLevels=(/-0.4,-0.2,0,0.4,0.8,1.0,1.4,1.8,2.0,2.4/)
res at cnLevels=(/-0.2,0,0.4,0.6,0.8,1.0,1.2,1.4,1.6,1.8,2.0,2.2,2.4/)
;res at cnLevelSelectionMode = "ExplicitLevels"
;res at cnLevels = (/-50,-45,-40,-35,-30,-25,-20,-15,-10,-5,0/)
;res at cnFillColors = (/12,11,10,9,8,7,6,5,4,3,2,2/)
;res at MainTitle = "REAL-TIME WRF"
;res at cnLevels = (/0,5,10,15,20,25,30,35,40,45,50/)
;res at cnLevels = (/0,10,20,30,40,50,60,70,80,90,100/)
;res at cnFillColors = (/2,2,3,4,5,6,7,8,9,10,11,12,13/)
;res at mpOutlineOn = False
res at InitTime = False
res at Footer = False
res at gsnMaximize = True
res at gsnPaperOrientation = "landscape"
pltres = True
pltres at NoTitles = True
pltres at gsnMaximize = True
pltres at FramePlot = False
pltres at PanelPlot = True
mpres = True
;mpres at mpFillOn = False
;mpres at mpGeophysicalLineColor = "black"
;mpres at mpGeophysicalLineThicknessF = 0.0 ; for type = "png"
; mpres at mpDataBaseVersion="Ncarg4_1"
; mpres at mpDataSetName="Earth..4"
; mpres at mpDataResolution="Finest"
; mpres at mpOutlineOn=False
lh_avg at description = "UHI"
;lh_avg at description = "Sensible Heat Flux"
lh_avg at units = "C"
; Plotting options for LH
opts = res
opts at cnFillOn = True
;opts at ContourParameters = (/ -50., 0., 5. /)
;opts at gsnSpreadColorEnd = -3 ; End third from the last color in color map
contour_lh = wrf_contour(a[0],wks,lh_avg,opts)
delete(opts)
; MAKE PLOTS
plot = wrf_map_overlays(a[0],wks,contour_lh,pltres,mpres)
shape_files = (/"/data/muhdomer/NSCC/scratch/Shapefiles/SGP_adm0.shp","/data/muhdomer/NSCC/scratch/Shapefiles/MYS_adm0.shp","/data/muhdomer/NSCC/scratch/Shapefiles/IDN_adm0.shp"/)
plot = add_coastline_sgmyid(shape_files,wks,plot,"black")
draw(plot)
frame(wks)
Muhammad Omer Mughal
MSc BSc Mechanical Engineering
PhD Research Scholar
Remote Sensing and Satellite Research Group
Department of Imaging and Applied Physics
Curtin University
Curtin University
Tel | +61 8 9266 7962
Fax | +61 8 9266 2377
Mobile | 0470 237 525
Email | m.mughal1 at postgrad.curtin.edu.au<mailto:m.lynch at curtin.edu.au>
Web | http://curtin.edu.au<http://curtin.edu.au/>
Curtin University is a trademark of Curtin University of Technology.
CRICOS Provider Code 00301J (WA), 02637B (NSW)
________________________________
From: Mary Haley <haley at ucar.edu>
Sent: Saturday, 14 October 2017 4:16:59 AM
To: Muhammad Omer Mughal
Cc: ncl-talk at ucar.edu
Subject: Re: [ncl-talk] Help Required in plotting average 2m temperature for 24 hours for 30 days in NCL
Muhammad,
I think this script is slightly more complicated than it needs to be.
It looks like your tc2 array is dimensioned ((/30,129,210/), and you want to skip indexes 0-11 of the leftmost dimension so that you are only averaging across indexes 12-29
? If so, then:
tc2_avg = dim_avg_n(tc2(12:,:,:),0)
This will give you back a 129 x 210 array, where each element is an average of the last 18 timesteps.
--Mary
On Wed, Oct 11, 2017 at 2:53 AM, Muhammad Omer Mughal <m.mughal1 at postgrad.curtin.edu.au<mailto:m.mughal1 at postgrad.curtin.edu.au>> wrote:
Hi Mary
Kindly see the script below which I use with the bash script to determine the average of 2M TEMPERATURE contained in the whole of WRF out directory. I intend to skip the first 12 hours (spin up time) and I also tend to plot the final average in one plot. Kindly let me know if you can help me out
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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "/data/muhdomer/NSCC/scratch/WRF_myfunc.ncl"
begin
DATADir = "./"
FILES = systemfunc (" ls -1 " + DATADir + "wrfout_d05* ")
a = addfiles(FILES+".nc","r")
tc2 = wrf_user_getvar(a,"T2",-1)
tc2 = tc2-273.16
times = a[:]->Times
ntimes = dimsizes(times) ; number of times in the file
ntimes_plot = ntimes(0) ; abandon the first 12 hours
lh3d = new((/30,129,210/), float)
lh2d = new((/129,210/), float)
c = 0
do i=time_serie, ntimes_plot-2,24
lh2d(:,:) = tc2(i,:,:)
lh3d(c,:,:) = lh2d(:,:)
c = c + 1
end do
lh_avg = dim_avg_n(lh3d,0)
type = "x11"
wks= gsn_open_wks(type,"temp_"+hour)
gsn_define_colormap(wks,"spread_15lev")
res = True
res at cnLevelSelectionMode = "ExplicitLevels"
res at cnLevels=(/24,25,26,27,28,29,30,31,32,33,34,35,36/)
res at InitTime = False
res at Footer = False
res at gsnMaximize = True
res at gsnPaperOrientation = "landscape"
pltres = True
pltres at NoTitles = True
pltres at gsnMaximize = True
pltres at FramePlot = False
pltres at PanelPlot = True
mpres = True
mpres at mpFillOn = False
mpres at mpGeophysicalLineColor = "black"
mpres at mpGeophysicalLineThicknessF = 0.0 ; for type = "png"
mpres at mpDataBaseVersion="Ncarg4_1"
mpres at mpDataSetName="Earth..4"
mpres at mpDataResolution="Finest"
mpres at mpOutlineOn=False
lh_avg at description = "temperature averaged"
lh_avg at units = "degree C"
opts = res
opts at cnFillOn = True
contour_lh = wrf_contour(a[0],wks,lh_avg,opts)
delete(opts)
plot = wrf_map_overlays(a[0],wks,contour_lh,pltres,mpres)
shape_files = (/"/data/muhdomer/NSCC/scratch/Shapefiles/MYS_adm0.shp","/data/muhdomer/NSCC/scratch/Shapefiles/IDN_adm0.shp","/data/muhdomer/NSCC/scratch/Shapefiles/MP14_REGION_NO_SEA_PL.shp"/)
plot = add_coastline_sgmyid(shape_files,wks,plot,"black")
draw(plot)
frame(wks)
end
Regards
Muhammad Omer
Muhammad Omer Mughal
MSc BSc Mechanical Engineering
PhD Research Scholar
Remote Sensing and Satellite Research Group
Department of Imaging and Applied Physics
Curtin University
Curtin University
Tel | +61 8 9266 7962
Fax | +61 8 9266 2377
Mobile | 0470 237 525
Email | m.mughal1 at postgrad.curtin.edu.au<mailto:m.lynch at curtin.edu.au>
Web | http://curtin.edu.au<http://curtin.edu.au/>
Curtin University is a trademark of Curtin University of Technology.
CRICOS Provider Code 00301J (WA), 02637B (NSW)
_______________________________________________
ncl-talk mailing list
ncl-talk at ucar.edu<mailto:ncl-talk at ucar.edu>
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171016/790f4fac/attachment.html>
More information about the ncl-talk
mailing list