Commit 977eb64b authored by unknown's avatar unknown
Browse files

v1.0.5.21 function TransfoParam_GR2M() cleaned

parent ef08ae63
Package: airGR
Type: Package
Title: Suite of GR Hydrological Models for Precipitation-Runoff Modelling
Version: 1.0.5.20
Date: 2017-01-23
Version: 1.0.5.21
Date: 2017-02-16
Authors@R: c(
person("Laurent", "Coron", role = c("aut", "trl")),
person("Charles", "Perrin", role = c("aut", "ths")),
......
TransfoParam_GR2M <- function(ParamIn,Direction){
NParam <- 2;
Bool <- is.matrix(ParamIn);
if(Bool==FALSE){ ParamIn <- rbind(ParamIn); }
if(ncol(ParamIn)!=NParam){ stop(paste("the GR2M model requires ",NParam," parameters \n",sep="")); return(NULL); }
if(Direction=="TR"){
ParamOut <- ParamIn;
ParamOut[,1] <- exp(ParamIn[,1]);
ParamOut[,2] <- ParamIn[, 2] / 4 + 2.5;
}
if(Direction=="RT"){
ParamOut <- ParamIn;
ParamOut[,1] <- log(ParamIn[,1]);
ParamOut[,2] <- (ParamIn[,2] - 2.5) * 4;
}
if(Bool==FALSE){ ParamOut <- ParamOut[1,]; }
return(ParamOut);
TransfoParam_GR2M <- function(ParamIn, Direction) {
NParam <- 2
Bool <- is.matrix(ParamIn)
if (Bool == FALSE) {
ParamIn <- rbind(ParamIn)
}
if (ncol(ParamIn) != NParam) {
stop(paste("the GR2M model requires ", NParam, " parameters \n", sep = ""))
return(NULL)
}
if (Direction == "TR") {
ParamOut <- ParamIn
ParamOut[, 1] <- exp(ParamIn[, 1])
ParamOut[, 2] <- ParamIn[, 2] / 4 + 2.5
}
if (Direction == "RT") {
ParamOut <- ParamIn
ParamOut[, 1] <- log(ParamIn[, 1])
ParamOut[, 2] <- (ParamIn[, 2] - 2.5) * 4
}
if (Bool == FALSE) {
ParamOut <- ParamOut[1, ]
}
return(ParamOut)
}
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment