modflowpy / flopy

A Python package to create, run, and post-process MODFLOW-based models.
https://flopy.readthedocs.io
Other
517 stars 313 forks source link

writing .cbc files #200

Closed porterma closed 7 years ago

porterma commented 7 years ago

I recently installed the development version of flopy. It seems that now I am not getting updated .cbc files when I rerun models that were previously writing this file. Any thoughts on your end?

langevin-usgs commented 7 years ago

Can you paste into this issue the contents of your name file and output control file?

langevin-usgs commented 7 years ago

And the syntax of your flopy output control instantiation?

porterma commented 7 years ago

Here they are.

Name file for MODFLOW-2005, generated by Flopy version 3.2.7.

xul:0; yul:1000; rotation:0; proj4_str:+init=EPSG:4326; units:meters; lenuni:2; length_multiplier:1.0 ;start_datetime:1-1-1970

LIST 2 tutorial2.list DIS 11 tutorial2.dis BAS6 13 tutorial2.bas LPF 15 tutorial2.lpf PCG 27 tutorial2.pcg GHB 23 tutorial2.ghb WEL 20 tutorial2.wel OC 14 tutorial2.oc DATA(BINARY) 51 tutorial2.hds REPLACE DATA(BINARY) 52 tutorial2.ddn REPLACE

oc = flopy.modflow.ModflowOc(mf, stress_period_data=spd, compact=True)

OC package for MODFLOW-2005, generated by Flopy.

HEAD PRINT FORMAT 0 HEAD SAVE UNIT 51 DRAWDOWN PRINT FORMAT 0 DRAWDOWN SAVE UNIT 52 COMPACT BUDGET AUX

period 1 step 1 print head print budget save head save budget save drawdown

period 2 step 1 print head print budget save head save budget save drawdown

... period 3 step 100 print head print budget save head save budget save drawdown

langevin-usgs commented 7 years ago

Thanks, that all looks good. Need two other things:

  1. Top few lines of LPF or BCF MODFLOW file
  2. flopy code for instantiating LPF or BCF object
jdhughes-usgs commented 7 years ago

You need to pass in a non-zero ipakcb during instantiation for packages you want to save budget information. So in your lpf pass a ipakcb=53 and you will get cell-by-cell flows. We eliminated default head, drawdown, and budget output in version 3.2.6

porterma commented 7 years ago

Thank you. I apologize for not catching that.

bdestombe commented 7 years ago

Thanks, I was having some difficulties with this as well.

Currently, tutorial 1 and 2 in Examples/Tutorials don't run because of it.