;*******************************************************************

 ; procedure to draw the Taylor Diagram Metrics Table

 ;

 ; AWWG METRICS 0.3 ;;;

 ; 27 Jun 06   ASP

 ;  5 Jul 06   DJS 

 ;

 ; Arguments:

 ;     mfname     : name of the output table

 ;     varNames   : variable names [metrics]

 ;     cases      : case (model) names

 ;     seasons    : season names

 ;     values     : array containing the values to be plotted

 ;     opt        : used to pass optional arguments

 procedure taylor_metrics_table(mfname[1]:string        \       ; plot name

                               ,varNames[*]:string      \       ; variables

                               ,cases[*]:string         \

                               ,seasons[*]:string       \

                               ,values[*][*][*]:numeric  \

                               ,opt:logical              )

 ;;;

 ;;; Output of metrics to a table [grid]

 begin  

     if (.not.isatt(opt,"pltType") .or. \

        (opt@pltType.eq."eps" .or. opt@pltType.eq."png" .or. opt@pltType.eq."gif")) then

         wks = gsn_open_wks("eps",mfname)    ; default 

     else

         wks = gsn_open_wks(opt@pltType,mfname)   

     end if

     xbox0  = (/0.0,1.0,1.0,0.0,0.0/) ; Box template.

     ybox0  = (/0.0,0.0,1.0,1.0,0.0/)

     nvar   = dimsizes(varNames)

     ncases = dimsizes(cases)

     nseas  = dimsizes(seasons)

     dimt   = dimsizes(values)        ; table dimensions

     if (ncases.ne.dimt(0)) then

        print("taylor_metrics_table fatal error: Number of case names ("+ncases+ \

              ") does not match the number of cases passed in ("+dimt(0)+")")

        exit

     end if

     if (nseas.ne.dimt(1)) then

        print("taylor_metrics_table fatal error: Number of season names ("+nseas+ \

              ") does not match the number of seasons passed in ("+dimt(1)+")")

        exit

     end if

     if (dimsizes(varNames).ne.dimt(2)) then

        print("taylor_metrics_table fatal error: Number of metric variables names ("+dimsizes(varNames)+ \

              ") does not match the number of metric variables passed in ("+dimt(2)+")")

        exit

     end if

     ncolm1 = dimt(0)*dimt(1)

     if (ncolm1.le.12) then

        if (ncolm1.le.6) then

           tt_width = .4

        end if

        if (ncolm1.ge.7.and.ncolm1.le.12) then

           tt_width = .2

        end if

        tt_height     = 0.1

        tt_theight    = 0.02

        mn_height     = 0.03  

        mn_theight    = 0.012   

     else

        if (ncolm1.le.20) then

           tt_width   = .15

           tt_height  = 0.06

           tt_theight = 0.012

           mn_height  = 0.018  

           mn_theight = 0.0072

        else

           tt_width   = .1

           tt_height  = 0.05

           tt_theight = 0.01

           mn_height  = 0.015  

           mn_theight = 0.006

        end if

     end if

     if (ncolm1.gt.26.or.dimt(2).gt.63) then

        print("Warning: Recommended maximium size of array passed into taylor_metrics_table is "+ \

              "26 columns (#cases * #seasons) by 63 rows (#variables), continuing")

     end if

     mn_width   = tt_width

     mv_height  = mn_height

     mv_width   = 0.1

     mv_theight = mn_theight

     tt_pRes = True

     tt_pRes@gsFillColor   = "CornflowerBlue"      ; background color for + values

     tt_tRes  = True

     tt_tRes@txFontHeightF = tt_theight

     mn_pRes = True

     mn_pRes@gsFillColor   = "Gray70"              ; background color for variable names

     mn_tRes = True

     mn_tRes@txFontHeightF = mn_theight

     mv_pRes = True

     mv_pRes@gsFillColor   = "White"               ; background color for reference values

     mv_tRes = True

     mv_tRes@txFontHeightF = mn_theight     

 ; Title

     xbox   = tt_width*xbox0

     ybox   = 1.-tt_height*ybox0

     ixtbox = xbox(0)+0.5*(xbox(1)-xbox(0))

     iytbox = ybox(0)+0.5*(ybox(2)-ybox(0))

     gsn_polygon_ndc(wks,xbox,ybox,tt_pRes)

     if (isatt(opt,"tableTitle") ) then

         gsn_text_ndc(wks, opt@tableTitle ,ixtbox,iytbox, tt_tRes)

     else

         gsn_text_ndc(wks,"CAM METRICS",ixtbox,iytbox, tt_tRes)

     end if

     gsn_polyline_ndc(wks,xbox,ybox,False)

     do im = 0,nvar-1

       ybox   = min(ybox)-(ybox0*mn_height)

       ixtbox = xbox(0)+0.5*(xbox(1)-xbox(0))

       iytbox = ybox(0)+0.5*(ybox(2)-ybox(0))

       gsn_polygon_ndc(wks,xbox,ybox,mn_pRes)

       gsn_polyline_ndc(wks,xbox,ybox,False)

       gsn_text_ndc(wks,varNames(im),ixtbox,iytbox, mn_tRes)

     end do

     do icase = 0, ncases-1

       ybox   = 1.-0.5*tt_height*ybox0

       xbox   = ((1.-tt_width)/ncases)*(xbox0+icase)+tt_width

       xboxi  = min(xbox)

       ixtbox = xbox(0)+0.5*(xbox(1)-xbox(0))

       iytbox = ybox(0)+0.5*(ybox(2)-ybox(0))

       gsn_polygon_ndc(wks,xbox,ybox,mn_pRes)

       gsn_polyline_ndc(wks,xbox,ybox,False)

       gsn_text_ndc(wks,cases(icase),ixtbox,iytbox, mn_tRes)

       do iseas = 0, nseas-1

       ybox     = 1.-0.5*tt_height-0.5*tt_height*ybox0

       xbox     = xboxi +  ((1.-tt_width)/(ncases*nseas))*(xbox0+iseas) 

       ixtbox   = xbox(0)+0.5*(xbox(1)-xbox(0))

       iytbox   = ybox(0)+0.5*(ybox(2)-ybox(0))

       gsn_polygon_ndc(wks,xbox,ybox,mn_pRes)

       gsn_polyline_ndc(wks,xbox,ybox,False)

       gsn_text_ndc(wks,seasons(iseas),ixtbox,iytbox, mn_tRes) 

       do im=0,nvar-1

           if (im .eq. 0) then

             ybox = 1.-tt_height-mv_height*ybox0

           else

             ybox = ybox - mv_height

           end if

           iytbox = ybox(0)+0.5*(ybox(2)-ybox(0))

           mv_pRes@gsFillColor = "White" 

           if (icase .gt. 0) then

             if (ismissing(values(icase,iseas,im)) .or. \    ; ???

                 ismissing(values(icase,iseas,im)))then 

     ;          print ("Missing values skipped")

             else

               if (values(icase,iseas,im).le.values(0,iseas,im)) then

     		  if (isatt(opt,"color0")) then

                    mv_pRes@gsFillColor = opt@color0

                 else

                    mv_pRes@gsFillColor = "DarkOliveGreen3"

     		  end if

               else

     		  if (isatt(opt,"color1")) then

                    mv_pRes@gsFillColor = opt@color1

     		  else	

     		   mv_pRes@gsFillColor = "IndianRed1"

     		  end if

               end if

             end if

           end if

           gsn_polygon_ndc(wks,xbox,ybox,mv_pRes)

           gsn_polyline_ndc(wks,xbox,ybox,False)

           gsn_text_ndc(wks,sprintf("%4.3f",values(icase,iseas,im)),ixtbox,iytbox, mv_tRes)     

         end do

       end do

     end do

     draw(wks)

   ;; activate the ;; lines if this is made a function

   ;;if (.not.isatt(opt,"gsnFrame") .or. opt@gsnFrame) then

         frame(wks)

                              ; if png or gif then use convert

         if (isatt(opt,"pltType")    .and. \ 

            (opt@pltType.eq."png" .or. opt@pltType.eq."gif")) then

             system("convert -trim +repage "+mfname+".eps "+mfname+"."+opt@pltType)

         end if

   ;;end if

 end