expertsArray = c(
	"CLM", #1
	"GAM", #2
	"CLM1", #3
	"CLM2", #4
	"S_AIRPARIF", #5
	"S_INERIS", #6
	"D_ESMERALDA", #7
	"D_PREVAIR", #8
	"D_PREVAIR2", #9
	"PERSIST", #10
	#new additions only for PQV_2014
	#TODO: default behavior on station != PQV_2014 ?
	"GAM_sepMar", #11
	"GAM_aprAug", #12
	"GAM_highPollution", #13
	"GAM_lowPollution", #14
	"GAM_hotTemperature", #15
	"GAM_coldTemperature", #16
	"GAM_eastWind", #17
	"GAM_westWind", #18
	"GAM_noRain", #19
	"GAM_rain" #20
)

stationsArray = c(
	"AIL", #1
	"ALE", #2
	"CAE", #3
	"CHD", #4
	"EVT", #5
	"HRI", #6
	"IFS", #7
	"JUS", #8
	"LIS", #9
	"MAS", #10
	"MRA", #11
	"NEI", #12
	"POS", #13
	"PQV", #14
	"SLO", #15
	"HRI_2014", #16
	"LIS_2014", #17
	"PQV_2014", #18
	"PQV2" #19
)

#' @title Get forecasts + observations
#'
#' @description Get forecasts of all specified experts for all specified stations, also with (ordered) dates and (unordered) stations indices.
#'
#' @param experts Names of the experts. Default: all
#' @param station Names of the stations. Default: all
#'
#' @export
getData = function(experts=expertsArray, stations=stationsArray)
{
	#no need because of "LazyData: true" in DESCRIPTION
	#data(list=stations, package="aggexp")
	data = as.data.frame(matrix(nrow=0, ncol=1 + length(experts) + 2))
	names(data) = c("Date", experts, "Measure", "Station")
	for (i in 1:length(stations))
	{
		stationInfo = get(stations[i])
		#date index is sufficient; also add station index
		stationInfo = cbind(Date = 1:nrow(stationInfo), stationInfo[,names(stationInfo) %in% experts], Measure = stationInfo[,"Measure"], Station = i)
		data = rbind(data, stationInfo)
	}

	#extra step: order by date (would be a DB request)
	data = data[order(data[,"Date"]),]

	return (data)
}