[ncl-talk] Overlay Wind plot issues.
Kwesi Quagraine
starskykwesi at gmail.com
Thu Jan 19 03:27:31 MST 2017
Hello NCLers, I am trying to make a wind overlay plot. I have been able to
make the plot from the overlay template.
[1] first plot in panel seem to only appear with some empty line through it
when gsnAddCyclic = False (please see attached) whereas vice-versa distorts
the other plots.
[2] Also I want to increase the font size of the contour labels.
Any help will be much appreciated.
Kwesi.
Here's my script;
begin
fili = "700_ncepJULY.nc"
a = addfile (diri+fili , "r")
uwnd = a->uwnd(0,0,:,:)
printVarSummary(uwnd)
fili = "700_eraJULY.nc"
b = addfile (diri+fili , "r")
ua = b->ua(0,0,:,:)
fili = "700_ccmaJULY.nc"
c = addfile (diri+fili , "r")
ua1 = c->ua(0,0,:,:)
printVarSummary(ua1)
fili = "700_cnrmJULY.nc"
d = addfile (diri+fili , "r")
ua2 = d->ua(0,0,:,:)
fili = "700_ichecJULY.nc"
e = addfile (diri+fili , "r")
ua3 = e->ua(0,0,:,:)
fili = "700_ipslJULY.nc"
f = addfile (diri+fili , "r")
ua4 = f->ua(0,0,:,:)
fili = "700_mirocJULY.nc"
g = addfile (diri+fili , "r")
ua5 = g->ua(0,0,:,:)
fili = "700_mohcJULY.nc"
h = addfile (diri+fili , "r")
ua6 = h->ua(0,0,:,:)
fili = "700_mpiJULY.nc"
i = addfile (diri+fili , "r")
ua7 = i->ua(0,0,:,:)
fili = "700_nccJULY.nc"
j = addfile (diri+fili , "r")
ua8 = j->ua(0,0,:,:)
wks = gsn_open_wks("pdf","julyuwnd at 700hPa") ; send graphics to
PNG file
;set resources
res = True
res at mpFillOn = False
res at mpMaxLatF = 25. ; specify the plot domain
res at mpMinLatF = 0. ;
res at mpMinLonF = -20. ;
res at mpMaxLonF = 20. ;
res at mpOutlineOn = True ; turn the map outline on
res at gsnDraw = False ; do not draw the plot
res at gsnFrame = False ; do not advance the frame
res at gsnAddCyclic =False
res at pmTickMarkDisplayMode = "Always" ; nicer tickmarks
res at mpCenterLonF = 0. ; Centers the plot at 180
res at cnLevelSelectionMode = "ExplicitLevels" ; use explicit levels
res at cnLevels = ispan(-16,4,2) ; set the contour levels
res at cnLineLabelsOn = True ; do not use line labels
res at cnFillOn = True ; color fill
res at cnLinesOn = True ; do not draw contour lines
res at cnFillPalette = "MPL_RdYlBu"
res at cnInfoLabelOn = False ; turn off contour info
label
;res at cnFillDrawOrder = "Predraw"
res at lbLabelBarOn = False ; No single label bar
sres = True ; set up a second resource
list
sres at cnLevelSelectionMode = "ExplicitLevels" ; use explicit levels
sres at cnLevels = ispan(-16,4,2) ; set the contour levels
sres at gsnAddCyclic = False
sres at pmTickMarkDisplayMode = "Always" ; nicer tickmarks
sres at cnInfoLabelOn = False ; turn off contour info
label
sres at tiXAxisOn = False
sres at tiYAxisOn = False
sres at tiMainOn = False
;define plots
plot = new( 10, graphic)
;make plots
res at gsnLeftString = "NCEP Wind"
plot(0) = gsn_csm_contour_map_overlay(wks,uwnd,uwnd,res,sres) ;
create the temperature plot
res at gsnLeftString = "ERA-INT Wind"
plot(1) = gsn_csm_contour_map_overlay(wks,ua,ua,res,sres)
res at gsnLeftString = "CCCMA Wind"
plot(2) = gsn_csm_contour_map_overlay(wks,ua1,ua1,res,sres)
res at gsnLeftString = "CNRM Wind"
plot(3) = gsn_csm_contour_map_overlay(wks,ua2,ua2,res,sres)
res at gsnLeftString = "ICHEC Wind"
plot(4) = gsn_csm_contour_map_overlay(wks,ua3,ua3,res,sres)
res at gsnLeftString = "IPSL Wind"
plot(5) = gsn_csm_contour_map_overlay(wks,ua4,ua4,res,sres)
res at gsnLeftString = "MIROC Wind"
plot(6) = gsn_csm_contour_map_overlay(wks,ua5,ua5,res,sres)
res at gsnLeftString = "MOHC Wind"
plot(7) = gsn_csm_contour_map_overlay(wks,ua6,ua6,res,sres)
res at gsnLeftString = "MPI Wind"
plot(8) = gsn_csm_contour_map_overlay(wks,ua7,ua7,res,sres)
res at gsnLeftString = "NCC-NOR Wind"
plot(9) = gsn_csm_contour_map_overlay(wks,ua8,ua8,res,sres)
resP = True ; panel options
resP at gsnMaximize = True ; maximize image
resP at gsnPanelLabelBar = True ; Add common label bar
resP at txString = "AEJ July"
gsn_panel(wks,plot,(/5,2/),resP)
end
var summary;
Variable: uwnd
Type: float
Total Size: 42048 bytes
10512 values
Number of Dimensions: 2
Dimensions and sizes: [lat | 73] x [lon | 144]
Coordinates:
lat: [90..-90]
lon: [ 0..357.5]
Number Of Attributes: 14
level : 700
time : 1678608
long_name : Monthly mean u wind
units : m/s
_FillValue : -9.96921e+36
missing_value : -9.96921e+36
precision : 2
least_significant_digit : 1
var_desc : u-wind
level_desc : Multiple levels
statistic : Mean
parent_stat : Other
dataset : NCEP Reanalysis Derived Products
actual_range : ( -68.19482, 124.4 )
Variable: ua1 (same for the rest)
Type: float
Total Size: 964004 bytes
241001 values
Number of Dimensions: 2
Dimensions and sizes: [lat | 401] x [lon | 601]
Coordinates:
lat: [ -15.. 25]
lon: [ -30.. 30]
Number Of Attributes: 8
plev : 70000
time : 15192.5
standard_name : eastward_wind
long_name : Eastward Wind
units : m s-1
_FillValue : 1e+20
missing_value : 1e+20
cell_methods : time: mean
--
Try not to become a man of success but rather a man of value-Albert Einstein
University of Cape Coast|College of Agriculture and Natural Sciences|Department
of Physics|
Team Leader|Recycle Up! Ghana|Technology Without Borders|
Other emails: kwesi.quagraine at ucc.edu.gh|kwesi.quagraine at teog.de|
Mobile: +233266173582
Skype: quagraine_cwasi
Twitter: @Pkdilly
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170119/f1cc8fb0/attachment-0001.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: julyuwnd at 700hPa.pdf
Type: application/pdf
Size: 952551 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170119/f1cc8fb0/attachment-0001.pdf
More information about the ncl-talk
mailing list