diff --git a/DESCRIPTION b/DESCRIPTION
index 329ad2511772ae0cb0cf466e1c8044ff39f31ef0..154ab031a52779d6dabdf328a56c857171bb1d1f 100755
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,8 +1,8 @@
 Package: hydroPSO
 Type: Package
 Title: Model-Independent Particle Swarm Optimisation for Environmental Models
-Version: 0.1-56-2
-Date: 2012-06-25
+Version: 0.1-56-3
+Date: 2012-06-28
 Author: Mauricio Zambrano-Bigiarini [aut, cre] and Rodrigo Rojas [ctb]
 Author@R: c(person("Mauricio", "Zambrano-Bigiarini", email = "mzb.devel@gmail.com", role=c("aut", "cre") ), person("Rodrigo", "Rojas", email = "Rodrigo.RojasMujica@gmail.com", role=c("ctb") ) )
 Maintainer: Mauricio Zambrano-Bigiarini <mzb.devel@gmail.com>
diff --git a/NAMESPACE b/NAMESPACE
index 8b23decadf2cf231051ac93c69cc0b2ffb98512a..ed5594c89c3607bfcd6d07d2611027e5bba65c95 100755
--- a/NAMESPACE
+++ b/NAMESPACE
@@ -1,8 +1,9 @@
-#importFrom("hydroTSM", hydroplot, hydropairs, dip, mip, yip, vector2zoo)
 importFrom("Hmisc", wtd.quantile, wtd.Ecdf)
 importFrom("sp", coordinates, spplot)
 importFrom("lattice", xyplot, axis.default)
 import(grid)
+#importFrom("hydroTSM", hydroplot, hydropairs, dip, mip, yip, vector2zoo)
+#importFrom("hydroGOF", ggof)
 #importFrom("scatterplot3d", scatterplot3d)
 #exportPattern("^[^\\.]")
 export(hydroPSO,
diff --git a/NEWS b/NEWS
index 651ad33db870f59f094ead417abdc060ec975bd7..f44cd310ef8c505f23cf8b32f0c403b3b36b2502 100755
--- a/NEWS
+++ b/NEWS
@@ -16,6 +16,13 @@ NEWS/ChangeLog for hydroPSO
                           -) default value for the exponent of the time-varying lambda was changed from 'TVlambda.exp= 1.7' to 'TVlambda.exp= 1'.
                              However, by default 'use.TVlambda= FALSE' still.
                           -) source code was tidy up
+        o 'plot_out'    : -) Fixed error that appeared when the 'plot_out' function was called from an R script running directly on the OS console 
+                             (not within the R editor), with the following error message:
+                             "Error in plot2(x = sim, y = obs, plot.type = "single", main = paste("Daily",  : 
+                              could not find function "hasArg"
+                             Calls: plot_results -> plot_out -> ggof -> plot2"
+        o 'plot_results': -) Same as immediately above (because 'plot_results' internally calls 'plot_out')
+                           
                              
 0.1-56	14-Jun-2012
         o 'hydroPSO'    : -) much less memory consumption for large values of 'maxit' and/or 'npart' (thanks to  Paul Smith for reporting it !)
diff --git a/R/plot_out.R b/R/plot_out.R
index fa698c1e1fa3c28fd98d2b851bacfdd545c23242..ef1ea4853f88c22f28583ee0fc6ad388f7752296 100755
--- a/R/plot_out.R
+++ b/R/plot_out.R
@@ -17,7 +17,7 @@
 ################################################################################
 # Author : Mauricio Zambrano-Bigiarini                                         #  
 # Started: 03-Feb-2012                                                         #
-# Updates: 15-Feb-2012 ; 22-Feb-2012                                           #        
+# Updates: 15-Feb-2012 ; 22-Feb-2012 ; 28-Jun-2012                             #        
 ################################################################################
 
 plot_out <- function(sim, obs, 
@@ -82,6 +82,12 @@ plot_out <- function(sim, obs,
     } else if ( !(MinMax %in% c("min", "max")) )
              stop("Invalid argument: 'MinMax' must be in c('min', 'max')")
   } # IF end
+  
+  # Checking 'hydroGOF' pacakge when ptype=="ts"
+  if (ptype=="ts") {
+   if (!require(hydroGOF))
+     stop("Invalid argument: You don't have the 'hydroGOF' package => You can not use 'ptype='ts' !!")
+  } # IF end
            
   # Checking 'class(sim)'    
   if ( (ptype=="corr") | (ptype=="ts") ) {
@@ -125,9 +131,6 @@ plot_out <- function(sim, obs,
     
   } else if (ptype=="ts") {
       main="Observed vs 'best' Simulation"
-             
-      if (!require(hydroGOF))
-        stop("Invalid argument: You don't have the 'hydroGOF' package => You can not use 'ptype='ts' !!")
        
       if (!is.null(dates)) {
         if (!is.zoo(obs)) obs <- zoo(obs, dates) 
@@ -135,6 +138,7 @@ plot_out <- function(sim, obs,
       } # IF end
       
       # Plotting Sim vs Obs
+      library(hydroGOF)
       ggof(sim=sim, obs=obs, main=main, ftype=ftype, FUN=FUN, 
            cex.main=cex.main, cex.axis=cex.axis, cex.lab=cex.lab, leg.cex=leg.cex) # dates.fmt="%Y-%m-%d")