Just saw the last two postings on this thread. I agree with Laura on the preference of using fc...+003 than fc...+000 for non-synoptic hours. In my view, if the purpose is to compare model to observations every three, or even every 1 hour, it is not a good idea, given the current situation, to perform 3 hourly or one hourly analysis.
Currently for the historical episode we are unable to find higher frequency lateral boundary data from ECMWF than 6-hourly analyses. Accordingly, you have received 4 observation/day data package centered around synoptic hour (00/06/12/18 UTC). Changing cycling interval will require quite some additional experiment runs with need for numerous modifications of scripts in order to fetch available observation and lateral boundary data.
On top of these, you are not likely to get better model data if that's what you are interested. I mean,e.g.
fc2006091503+000
will unlikely be of higher quality than
fc2006091500+003
In the above example, obtaining first model data requires quite some extra effort, but it may not be of better quality than the second one, which you already have now. This is due to typical shortcomings in assimilation technique and in moisture spin-up. Furthermore, if the quantity you want to compare is precipitation or other things that are closely related to moisture process, I'd actually recommend use of forecasts with longer forecast lead time than 6 h.