*   initialization:
    "reinit"
    "set display color white"
    "c"
*   Open file:
    "open Cor_dr_clim_spe_mod.ctl"
    "set gxout grfill"; "set frame on"; "set mpdraw on"; "set mpdset hires"
    "set ylab on"; "set xlab on"; "set grid on"; "set grads off"

    "set lon 102 112"
    "set lat 8 24"
    Nvar=24
    Nmodels=3
    imod=0
    While (imod<Nmodels)
      imod=imod+1
      ivar=0
      While (ivar<Nvar)
        ivar=ivar+1
        Co_clor()
        d_stn()
        if (ivar<10);  "d spe0"ivar; endif
        if (ivar>=10); "d spe" ivar; endif
*       Sub_Reg2
       Main_tit(4.2, 10.5, "Correlation between drought index and climate index ("imod")")
*       "draw title Correlation between drought index and climate index ("imod")"
*       x1=7.0; y1=1; y2=10
*       "run xcbar "x1" "x1+0.2" "y1" "y2" -fs 2 -xo 0.1 -fw 0.08 -fh 0.1 -fx -0.0 -line on -edge triangle ";
       pull dd
       "c"
      Endwhile
    Endwhile
    
*******************************************************
function Main_tit(x, y, title)
  "set string 1 c 5.5 0"; "set strsiz 0.12 0.13"
  "draw string "x" "y" "title
  "set string 1 c 5.5 0"
return

********************************************
function d_stn()
    "set gxout stnmark"
    "set digsiz 0.3"
    "set cmark 3"
    "set cthick 1"
*    "d sensl.1*120"
    "set grads off"
return
*******************************************************
    function Co_clor()
      levels = " -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 -0.05 0 0.05 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 "
      "run color.gs -gxout grfill -kind darkgreen->green-(0)->orange->red -levs" levels
      "set grads off"
    return

*******************************************************
    function P_clor()
      levels = " -300 -200 -100 -50 -20 -5 0 5 20 50 100 200 300 "
      "run color.gs -gxout grfill -kind red->orange->white-(0)->white->lightblue->blue -levs" levels
      "set grads off"
    return

*******************************************************
    function T_clor()
      levels = " -4 -3 -2 -1 -0.5 0 0.5 1 2 3 4 "
      "run color.gs -gxout grfill -kind darkgreen->green->white-(0)->white->orange->red -levs" levels
      "set grads off"
    return
*******************************************************