begin
;;;; Define Colors 

  greencol = (/(/187,255,187/),(/134,255,134/), \
               (/100,250,100/),(/ 50,221, 50/),(/0,187,0/), \
               (/0  ,134,  0/),(/0,80,0/) /)/255.0

  x = (/"red","yellow" /)
  opt                  = True
  opt@NumColorsInRange = (/ 11/)
  yr = span_named_colors(x,opt)
  yr1 = yr(2::,:)
  dim = dimsizes(yr1)
  yred = new((/dim(0)+1,dim(1)/), typeof(yr1))
  yred(0:dim(0)-1,:) = yr1
    
  yred(dim(0),0) = 1.0
  yred(dim(0),1) = 1.0
  yred(dim(0),2) = 0.5

  x = (/"greenyellow","darkgreen" /)
  opt                  = True
  opt@NumColorsInRange = (/ 12/)
  yg = span_named_colors(x,opt)
  dim = dimsizes(yg)
  xdim = dim(0) - 1
  ygreen = yg(2:xdim,:)
  ygreen(0,:) = 1.0
  dim = dimsizes(yred)
  y1d = dim(0)
  dim = dimsizes(ygreen)
  y3d =  dim(0)
  y = y1d+y3d
  coloranomaly = new((/y,dim(1)/),typeof(yred))
  coloranomaly(0:y1d-1,:) = yred
  coloranomaly(y1d:(y1d+y3d-1),:) = ygreen
  coloranomaly(13:16,:) = greencol(0:3,:)
  coloranomaly(y-1,:)   = greencol(6,:)
  delete(greencol)
  greenanomaly = coloranomaly(12:y-1,:) 
  redanomaly = coloranomaly(0:10,:)
  redanomaly = redanomaly(::-1,:)
  delete([/yr,yr1,y3d,dim,ygreen,yred,x,opt,y,y1d/])

;;;;;;;;;;;;;;;;;;;;;
;; precip above normal color

  print(coloranomaly)

   opt  = True
   opt@fout  = "colAnomalyPrecWRF.txt"
   fmtx   = "3(f7.5,1x) "
   write_matrix(coloranomaly,fmtx,opt)
   
  
  end