[ncl-talk] Changing x-axis and y-axis range

Dipti Sharma sdipti596 at gmail.com
Tue May 16 13:37:22 MDT 2017


Thank you Adam.


Dipti

On Tue, May 16, 2017 at 3:08 PM, Adam Phillips <asphilli at ucar.edu> wrote:

> Hi Dipti,
> Yes. You would need to set a couple of resources that control the range of
> the x-axis:
> r at trXMinF = 0.0
> r at trXMaxF = 0.001.
>
> If that does not fix the issue please let the ncl-talk email list know.
> Adam
>
> On Tue, May 16, 2017 at 10:39 AM, Dipti Sharma <sdipti596 at gmail.com>
> wrote:
>
>> All,
>>
>> The below script works fine but I can not change the x-axix range/y-axis range (values).
>> it is 0. to 0.5 now but I want to make it going from 0 to 0.001 for instance. Is it possible to modify it?
>>
>> Thank you in advance for your help on this.
>>
>> Dipti
>>
>> ;Here is the script
>>
>> ;spec_4.ncl
>> ;
>> ; Concepts illustrated:
>> ;   - Calculating and plotting spectra
>> ;   - Calculating confidence intervals
>> ;   - Making an axis logarithmic in an XY plot
>> ;************************************************
>> 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/csm/shea_util.ncl"
>> ;************************************************
>> begin
>> ;************************************************
>> ; variable and file handling
>> ;************************************************
>>    fn  = "SOI_Darwin.nc" ; define filename
>>    in  = addfile(fn,"r")                                  ; open netcdf file
>>    soi  = in->DSOI                                        ; get data
>> ;************************************************
>> ; set function arguments
>> ;************************************************
>>   d   = 0    ; detrending opt: 0=>remove mean 1=>remove mean + detrend
>>   sm  = 21   ; smooth: should be at least 3 and odd
>>   pct = 0.10 ; percent taper: (0.0 <= pct <= 1.0) 0.10 common.
>> ;************************************************
>> ; calculate spectrum
>> ;************************************************
>>   sdof = specx_anal(soi,d,sm,pct)
>> ;************************************************
>> ; calculate confidence interval [here 5 and 95%]
>> ; return 4 curves to be plotted
>> ;************************************************
>>   splt = specx_ci (sdof, 0.05, 0.95)
>> ;************************************************
>> ; plotting
>> ;************************************************
>>    wks  = gsn_open_wks("ps","spec")              ; Opens a ps file
>>
>>    r = True                                      ; plot mods desired
>>    r at tiMainString = "SOI"                        ; title
>>    r at tiXAxisString = "Frequency (cycles/month)"  ; xaxis
>>    r at tiYAxisString = "Variance"                  ; yaxis
>> ;***********************************************
>> ; Generate log plot showing "red noise" confidence bounds
>> ; (a) log scaling and (b) the Band Width
>> ;***********************************************
>>    r at trYLog              = True                 ; log scaling
>>    r at trYMinF             = 0.10                 ; manually set lower limit
>>    r at trYMaxF             = 30.0                 ;   "          upper
>>    r at gsnFrame            = False                ; do not advance frame
>>    plot  = gsn_csm_xy(wks,sdof at frq, splt,r)
>>
>>    xf   = (/0.40, 0.40+sdof at bw/)                ; Create band width line
>>    ys   = 0.75*max(sdof at spcx)                   ; 75% up Y axis
>>    yv   = (/ys,ys/)
>>    rpl  = True                                  ; resources for polyline
>>    rpl at gsLineThicknessF  = 2                    ; Define line thickness
>>    gsn_polyline(wks,plot,xf,yv,rpl)             ; Draw BandWidth
>>
>>    txres= True                                  ; label BW line
>>    txres at txFontHeightF = 0.015                  ; font height
>>    txres at txJust        = "CenterLeft"           ; Set lable location
>>    gsn_text(wks,plot,"BW",0.41+sdof at bw,ys,txres); Label
>>    frame (wks)
>> end
>>
>>
>>
>> On Sat, May 13, 2017 at 2:00 PM, <ncl-talk-request at ucar.edu> wrote:
>>
>>> Send ncl-talk mailing list submissions to
>>>         ncl-talk at ucar.edu
>>>
>>> To subscribe or unsubscribe via the World Wide Web, visit
>>>         http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>> or, via email, send a message with subject or body 'help' to
>>>         ncl-talk-request at ucar.edu
>>>
>>> You can reach the person managing the list at
>>>         ncl-talk-owner at ucar.edu
>>>
>>> When replying, please edit your Subject line so it is more specific
>>> than "Re: Contents of ncl-talk digest..."
>>>
>>>
>>> Today's Topics:
>>>
>>>    1. Re: Overlay profile of topography (Appo derbetini)
>>>
>>>
>>> ----------------------------------------------------------------------
>>>
>>> Message: 1
>>> Date: Fri, 12 May 2017 19:17:56 +0100
>>> From: Appo derbetini <appopson4 at gmail.com>
>>> Subject: Re: [ncl-talk] Overlay profile of topography
>>> To: Mary Haley <haley at ucar.edu>
>>> Cc: NCL <ncl-talk at ucar.edu>
>>> Message-ID:
>>>         <CAHVpkHj4MrBrGt1Yrcq1RGkVaF1Fy-fG5e2x6zcg9LOvT77S9Q at mail.gm
>>> ail.com>
>>> Content-Type: text/plain; charset="utf-8"
>>>
>>> Thank you Mary
>>>
>>>
>>>
>>> 2017-05-12 16:07 GMT+01:00 Mary Haley <haley at ucar.edu>:
>>>
>>> > Hi Appo,
>>> >
>>> > The first thing I tell people when they are trying to overlay plots and
>>> > not seeing both plots is to draw them individually before you call
>>> > overlay.  This way you can make sure the individual plots look okay
>>> before
>>> > you get to the overlay call.
>>> >
>>> > All this involves in your case is commenting out the gsnDraw and
>>> gsnFrame
>>> > settings.
>>> >
>>> > When you do this, you'll see that the Y axis on the wind profile ranges
>>> > from 20000 to 100000, while on the humidity plot it ranges from 200 to
>>> 1000.
>>> >
>>> > It's not well-advertised, but under the hood, gsn_csm_pres_hgt converts
>>> > the pressure axis to mb before plotting. (I've fixed the documentation
>>> to
>>> > make this more clear.)
>>> >
>>> > gsn_csm_vector doesn't do this, so you will need to scale the pressure
>>> > values yourself before plotting.
>>> >
>>> > I suggest scaling all the arrays, just to be consistent, even for r00:
>>> >
>>> > ;---Convert levels from Pa to hPa for consistent plotting
>>> >
>>> >    r00&plev = r00&plev * 0.01
>>> >    udiv00&plev = udiv00&plev * 0.01
>>> >    omega00&plev = omega00&plev * 0.01
>>> >    r00&plev at units  = "mb"
>>> >    udiv00&plev at units  = "mb"
>>> >    omega00&plev at units = "mb"
>>> >
>>> > --Mary
>>> >
>>> >
>>> > On Wed, May 10, 2017 at 3:02 AM, Appo derbetini <appopson4 at gmail.com>
>>> > wrote:
>>> >
>>> >> Dear Mary,
>>> >>
>>> >> I am very surprise.
>>> >> I try the script with another dataset.
>>> >> As you can see in the attached plot, only humidity is displayed.
>>> >> Vertical profile of wind is not appearing.
>>> >>
>>> >> What i going wrong?
>>> >>
>>> >> Best regards,
>>> >>
>>> >> Appo
>>> >>
>>> >>
>>> >>
>>> >>
>>> >>
>>> >> 2017-05-08 20:22 GMT+01:00 Appo derbetini <appopson4 at gmail.com>:
>>> >>
>>> >>> Thank you very much Mary,
>>> >>> Now, I have a perfect plot.
>>> >>> Cheers,
>>> >>> Appo
>>> >>>
>>> >>> 2017-05-07 16:08 GMT+01:00 Mary Haley <haley at ucar.edu>:
>>> >>>
>>> >>>> Hi,
>>> >>>>
>>> >>>> If by "optimized for this plot" you mean is the plotting being done
>>> >>>> correctly, then I would say "yes".
>>> >>>>
>>> >>>> However, I did notice that you have two slightly different scales on
>>> >>>> your color bar. The top one goes from 30 to 90 in steps of 5, while
>>> the
>>> >>>> bottom one goes from 20 to 90 in steps of 5.
>>> >>>>
>>> >>>> If somebody is not paying attention when looking at these images,
>>> they
>>> >>>> may not realize that the bottom purple color represents values < 30
>>> in the
>>> >>>> top plot, while it represents values < 20 in the bottom plot.
>>> >>>>
>>> >>>> I recommend using the same color scale for both plots, and then
>>> using a
>>> >>>> single panel labelbar so that it's more clear what the colors mean.
>>> >>>>
>>> >>>> To make sure you use the same contour levels for both plots, set
>>> these
>>> >>>> resources:
>>> >>>>
>>> >>>>  res at cnLevelSelectionMode = "ManualLevels" ; manually set contour
>>> levels
>>> >>>>  res at cnMinLevelValF       = 20.            ; minimum contour level
>>> >>>>  res at cnMaxLevelValF       = 90.            ; maximum contour level
>>> >>>>  res at cnLevelSpacingF        = 5.
>>> >>>>
>>> >>>> You then you need to turn off the individual labelbars for both
>>> plots, so you can set a common labelbar in the panel:
>>> >>>>
>>> >>>>  res at lbLabelBarOn = False
>>> >>>>
>>> >>>> Now, in the panel resources, set:
>>> >>>>
>>> >>>> pres at gsnPanelLabelBar = True
>>> >>>> pres2 at lbOrientation   = "vertical" ; set this if you want a
>>> vertical label bar
>>> >>>>
>>> >>>> If you need more customization of your panel plot, then we have
>>> some examples on our panel page:
>>> >>>>
>>> >>>> http://www.ncl.ucar.edu/Applications/panel.shtml
>>> >>>>
>>> >>>> --Mary
>>> >>>>
>>> >>>>
>>> >>>>
>>> >>>>
>>> >>>> On Fri, May 5, 2017 at 4:07 AM, Appo derbetini <appopson4 at gmail.com
>>> >
>>> >>>> wrote:
>>> >>>>
>>> >>>>> Dear Mary,
>>> >>>>> I applied your suggestions.
>>> >>>>> Now, things are going like a charm.
>>> >>>>> Here attached figure obtained and script.
>>> >>>>> I guess that this script is optimized  for this plot?
>>> >>>>>
>>> >>>>> Thank you very much
>>> >>>>>
>>> >>>>> begin
>>> >>>>>
>>> >>>>>     in1  = addfile("uvwr.mean.JJAS.nc", "r")
>>> >>>>> ; open netcdf file
>>> >>>>>     input4  = addfile("topo.nc", "r")
>>> >>>>>
>>> >>>>>    u00     =  in1->u
>>> >>>>>    v00     =  in1->v
>>> >>>>>    w00     =  in1->w
>>> >>>>>    rhum00  =  in1->r
>>> >>>>>
>>> >>>>>    topo1 = input4->HT
>>> >>>>>
>>> >>>>>
>>> >>>>>    lat = in1->latitude
>>> >>>>>    lon = in1->longitude
>>> >>>>>
>>> >>>>>
>>> >>>>>    dv00    = uv2dvF (u00, v00)
>>> >>>>>    uvd00   = dv2uvF (dv00)
>>> >>>>>    ucomp_div00 = uvd00(0, :, : , :, :)
>>> >>>>>    copy_VarCoords(u00, ucomp_div00)
>>> >>>>>
>>> >>>>>
>>> >>>>>    sellat = 3.0
>>> >>>>>    latmin = -5.0
>>> >>>>>    latmax = 3.0
>>> >>>>>    lonmin = min(lon)
>>> >>>>>    lonmax = 20.0
>>> >>>>>
>>> >>>>>    sellat2 = -1.0
>>> >>>>>    latmin2 = -4.0
>>> >>>>>    latmax2 = -1.0
>>> >>>>>    lonmin = min(lon)
>>> >>>>>    lonmax = 20.0
>>> >>>>>
>>> >>>>>    udiv00 = dim_avg_n_Wrap(ucomp_div00(:, :, {latmin:latmax},
>>> >>>>> {lonmin:lonmax}), 2)
>>> >>>>>    omega00 = dim_avg_n_Wrap(w00(:, :, {latmin:latmax},
>>> >>>>> {lonmin:lonmax}), 2)
>>> >>>>>    r00 = dim_avg_n_Wrap(rhum00(:, :, {latmin:latmax},{lonmin:lonmax
>>> }),
>>> >>>>> 2)
>>> >>>>>
>>> >>>>>    udiv002 = dim_avg_n_Wrap(ucomp_div00(:, :, {latmin2:latmax2},
>>> >>>>> {lonmin:lonmax}), 2)
>>> >>>>>    omega002 = dim_avg_n_Wrap(w00(:, :, {latmin2:latmax2},
>>> >>>>> {lonmin:lonmax}), 2)
>>> >>>>>    r002 = dim_avg_n_Wrap(rhum00(:, :,
>>> {latmin2:latmax2},{lonmin:lonmax}),
>>> >>>>> 2)
>>> >>>>>
>>> >>>>>
>>> >>>>>
>>> >>>>>    elev = 1013.25*(1 - topo1*0.0065/288.15)^5.25145
>>> >>>>>
>>> >>>>>    copy_VarCoords(topo1, elev)
>>> >>>>>
>>> >>>>> ;  topo = dim_avg_n_Wrap(elev({latmin:latmax},{lonmin:lonmax}), 0)
>>> >>>>>
>>> >>>>>    topo = elev({sellat}, {lonmin:lonmax})
>>> >>>>>    topo2 = elev({sellat2}, {lonmin:lonmax})
>>> >>>>>
>>> >>>>> ;---create plot
>>> >>>>>    wks   = gsn_open_wks ("eps", "Overlay_Rhum_udiv_JJAS_Timmea
>>> n_5S3N")
>>> >>>>> ; open workstation
>>> >>>>>    plot = new(2, graphic)
>>> >>>>>    res                     = True                  ; plot mods
>>> desired
>>> >>>>>    res at gsnDraw             = False    ; turn off draw
>>> >>>>>    res at gsnFrame            = False    ; turn off frame
>>> >>>>>    res at gsnMaximize         = True
>>> >>>>>
>>> >>>>>    res at gsnRightString      = ""
>>> >>>>>    res at gsnLeftString       = ""
>>> >>>>>
>>> >>>>>    cnres = res
>>> >>>>>    cnres at cnFillOn          = True               ; turn on color
>>> fill
>>> >>>>>    cnres at cnFillPalette     = "MPL_rainbow"
>>> >>>>>    cnres at lbOrientation     = "Vertical"
>>> >>>>>    cnres at pmLabelBarOrthogonalPosF = 0.08
>>> >>>>>    cnres at tiYAxisString      = "Pressure (hPa)"
>>> >>>>>
>>> >>>>>    vcres = res
>>> >>>>>    vcres at vcRefMagnitudeF = 10.0                ; define vector
>>> ref mag
>>> >>>>>    vcres at vcRefLengthF    = 0.1              ; define length of
>>> vec ref
>>> >>>>>    vcres at vcGlyphStyle    = "CurlyVector"      ; turn on curley
>>> vectors
>>> >>>>>    vcres at vcMinDistanceF  = 0.04               ; thin out vectors
>>> >>>>>    vcres at vcMapDirection  = False
>>> >>>>>    vcres at vcLineArrowThicknessF   = 3.0
>>> >>>>>    vcres at vcVectorDrawOrder  = "Draw"        ; draw vectors last
>>> >>>>>
>>> >>>>>
>>> >>>>>
>>> >>>>>    contour_plota1  =  gsn_csm_pres_hgt(wks, r00(0, :,
>>> >>>>> {lonmin:lonmax}), cnres)
>>> >>>>>    vector_plota1   = gsn_csm_vector(wks, udiv00(0, :,
>>> >>>>> {lonmin:lonmax}), -200.0*omega00(0, :, {lonmin:lonmax}), vcres )
>>> >>>>>
>>> >>>>>    contour_plota2  =  gsn_csm_pres_hgt(wks, r002(0, :,
>>> >>>>> {lonmin:lonmax}), cnres)
>>> >>>>>    vector_plota2   = gsn_csm_vector(wks, udiv002(0, :,
>>> >>>>> {lonmin:lonmax}), -200.0*omega002(0, :, {lonmin:lonmax}), vcres )
>>> >>>>>
>>> >>>>>
>>> >>>>> ;---Add topo field using a filled polygon.
>>> >>>>>    getvalues contour_plota1
>>> >>>>>       "trYMinF" : ymin
>>> >>>>>       "trYMaxF" : ymax
>>> >>>>>    end getvalues
>>> >>>>>
>>> >>>>> ;---Create new X,Y arrays that form a closed polygon.
>>> >>>>>    nlon  = dimsizes(topo&lon)
>>> >>>>>    xtopo = new(nlon+3,typeof(topo&lon))
>>> >>>>>    ytopo = new(nlon+3,typeof(topo))
>>> >>>>>
>>> >>>>>    xtopo(0:nlon-1) = topo&lon
>>> >>>>>    ytopo(0:nlon-1) = topo
>>> >>>>>    xtopo(nlon)     = topo&lon(nlon-1)
>>> >>>>>    ytopo(nlon)     = ymax          ; Use actual Y max of contour
>>> plot
>>> >>>>>    xtopo(nlon+1)   = topo&lon(0)
>>> >>>>>    ytopo(nlon+1)   = ymax
>>> >>>>>    xtopo(nlon+2)   = topo&lon(0)   ; This last point closes
>>> >>>>>    ytopo(nlon+2)   = topo(0)       ; the polygon.
>>> >>>>>
>>> >>>>>
>>> >>>>>
>>> >>>>>
>>> >>>>> ;---Add topo field using a filled polygon.
>>> >>>>>    getvalues contour_plota2
>>> >>>>>       "trYMinF" : ymin2
>>> >>>>>       "trYMaxF" : ymax2
>>> >>>>>    end getvalues
>>> >>>>>
>>> >>>>> ;---Create new X,Y arrays that form a closed polygon.
>>> >>>>>    nlon  = dimsizes(topo2&lon)
>>> >>>>>    xtopo2 = new(nlon+3,typeof(topo2&lon))
>>> >>>>>    ytopo2 = new(nlon+3,typeof(topo2))
>>> >>>>>
>>> >>>>>    xtopo2(0:nlon-1) = topo2&lon
>>> >>>>>    ytopo2(0:nlon-1) = topo2
>>> >>>>>    xtopo2(nlon)     = topo2&lon(nlon-1)
>>> >>>>>    ytopo2(nlon)     = ymax2          ; Use actual Y max of contour
>>> plot
>>> >>>>>    xtopo2(nlon+1)   = topo2&lon(0)
>>> >>>>>    ytopo2(nlon+1)   = ymax2
>>> >>>>>    xtopo2(nlon+2)   = topo2&lon(0)   ; This last point closes
>>> >>>>>    ytopo2(nlon+2)   = topo2(0)       ; the polygon.
>>> >>>>>
>>> >>>>> ;---Add the polygon to the contour plot.
>>> >>>>>    gnres = True
>>> >>>>>    gnres at gsFillColor = "gray20"
>>> >>>>>
>>> >>>>>    id = gsn_add_polygon(wks, contour_plota1, xtopo, ytopo, gnres)
>>> >>>>>    i2 = gsn_add_polygon(wks, contour_plota2, xtopo2, ytopo2, gnres)
>>> >>>>>
>>> >>>>>     overlay(contour_plota1, vector_plota1)
>>> >>>>>     overlay(contour_plota2, vector_plota2)
>>> >>>>>
>>> >>>>>      plot(0)  = contour_plota1
>>> >>>>>      plot(1)  = contour_plota2
>>> >>>>>
>>> >>>>>    ;draw(contour_plot)     ; This draws everything
>>> >>>>>
>>> >>>>>
>>> >>>>>
>>> >>>>>
>>> >>>>>    pres = True
>>> >>>>>    pres at gsnMaximize = True
>>> >>>>>    gsn_panel(wks, plot, (/2, 1/), pres)
>>> >>>>>
>>> >>>>>
>>> >>>>>
>>> >>>>> end
>>> >>>>>
>>> >>>>>
>>> >>>>>
>>> >>>>> 2017-05-05 0:39 GMT+01:00 Mary Haley <haley at ucar.edu>:
>>> >>>>>
>>> >>>>>> I'm so sorry I didn't get back to you sooner on this.  I've been
>>> >>>>>> swamped and then was out-of-town for a few days.
>>> >>>>>>
>>> >>>>>> It's a little hard to follow what's going on in this script
>>> without
>>> >>>>>> being able to run it.
>>> >>>>>>
>>> >>>>>> I think the issue is that you are using the "res" variable in your
>>> >>>>>> call to gsn_panel, and "res" has both gsnDraw and gsnFrame set to
>>> False.
>>> >>>>>> This will cause *no* plot to appear when you call gsn_panel.  I'm
>>> surprised
>>> >>>>>> you are getting any plot at all.
>>> >>>>>>
>>> >>>>>> I recommend creating a new resource variable for gsn_panel,
>>> because
>>> >>>>>> you really don't want to use the one that you used for the plots
>>> themselves:
>>> >>>>>>
>>> >>>>>> pres = True
>>> >>>>>> pres at gsnMaximize = True
>>> >>>>>> gsn_panel(wks, plot, (/2, 1/), pres)
>>> >>>>>>
>>> >>>>>> If you
>>> >>>>>> ? continue to have problems, then it would help if you send me the
>>> >>>>>> two data files. If you can't send the files, then send me the
>>> image you're
>>> >>>>>> getting, and any error messages ?that the script produces.
>>> >>>>>>
>>> >>>>>> ?Thanks,
>>> >>>>>>
>>> >>>>>> --Mary
>>> >>>>>> ?
>>> >>>>>>
>>> >>>>>>
>>> >>>>>>
>>> >>>>>> On Tue, Apr 25, 2017 at 12:56 AM, Appo derbetini <
>>> appopson4 at gmail.com
>>> >>>>>> > wrote:
>>> >>>>>>
>>> >>>>>>> Good morning Mary,
>>> >>>>>>> Sorry , i forgot to add ncl-users email to the reply.
>>> >>>>>>>
>>> >>>>>>> I am attaching the modified version of the script were I'm
>>> trying to
>>> >>>>>>> create a panel plot for two transects at latitude sellat and
>>> sellat2.
>>> >>>>>>> Unfortunately, only one plot appears.
>>> >>>>>>> Really, I don't know what is wrong this script.
>>> >>>>>>>
>>> >>>>>>> Your help will be appreciated.
>>> >>>>>>>
>>> >>>>>>>
>>> >>>>>>>
>>> >>>>>>> begin
>>> >>>>>>>
>>> >>>>>>>     in1  = addfile("uvwr.mean.JJAS.nc",
>>> >>>>>>> "r")                         ; open netcdf file
>>> >>>>>>>     input4  = addfile("atanas_topo.nc", "r")
>>> >>>>>>>
>>> >>>>>>>    u00     =  in1->u
>>> >>>>>>>    v00     =  in1->v
>>> >>>>>>>    w00     =  in1->w
>>> >>>>>>>    rhum00  =  in1->r
>>> >>>>>>>
>>> >>>>>>>    topo1 = input4->HT
>>> >>>>>>>
>>> >>>>>>>
>>> >>>>>>>    lat = in1->latitude
>>> >>>>>>>    lon = in1->longitude
>>> >>>>>>>
>>> >>>>>>>
>>> >>>>>>>    dv00    = uv2dvF (u00, v00)
>>> >>>>>>>    uvd00   = dv2uvF (dv00)
>>> >>>>>>>    ucomp_div00 = uvd00(0, :, : , :, :)
>>> >>>>>>>    copy_VarCoords(u00, ucomp_div00)
>>> >>>>>>>
>>> >>>>>>>
>>> >>>>>>>    sellat = 3.0
>>> >>>>>>>    sellat2 = -2.5
>>> >>>>>>>    latmin = -5.0
>>> >>>>>>>    latmax = 3.0
>>> >>>>>>>    lonmin = min(lon)
>>> >>>>>>>    lonmax = 20.0
>>> >>>>>>>
>>> >>>>>>>    udiv00 = dim_avg_n_Wrap(ucomp_div00(:, :, {latmin:latmax},
>>> >>>>>>> {lonmin:lonmax}), 2)
>>> >>>>>>>    omega00 = dim_avg_n_Wrap(w00(:, :, {latmin:latmax},
>>> >>>>>>> {lonmin:lonmax}), 2)
>>> >>>>>>>    r00 = dim_avg_n_Wrap(rhum00(:, :,
>>> {latmin:latmax},{lonmin:lonmax}),
>>> >>>>>>> 2)
>>> >>>>>>>
>>> >>>>>>>    ;udiv00  = ucomp_div00(:, :, {sellat}, {lonmin:lonmax})
>>> >>>>>>>    ;omega00 = w00(:, :, {sellat}, {lonmin:lonmax})
>>> >>>>>>>    ;r00     = rhum00(:, :, {sellat},{lonmin:lonmax})
>>> >>>>>>>
>>> >>>>>>>
>>> >>>>>>>    elev = 1013.25*(1 - topo1*0.0065/288.15)^5.25145
>>> >>>>>>>
>>> >>>>>>>    copy_VarCoords(topo1, elev)
>>> >>>>>>>
>>> >>>>>>> ;  topo = dim_avg_n_Wrap(elev({latmin:latmax},{lonmin:lonmax}),
>>> 0)
>>> >>>>>>>
>>> >>>>>>>    topo  = elev({sellat}, {lonmin:lonmax})
>>> >>>>>>>    topo2 = elev({sellat2}, {lonmin:lonmax})
>>> >>>>>>>
>>> >>>>>>> ;---create plot
>>> >>>>>>>    wks   = gsn_open_wks ("eps", "Overlay_Rhum_udiv_topo")
>>> ;
>>> >>>>>>> open workstation
>>> >>>>>>>     plot = new(2, graphic)
>>> >>>>>>>    res                     = True                  ; plot mods
>>> >>>>>>> desired
>>> >>>>>>>    res at gsnDraw             = False    ; turn off draw
>>> >>>>>>>    res at gsnFrame            = False    ; turn off frame
>>> >>>>>>>    res at gsnMaximize         = True
>>> >>>>>>>
>>> >>>>>>>    res at gsnRightString      = ""
>>> >>>>>>>    res at gsnLeftString       = ""
>>> >>>>>>>
>>> >>>>>>>    cnres = res
>>> >>>>>>>    cnres at cnFillOn          = True               ; turn on color
>>> fill
>>> >>>>>>>    cnres at cnFillPalette     = "MPL_rainbow"
>>> >>>>>>>    cnres at lbOrientation     = "Vertical"
>>> >>>>>>>    cnres at pmLabelBarOrthogonalPosF = 0.08
>>> >>>>>>>    cnres at tiYAxisString      = "Pressure (hPa)"
>>> >>>>>>>
>>> >>>>>>>    vcres = res
>>> >>>>>>>    vcres at vcRefMagnitudeF = 10.0                ; define vector
>>> ref
>>> >>>>>>> mag
>>> >>>>>>>    vcres at vcRefLengthF    = 0.1              ; define length of
>>> vec
>>> >>>>>>> ref
>>> >>>>>>>    vcres at vcGlyphStyle    = "CurlyVector"      ; turn on curley
>>> >>>>>>> vectors
>>> >>>>>>>    vcres at vcMinDistanceF  = 0.04               ; thin out vectors
>>> >>>>>>>    vcres at vcMapDirection  = False
>>> >>>>>>>    vcres at vcLineArrowThicknessF   = 3.0
>>> >>>>>>>    vcres at vcVectorDrawOrder  = "Draw"        ; draw vectors last
>>> >>>>>>>
>>> >>>>>>>    contour_plot_1 = gsn_csm_pres_hgt(wks, r00(0, :,
>>> >>>>>>> {lonmin:lonmax}), cnres)
>>> >>>>>>>    vector_plot_1  = gsn_csm_vector(wks, udiv00(0, :,
>>> >>>>>>> {lonmin:lonmax}), -200.0*omega00(0, :, {lonmin:lonmax}), vcres )
>>> >>>>>>>
>>> >>>>>>>    contour_plot_2 = gsn_csm_pres_hgt(wks, r00(0, :,
>>> >>>>>>> {lonmin:lonmax}), cnres)
>>> >>>>>>>    vector_plot_2  = gsn_csm_vector(wks, udiv00(0, :,
>>> >>>>>>> {lonmin:lonmax}), -200.0*omega00(0, :, {lonmin:lonmax}), vcres )
>>> >>>>>>>
>>> >>>>>>>
>>> >>>>>>> ;---Add topo field using a filled polygon.
>>> >>>>>>>    getvalues contour_plot_1
>>> >>>>>>>       "trYMinF" : ymin
>>> >>>>>>>       "trYMaxF" : ymax
>>> >>>>>>>    end getvalues
>>> >>>>>>>
>>> >>>>>>>
>>> >>>>>>> ;---Add topo field using a filled polygon.
>>> >>>>>>>    getvalues contour_plot_2
>>> >>>>>>>       "trYMinF" : ymin2
>>> >>>>>>>       "trYMaxF" : ymax2
>>> >>>>>>>    end getvalues
>>> >>>>>>>
>>> >>>>>>>
>>> >>>>>>>
>>> >>>>>>> ;---Create new X,Y arrays that form a closed polygon.
>>> >>>>>>>    nlon  = dimsizes(topo&lon)
>>> >>>>>>>    xtopo = new(nlon+3,typeof(topo&lon))
>>> >>>>>>>    ytopo = new(nlon+3,typeof(topo))
>>> >>>>>>>
>>> >>>>>>>    xtopo(0:nlon-1) = topo&lon
>>> >>>>>>>    ytopo(0:nlon-1) = topo
>>> >>>>>>>    xtopo(nlon)     = topo&lon(nlon-1)
>>> >>>>>>>    ytopo(nlon)     = ymax          ; Use actual Y max of contour
>>> plot
>>> >>>>>>>    xtopo(nlon+1)   = topo&lon(0)
>>> >>>>>>>    ytopo(nlon+1)   = ymax
>>> >>>>>>>    xtopo(nlon+2)   = topo&lon(0)   ; This last point closes
>>> >>>>>>>    ytopo(nlon+2)   = topo(0)       ; the polygon.
>>> >>>>>>>
>>> >>>>>>>
>>> >>>>>>>    nlon2  = dimsizes(topo2&lon)
>>> >>>>>>>    xtopo2 = new(nlon2+3,typeof(topo2&lon))
>>> >>>>>>>    ytopo2 = new(nlon2+3,typeof(topo2))
>>> >>>>>>>
>>> >>>>>>>    xtopo2(0:nlon-1) = topo2&lon
>>> >>>>>>>    ytopo2(0:nlon-1) = topo2
>>> >>>>>>>    xtopo2(nlon)     = topo2&lon(nlon-1)
>>> >>>>>>>    ytopo2(nlon)     = ymax2         ; Use actual Y max of contour
>>> >>>>>>> plot
>>> >>>>>>>    xtopo2(nlon+1)   = topo2&lon(0)
>>> >>>>>>>    ytopo2(nlon+1)   = ymax2
>>> >>>>>>>    xtopo2(nlon+2)   = topo2&lon(0)   ; This last point closes
>>> >>>>>>>    ytopo2(nlon+2)   = topo2(0)       ; the polygon.
>>> >>>>>>>
>>> >>>>>>>
>>> >>>>>>> ;---Add the polygon to the contour plot.
>>> >>>>>>>    gnres = True
>>> >>>>>>>    gnres at gsFillColor = "gray20"
>>> >>>>>>>
>>> >>>>>>>    id = gsn_add_polygon(wks, contour_plot_1, xtopo, ytopo, gnres)
>>> >>>>>>>    id2 = gsn_add_polygon(wks, contour_plot_2, xtopo2, ytopo2,
>>> gnres)
>>> >>>>>>>    overlay(contour_plot_1, vector_plot_1)
>>> >>>>>>>    overlay(contour_plot_2, vector_plot_2)
>>> >>>>>>>
>>> >>>>>>>
>>> >>>>>>>      draw(contour_plot_1)     ; This draws everything
>>> >>>>>>>      draw(contour_plot_2)     ; This draws everything
>>> >>>>>>>
>>> >>>>>>>    plot(0) = contour_plot_1
>>> >>>>>>>
>>> >>>>>>>    plot(1) = contour_plot_2
>>> >>>>>>>
>>> >>>>>>>     gsn_panel(wks,plot, (/2, 1/), res)               ; now draw
>>> as
>>> >>>>>>> one plot
>>> >>>>>>>
>>> >>>>>>>
>>> >>>>>>> end
>>> >>>>>>>
>>> >>>>>>>
>>> >>>>>>> 2017-04-24 21:58 GMT+01:00 Mary Haley <haley at ucar.edu>:
>>> >>>>>>>
>>> >>>>>>>> Hi,
>>> >>>>>>>>
>>> >>>>>>>> In the future, please email ncl-talk with follow-up questions so
>>> >>>>>>>> everybody can benefit from the answers.
>>> >>>>>>>>
>>> >>>>>>>> What is the exact error message you're getting?
>>> >>>>>>>>
>>> >>>>>>>> --Mary
>>> >>>>>>>>
>>> >>>>>>>>
>>> >>>>>>>> On Mon, Apr 24, 2017 at 4:08 AM, Appo derbetini <
>>> >>>>>>>> appopson4 at gmail.com> wrote:
>>> >>>>>>>>
>>> >>>>>>>>> Hi Mary,
>>> >>>>>>>>>
>>> >>>>>>>>> Since you send me the script that plot masked fields, I am
>>> working
>>> >>>>>>>>> to do it for many latitude in a panel plot.
>>> >>>>>>>>> Always i got many errors saying that i cannot overlay.
>>> >>>>>>>>> How to solve it?
>>> >>>>>>>>> Thank you in advance.
>>> >>>>>>>>>
>>> >>>>>>>>> Appo
>>> >>>>>>>>>
>>> >>>>>>>>> 2017-04-10 8:23 GMT+01:00 Appo derbetini <appopson4 at gmail.com
>>> >:
>>> >>>>>>>>>
>>> >>>>>>>>>> Thank you very much.
>>> >>>>>>>>>> It's look like what I'm trying to do.
>>> >>>>>>>>>> regards
>>> >>>>>>>>>>
>>> >>>>>>>>>> 2017-04-09 22:11 GMT+01:00 Mary Haley <haley at ucar.edu>:
>>> >>>>>>>>>>
>>> >>>>>>>>>>> I'm sorry I didn't get back to this sooner, but I was having
>>> >>>>>>>>>>> trouble trying to understand what you are doing in the
>>> script.
>>> >>>>>>>>>>>
>>> >>>>>>>>>>> First, I'm not sure it makes sense to take an average of
>>> >>>>>>>>>>> topographical data across a set of latitudes and plot that
>>> as a topo layer.
>>> >>>>>>>>>>>
>>> >>>>>>>>>>> But, that aside, I think you just want to plot the topo line
>>> as
>>> >>>>>>>>>>> a filled polygon instead of a separate XY plot.
>>> >>>>>>>>>>>
>>> >>>>>>>>>>> I'm not sure what I've attached is correct, but hopefully it
>>> >>>>>>>>>>> gives you an idea of what to do.
>>> >>>>>>>>>>>
>>> >>>>>>>>>>> Note that I'm using "sellat" to select a single latitude
>>> value,
>>> >>>>>>>>>>> rather than averaging across a range of latitudes.
>>> >>>>>>>>>>>
>>> >>>>>>>>>>> --Mary
>>> >>>>>>>>>>>
>>> >>>>>>>>>>>
>>> >>>>>>>>>>> On Sun, Apr 9, 2017 at 2:10 AM, Appo derbetini <
>>> >>>>>>>>>>> appopson4 at gmail.com> wrote:
>>> >>>>>>>>>>>
>>> >>>>>>>>>>>> Dear Mary,
>>> >>>>>>>>>>>>
>>> >>>>>>>>>>>> Despite applying what you suggest the problem remain the
>>> same
>>> >>>>>>>>>>>> Thank you very much
>>> >>>>>>>>>>>>
>>> >>>>>>>>>>>> 2017-04-06 8:53 GMT+01:00 Appo derbetini <
>>> appopson4 at gmail.com>:
>>> >>>>>>>>>>>>
>>> >>>>>>>>>>>>> Dear Mary,
>>> >>>>>>>>>>>>>
>>> >>>>>>>>>>>>> Datasets are uploaded as indicated.
>>> >>>>>>>>>>>>>
>>> >>>>>>>>>>>>> What I really trying to do is something like Mask Example
>>> 14.
>>> >>>>>>>>>>>>>
>>> >>>>>>>>>>>>> But I don't want to interpolate topography to grid of wind.
>>> >>>>>>>>>>>>>
>>> >>>>>>>>>>>>> Here attached figure produced by the script
>>> >>>>>>>>>>>>>
>>> >>>>>>>>>>>>> Thank you
>>> >>>>>>>>>>>>>
>>> >>>>>>>>>>>>>
>>> >>>>>>>>>>>>> 2017-04-05 15:36 GMT+01:00 Mary Haley <haley at ucar.edu>:
>>> >>>>>>>>>>>>>
>>> >>>>>>>>>>>>>> I don't think I've ever seen this error before.
>>> >>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>> It's hard to debug this particular script without having
>>> the
>>> >>>>>>>>>>>>>> data so I can run it. But, I'm wondering if this line
>>> might be part of the
>>> >>>>>>>>>>>>>> problem:
>>> >>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>>   plot  = gsn_csm_pres_hgt_vector(wks, r00(0, :,
>>> >>>>>>>>>>>>>> {lonmin:lonmax}), udiv00(0, :, {lonmin:lonmax}),
>>> -200.0*omega00(0, :,
>>> >>>>>>>>>>>>>> {lonmin:lonmax}), res )
>>> >>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>> When you do an arithmetic operation on an array, like:
>>> >>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>>   -200. * omega00
>>> >>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>> you end up stripping all of the metadata off omega00.
>>> What's
>>> >>>>>>>>>>>>>> getting passed into gsn_csm_pres_hgt has no metadata.
>>> This might be okay,
>>> >>>>>>>>>>>>>> however, because I think this routine might just use the
>>> metadata from
>>> >>>>>>>>>>>>>> "udiv00" instead.
>>> >>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>> Still, it's worth trying this:
>>> >>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>>  omega00_scale = omega00(0, :, {lonmin:lonmax})   ;
>>> subsets
>>> >>>>>>>>>>>>>> omega00 *and* copies all metadata
>>> >>>>>>>>>>>>>>  omega00_scale = -200.0 * omega00_scale
>>> >>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>>   plot  = gsn_csm_pres_hgt_vector(wks, r00(0, :,
>>> >>>>>>>>>>>>>> {lonmin:lonmax}), udiv00(0, :, {lonmin:lonmax}),
>>> omega00_scale, res )
>>> >>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>> If you continue to have problems with this script, could
>>> you
>>> >>>>>>>>>>>>>> upload your data (if it's not too large) to our ftp:
>>> >>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>> ftp ftp.cgd.ucar.edu
>>> >>>>>>>>>>>>>> anonymous
>>> >>>>>>>>>>>>>> <use your email address for the password>
>>> >>>>>>>>>>>>>> cd incoming
>>> >>>>>>>>>>>>>> put uvwr.mean.JJAS.nc
>>> >>>>>>>>>>>>>> put atanas_topo.nc
>>> >>>>>>>>>>>>>> quit
>>> >>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>> Thanks,
>>> >>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>> --Mary
>>> >>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>> On Tue, Apr 4, 2017 at 9:18 AM, Appo derbetini <
>>> >>>>>>>>>>>>>> appopson4 at gmail.com> wrote:
>>> >>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>>> Dear all,
>>> >>>>>>>>>>>>>>> I am trying to overlay topography profil on a profil of
>>> wind
>>> >>>>>>>>>>>>>>> and relative humidity.
>>> >>>>>>>>>>>>>>> I don't want to interpolate topography to show some
>>> details.
>>> >>>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>>> Atached are script and datasets used for the plot.
>>> >>>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>>> But I having errors.
>>> >>>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>>> warning:LlDataPolygon: point 6.150009,0.000000 outside
>>> data
>>> >>>>>>>>>>>>>>> domain
>>> >>>>>>>>>>>>>>> warning:LlDataPolygon: point 19.983337,0.000000 outside
>>> data
>>> >>>>>>>>>>>>>>> domain
>>> >>>>>>>>>>>>>>> warning:LlDataLineTo: point 6.150009,0.000000 outside
>>> data
>>> >>>>>>>>>>>>>>> domain
>>> >>>>>>>>>>>>>>> warning:LlDataLineTo: point 19.983337,0.000000 outside
>>> data
>>> >>>>>>>>>>>>>>> domain
>>> >>>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>>> Any help will be appreciated.
>>> >>>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>>> Best regards.
>>> >>>>>>>>>>>>>>>
>>> >>>>>>>>>>>>>>> _______________________________________________
>>> >>>>>>>>>>>>>>> ncl-talk mailing list
>>> >>>>>>>>>>>>>>> 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/20170
>>> 512/89245831/attachment-0001.html
>>>
>>> ------------------------------
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>> End of ncl-talk Digest, Vol 162, Issue 26
>>> *****************************************
>>>
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
>
> --
> Adam Phillips
> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
> www.cgd.ucar.edu/staff/asphilli/   303-497-1726 <(303)%20497-1726>
>
> <http://www.cgd.ucar.edu/staff/asphilli>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170516/8dc34afa/attachment.html 


More information about the ncl-talk mailing list