integer ora,giorno,mese,anno,rec
	character data*60,file*80,title*60,source(3)*5
	parameter (nlon=480,nlat=480)
	real rain(nlon,nlat),slon,slat,elon,elat
	data source /'Radar','  WRF','MICRA'/

	call montemidiagrid(1,1,slon,slat,map)
	call montemidiagrid(nlon,nlat,elon,elat,map)
	call mvsetflags('Numero di Colori',20.0)
	call mvsetflags('Palette di colori',4.0)
	call mvsetflags('Colore Confini',1.0)
	call mvsetflags('data style',3.0)
	ora=10 ; giorno=20 ; mese=6 ; anno=2010
	rec=index1h(ora,giorno,mese,anno)
	call datafromhour (ora,giorno,mese,anno,data)
	do id=1,3
	   if (id.eq.1) then
	      call montemidiamap(rec,anno,rain,ngood)
	   else if (id.eq.2) then
	      call wrfmap(rec,anno,rain,ngood)
	   else
	      call micramap(rec,anno,rain,ngood)
	   endif
	   call settaiconfini(slon,elon,slat,elat)
	   xmin=1.0 ; xmax=vecmax(rain,nlon*nlat)
	   call boxplot(rain,nlon,nlat,xmin,xmax,2,20,' ')
	   call boxlegend(xmin,xmax,2,20,'f5.1')
	   call settaiconfini(slon,elon,slat,elat)
	   call confinidelleprovincie(slon,elon,slat,elat,1)
	   title=source(id)//' hourly precipitation estimation (mm)'
	   call displayexample('example09',data,title)
	enddo
	end