R/thermodynamics.LI.R

In metno/wgen: Stochastic Weather Generator

#' thermodynamics.LI
#'
#' LI: Lifted Index: T500-Tmoistadiab500
#' @param T
#' @param Twbp
#' @param p
#' @keywords thermodynamics
#' @export
#' @examples
#' thermodynamics.LI()
thermodynamics.LI<-function(T,Twbp,p){
   T500<-T[which(p==500)]
   Tmoistadiab500<-thermodynamics.moistadiab(Twbp,500)
  
   results<-T500-Tmoistadiab500
   return(results)
}