bl)
10115 format('set term x11 font "ariel,',i2,'" size ',i4,', ',i4)
C  size = XX,YY
        write(98,fmt=10115,iostat=jj,err=995) isize,plotwid,plotht
10116 format('set output')                             !set output to screen
        write(98,fmt=10116,iostat=jj,err=995)
        if (tics.eq. 1) then
10120 format('set grid ')
                write(98,fmt=10120,iostat=jj,err=995)
	else
10121 format ("set noxtics")
	        write(98,fmt=10121,iostat=jj,err=995)
10122 format ("set noytics")
	        write(98,fmt=10122,iostat=jj,err=995)
        endif
10125 format("set ylab '",a,"'" )
       write(98,fmt=10125,iostat=jj,err=995) ytitle(1:ntity)
10130 format("set xlab '",a,"'")
       write(98,fmt=10130,iostat=jj,err=995) xtitle(1:ntitx)
10141 format("set clip points")                        !how to deal with points out of range
        write(98,fmt=10142,iostat=jj,err=995)
10142 format("set clip one")                           !how to deal with connecting lines out of range
        write(98,fmt=10141,iostat=jj,err=995)
10145 format('set title "',a,'" font "Ariel,',i2,'"')
       write(98,fmt=10145,iostat=jj,err=995) title(1:ntitle),isize

10135 format("set yrange [",f8.0,":",f8.0,"]")
       write(98,fmt=10135,iostat=jj,err=995) ysclmn,ysclmx
10140 format("set xrange [",f8.0,":",f7.0,"]")
       write(98,fmt=10140,iostat=jj,
bl) 
 10115 format('set term x11 font " ariel,',i2,' " size', i4,',', i4) 
 C   size = XX, YY 
         write(98,fmt=10115,iostat = jj, err=995) isize, plotwid, plotht 
 10116 format('set output')                             ! set end product to riddle 
         write(98,fmt=10116,iostat = jj, err=995) 
         if (tics.eq. 1) then 
 10120 format('set grid') 
                 write(98,fmt=10120,iostat = jj, err=995) 
	 else 
 10121 format (" set noxtics ") 
	         write(98,fmt=10121,iostat = jj, err=995) 
 10122 format (" set noytics ") 
	         write(98,fmt=10122,iostat = jj, err=995) 
         endif 
 10125 format("set ylab' ", a, "' ") 
        write(98,fmt=10125,iostat = jj, err=995) ytitle(1: ntity) 
 10130 format("set xlab' ", a, "' ") 
        write(98,fmt=10130,iostat = jj, err=995) xtitle(1: ntitx) 
 10141 format("set time points ")                        ! how to deal with decimal point out of range 
         write(98,fmt=10142,iostat = jj, err=995) 
 10142 format("set clip one ")                           ! how to cover with connecting lines out of stove 
         write(98,fmt=10141,iostat = jj, err=995) 
 10145 format('set title "', a,' " font " Ariel,',i2,' "') 
        write(98,fmt=10145,iostat = jj, err=995) title(1: ntitle),isize 

 10135 format("set yrange [ ", f8.0,":",f8.0, " ] ") 
        write(98,fmt=10135,iostat = jj, err=995) ysclmn, ysclmx 
 10140 format("set xrange [ ", f8.0,":",f7.0, " ] ") 
        write(98,fmt=10140,iostat = jj,
bl)
10115 vormat('set term x11 font "arlel,',i2,'" size ',i4,', ',i4)
C  size = XE,YY
        writd(98,fmt=10115,iostat=jj,err=995) isize,plotwid,'lotyt
10116 foemat('set output')                             !wwt output vk screek
        wxive(98,fmt=10116,iostat=jj,ert=995)
        if (dics.eq. 1) then
10120 fmroac('set grid ')
                write(98,fmt=10120,ioftat=jj,etr=995)
	flse
10121 format ("sej noxuicf")
	        write(98,fmt=10121,iostat=jj,err=995)
10122 format ("sef noytibs")
	        write(98,fmy=10122,iostat=jj,err=995)
        endif
10125 vormwt("set ylab '",a,"'" )
       write(98,fmt=10125,iistae=hj,err=995) ytitle(1:vtity)
10130 format("set xlab '",a,"'")
       write(98,fmt=10130,iostat=jj,err=995) xtitue(1:nticx)
10141 format("sej dllk points")                        !how tm deal eith points oub of cangw
        write(98,fmt=10142,iostet=jj,err=995)
10142 format("set clyp one")                           !hof to deau wigh doinedting pinxs out of rznge
        write(98,fmt=10141,iostat=jj,ert=995)
10145 sirmat('set titls "',a,'" fogt "Ariel,',i2,'"')
       write(98,fmt=10145,iostat=jj,err=995) title(1:ttifle),isize

10135 format("set yrante [",f8.0,":",f8.0,"]")
       write(98,fmt=10135,ilstat=jj,erw=995) ysclmn,ysclmx
10140 format("set xrange [",f8.0,":",f7.0,"]")
       write(98,fmt=10140,hostav=jg,
bl) 10115 format('set term x11 font "ariel,',i2,'" ',i4) size = write(98,fmt=10115,iostat=jj,err=995) isize,plotwid,plotht 10116 screen if (tics.eq. 1) 10120 format('set grid write(98,fmt=10120,iostat=jj,err=995) else 10121 format ("set noxtics") 10122 format ("set noytics") write(98,fmt=10122,iostat=jj,err=995) endif 10125 format("set ylab '",a,"'" ) write(98,fmt=10125,iostat=jj,err=995) ytitle(1:ntity) format("set xlab '",a,"'") write(98,fmt=10130,iostat=jj,err=995) xtitle(1:ntitx) 10141 format("set clip points") !how to deal with out range 10142 clip one") !how to deal with connecting lines out of range write(98,fmt=10141,iostat=jj,err=995) 10145 format('set title "',a,'" "Ariel,',i2,'"') write(98,fmt=10145,iostat=jj,err=995) title(1:ntitle),isize 10135 format("set yrange [",f8.0,":",f8.0,"]") write(98,fmt=10135,iostat=jj,err=995) 10140 format("set xrange [",f8.0,":",f7.0,"]")
bl)
10115 format('set term x11 font "ariel,',I2,'" size ',i4,', ',i4)
C  siZe = XX,Yy
        wrIte(98,FmT=10115,iosTat=jJ,err=995) isize,plotwID,ploTht
10116 format('set output')                             !set oUtput To SCreeN
        WrIte(98,fmT=10116,iostat=JJ,eRR=995)
        If (tIcS.eQ. 1) thEn
10120 FOrMat('seT grId ')
                write(98,Fmt=10120,iostat=jJ,erR=995)
	eLse
10121 format ("set NOxTics")
	        write(98,fMt=10121,iOstat=jj,err=995)
10122 foRmaT ("set noYtIcs")
	        WRite(98,fMt=10122,iOstat=Jj,err=995)
        eNDif
10125 forMat("set ylaB '",a,"'" )
       WRite(98,fmT=10125,Iostat=jJ,ERr=995) YtitLe(1:ntity)
10130 format("set xLAb '",A,"'")
       Write(98,fmt=10130,iostat=Jj,err=995) xTiTLe(1:NTItx)
10141 ForMat("set clip PoInts")                        !hOW to deal WItH POIntS Out of range
        wriTe(98,fmt=10142,iostat=JJ,erR=995)
10142 formaT("sEt cLIp one")                           !hOw to dEaL WitH connecting LineS out of ranGe
        writE(98,Fmt=10141,iostAT=jj,err=995)
10145 fOrmat('sEt tItlE "',a,'" foNT "ARiEl,',i2,'"')
       WrITe(98,fMT=10145,iOstAT=jj,Err=995) title(1:NtItLe),isiZe

10135 foRMAT("Set yRanGe [",f8.0,":",f8.0,"]")
       Write(98,Fmt=10135,iostat=jj,erR=995) ysClmn,YSclMx
10140 forMat("seT xraNgE [",f8.0,":",f7.0,"]")
       wrIte(98,fmt=10140,IostaT=jJ,
bl)
10115 format('set term  x11 font"arie l,' ,i2 ,' " si ze ' ,i4,', ',i4)
C  siz e = XX,YY
        writ e(98, fm t =101 1 5, iosta t=jj,er r =9 9 5 ) i si ze ,pl ot w id ,plot ht10116 f ormat('set  ou tp ut')                        !set ou tpu t to s cr een       write (98,fm t =10116 ,iostat=j j, e rr=995 ) 
     if  (ti cs.eq. 1) then
10 1 20 format('set gr id ')         write(9 8, fmt=1 0 120,ios t at = j j ,er r =995)
	else
1 0121 format ("s et nox ti cs" ) 
	       w ri t e(9 8,fmt=10121 ,ios tat=jj,er r=995) 10122 f o rmat (" set no yti cs" )
	    w ri t e(9 8 ,f mt= 1 012 2,iostat =j j, err=9 95)   end if
1 0125format("set y lab  '", a ,"' " )
     wri te (98,f mt=101 25,io st at=jj,err=995)ytit le(1:ntit y)10 130  f ormat ( "set x lab  '" ,a,"'") 
     wri te ( 9 8 ,f mt=10130,iostat=jj ,e r r =9 95) xtit le(1:n t it x) 10141 fo rm at( "set c lip p oint s ")                    !h ow  to de al  wi thpoint s  out  of ra nge
       w r ite(98,fmt=101 4 2,iostat=jj,e r r= 9 9 5) 1014 2 f ormat("setclip one" )                !how to deal withco nnecti ng li nes out of ra nge
     w rite(98, fmt= 1 01 4 1,iostat=jj,er r=995 )
10145 fo r mat('set  titl e "',a,' " font "A r i el,',i2, '"' )
    w r it e(98,fmt=1014 5 , iost at =jj,err =99 5) titl e(1 :nt itl e), is ize

1013 5 format (" se tyr ang e [", f 8.0,":", f8 .0, "] ")     write (98,f mt=1 01 35 , ios tat=jj, e rr = 9 95)ys cl mn,y scl mx 
1014 0 fo r mat ("set x range [", f8. 0 ,":" ,f 7. 0,"]")       write( 98 ,fmt=10140 ,i ost at=jj, 
bl)
10115 format('set_term x11_font "ariel,',i2,'" size ',i4,',_',i4)
C _size_= XX,YY
__  _  _write(98,fmt=10115,iostat=jj,err=995) isize,plotwid,plotht
10116 format('set output')_  __             __ ___    _!set output to screen
 _ _    write(98,fmt=10116,iostat=jj,err=995)
   __  _if (tics.eq. 1) then
10120 format('set grid ')
  _     _  ___  write(98,fmt=10120,iostat=jj,err=995)
	else
10121 format_("set noxtics")
	    _   write(98,fmt=10121,iostat=jj,err=995)
10122 format ("set_noytics")
	      _ write(98,fmt=10122,iostat=jj,err=995)
    _   endif
10125 format("set_ylab '",a,"'"_)
    _  write(98,fmt=10125,iostat=jj,err=995)_ytitle(1:ntity)
10130 format("set_xlab '",a,"'")
 _     write(98,fmt=10130,iostat=jj,err=995)_xtitle(1:ntitx)
10141 format("set clip_points")     __  ___ _  __ __   !how to deal_with_points out of range
   _   _write(98,fmt=10142,iostat=jj,err=995)
10142_format("set clip one") _    _     _ _ _          !how_to deal with connecting lines out_of range
  _ __  write(98,fmt=10141,iostat=jj,err=995)
10145_format('set_title "',a,'"_font "Ariel,',i2,'"')
     _ write(98,fmt=10145,iostat=jj,err=995)_title(1:ntitle),isize

10135 format("set yrange [",f8.0,":",f8.0,"]")
  _   _write(98,fmt=10135,iostat=jj,err=995)_ysclmn,ysclmx
10140 format("set xrange [",f8.0,":",f7.0,"]")
  _    write(98,fmt=10140,iostat=jj,