dados_dia <- data.frame(dados_dia)

    ICL <- c(0,0.3,0.6,0.9)
    VEL <- c(0,2,4,8)
    URE <- c(0,70,80,90,100)
    TPR <- c(0,14,19,24,40)



    for(i in 2:4){
        for(j in 2:4){
            for(k in 2:5){
                for(l in 2:5){

print(c('ICL > ',ICL[i-1],' e ICL <= ',ICL[i]))
  print(c('VEL > ', VEL[j-1],' e VEL <= ',VEL[j]))
    print(c('URE > ',URE[k-1],' e URE <= ',URE[k]))
      print(c('TPR > ',TPR[l-1],' e TPR <= ',TPR[l])) ## essas são as combinações


print(c(dados_dia$date[dados_dia$I.Cla > ICL[i-1] &                                    dados_dia$I.Cla <= ICL[i] &
                                                      dados_dia$speed > VEL[j-1] &
                                                      dados_dia$speed <= VEL[j] &
                                                      dados_dia$rH > URE[k-1] &
                                                      dados_dia$rH <= URE[k] &
                                                      dados_dia$Tair > TPR[l-1] &
                                                      dados_dia$Tair <= TPR[l]] )) ## Aqui sai as datas em que ocorrem cada combinação acima
                print('')
                print('')
                print('')
                }
            }
        }
    }