[ncl-talk] Plotting spatial map of SPI

Dennis Shea shea at ucar.edu
Wed Apr 19 09:25:39 MDT 2017


Please do not use ncl-talk digest to post questions. Other, non related
questions are appended to your email.

---
+Replace

  yyyymm = f->time
  yyyymm = yyyymm/100

  year    = yyyymm/100
  yrStrt  = 1901
  yrLast  = 2014
  nyear   = yrLast-yrStrt+1
  yyyymm  = yyyymm_time(yrStrt, yrLast, "integer")
  yrfrac  = (/ yyyymm_to_yyyyfrac(yyyymm, 0.0) /)

++With

  yyyymm = cd_calendar(prc&time, -1)
  printVarSummary(yyyymm)

  prc&time = (/ yyyymm /)
  printVarSummary(prc&time)
  yrfrac   = (/ yyyymm_to_yyyyfrac(yyyymm, 0.0) /)  ; (/.../) no meta data

Please see:
   http://www.ncl.ucar.edu/Document/Functions/date.shtml

Specifically:
   http://www.ncl.ucar.edu/Document/Functions/Built-in/cd_calendar.shtml

http://www.ncl.ucar.edu/Document/Functions/Contributed/yyyymm_to_yyyyfrac.shtml

====
By default, will plot the entire grid. You have only a small area. You must
limit the region.


++Add
  latS = min(prc&lat)
  latN = max(prc&lat)
  lonL = min(prc&lon)
  lonR = max(prc&lon)

  res at mpMinLatF = latS
  res at mpMaxLatF = latN
  res at mpMinLonF = lonL
  res at mpMaxLonF = lonR

===
++ Change
  res at gsnCenterString  = yyyymm({192112})

++ To
  res at gsnCenterString  = spi&time({192112})

====
LOOK at:

Variable: spi
Number of Dimensions: 4
Dimensions and sizes:   [run | 1] x [time | 1368] x [lat | 4] x [lon | 9]
Coordinates:
            run: [12..12]
            time: [190101..201412]
            lat: [3.75..5.25]
            lon: [36.75..40.75]

You use:

    spi({192112},:,:,:)

It should be

    spi(0,{192112},:,:)   ; run=0, {time}, lat, lon



On Tue, Apr 18, 2017 at 10:07 AM, Dawit Abebe <dawit.kan at gmail.com> wrote:

>
>
> Dear All,
>
> I was trying to plot a spatial map of SPI but having a challenge indexing
> coordinates. I pasted my script as well as the point a stuck (red text)
> from NCL output below.
>
> Thank you so much for your support as usual.
>
> Dawit
>
> *1. My script*
>
> ;===========================================================
> ==============================
> ;spi spatial map 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"
>
> begin
>
>     diri  = "./"
>     fili  = "pre_cru_1901_2014.nc"
>     f     = addfile(diri+fili, "r")
>
>     prc  = f->pre
>     prc at _FillValue = -999.0
>     pmsg = prc at _FillValue
>     printVarSummary(prc)
>
>     dimp = dimsizes(prc)
>     ntim = dimp(0)
>     nlat = dimp(1)
>     mlon = dimp(2)
>
> run  = (/12/)
> nrun = dimsizes(run)
>
> spi = new((/nrun,ntim,nlat,mlon/), typeof(prc), pmsg)
>
>     opt = True
>
>     do nr=0,nrun-1
>    spi(nr,:,:,:)= dim_spi_n(prc, run(nr), True, 0)
> end do
>
> copy_VarCoords(prc,spi(0,:,:,:))
> spi at long_name = "SPI"
> spi!0         = "run"
> spi&run       =  run
> printVarSummary(spi)
>
>  ;==========================================================
> ==============================
>  ; Plot parameter
>  ;==========================================================
> ==============================
>     dimp = dimsizes(prc)
>     ntim = dimp(0)
>     nlat = dimp(1)
>     mlon = dimp(2)
>
>   yyyymm = f->time
>   yyyymm = yyyymm/100
>
>   year    = yyyymm/100
>   yrStrt  = 1901
>   yrLast  = 2014
>   nyear   = yrLast-yrStrt+1
>   yyyymm  = yyyymm_time(yrStrt, yrLast, "integer")
>   yrfrac  = (/ yyyymm_to_yyyyfrac(yyyymm, 0.0) /)
>
>   prc&time = (/ yyyymm /)
>
> ;===========================================================
> =============================
>
>   wks          = gsn_open_wks ("ps","spi")
>   gsn_define_colormap(wks,"BlueWhiteOrangeRed")
>   gsn_reverse_colormap(wks)                                    ; Reverse
> the color map.
>
>   res                  = True
>   res at gsnDraw          = False                                 ; don't
> draw
>   res at gsnFrame         = False                                 ; don't
> advance frame
>   res at gsnSpreadColors  = True
>
>   res at cnFillOn         = True                                   ; color
> Fill
>   res at cnFillMode       = "RasterFill"                           ; Raster
> Mode
>   res at cnLinesOn        =  False                                 ; Turn
> off contour lines
>   res at cnLevelSelectionMode = "ManualLevels"                     ; set
> manual contour levels
>   res at cnMinLevelValF   = -3.0                                   ; set min
> contour level
>   res at cnMaxLevelValF   =  3.0                                   ; set max
> contour level
>   res at cnLevelSpacingF  =  0.5                                   ; set
> contour spacing
>   res at lbLabelBarOn     = False                                  ; turn
> off individual cb's
>
>   resP                 = True                                   ; panel
> resources
>   resP at txString        = False
>   resP at gsnPanelLabelBar= True                                   ; add
> common colorbar
>
>   plot = new ( 5, "graphic")
>
>   do nl=0,nlat-1
>      do ml=0,mlon-1
>         do nr=0,nrun-1
>            spi(nr,:,:,:)= dim_spi_n(prc, run(nr), True, 0)
>            spi at long_name = "SPI"
>            spi at units = "nrun="+run(nr)
>            copy_VarCoords(prc,spi(0,:,:,:))
>       spi!0         = "run"
>       spi&run       =  run
>       printVarSummary(spi)
>
>            res at gsnCenterString  = yyyymm({192112})
>            plot(0) = gsn_csm_contour_map_ce(wks,spi({192112},:,:,:), res)
>
>            res at gsnCenterString  = yyyymm({192212})
>            plot(1) = gsn_csm_contour_map_ce(wks,spi({192212},:,:,:), res)
>
>            res at gsnCenterString  = yyyymm({192712})
>            plot(2) = gsn_csm_contour_map_ce(wks,spi({192712},:,:,:), res)
>
>            res at gsnCenterString  = yyyymm({197312})
>            plot(3) = gsn_csm_contour_map_ce(wks,spi({197312},:,:,:), res)
>
>            res at gsnCenterString  = yyyymm({198412})
>            plot(4) = gsn_csm_contour_map_ce(wks,spi({198412},:,:,:), res)
>         end do
>       end do
>      end do
>
>      gsn_panel(wks,plot,(/5,1/),resP)          ; now draw as one plot
>
> end
>
> *2. NCL output from exciting the script above*
>
> (my_ncl) Dagemawes-MacBook-Pro:~ dawitabebe$ ncl -Q
> spi12_plotmap_script.ncl
>
> Variable: prc
> Type: double
> Total Size: 393984 bytes
>             49248 values
> Number of Dimensions: 3
> Dimensions and sizes: [time | 1368] x [lat | 4] x [lon | 9]
> Coordinates:
>             time: [ 380..41987]
>             lat: [3.75..5.25]
>             lon: [36.75..40.75]
> Number Of Attributes: 5
>   correlation_decay_distance : 450
>   long_name : precipitation
>   units : mm/month
>   missing_value : 9.96920996839e+36
>   _FillValue : -999
>
> Variable: spi
> Type: double
> Total Size: 393984 bytes
>             49248 values
> Number of Dimensions: 4
> Dimensions and sizes: [run | 1] x [time | 1368] x [lat | 4] x [lon | 9]
> Coordinates:
>             run: [12..12]
>             time: [ 380..41987]
>             lat: [3.75..5.25]
>             lon: [36.75..40.75]
> Number Of Attributes: 2
>   long_name : SPI
>   _FillValue : -999
>
> Variable: spi
> Type: double
> Total Size: 393984 bytes
>             49248 values
> Number of Dimensions: 4
> Dimensions and sizes: [run | 1] x [time | 1368] x [lat | 4] x [lon | 9]
> Coordinates:
>             run: [12..12]
>             time: [190101..201412]
>             lat: [3.75..5.25]
>             lon: [36.75..40.75]
> Number Of Attributes: 4
>   units : nrun=12
>   long_name : SPI
>   run : 12
>   _FillValue : -999
> fatal:NclOneDValGetClosestIndex: finish coordinate index out of range,
> can't continue
> fatal:Could not obtain coordinate indexes, unable to perform subscript
> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 98 in
> file spi12_plotmap_script.ncl
>
> (my_ncl) Dagemawes-MacBook-Pro:~ dawitabebe$
> ……………………………………...
>
>
> On Apr 18, 2017, at 10:29 AM, 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: Problem with gsn_add_text (Stavros Dafis)
>   2. Problem with contour plot (Guilherme Martins)
>   3. Re: Problem with gsn_add_text (Rick Brownrigg)
>   4. Re: Problem with contour plot (Mary Haley)
>   5. how to draw this kind of plot (wen)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Mon, 17 Apr 2017 21:48:31 +0300
> From: Stavros Dafis <sdafis at cc.uoi.gr>
> Subject: Re: [ncl-talk] Problem with gsn_add_text
> To: Rick Brownrigg <brownrig at ucar.edu>
> Cc: "ncl-talk at ucar.edu" <ncl-talk at ucar.edu>
> Message-ID: <1492454911.58f50dffc9203 at webmail.uoi.gr>
> Content-Type: text/plain; charset=ISO-8859-7
>
> Dear Rick, thank you for your suggestion, I have tried to assign both the
> FillValue and the missing value (@missing_value) but nothing seems to work.
>
> I forgot to mention that I am using NCL 6.2.0, 6.3.0 and 6.4.0, on
> different
> machines.
>
>
> --
> Stavros NTAFIS (DAFIS)
> -----------------------------------------------
> Physicist - Meteorologist, M.Sc.
> Research Associate, National Observatory of Athens
> (+33)9 81 94 22 12 <+33%209%2081%2094%2022%2012>, Mobile: +30 697 04 20
> 242
> ---------------
> Weather charts:
> http://www.meteo.gr
> http://www.meteo.gr/meteomaps/
>
>
> Quoting Rick Brownrigg <brownrig at ucar.edu>:
>
> Hi Stavros,
>
> I don't know for certain, but I suspect you need to assign the _FillValue
> to the string array "points" that you constructed from "t2" (I don't see
> how the attributes from t2 would carry to points via the sprintf() ).
>
> Hope that helps...
> Rick
>
> On Sun, Apr 16, 2017 at 12:56 PM, Stavros Dafis <sdafis at cc.uoi.gr> wrote:
>
> Hello everyone, I am facing a problem with the function gsn_add_text. In
> the
> description of the function, it is mentioned that no missing values will be
> plotted:
>
> &#8220;This function attaches a list of text strings to the given plot.
> The text
> strings will be drawn only when the plot is drawn. If a missing value is
> encountered in the text strings or x/y values, then no text string will be
> drawn at this pair.&#8221;
>
> My script reads a netCDF file with 3 variables t2(lat,lon), lat(lat),
> lon(lon),
> so in order to attach values on the plot using this function, I have to
> make
> new 2-D variables for lat/lon and define the string to be plotted on the
> map.
> Values are assigned just fine on the plot, over the land, but missing
> values
> are also printed over the sea. I tried to use landsea masking but my data
> have
> high horizontal resolution (5km). Any suggestions?
>
> I have the same problem with WRF data, with already 2-D lat/lon arrays.
>
> This is the part of the code which prints the values:
> .
> .
> .
> contour = gsn_csm_contour_map(wks,t2,res)
>
> t2 at missing_value = -1
> t2 at _FillValue = -1
> (default missing values = -999)
>
> lat2 = new(dimsizes(t2),typeof(lat),-1)
> lon2 = new(dimsizes(t2),typeof(lon),-1)
>
> nb = dimsizes(t2)
>
> do nl=0,nb(0)-1,1
> do ml=0,nb(1)-1,1
>
> lat2(nl,ml) = lat(nl)
> lon2(nl,ml) = lon(ml)
>
> end do
> end do
>
> points = sprintf("%3.0f", t2)
> ;points = ""+t2(:,:)+""
>
> values = gsn_add_text(wks,contour,points,lon2(:,:),lat2(:,:),True)
>
>
>
> --
> Stavros NTAFIS (DAFIS)
> -----------------------------------------------
> Physicist - Meteorologist, M.Sc.
> Research Associate, National Observatory of Athens
> (+33)9 81 94 22 12 <+33%209%2081%2094%2022%2012>, Mobile: +30 697 04 20
> 242
> ---------------
> Weather charts:
> http://www.meteo.gr
> http://www.meteo.gr/meteomaps/
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
>
> ------------------------------
>
> Message: 2
> Date: Mon, 17 Apr 2017 19:00:02 +0000
> From: Guilherme Martins <jgmsantos at gmail.com>
> Subject: [ncl-talk] Problem with contour plot
> To: NCL <ncl-talk at ucar.edu>
> Message-ID:
> <CAH+t_hMaRUDHM0wMkLC9cz2pFmMYGXgGC+X+bcxDLo=+8fgWQg at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Dear users,
>
> After searching extensively on the NCL website I was not able to find a
> solution to my problem. I only want values within the contour of the
> continent, but some edges appear in the figure that are stemming from the
> mask that I applied and I do not know how to remove this. The figure is
> attached for a better view of my problem.
>
> Best regards,
>
> Guilherme.
> --
> ------------------------------------------------------------
> ------------------------------------
> Instituto Nacional de Pesquisas Espaciais (INPE)
> Centro de Previs?o de Tempo e Estudos Clim?ticos (CPTEC)
> Divis?o de Sat?lites e Sistemas Ambientais (DSA)
> Programa de Monitoramento de Queimadas
> Telefone (INPE/CP): +55 12 3186-9205 <+55%2012%203186-9205> || Celular
> (TIM): +55 12 98111-4292 <+55%2012%2098111-4292>
> E-mail: guilherme.martins at inpe.br || jgmsantos at gmail.com
> Skype: guilherme.martins.
> Homepage: https://sites.google.com/site/jgmsantos
> Curr?culo Lattes: http://lattes.cnpq.br/5997657584785803
> ------------------------------------------------------------
> ------------------------------------
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20170417/c3215b0e/attachment-0001.html
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: figure.png
> Type: image/png
> Size: 58495 bytes
> Desc: not available
> Url : http://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20170417/c3215b0e/attachment-0001.png
>
> ------------------------------
>
> Message: 3
> Date: Mon, 17 Apr 2017 14:32:43 -0600
> From: Rick Brownrigg <brownrig at ucar.edu>
> Subject: Re: [ncl-talk] Problem with gsn_add_text
> To: Stavros Dafis <sdafis at cc.uoi.gr>
> Cc: "ncl-talk at ucar.edu" <ncl-talk at ucar.edu>
> Message-ID:
> <CAGKRhbHX1v-1sJ3G=QrXGA+0wSJ2KKsBB_7GqbCnp9krh9fkkQ at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Hi,
>
> You used -1 as _FillValue for variable "t2". Note it is an integer value --
> what type is the array t2 -- float or int?  In any case, you created the
> formatted string array "points" from t2 as "%3.0f", so that any data values
> in t2 that were -1 are now " -1", and thus it seems your _FillValue should
> be " -1". Did you do this, and do the types between data and _FillValue
> agree everywhere?
>
> If that's not the issue, I think we may need to see script and dataset to
> figure out what the issue is.
>
> HTH...
> Rick
>
> On Mon, Apr 17, 2017 at 12:48 PM, Stavros Dafis <sdafis at cc.uoi.gr> wrote:
>
> Dear Rick, thank you for your suggestion, I have tried to assign both the
> FillValue and the missing value (@missing_value) but nothing seems to work.
>
> I forgot to mention that I am using NCL 6.2.0, 6.3.0 and 6.4.0, on
> different
> machines.
>
>
> --
> Stavros NTAFIS (DAFIS)
> -----------------------------------------------
> Physicist - Meteorologist, M.Sc.
> Research Associate, National Observatory of Athens
> (+33)9 81 94 22 12 <+33%209%2081%2094%2022%2012>, Mobile: +30 697 04 20
> 242
> ---------------
> Weather charts:
> http://www.meteo.gr
> http://www.meteo.gr/meteomaps/
>
>
> Quoting Rick Brownrigg <brownrig at ucar.edu>:
>
> Hi Stavros,
>
> I don't know for certain, but I suspect you need to assign the _FillValue
> to the string array "points" that you constructed from "t2" (I don't see
> how the attributes from t2 would carry to points via the sprintf() ).
>
> Hope that helps...
> Rick
>
> On Sun, Apr 16, 2017 at 12:56 PM, Stavros Dafis <sdafis at cc.uoi.gr>
>
> wrote:
>
>
> Hello everyone, I am facing a problem with the function gsn_add_text.
>
> In
>
> the
> description of the function, it is mentioned that no missing values
>
> will be
>
> plotted:
>
> &#8220;This function attaches a list of text strings to the given plot.
> The text
> strings will be drawn only when the plot is drawn. If a missing value
>
> is
>
> encountered in the text strings or x/y values, then no text string
>
> will be
>
> drawn at this pair.&#8221;
>
> My script reads a netCDF file with 3 variables t2(lat,lon), lat(lat),
> lon(lon),
> so in order to attach values on the plot using this function, I have to
> make
> new 2-D variables for lat/lon and define the string to be plotted on
>
> the
>
> map.
> Values are assigned just fine on the plot, over the land, but missing
> values
> are also printed over the sea. I tried to use landsea masking but my
>
> data
>
> have
> high horizontal resolution (5km). Any suggestions?
>
> I have the same problem with WRF data, with already 2-D lat/lon arrays.
>
> This is the part of the code which prints the values:
> .
> .
> .
> contour = gsn_csm_contour_map(wks,t2,res)
>
> t2 at missing_value = -1
> t2 at _FillValue = -1
> (default missing values = -999)
>
> lat2 = new(dimsizes(t2),typeof(lat),-1)
> lon2 = new(dimsizes(t2),typeof(lon),-1)
>
> nb = dimsizes(t2)
>
> do nl=0,nb(0)-1,1
> do ml=0,nb(1)-1,1
>
> lat2(nl,ml) = lat(nl)
> lon2(nl,ml) = lon(ml)
>
> end do
> end do
>
> points = sprintf("%3.0f", t2)
> ;points = ""+t2(:,:)+""
>
> values = gsn_add_text(wks,contour,points,lon2(:,:),lat2(:,:),True)
>
>
>
> --
> Stavros NTAFIS (DAFIS)
> -----------------------------------------------
> Physicist - Meteorologist, M.Sc.
> Research Associate, National Observatory of Athens
> (+33)9 81 94 22 12 <+33%209%2081%2094%2022%2012>, Mobile: +30 697 04 20
> 242
> ---------------
> Weather charts:
> http://www.meteo.gr
> http://www.meteo.gr/meteomaps/
>
>
>
> _______________________________________________
> 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/20170417/01bd2242/attachment-0001.html
>
> ------------------------------
>
> Message: 4
> Date: Mon, 17 Apr 2017 15:02:08 -0600
> From: Mary Haley <haley at ucar.edu>
> Subject: Re: [ncl-talk] Problem with contour plot
> To: Guilherme Martins <jgmsantos at gmail.com>
> Cc: NCL <ncl-talk at ucar.edu>
> Message-ID:
> <CACNN_CJz1yvpMXMABZ8FhSEov1T0qN63rzsNBBVXBhLos3UKjg at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Guilherme,
>
> Please see examples shapefiles_21.ncl at:
>
> http://www.ncl.ucar.edu/Applications/shapefiles.shtml#ex21
>
> It explains why you see these "blocky" contours in NCL.
>
> Simply put, if you have four points that represent a box at the edge of a
> domain, and three of them have data one them and one is missing, then NCL
> will only draw contours in half of the box, because it won't extrapolate to
> that missing point.
>
> The real question is: do you simply want the graphics not to have gaps, or
> do you actually want there to be data in locations that fall just outside
> the geographical boundary?  You can do things to force values just outside
> the geographical boundary, but this depends on how you are masking in the
> first place.
>
> If you are masking based on a shapefile, then see example shapefiles_18.ncl
> on the same page:
>
> http://www.ncl.ucar.edu/Applications/shapefiles.shtml#ex18
>
> Note, however, that this mask is now returning values outside the range of
> interest.
>
> If you are simply trying to mask your data *graphically* by some
> geographical area, then instead of trying to use a mask function, you can
> draw the full data, and then draw polygons filled in with white in the
> areas you don't want to see graphically.  That's what example 21 is
> showing.
>
> --Mary
>
>
> On Mon, Apr 17, 2017 at 1:00 PM, Guilherme Martins <jgmsantos at gmail.com>
> wrote:
>
> Dear users,
>
> After searching extensively on the NCL website I was not able to find a
> solution to my problem. I only want values within the contour of the
> continent, but some edges appear in the figure that are stemming from the
> mask that I applied and I do not know how to remove this. The figure is
> attached for a better view of my problem.
>
> Best regards,
>
> Guilherme.
> --
> ------------------------------------------------------------
> ------------------------------------
> Instituto Nacional de Pesquisas Espaciais (INPE)
> Centro de Previs?o de Tempo e Estudos Clim?ticos (CPTEC)
> Divis?o de Sat?lites e Sistemas Ambientais (DSA)
> Programa de Monitoramento de Queimadas
> Telefone (INPE/CP): +55 12 3186-9205 <+55%2012%203186-9205>
> <+55%2012%203186-9205> || Celular
> (TIM): +55 12 98111-4292 <+55%2012%2098111-4292> <+55%2012%2098111-4292>
> E-mail: guilherme.martins at inpe.br || jgmsantos at gmail.com
> Skype: guilherme.martins.
> Homepage: https://sites.google.com/site/jgmsantos
> Curr?culo Lattes: http://lattes.cnpq.br/5997657584785803
> ------------------------------------------------------------
> ------------------------------------
>
> _______________________________________________
> 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/20170417/96b2b947/attachment-0001.html
>
> ------------------------------
>
> Message: 5
> Date: Tue, 18 Apr 2017 15:29:04 +0800 (CST)
> From: wen <wenguanhuan at 163.com>
> Subject: [ncl-talk] how to draw this kind of plot
> To: ncl-list <ncl-talk at ucar.edu>
> Message-ID: <5de3ed86.9f7f.15b7ff6fb6f.Coremail.wenguanhuan at 163.com>
> Content-Type: text/plain; charset="gbk"
>
> Hi all,
>  I am trying to use TRMM 3B42 data  to draw this kind of plot. I haven't
> found a example on the website to draw it.  Do anyone know  how to make it.
> Which function should I
> use. Thank you!
>
>
>
>
>
>
>
>
> --
>
> ??
> ???
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20170418/b09be439/attachment.html
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: ??1.png
> Type: image/png
> Size: 249741 bytes
> Desc: not available
> Url : http://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20170418/b09be439/attachment.png
>
> ------------------------------
>
> _______________________________________________
> 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 161, Issue 23
> *****************************************
>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170419/561025ae/attachment.html 


More information about the ncl-talk mailing list