[ncl-talk] Overlay plots

Amadou Coulibaly mpapin24 at gmail.com
Wed May 18 03:32:52 MDT 2016


Hello NCL users,

I want to have an overlay plot aver time and levels. Let say I have:
height_levels = ispan(100, 25900, 500) ; height levels to plot - in meter
nlevels = dimsizes(height_levels)
times = wrf_user_getvar(a,"times",-1)  ; get all times in the file
ntimes = dimsizes(times)    ; number of times in the file

I am working 3-hourly data of on one day, which means that ntimes=8.
I want to plot the contour map of the following parameter:
Rot_st = new((/ntimes, nlevels, nlat, nlon/), float)

Variable: Rot_st
Type: float
Total Size: 20653568 bytes
            5163392 values
Number of Dimensions: 4
Dimensions and sizes:    [8] x [52] x [107] x [116]
Coordinates:
Number Of Attributes: 1
  _FillValue :    9.96921e+36

​The following command lines plot that parameter at each time for each
level. But I want to have a plot averaged over times and levels​.

  do it = 0,ntimes-1             ; TIME LOOP
      print("Working on time: " + times(it) )
      res at TimeLabel = times(it)   ; Set Valid time to use on plots

      do level = 0, nlevels-1            ; LOOP OVER LEVELS
           pressure = height_levels(level)
           ; Add some level info to the plot
           res at PlotLevelID = pressure + "m"

    ; Plotting options for RH
        opts = res
        opts at cnFillOn = True
        opts at pmLabelBarOrthogonalPosF = -0.1
        opts at ContourParameters = (/0.6/)
        opts at cnFillColors = (/"White","White","White", \
                              "White","Chartreuse","Green",\
                              "Green3","Green4", \
                              "ForestGreen","PaleGreen4"/)
        contour_rot =         wrf_contour(a,wks,Rot_st(it,level,:,:),opts)
        delete(opts)

    ; MAKE PLOTS
      plot = wrf_map_overlays(a,wks,(/contour_rot/),pltres,mpres)

    end do      ; END OF LEVEL LOOP

  end do        ; END OF TIME LOOP

​Is there a way to transform the above command lines in such that to have a
plot averaged over times and levels?

Best regards​

-- 
* COULIBALY   AMADOU   *
PhD  Student  on  West  African  Climate  System (WACS)
FUTA - Federal  University  of  Technology of Akure, Nigeria

*Visiting Student - *
*University of Cologne, Germany**Institute of Geophysics and Meteorology*
Pohligstr. 3 / Office 3.102
D-50969 Köln

*Project*: WASCAL (West African Science Service Centre on Climate Change
and Adapted Land Use)
Phone:(+234) 812 965 3659 /(+223) 90 08 26 71 /(+226) 64 57 37 27 /+49
15218352574

E-mail: mpapin24 at gmail.com / coulibalya68 at yahoo.com

"*The time is always right to do right": Nelson Mandela*

*"Character is like a tree and reputation like a shadow.  The shadow is
what we think of it; the tree is the real thing" : Abraham Lincoln*

*"Do what you can, with what you have, where you are" Theodore Roosevelt*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160518/81775956/attachment.html 


More information about the ncl-talk mailing list