parameter (nmes=4000,nplot=3) ; real p(nmes),lat(nmes),lon(nmes)
	integer ora,giorno,mese,anno ; character data*60,str*60,title*60
	real slo(nplot),elo(nplot),sla(nplot),ela(nplot)
	data slo(1),elo(1),sla(1),ela(1) /07.8,13.3,43.5,46.0/
	data slo(2),elo(2),sla(2),ela(2) /11.5,14.5,41.2,43.6/
	data slo(3),elo(3),sla(3),ela(3) /13.5,17.5,38.2,42.6/
	ora=8 ; giorno=16 ; mese=12 ; anno=2011 
	call hourlyrain(ora,giorno,mese,anno,p,lat,lon,n)
	call datafromhour (ora,giorno,mese,anno,data)
	title='Precipitation (mm) - '//data(1:lenstr(data))
	do iplot=1,nplot
	   slon=slo(iplot) ; elon=elo(iplot) ; slat=sla(iplot) ; elat=ela(iplot)
	   call atlante (slon,elon,slat,elat,1,1,1,0,0)
c	   call gstxci(36)
	   do i=1,n 
	      write (str,'(f6.1)') p(i) ; call noinspace(str)
	      call wtstr(lon(i),lat(i),str(1:lenstr(str)),10,0,0)
	   enddo
	   call displayexample('example28',title,
	2	'CETEMPS - Integrated rain gauges database')
	enddo
	end