diff --git a/R/RunModel.GR.R b/R/RunModel.GR.R
index 8cb78f8c1dcb703471c796d9fa4f365526fff3cd..234fdfe32411b5f66048cbee3cb258ddd23cf120 100644
--- a/R/RunModel.GR.R
+++ b/R/RunModel.GR.R
@@ -22,6 +22,8 @@ RunModel.GR <- function(x, RunOptions, Param, ...) {
     # All parameters
     iFirstParamRunOffModel <- 1
   }
+  # Avoiding Error in `FUN_MOD(x, RunOptions, Param)`: NA/NaN/Inf in foreign function call (arg 7)
+  RunOptions$IniStates[is.na(RunOptions$IniStates)] <- 0
 
   FUN_MOD <- match.fun(x$FUN_MOD)
   OutputsModel <- FUN_MOD(x, RunOptions = RunOptions,