skipZero option in sensitivity

RBBGCMuso/R/musoMonte.R

@@ -31,6 +31,7 @@ musoMonte <- function(settings=NULL,
                      debugging = FALSE,
                      keepEpc = FALSE,
                      constrains = NULL,
+                     skipZero = TRUE,
                      ...){
 
 
@@ -150,6 +151,11 @@ musoMonte <- function(settings=NULL,
                 for(j in 1:numVars){
                     tmp2[j]<-funct[[j]](tmp[,j])
                 }
+                if(skipZero){
+                    if(tmp2[j]==0){
+                        tmp2[j] <- NA
+                    }
+                }
             } else {
                 for(j in 1:numVars){
                     tmp2[j]<-NA

RBBGCMuso/R/musoSensi.R

@@ -30,6 +30,7 @@ musoSensi <- function(monteCarloFile = NULL,
                      plotName = "sensitivity.png",
                      plotTitle = "Sensitivity",
                      skipSpinup = TRUE,
+                     skipZero = TRUE,
                      dpi=300){
 
     if(is.null(parameters)){
@@ -91,7 +92,8 @@ musoSensi <- function(monteCarloFile = NULL,
                       outVars = outVars,
                       fun = fun,
                       varIndex = varIndex,
-                      skipSpinup = skipSpinup
+                      skipSpinup = skipSpinup,
+                      skipZero=skipZero
                       )
         M <- cbind(seq_along(M[,1]),M)
         yInd <-  grep("mod.", colnames(M))[varIndex]

RBBGCMuso/man/getAnnualOutputList.Rd

@@ -0,0 +1,14 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/getOutPutList.R
+\name{getAnnualOutputList}
+\alias{getAnnualOutputList}
+\title{getAnnualOutputList}
+\usage{
+getAnnualOutputList(settings = NULL)
+}
+\arguments{
+\item{settings}{bla}
+}
+\description{
+bla bla
+}