# NOT RUN {
#--------------------------------------------------------
# Example 1: Figure 6 (right) from S<U+00E1>enz et al. (2020)
#--------------------------------------------------------
# Load the data
data("WRF")
# Parameters
Uxlim=c(60,110)
Uylim=c(0,50)
Uxlab<-"Qx (kg/m/s)"
Uylab<-"Qy (kg/m/s)"
plotmain<-"Water vapour transport"
sfactor<-0.1
ref<-WRF[["ref"]][,c("mod","U","V")]
mod<-WRF[["mod"]][,c("mod","U","V")]
isBad=((is.na(ref$U)) |(is.na(ref$V)))
isOK=(!isBad)
ref<-ref[isOK,]
nmod1<-which(mod$mod=="WrfN")
mod1<-mod[nmod1,][isOK,]
nmod2<-which(mod$mod=="WrfD")
mod2<-mod[nmod2,][isOK,]
nmod3<-which(mod$mod=="ERAI")
mod3<-mod[nmod3,][isOK,]
mod<-rbind(mod1,mod2,mod3)
# Index
sIWRF=SailoR.Indices(ref,mod)
# Index table
sIWRF.table<-SailoR.Table(sIWRF)
# plot
SailoR.Plot(ref, mod, ColourList=NULL, sfactor, docenter=TRUE,
Uxlim, Uylim, Uxlab, Uylab, plotmain, plotlegend=TRUE,
plotRMSElegend=TRUE, plotscalelegend=TRUE,
RMSE_legend_Rounding=0, RMSE_legend_units = " kg/m/s",
referenceName="Sounding")
# }
Run the code above in your browser using DataLab