rm(list = ls())
a = 21.8745584  # 冰面
b = 7.66        #  冰面
T_d <- data.frame() 
T <- data$T2 
rh <- data$rh2 
et <- c()
# 初始值
for (i in 1:39) {
  e_s <- 6.1078 * exp((a * (T[i] - 273.15)/(T[i] - b)))
  e <- e_s * rh[i] / 100
  if(e < e_s){
    while(e < e_s){
         et[i] <- T[i] - 0.05
         e_s1 <- 6.1078 * exp((a * (et[i] - 273.15)/(et[i] - b)))
         e_s <- e_s1
         T[i] <- et[i]
         message('正在迭代温度:',i)
         if(e > e_s){
         T_d[i,1] <- T[i]}}
}else{
         T_d[i,1] <- T[i]}
  message(paste('正在迭代',i,sep=':'))
}
write.csv(T_d,'d:/T_d.csv') 
                                                                              2019-12-12