|||
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
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-20 00:12
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社