******* Phan loai cap TC cua JMA:
* 2 : Tropical Depression (TD)
* 3 : Tropical Storm (TS)
* 4 : Severe Tropical Storm (STS)
* 5 : Typhoon (TY)
* 6 : Extra-tropical Cyclone (L)
* 7 : Just entering into the responsible area of Japan Meteorological Agency (JMA)
* 8 : Not used
* 9 : Tropical Cyclone of TS intensity or higher
***********************

   "reinit"
   "set display white"; "c"
   "set mpdraw on"; "set frame on"; "set grid on";"set xlab on";"set ylab on";

    mark.2=0; mark.3=0; mark.4=0; mark.5=0; mark.6=0; mark.7=0; mark.9=0
*    colo.2=4; colo.3=10; colo.4=12; colo.5=2; colo.6=1; colo.7=12; colo.9=12
    colo.2=5; colo.3=7; colo.4=8; colo.5=6; colo.6=0; colo.7=12; colo.9=12
    thic.2=4; thic.3=4; thic.4=4; thic.5=4; thic.6=4; thic.7=4; thic.9=4
    styl.2=1; styl.3=1; styl.4=1; styl.5=1; styl.6=1; styl.7=1; styl.9=1;
********************************************************
    MON.1="Jan"; MON.2="Feb"; MON.3="Mar"; MON.4="Apr"; MON.5="May"; MON.6="Jun";
    MON.7="Jul"; MON.8="Aug"; MON.9="Sep"; MON.10="Oct"; MON.11="Nov"; MON.12="Dec";
********************************************************
*  Years to be plotted: Case=3
    yr0=1997; yr1=1997
    Nyear=yr1-yr0+1
*  ENSO years: Case=1 for El Nino, Case=2 for La Nina
    NEL=15; NLA=18
    E.1.1=1951; E.1.2=1957; E.1.3=1963; E.1.4=1965; E.1.5=1972; E.1.6=1977; E.1.7=1982; E.1.8=1987;
    E.1.9=1991; E.1.10=1994; E.1.11=1997; E.1.12=2002; E.1.13=2004; E.1.14=2009; E.1.15=2015

    E.2.1=1954; E.2.2=1955; E.2.3=1956; E.2.4=1964; E.2.5=1970; E.2.6=1971; E.2.7=1973; E.2.8=1974;
    E.2.9=1975; E.2.10=1984; E.2.11=1985; E.2.12=1988; E.2.13=1998; E.2.14=1999; E.2.15=2000;
    E.2.16=2007; E.2.17=2010; E.2.18=2011

    Tit.1="Tropical cyclone tracks during El Nino events"
    Tit.2="Tropical cyclone tracks during La Nina events"
    Tit.3="Tropical cyclone tracks during "yr0"-"yr1

*  ------ Chose Case for plotting: -------
    Case=3
    if (Case=1); Ncase=NEL; ica=1; Fout="EL"; endif
    if (Case=2); Ncase=NLA; ica=2; Fout="LA"; endif
    if (Case=3); Ncase=Nyear; Fout=yy0"-"yy1; endif
********************************************************
    OutDir="/home/tanpv/Progs/TCs/TC_Statics/PDF/"

    NInt_tc=4
    inten.1="2"; inten.2="3"; inten.3="4"; inten.4="5"; inten.5="6"; inten.6="7"; inten.7="9"
    cdn.1="(TD)"; cdn.2="(TS)"; cdn.3="(STS)"; cdn.4="(TY)"

    Nrow=3; Ncol=4
    x00=0.5; x01=10.5; y00=1.0; y01=8.0
    dx=(x01-x00)/Ncol; dy=(y01-y00)/Nrow
    ddx=0.1; ddy=-0.1

   "sdfopen 4TC_paper.nc"
   lat00=0; lat01=35; lon00=100; lon01=145; dlat=lat01-lat00; dlon=lon01-lon00
   "set lon "lon00" "lon01; "set lat "lat00" "lat01
*   "set map auto"; "draw map"
   "set ccolor 0"; "set clab off"; "set xlab off"; "set ylab off"; "set grads off"
   "set mpdset mres"
   irow=0; imon=0
   While (irow<Nrow)
     irow=irow+1
     if (irow=1); y2=y01; y1=y2-dy; y11=y1; endif
     if (irow>1); y2=y11+ddy; y1=y2-dy; y11=y1; endif
     icol=0
     While (icol<Ncol)
       icol=icol+1
       if (icol=1); x1=x00; x2=x1+dx; x11=x2; endif
       if (icol>1); x1=x11+ddx; x2=x1+dx; x11=x2; endif
       "set parea "x1" "x2" "y1" "y2
       "set xlab off"; if (irow=Nrow); "set xlab on"; "set xlint 10"; endif
       "set ylab off"; if (icol=1); "set ylab on"; "set ylint 10"; endif
       
       "set ccolor 0"; "d tc"
       HS_TS
*       "Title_lonlat -lon "lon " -lat "lat " -str "strsize " -ang "angle " - tit String"
       imon=imon+1
       draw_tit(x1+0.3, y2-0.2,MON.imon)

       if (imon<10); filename="/home/tanpv/Progs/TCs/TC_Dat/TC_tr_for_GrADS_0"imon".txt"; endif
       if (imon>=10); filename="/home/tanpv/Progs/TCs/TC_Dat/TC_tr_for_GrADS_"imon".txt"; endif

       code=0
       res=read (filename)
       code=sublin(res,1)
       Ntc=sublin(res,2)
       While (code=0)
         res=read (filename)
         code=sublin(res,1)
*         if (code=2); "draw map"; "gxprint "OutDir"Tracks_"Fout"_"imon".pdf pdf"; "close 1"; break; endif
*         if (code=2); "q files"; "close 2"; break; endif
         if (code=2); break; endif
         if (code=9); say "something wrong !! "code; "close 2"; break; endif
         xx=sublin(res,2)
         Nlf=subwrd(xx,1)
         res=read (filename)
         code=sublin(res,1)
         xx = sublin(res,2)
         id = subwrd(xx,1); yr = subwrd(xx,2); mm = subwrd(xx,3); dd = subwrd(xx,4); hh = subwrd(xx,5);
         lon1 = subwrd(xx,6); lat1 = subwrd(xx,7); CD1 = subwrd(xx,8)
         if (CD1=9); CD1=3; endif
         lon1 = lon1*0.10; lat1 = lat1*0.10
         "q w2xy "lon1" "lat1
         xx1 = subwrd(result,3)
         yy1 = subwrd(result,6)
         ilf=1
         While (ilf<Nlf)
            ilf=ilf+1
            res=read (filename)
            code=sublin(res,1)
            xx = sublin(res,2)
            id = subwrd(xx,1)
            yr = subwrd(xx,2); mm = subwrd(xx,3); dd = subwrd(xx,4); hh = subwrd(xx,5);
            lon2 = subwrd(xx,6); lat2 = subwrd(xx,7)
            lon2 = lon2/10; lat2 = lat2/10; CD2 = subwrd(xx,8); if (CD2=9); CD2=3; endif
            "q w2xy "lon2" "lat2
            xx2 = subwrd(result,3)
            yy2 = subwrd(result,6)

            yyyy=0
            While (yyyy<Ncase)
              yyyy=yyyy+1;
              if (Case=3);
                 if (yr>=yr0 & yr<=yr1);
                    "set line "colo.CD2" 1 4"
                    if (xx1>x1 & xx1<x2 & xx2>x1 & xx2<x2 & yy1>y1 & yy1<y2 & yy2<y2);
                       "draw line "xx1" "yy1" "xx2" "yy2;
                    endif
                    xx1 = xx2; yy1 = yy2;
                 endif
                 yyyy=Ncase
              endif

              if (Case=1 | Case=2);
                 if (yr=E.ica.yyyy);
*                    say yyyy" "yr" "E.ica.yyyy" "CD2
                    "set line "colo.CD2" 1 4"
                    if (xx1>x1 & xx1<x2 & xx2>x1 & xx2<x2 & yy1>y1 & yy1<y2 & yy2<y2);
                       "draw line "xx1" "yy1" "xx2" "yy2;
                    endif
                    xx1 = xx2; yy1 = yy2;
                 endif
              endif
            Endwhile

*            if (code=2); "q files"; "close 2"; break; endif
            if (code=2); break; endif
            if (code=9); say "something wrong !! "code; "q files"; "close 2"; break; endif
         Endwhile
*         "q files"; 
*         File=sublin(result,2); say "File=" File
*         pull dd
         "draw map";
       Endwhile
       rc = close(filename)
     Endwhile
   Endwhile
   icd=0
   While (icd<NInt_tc)
      icd=icd+1; x0=x00+(icd-1)*dx+0.3; xx0=x0+0.2
      CDtc=inten.icd
      draw_tit(x0+0.9, y00-0.65, inten.icd " "cdn.icd)
      "set line "colo.CDtc" 1 4"
      "draw line "x0" "y00-0.75" "x0+dx-0.7" "y00-0.75
      draw_tit((x00+x01)/2, y00-0.55,"TC intensity")
   Endwhile
   Main_title ((x00+x01)/2, y01+0.2, Tit.Case)
   "gxprint " OutDir "Tracks_" Fout ".pdf pdf"
   "close 1"
   pull dd
*   "draw map";"gxprint "OutDir"Tracks_"Fout"_"imon".pdf pdf"
*   say "next mon: "imon
*   "c"
* Endwhile
* say "Done !!!"
*   pull dd
*  ----------------------------------------------------
*******************************************************
   function Main_title (x, y, title)
     "set string 1 c 5.5 0"
     "set strsiz 0.11 0.13"
     "draw string "x" "y" "title
   return
*******************************************************

   function draw_tit(x, y, title)
      "set string 1 c 5.5"; "set strsiz 0.12 0.12"
      "draw string "x" "y" "title
   return
*******************************************************

   function draw_hsts(x, y, title)
      "set string 1 c 4.0 45"; "set strsiz 0.10 0.10"
      "draw string "x" "y" "title
   return
*******************************************************