"reinit" "set display color white"; "c" 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" Dir_Inp="/home/tanpv/Data/562_VNU/0.25_VN/Dr_Chars/" Dat_out="/home/tanpv/Progs/562_VNU/PDF_Ope/" Yr_Ope=2024 Mn_Ope=01 Nlead=6; Nvar=2 Varnm.1.1="Monthly Rainfall Amount (mm/month)"; Varnm.2.1="Monthly 2-m Temperature (C degree)" Varnm.1.2="Monthly Relative Rainfall Anomaly (%)"; Varnm.2.2="Monthly 2-m Temperature Anomaly (C degree)" V.1="Pr"; V.2="T2" Var.1.1="pr"; Var.1.2="pr_ano" Var.2.1="t2"; V.2.2="t2_ano" Ncol=Nlead; Nrow=2 ;* 5 cases, 4 months (1,4,7,10) y00=8.1; y01=0.3 dy=(y00-y01)/Nrow x00=0.5; x01=11.0 dx=(x01-x00)/Ncol-0.15 "sdfopen "Dir_Inp "Forecast_Ope_Pr_T2_and_Anomaly.nc" ivar=0 While (ivar1); y2=y11-0.1; y1=y2-dy; y11=y1; endif While (icol12); Month=Month-12; Year=Year+1; endif if (icol=1); x1=x00; x2=x1+dx; x11=x2; endif if (icol>1); x1=x11+0.02; x2=x1+dx; x11=x2; endif "set parea "x1" "x2" "y1" "y2 "set mpdraw off"; "set frame on"; "set xlab on"; "set grid off"; "set ylab off";"set xlint 4"; "set ylint 2" "set lon 101 111"; "set lat 7.5 24.0" if (icol=1); "set ylab on"; "set grid off"; endif ;* if (irow=Nrow); "set xlab on"; "set grid off"; endif if (ivar=1 & irow=1) ;* pr Pr_Clor(); "d pr_fcs(t="icol")" ;* "d "Var.ivar.irow"(t="icol")" Con_pr(); "d pr_fcs(t="icol")" ;* "d "Var.ivar.irow"(t="icol")" endif if (ivar=1 & irow=2) ;* pr_ano Pr_ano(); "d pr_ano(t="icol")*100" ;* "d "Var.ivar.irow"(t="icol")" Con_pr_ano(); "d pr_ano(t="icol")*100" ;* "d "Var.ivar.irow"(t="icol")" endif if (ivar=2 & irow=1) ;* t2 t2_Clor(); "d t2_fcs(t="icol")" ;* "d "Var.ivar.irow"(t="icol")" Con_t2(); "d t2_fcs(t="icol")" ;* "d "Var.ivar.irow"(t="icol")" endif if (ivar=2 & irow=2) ;* t2_ano t2_ano(); "d t2_ano(t="icol")" ;* "d "Var.ivar.irow"(t="icol")" Con_t2_ano(); "d t2_ano(t="icol")" ;* "d "Var.ivar.irow"(t="icol")" endif draw_tit((x1+x2)/2, y2-0.2, Month"-"Year); if (icol=Ncol);"xcbar "x2+0.05" "x2+0.25" "y1+0.4" "y2-0.4" -fs 1 -xo 0.1 -fw 0.08 -fh 0.1 -fx -0.0 -line on ";endif Sub_Reg2 Endwhile Main_tit((x00+x01)/2,y2+0.05,Varnm.ivar.irow" forecast initialed at "Mn_Ope"-"Yr_Ope) Endwhile "gxprint " Dat_out V.ivar"_forecast.png png" pull dd "c" Endwhile "close 1" ******************************************************* function Main_tit(x, y, title) "set string 1 c 5.5 0"; "set strsiz 0.11 0.13" "draw string "x" "y" "title "set string 1 c 5.5 0" return ******************************************************* function draw_tit(x, y, title) "set string 1 c 5.5 0"; "set strsiz 0.09 0.10" "draw string "x" "y" "title "set string 1 c 5.5 0" return ******************************************************* function draw_ver_tit(x, y, title) "set string 1 c 5.5 90"; "set strsiz 0.09 0.10" "draw string "x" "y" "title "set string 1 c 5.5 0" return ******************************************************* function Line_pl(mark,colo,thick,styl) "set gxout line" *"set gxout stnmark" "set cmark "mark "set ccolor "colo "set cthick "thick "set cstyle "styl return ******************************************** function Bar_pl(clor,gap) "set gxout bar" "set ccolor "clor "set bargap "gap "set barbase 0" "set baropts outline" "set cthick 2" "set baropts filled" *"set ccolor 1" return ******************************************** function Con_pr() levels = " 5 10 20 30 50 70 100 150 200 250 300 350 400 450 500 " "set clevs " levels "set ccolor 4" "set clab off" "set gxout contour" "set grads off" return ******************************************** function Pr_Clor() * "run color60_4rainfall.gs" * levels = " 5 10 20 30 40 50 60 70 80 90 " levels = " 5 10 20 30 50 70 100 120 140 160 180 200 240 280 320 350 400 450 500 " * "run color.gs -gxout grfill -kind green->lightgreen->yellow->orange->red -levs " levels * "run color.gs -gxout grfill -kind lightblue->blue->darkblue -levs " levels colors = " 32 36 40 48 50 52 56 58 60 62 64 66 68 70 74 76 80 84 88 90 " "set clevs " levels "set ccols " colors "set gxout grfill" "run color60_4rainfall.gs" "set grads off" return ******************************************** function Con_pr_ano() levels = " -100 -80 -60 -40 -30 -20 -10 -5 0 5 10 20 30 40 60 80 100 " "set clevs " levels "set ccolor 4" "set clab off" "set gxout contour" "set grads off" return ******************************************** function Pr_ano() levels = "-100 -80 -60 -40 -30 -20 -10 -5 0 5 10 20 30 40 60 80 100 " "run color.gs -gxout grfill -kind red->orange->yellow->white-(0)->white->lightblue->blue->darkblue -levs " levels return ******************************************** function Con_t2() levels = " 10 12 14 16 18 20 22 24 26 28 30 32 " "set clevs " levels "set ccolor 4" "set clab off" "set gxout contour" "set grads off" return ******************************************** function t2_Clor() levels = " 10 12 14 16 18 20 22 24 26 28 30 32 " "run color.gs -gxout grfill -kind darkgreen->green->lightgreen->yellow->orange->red->darkred -levs " levels "set grads off" return ******************************************** function Con_t2_ano() levels = " -2 -1.5 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 1.5 2 " "set clevs " levels "set ccolor 4" "set clab off" "set gxout contour" "set grads off" return ******************************************** function t2_ano() levels = " -2 -1.5 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 1.5 2 " "run color.gs -gxout grfill -kind darkgreen->green->lightgreen-(0)->yellow->orange->red -levs " levels "set grads off" return ********************************************