*ifile='enspctl_q90_TC_freq_201811'
*lead.1='leadtime1'
*lead.2='leadtime2'
*lead.3='leadtime3'
*lead.4='leadtime4'
dir='/work/users/hapt/for_Jack_sfc/control/TC_prob/leadtime'
odir='/work/users/hapt/for_Jack_sfc/control/TC_prob/image'
*nlead=4
nlead=1
ilead=1
nfile=3
file.1='percentage_enssum_prob_0tc'
file.2='percentage_enssum_prob_1tc_2tc'
file.3='percentage_enssum_prob_ge_3tc'

while (ilead <= nlead)
      ifile=1
      while(ifile <= nfile)
            'reinit'
            'sdfopen ' dir ilead '/' file.ifile '_lt' ilead '.nc'
            'set lon 100 140'
            'set lat 0 30'
            'set xlint 5'
            'set xlopts 1 5 0.15'
            'set ylopts 1 5 0.15'
            'set ylint 5'
*            'set parea 1 10.5 0.4 8.5'
            'set display color white'
            'set grads off'
            'c'
            'set gxout grfill'
*            'color 0 100 5 -kind white->cyan->pink->red->yellow->green->lightblue->blue->gray'
*            'color 5 30 2.5 -kind white->cyan->pink->red->green->yellow'
            'set gxout shaded'
            'set mpdset hires'
            'd tc_freq'
            'cbar'
            'draw title SEP-2018_' file.ifile '_lt' ilead
            'printim ' odir '/'file.ifile '_lt' ilead '.png'   
*            'c'
            ifile=ifile+1
      endwhile
      ilead=ilead+1
endwhile