function main(arg)
*datum=subwrd(arg,1)
rc = gsfallow("on")
*'open ulrichobs.ctl'
'open linterpall_ulrich_new.ctl'
'set display color white'
'set gxout shaded'
'c'
'set grads off'
*----- odredjivanje start date ------
*'set t 1 120'
'set t 1 89'
'c'
'set grads off'
'set z 1 70'
'set x 75'

*'set clab forced'
'set gxout shaded'
*'set clopts 1 1 0.25'
'run col_prec.gs'
'set clevs 0.01 0.1 0.5 1 2 3 4 8' 
'set clevs 0.0001 0.001 0.01 0.1 0.5 1 2 3 4 8' 
'set clevs 0.001 0.01 0.1 0.5 1 2 3 4 8' 
'set clevs 0.001 0.01 0.1 1 2 4 6 8 10 ' 
'set clevs  0.01 0.1 0.5 1 2 4 6 8 10 ' 
*'d log10(fimmint+1)'
'd log10(fimmui+1)'
'set gxout contour'
'set clab forced'
'set ccolor 2'
'd tint-273'

 'set strsiz 0.12 0.13'
 'set string 1 bl 5'
'q gxinfo'
   xline=sublin(result,3)
   yline=sublin(result,4)
   rec3 = sublin(result,3)
   rec4 = sublin(result,4)
   xlb=subwrd(xline,4)
   ylb=subwrd(yline,4)
   xur=subwrd(xline,6)
   yur=subwrd(yline,6)
   rec3 = sublin(result,3)
   rec4 = sublin(result,4)

   naslov1='DREAM8-assim: LOG10(number of ice nuclei due to dust) and temperature (`3.`0C)'
   naslov2='Station: Soverato'

'draw string 'xlb' 'yur+0.3'  'naslov1
'draw string 'xlb' 'yur+0.1'  'naslov2

*'run logobw.gs 'xur-0.7' 'yur+0.1' 0.7'
*'run bwfulllogo.gs 'xur-1.2' 'yur+0.35' 1.2'
***'run bwfulllogo.gs 'xur-1.' 'yur+0.25' 1.'

*'run bwfulllogo_bojan.gs 'xur-1.' 'yur+0.1' 1.'
*'run bwfulllogo_bojan.gs 'xur-1.45' 'yur+0.1' 1.35'
*'run bwfulllogo_bojan.gs 'xlb' 'ylb-0.33' 1.35'

   ylo = subwrd(yline,4)-0.55
   ymx=ylo+0.07
   xlo = subwrd(xline,4)
   xmx = subwrd(xline,6)
   fh=0.1  
   fw=0.1
   'xcbar 'xlo' 'xmx' 'ylo-0.1' 'ymx-0.1' -fh 'fh' -fw 'fw   

*   rc=arrow(xlo,ylo,0.5,20)

  'set line 1 1 1'

*  'draw shp /clima/grads/data/RWDB_Ad1-Ln.shp'

 'print q.eps'
 '!convert -rotate 90 -density 300x300 +antialias -flatten q.eps p.gif'
 '!convert p.gif -resize 40% +antialias q.gif'
 '!cp q.gif fimm_dust_alldays.gif'
 '!rm -v q.eps'
 '!rm -v p.gif'

'close 1'

*'quit'

function arrow(x,y,len,scale)
'set line 1 1 4'
'draw line 'x-len/2.' 'y' 'x+len/2.' 'y
'draw line 'x+len/2.-0.05' 'y+0.025' 'x+len/2.' 'y
'draw line 'x+len/2.-0.05' 'y-0.025' 'x+len/2.' 'y
'set string 1 c'
'set strsiz 0.1'
'draw string 'x' 'y-0.1' 'scale
return