[ncl-talk] Overlay Wind plot issues.

Kwesi Quagraine starskykwesi at gmail.com
Thu Jan 19 09:44:41 MST 2017


Thank you Mary.



On Jan 19, 2017 6:43 PM, "Mary Haley" <haley at ucar.edu> wrote:

> You can simply change the value after you create the first plot. Right now
> you have:
>
>
> ;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)
>    . . .
>
> So right before you create plot(0), set the resource to True:
>
>   ;make plots
>   res at gsnLeftString           = "NCEP Wind"
>   res at gsnAddCyclic            = True
>   plot(0)   = gsn_csm_contour_map_overlay(wks,uwnd,uwnd,res,sres)   ;
> create the temperature plot
>
> and then before you create the rest of the plots. set it to False.  It
> will then stay equal to False for the rest of the code:
>
>   res at gsnLeftString           = "ERA-INT Wind"
>   res at gsnAddCyclic            = False
>   plot(1)   = gsn_csm_contour_map_overlay(wks,ua,ua,res,sres)
>
>
> On Thu, Jan 19, 2017 at 8:40 AM, Kwesi Quagraine <starskykwesi at gmail.com>
> wrote:
>
>> Thanks Mary, any idea how I can set res at gsnAddCyclic to False for only
>> that plot? since the res is defined for all the other plots in my script?
>>
>> Thanks in advance.
>>
>> On Thu, Jan 19, 2017 at 5:37 PM, Mary Haley <haley at ucar.edu> wrote:
>>
>>> Dear Kwesi,
>>>
>>> Your uwnd data is global (goes from 0 to 357.5), but it doesn't quite go
>>> all the way to to 360, so you need to set res at gsnAddCyclic to True for
>>> that data so the longitude cyclic point gets added.
>>>
>>> For the other plots, which appear to all be regional, set
>>> res at gsnAddCyclic to False.
>>>
>>> --Mary
>>>
>>>
>>> On Thu, Jan 19, 2017 at 3:27 AM, Kwesi Quagraine <starskykwesi at gmail.com
>>> > wrote:
>>>
>>>> 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 <+233%2026%20617%203582>
>>>> Skype: quagraine_cwasi
>>>> Twitter: @Pkdilly
>>>>
>>>>
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at ucar.edu
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>>
>>>
>>
>>
>> --
>> 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 <+233%2026%20617%203582>
>> Skype: quagraine_cwasi
>> Twitter: @Pkdilly
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170119/3f1c2415/attachment.html 


More information about the ncl-talk mailing list