Commit f8e091fa authored by unknown's avatar unknown
Browse files

v1.0.5.19 function TransfoParam_CemaNeige() cleaned

Showing with 30 additions and 20 deletions
+30 -20
Package: airGR Package: airGR
Type: Package Type: Package
Title: Suite of GR Hydrological Models for Precipitation-Runoff Modelling Title: Suite of GR Hydrological Models for Precipitation-Runoff Modelling
Version: 1.0.5.18 Version: 1.0.5.19
Date: 2017-01-23 Date: 2017-01-23
Authors@R: c( Authors@R: c(
person("Laurent", "Coron", role = c("aut", "trl")), person("Laurent", "Coron", role = c("aut", "trl")),
......
TransfoParam_CemaNeige <- function(ParamIn,Direction){ TransfoParam_CemaNeige <- function(ParamIn, Direction) {
NParam <- 2 NParam <- 2
Bool <- is.matrix(ParamIn) Bool <- is.matrix(ParamIn)
if(Bool==FALSE){ ParamIn <- rbind(ParamIn) }
if(ncol(ParamIn)!=NParam){ stop(paste("the CemaNeige module requires ",NParam," parameters \n",sep="")); return(NULL) } if (Bool == FALSE) {
ParamIn <- rbind(ParamIn)
if(Direction=="TR"){ }
ParamOut <- ParamIn;
ParamOut[,1] <- (ParamIn[,1]+9.99)/19.98 ### CemaNeige X1 (weighting coefficient for snow pack thermal state) if (ncol(ParamIn) != NParam) {
ParamOut[,2] <- exp(ParamIn[, 2]) / 200 ### CemaNeige X2 (degree-day melt coefficient) stop(paste( "the CemaNeige module requires ", NParam, " parameters \n", sep = "" ))
} return(NULL)
if(Direction=="RT"){ }
ParamOut <- ParamIn;
ParamOut[,1] <- ParamIn[,1]*19.98-9.99 ### CemaNeige X1 (weighting coefficient for snow pack thermal state) if (Direction == "TR") {
ParamOut[,2] <- log(ParamIn[,2] * 200) ### CemaNeige X2 (degree-day melt coefficient) ParamOut <- ParamIn
} ParamOut[, 1] <- (ParamIn[, 1] + 9.99) / 19.98 ### CemaNeige X1 (weighting coefficient for snow pack thermal state)
ParamOut[, 2] <- exp(ParamIn[, 2]) / 200 ### CemaNeige X2 (degree-day melt coefficient)
if(Bool==FALSE){ ParamOut <- ParamOut[1,] } }
if (Direction == "RT") {
ParamOut <- ParamIn
ParamOut[, 1] <- ParamIn[, 1] * 19.98 - 9.99 ### CemaNeige X1 (weighting coefficient for snow pack thermal state)
ParamOut[, 2] <- log(ParamIn[, 2] * 200) ### CemaNeige X2 (degree-day melt coefficient)
}
if (Bool == FALSE) {
ParamOut <- ParamOut[1, ]
}
return(ParamOut) return(ParamOut)
} }
Supports Markdown
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