use chymdata , only : nlon,nlat,wkm1,mchym,rchym
	character source*50 
	parameter (nmes=5000) ; real p(nmes),lat(nmes),lon(nmes)
	integer ora,giorno,mese,anno ; character data*60
	logical realrange
	data ora,giorno,mese,anno /10,20,06,2010/ 
	call datafromhour (ora,giorno,mese,anno,data)
	call chymreadstaticfields('tmp/Po.chym',0)
	call chymplot(1,'c2L')
	slon=rchym(1) ; slat=rchym(2) ; dij=rchym(5) ; dji=rchym(5)
	elon=slon+(nlon-1)*dij ; elat=slat+(nlat-1)*dji
	call hourlyrain(ora,giorno,mese,anno,p,lat,lon,n)
	wkm1=-9999
	do ii=1,n
	   if (realrange(lat(ii),slat,elat).and.realrange(lon(ii),slon,elon)) 
	2	then
	      call locateij(lat(ii),lon(ii),slat,slon,dij,dji,nlat,nlon,i,j)
	      wkm1(i,j)=10.5
	   endif
	enddo
	call mvsetflags('Boxplot style',2.0)
	call boxplot(wkm1,nlon,nlat,10.0,11.0,36,36,' ')
	call mvsetflags('Logo',2.0)
	call displayexample('example31','Location of DEWETRA sensors',data)
	end