PotEvap and Precip are reversed in RunModel_GR1A outputs
PotEvap and Precip are reversed in the RunModel_GR1A outputs.
The code of the GR1A is duplicated in the RunModel_GR1A()
function. The code is written in R and in Fortran. It is the R code that is used to run the model. All RunModel_*()
functions return first PotEvap and second Precip, but it is inversed in the GR1A R code (not in the Fortran code):
PEQ <- rbind(c(NA,NA,NA),cbind(P1,E1,Q1))
Maybe it is a good idea to remove this R code.