real x(366*24*4),y(366*24*4),x1(366*24*4),y1(366*24*4) character com*50,sens*50,str*60,title*128,latlon*60 integer month,year month=11 ; year=2010 i1=index15m(00,00,1,month,year) i2=index15m(45,23,monthlen(month,year),month,year) call monthofyear(month,year,str) do irec=571,575,2 call museoanagrafica('dewrain',irec,com,sens,xlat,xlon) call museotimeseries('dewrain',irec,year,x1,y1,n1) n=0 ; dx=1./(4*24) do i=i1,i2 if (nint(y1(i)).ne.-9999.and.nint(y1(i)).ne.-10) then n=n+1 x(n)=(i-i1)*dx ; y(n)=y1(i) endif enddo call plottailmese(x,y,366*24*4,1,n,sens(1:lenstr(sens))) call latlon2str(xlat,xlon,latlon) write (title,'(a,i4,a,i4)')'MuSEO dewetra data base, Rec. N.',irec, 2 ' - '//latlon(1:lenstr(latlon))//', year ',year call displayexample('example27',com(1:lenstr(com))//' - '//str,title) enddo end