Error reading .dis file

[tungttnguyen]

I'm getting the following errors when reading the attached .dis file. Can you help me? Thanks!

dis <- read_dis(paste0("example", ".dis"))
Error in gsub(paste0("([[:alnum:][:punct:][:space:]]{", nPerNum, "})"), :
invalid regular expression '([[:alnum:][:punct:][:space:]]{NA})', reason 'Invalid contents of {}'
In addition: Warning messages:
1: NAs introduced by coercion
2: In read_modflow_array(dis_lines, 1, dis$nrow, 1, ndim = 1) :
NAs introduced by coercion

Session Info

devtools::session_info()
Session info ----------------------------------------------------------------------------------------------------------------------------
setting value
version R version 3.4.0 (2017-04-21)
system x86_64, linux-gnu
ui RStudio (1.0.143)
language (EN)
collate en_US.UTF-8
tz America/Los_Angeles
date 2017-04-28
Packages --------------------------------------------------------------------------------------------------------------------------------
package * version date source
acepack 1.4.1 2016-10-29 cran (@1.4.1)
akima 0.6-2 2016-12-20 cran (@0.6-2)
animation 2.5 2017-03-30 cran (@2.5)
automap 1.0-14 2013-08-29 cran (@1.0-14)
backports 1.0.5 2017-01-18 cran (@1.0.5)
base64enc 0.1-3 2015-07-28 cran (@0.1-3)
checkmate 1.8.2 2016-11-02 cran (@1.8.2)
class 7.3-14 2015-08-30 CRAN (R 3.2.2)
cluster 2.0.6 2017-03-16 CRAN (R 3.3.3)
colorspace 1.3-2 2016-12-14 cran (@1.3-2)
curl 2.5 2017-04-14 CRAN (R 3.4.0)
data.table 1.10.4 2017-02-01 cran (@1.10.4)
DEoptim 2.2-4 2016-12-19 cran (@2.2-4)
devtools * 1.12.0 2016-12-05 CRAN (R 3.4.0)
digest 0.6.4 2013-12-03 CRAN (R 3.0.2)
directlabels 2017.03.31 2017-04-08 cran (@2017.03)
e1071 1.6-8 2017-02-02 cran (@1.6-8)
FNN 1.1 2013-07-31 cran (@1.1)
foreign 0.8-67 2016-09-13 CRAN (R 3.3.1)
Formula 1.2-1 2015-04-07 cran (@1.2-1)
ggplot2 2.2.1 2016-12-30 cran (@2.2.1)
git2r 0.18.0 2017-01-01 CRAN (R 3.4.0)
gridExtra 2.2.1 2016-02-29 cran (@2.2.1)
gstat 1.1-5 2017-03-12 cran (@1.1-5)
gtable 0.2.0 2016-02-26 cran (@0.2.0)
Hmisc 4.0-2 2016-12-31 cran (@4.0-2)
hms 0.3 2016-11-22 cran (@0.3)
htmlTable 1.9 2017-01-26 cran (@1.9)
htmltools 0.3.5 2016-03-21 CRAN (R 3.4.0)
htmlwidgets 0.8 2016-11-09 cran (@0.8)
httpuv 1.3.3 2015-08-04 CRAN (R 3.4.0)
httr 1.2.1 2016-07-03 CRAN (R 3.4.0)
hydroGOF 0.3-8 2014-02-04 cran (@0.3-8)
hydroPSO 0.3-4 2014-04-13 cran (@0.3-4)
hydroTSM 0.4-2-1 2014-01-23 cran (@0.4-2-1)
intervals 0.15.1 2015-08-27 cran (@0.15.1)
jsonlite 1.4 2017-04-08 CRAN (R 3.4.0)
knitr 1.15.1 2016-11-22 cran (@1.15.1)
lattice 0.20-35 2017-03-25 CRAN (R 3.3.3)
latticeExtra 0.6-28 2016-02-09 cran (@0.6-28)
lazyeval 0.2.0 2016-06-12 cran (@0.2.0)
lhs 0.14 2016-08-09 cran (@0.14)
magrittr 1.5 2014-11-22 cran (@1.5)
Matrix 1.2-8 2017-01-20 CRAN (R 3.3.2)
memoise 1.1.0 2017-04-21 CRAN (R 3.4.0)
mime 0.5 2016-07-07 CRAN (R 3.4.0)
munsell 0.4.3 2016-02-13 cran (@0.4.3)
nnet 7.3-12 2016-02-02 CRAN (R 3.2.3)
plyr 1.8.4 2016-06-08 cran (@1.8.4)
quadprog 1.5-5 2013-04-17 cran (@1.5-5)
R6 2.2.0 2016-10-05 CRAN (R 3.4.0)
RColorBrewer 1.1-2 2014-12-07 cran (@1.1-2)
Rcpp 0.12.10 2017-03-19 cran (@0.12.10)
readr 1.1.0 2017-03-22 cran (@1.1.0)
reshape 0.8.6 2016-10-21 cran (@0.8.6)
rgdal 1.2-7 2017-04-25 CRAN (R 3.4.0)
rgl 0.98.1 2017-03-08 cran (@0.98.1)
RMODFLOW * 0.4.0 2017-04-28 Github (8a9578d)
rpart 4.1-11 2017-04-21 CRAN (R 3.4.0)
scales 0.4.1 2016-11-09 cran (@0.4.1)
shiny 1.0.3 2017-04-26 CRAN (R 3.4.0)
sp 1.2-4 2016-12-22 cran (@1.2-4)
spacetime 1.2-0 2016-09-03 cran (@1.2-0)
stringi 1.1.5 2017-04-07 cran (@1.1.5)
stringr 1.2.0 2017-02-18 cran (@1.2.0)
survival 2.41-3 2017-04-04 CRAN (R 3.3.3)
tibble 1.3.0 2017-04-01 cran (@1.3.0)
withr 1.0.2 2016-06-20 CRAN (R 3.4.0)
xtable 1.8-2 2016-02-05 cran (@1.8-2)
xts 0.9-7 2013-08-29 CRAN (R 3.0.1)
zoo 1.8-0 2017-04-12 cran (@1.8-0)

example.dis.tar.gz

-- Tung

[rogiersbart]

Hi Tung,

fixed-format control records (see online guide to MODFLOW) cannot be dealt with at the moment. I'm planning a major update of RMODFLOW in the near future, and it's one of the things I will be implementing. If you think you know what changes are required, you are of course free to create a pull request.

Kind regards,
Bart

[cneyens]

Note that reading fixed-format array headers is supported in the current development version. Right now, you'll need to supply an additional rmf NAM object to rmf_read_dis in order to do this, which might require further testing.