[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